Treating Vacuum Regions

I am simulating a micrometric flow problem which generates vacuum regions near solid boundaries. But if initialise my lbm code with vacuum regions, it creates strong density gradients and flow to these regions. Is there any effective method to treat such kind of problems ? Please point me to proposed theories / papers in this area.

Thanks and Regards,
Omkar