The boundary layer equations for plane, incompressible, and steady flow are 

$$\matrix{ {u{{\partial u} \over {\partial x}} + v{{\partial u} \over {\partial y}} = - {1 \over \varrho }{{\partial p} \over {\partial x}} + v{{{\partial ^2}u} \over {\partial {y^2}}},} \cr {0 = {{\partial p} \over {\partial y}},} \cr {{{\partial u} \over {\partial x}} + {{\partial v} \over {\partial y}} = 0.} \cr }$$

Boundary Layer Theory

An efficient ghost-cell immersed boundary method (GCIBM) for simulating turbulent flows in complex geometries is presented. A boundary condition is enforced through a ghost cell method. The reconstruction procedure allows systematic development of numerical schemes for treating the immersed boundary while preserving the overall second-order accuracy of the base solver. Both Dirichlet and Neumann boundary conditions can be treated. The current ghost cell treatment is both suitable for staggered and non-staggered Cartesian grids. The accuracy of the current method is validated using flow past a circular cylinder and large eddy simulation of turbulent flow over a wavy surface. Numerical results are compared with experimental data and boundary-fitted grid results. The method is further extended to an existing ocean model (MITGCM) to simulate geophysical flow over a three-dimensional bump. The method is easily implemented as evidenced by our use of several existing codes.

A Ghost-Cell Immersed Boundary Method for Flow in Complex Geometry

The lattice Boltzmann method (LBM) has experienced tremendous advances and has been well accepted as a useful method to simulate various fluid behaviors. For computational microfluidics, LBM may present some advantages, including the physical representation of microscopic interactions, the uniform algorithm for multiphase flows, and the easiness in dealing with complex boundary. In addition, LBM-like algorithms have been developed to solve microfluidics-related processes and phenomena, such as heat transfer, electric/magnetic field, and diffusion. This article provides a practical overview of these LBM models and implementation details for external force, initial condition, and boundary condition. Moreover, recent LBM applications in various microfluidic situations have been reviewed, including microscopic gaseous flows, surface wettability and solid–liquid interfacial slip, multiphase flows in microchannels, electrokinetic flows, interface deformation in electric/magnetic field, flows through porous structures, and biological microflows. These simulations show some examples of the capability and efficiency of LBM in computational microfluidics.

Lattice Boltzmann method for microfluidics: models and applications

Accurate modeling of gas microflow is crucial for the microfluidic devices in MEMS. Gas microflows through these devices are often in the slip and transition flow regimes, characterized by the Knudsen number of the order of 10−2~100. An increasing number of researchers now dedicate great attention to the developments in the modeling of non-equilibrium boundary conditions in the gas microflows, concentrating on the slip model. In this review, we present various slip models obtained from different theoretical, computational and experimental studies for gas microflows. Correct descriptions of the Knudsen layer effect are of critical importance in modeling and designing of gas microflow systems and in predicting their performances. Theoretical descriptions of the gas-surface interaction and gas-surface molecular interaction models are introduced to describe the boundary conditions. Various methods and techniques for determination of the slip coefficients are reviewed. The review presents the considerable success in the implementation of various slip boundary conditions to extend the Navier–Stokes (N–S) equations into the slip and transition flow regimes. Comparisons of different values and formulations of the first- and second-order slip coefficients and models reveal the discrepancies arising from different definitions in the first-order slip coefficient and various approaches to determine the second-order slip coefficient. In addition, no consensus has been reached on the correct and generalized form of higher-order slip expression. The influences of specific effects, such as effective mean free path of the gas molecules and viscosity, surface roughness, gas composition and tangential momentum accommodation coefficient, on the hybrid slip models for gas microflows are analyzed and discussed. It shows that although the various hybrid slip models are proposed from different viewpoints, they can contribute to N–S equations for capturing the high Knudsen number effects in the slip and transition flow regimes. Future studies are also discussed for improving the understanding of gas microflows and enabling us to exactly predict and actively control gas slip.

A review on slip models for gas microflows

https://www.pure.ed.ac.uk/ws/files/17079048/ScanlonEtAlCandF2010.pdf

An open source, parallel DSMC code for rarefied gas flows in arbitrary geometries

The kinetic-theory-based numerical schemes, such as direct simulation Monte Carlo (DSMC) and information preservation (IP), can be readily used to solve transition flow regimes. However, their high computational cost still promotes the researchers to extend the Navier–Stokes (NS) equations beyond the slip flow and to the transition regime applications. Evidently, a suitable extension would accurately predict both the local velocity profiles and the mass flow rate magnitude as well as the streamwise pressure distribution. The second-order slip velocity model derived from kinetic theory can provide relatively accurate velocity profiles up to a Knudsen (Kn) number of around 0.5; however, its mass flow rate accuracy decreases as Knudsen number approaches the upper bound. One remedy is to consider the rarefaction effects in calculating the NS viscosity coefficient. In this work, we use the shear stress distribution derived from our IP simulations, extend an analytical expression for the viscosity coefficient, ...

/pdf/extending-the-navier-stokes-solutions-to-transition-regime-rmzny92g1t.pdf

Extending the Navier–Stokes solutions to transition regime in two-dimensional micro- and nanochannel flows using information preservation scheme

We use an extended direct simulation Monte Carlo (DSMC) method, applicable to unstructured meshes, to numerically simulate a wide range of rarefaction regimes from subsonic to supersonic flows through micro/nanoscale converging–diverging nozzles. Our unstructured DSMC method considers a uniform distribution of particles, employs proper subcell geometry, and follows an appropriate particle tracking algorithm. Using the unstructured DSMC, we study the effects of back pressure, gas/surface interactions (diffuse/specular reflections), and Knudsen number on the flow field in micro/nanoscale nozzles. If we apply the back pressure at the nozzle outlet, a boundary layer separation occurs before the outlet and a region with reverse flow appears inside the boundary layer. Meanwhile, the core region of inviscid flow experiences multiple shock-expansion waves. In order to accurately simulate the outflow, we extend a buffer zone at the nozzle outlet. We show that a high viscous force creation in the wall boundary layer prevents any supersonic flow formation in the divergent part of the nozzle if the Knudsen number exceeds a moderate magnitude. We also show that the wall boundary layer prevents forming any normal shock in the divergent part. In reality, Mach cores would appear at the nozzle center followed by bow shocks and expansion region. We compare the current DSMC results with the solution of the Navier–Stokes equations subject to the velocity slip and temperature jump boundary conditions. We use OpenFOAM as a compressible flow solver to treat the Navier–Stokes equations.

Study of subsonic–supersonic gas flow through micro/nanoscale nozzles using unstructured DSMC solver

We use a direct simulation Monte Carlo (DSMC) method to simulate gas heating/cooling and choked subsonic flows in micro/nanoscale channels subject to either constant wall temperature or constant/variable heat flux boundary conditions. We show the effects of applying various boundary conditions on the mass flow rate and the flow parameters. We also show that it is necessary to add a buffer zone at the end of the channel if we wish to simulate more realistic conditions at the channel outlet. We also discuss why applying equilibrium-based Maxwellian distribution on molecules coming from the channel outlet, where the flow is nonequilibrium, will not disturb the DSMC solution. The current velocity, pressure, and mass flow rate results are compared with different analytical solutions of the Navier―Stokes equations. Although there are good agreements between the DSMC results and the analytical solutions in low compressible flow, the analytical solutions yield incorrect velocity and mass flow rate values in short micro/nanochannel flows with high compressibility and/or choked flow conditions.

Direct Simulation Monte Carlo Solution of Subsonic Flow Through Micro/Nanoscale Channels

The main objective of this work is to develop a novel moving-mesh finite-volume method capable of solving the seepage problem in domains with arbitrary geometries. One major difficulty in analysing the seepage problem is the position of phreatic boundary which is unknown at the beginning of solution. In the current algorithm, we first choose an arbitrary solution domain with a hypothetical phreatic boundary and distribute the finite volumes therein. Then, we derive the conservative statement on a curvilinear co-ordinate system for each cell and implement the known boundary conditions all over the solution domain. Defining a consistency factor, the inconsistency between the hypothesis boundary and the known boundary conditions is measured at the phreatic boundary. Subsequently, the preceding mesh is suitably deformed so that its upper boundary matches the new location of the phreatic surface. This tactic results in a moving-mesh procedure which is continued until the nonlinear boundary conditions are fully satisfied at the phreatic boundary. To validate the developed algorithm, a number of seepage models, which have been previously targeted by the other investigators, are solved. Comparisons between the current results and those of other numerical methods as well as the experimental data show that the current moving-grid finite-volume method is highly robust and it provides sufficient accuracy and reliability. Copyright © 2007 John Wiley & Sons, Ltd.

Masoud Darbandi

Papers

An open source, parallel DSMC code for rarefied gas flows in arbitrary geometries

Extending the Navier–Stokes solutions to transition regime in two-dimensional micro- and nanochannel flows using information preservation scheme

Study of subsonic–supersonic gas flow through micro/nanoscale nozzles using unstructured DSMC solver

Direct Simulation Monte Carlo Solution of Subsonic Flow Through Micro/Nanoscale Channels

A moving-mesh finite-volume method to solve free-surface seepage problem in arbitrary geometries