Reduced-Order Modeling of Low Mach Number Unsteady Microchannel Flows

2014 ◽  
Vol 136 (5) ◽  
Author(s):  
Leila Issa ◽  
Issam Lakkis
Keyword(s):  
Mach Number ◽  
Boundary Conditions ◽  
Stokes Equations ◽  
Solid Boundary ◽  
Two Dimensional ◽  
Slip Boundary Conditions ◽  
Reduced Order ◽  
Low Mach Number ◽  
First Order ◽  
Slip Boundary

We present reduced-order models of unsteady low-Mach-number ideal gas flows in two-dimensional rectangular microchannels subject to first-order slip-boundary conditions. The pressure and density are related by a polytropic process, allowing for isothermal or isentropic flow assumptions. The Navier–Stokes equations are simplified using low-Mach-number expansions of the pressure and velocity fields. Up to first order, this approximation results in a system that is subject to no-slip condition at the solid boundary. The second-order system satisfies the slip-boundary conditions. The resulting equations and the subsequent pressure-flow-rate relationships enable modeling the flow using analog circuit components. The accuracy of the proposed models is investigated for steady and unsteady flows in a two-dimensional channel for different values of Mach and Knudsen numbers.

Transition and turbulence in horizontal convection: linear stability analysis

2017 ◽  
Vol 821 ◽  
pp. 31-58 ◽  
Author(s):  
Pierre-Yves Passaggia ◽  
Alberto Scotti ◽  
Brian White
Keyword(s):  
Stability Analysis ◽  
Rayleigh Number ◽  
Boundary Conditions ◽  
Stokes Equations ◽  
Critical Rayleigh Number ◽  
Base Flow ◽  
Two Dimensional ◽  
Slip Boundary Conditions ◽  
Slip Boundary ◽  
Horizontal Convection

The linear instability mechanisms of horizontal convection in a rectangular cavity forced by a horizontal buoyancy gradient along its surface are investigated using local and global stability analyses for a Prandtl number equal to unity. The results show that the stability of the base flow, a steady circulation characterized by a narrow descending plume and a broad upwelling region, depends on the Rayleigh number, $Ra$. For free-slip boundary conditions at a critical value of $Ra\approx 2\times 10^{7}$, the steady base flow becomes unstable to three-dimensional perturbations, characterized by counter-rotating vortices originating within the plume region. A Wentzel–Kramers–Brillouin (WKB) method applied along closed streamlines demonstrates that this instability is of a Rayleigh–Taylor type and can be used to accurately reconstruct the global instability mode. In the case of no-slip boundary conditions, the base flow also becomes unstable to a self-sustained two-dimensional instability whose critical Rayleigh number is $Ra\approx 1.7\times 10^{8}$. Beyond this critical $Ra$, two-dimensional equilibrium stationary states of the Navier–Stokes equations are computed using the selective frequency damping method. The two-dimensional onset of instability is shown to be characterized by a family of modes also originating within the plume. A local spatio-temporal stability analysis shows that the flow becomes absolutely unstable at the origin of the plume. Taken together, these results illustrate the mechanisms behind the onset of turbulence that has been observed in horizontal convection.

Numerical Investigation of Combined Effects of Rarefaction and Compressibility for Gas Flow in Microchannels and Microtubes

2009 ◽  
Vol 131 (10) ◽  
Author(s):  
Xiaohong Yan ◽  
Qiuwang Wang
Keyword(s):  
Boundary Conditions ◽  
Friction Factor ◽  
Knudsen Number ◽  
Velocity Gradient ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Slip Boundary Conditions ◽  
Navier Stokes Equations ◽  
First Order ◽  
Slip Boundary

In this paper, first, the Navier–Stokes equations for incompressible fully developed flow in microchannels and microtubes with the first-order and second-order slip boundary conditions are analytically solved. Then, the compressible Navier–Stokes equations are numerically solved with slip boundary conditions. The numerical methodology is based on the control volume scheme. Numerical results reveal that the compressibility effect increases the velocity gradient near the wall and the friction factor. On the other hand, the increment of velocity gradient near the wall leads to a much larger slip velocity than that for incompressible flow with the same value of Knudsen number and results in a corresponding decrement of friction factor. General correlations for the Poiseuille number (fRe), the Knudsen number (Kn), and the Mach number (Ma) containing the first-order and second-order slip coefficients are proposed. Correlations are validated with available experimental and numerical results.

scholarly journals Bifurcations in a quasi-two-dimensional Kolmogorov-like flow

2017 ◽  
Vol 828 ◽  
pp. 837-866 ◽  
Author(s):  
Jeffrey Tithof ◽  
Balachandra Suri ◽  
Ravi Kumar Pallantla ◽  
Roman O. Grigoriev ◽  
Michael F. Schatz
Keyword(s):  
Boundary Conditions ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Quantitative Agreement ◽  
Navier Stokes ◽  
Dielectric Fluid ◽  
Two Dimensional ◽  
Slip Boundary Conditions ◽  
Slip Boundary ◽  
2D Model

We present a combined experimental and theoretical study of the primary and secondary instabilities in a Kolmogorov-like flow. The experiment uses electromagnetic forcing with an approximately sinusoidal spatial profile to drive a quasi-two-dimensional (Q2D) shear flow in a thin layer of electrolyte suspended on a thin lubricating layer of a dielectric fluid. Theoretical analysis is based on a two-dimensional (2D) model (Suri et al., Phys. Fluids, vol. 26 (5), 2014, 053601), derived from first principles by depth-averaging the full three-dimensional Navier–Stokes equations. As the strength of the forcing is increased, the Q2D flow in the experiment undergoes a series of bifurcations, which is compared with results from direct numerical simulations of the 2D model. The effects of confinement and the forcing profile are studied by performing simulations that assume spatial periodicity and strictly sinusoidal forcing, as well as simulations with realistic no-slip boundary conditions and an experimentally validated forcing profile. We find that only the simulation subject to physical no-slip boundary conditions and a realistic forcing profile provides close, quantitative agreement with the experiment. Our analysis offers additional validation of the 2D model as well as a demonstration of the importance of properly modelling the forcing and boundary conditions.

Reduced Order Models of Low Mach Number Isothermal Flows in Microchannels

2013 ◽  
Author(s):  
Leila Issa ◽  
Issam Lakkis
Keyword(s):  
Mach Number ◽  
Analog Circuit ◽  
Stokes Equations ◽  
Momentum Equation ◽  
Ideal Gas ◽  
Reduced Order Models ◽  
Two Dimensional ◽  
Reduced Order ◽  
Low Mach Number ◽  
Circuit Components

We present reduced order models of unsteady low Mach number isothermal ideal gas flows in two-dimensional rectangular microchannels subject to first order slip boundary conditions. The Navier-Stokes equations are simplified using Low Mach Number expansions of the pressure and velocity fields. This approximation allows decoupling the density from spatial pressure variations, thus simplifying the momentum equation. The resulting diffusion equation and the subsequent pressure-flow-rate relationship enables modeling the flow using analog circuit components. The accuracy of the proposed models is investigated for steady and unsteady flows in a two-dimensional channel for different values of Reynolds and Knudsen numbers.

Hydrodynamic interaction of two unequal-sized spheres in a slightly rarefied gas: resistance and mobility functions

1989 ◽  
Vol 207 ◽  
pp. 353-378 ◽  
Author(s):  
Ruoxian Ying ◽  
Michael H. Peters
Keyword(s):  
Boundary Conditions ◽  
Hydrodynamic Interaction ◽  
Stokes Equations ◽  
Rarefied Gas ◽  
Zero Order ◽  
Knudsen Layer ◽  
Slip Boundary Conditions ◽  
First Order ◽  
Slip Boundary ◽  
Slightly Rarefied Gas

The problem of the hydrodynamic interaction of two unequal-sized spheres in a slightly rarefied gas is treated following the singular perturbation scheme of Sone & Onishi (1978), valid at small, but finite, particle Knudsen numbers. In this method the solution to the linearized BGKW transport equation governing the gas molecular motion consists of two parts: one describing a Knudsen layer where the actual microscopic boundary conditions are applied and the other describing a Hilbert region where the Stokes equations of continuum hydrodynamics hold. The Knudsen-layer solution establishes the ‘slip’ boundary conditions for the Stokes equations. Here we clearly distinguish between particle ‘slip’ due to the type of boundary conditions and particle ‘slip’ due to lengthscale effects as measured by the Knudsen number. The present analysis has been carried out to first order in particle Knudsen number for the case of diffuse reflective molecular boundary conditions. General relationships between the first- and zero-order velocity fields, both of which are written in the form of Lamb's (1932) solution to the Stokes equation, are established. It is illustrated how these general relationships can be used to determine the force and torque acting on a single sphere translating and rotating in a slightly rarefied gas. Finally, we have treated the two-sphere problem in a slightly rarefied gas using the twin multipole expansion method of Jeffrey & Onishi (1984). Here again, general relationships are established between the solutions of the first-order fluid velocity field and the zero-order velocity field, the latter being shown to recover Jeffrey & Onishi's results for stick boundary conditions. These general relationships are subsequently used to determine the complete resistance and mobility matrices of the two-sphere system. The symmetric properties of the resistance and mobility matrices are demonstrated for slip boundary conditions, in agreement with the general proof of Landau & Lifshitz (1980) and Bedeaux, Albano & Mazur (1977).

Parallel iterative finite-element algorithms for the Navier–Stokes equations with nonlinear slip boundary conditions

2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Kangrui Zhou ◽  
Yueqiang Shang
Keyword(s):  
Finite Element ◽  
Boundary Conditions ◽  
Parallel Algorithms ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Slip Boundary Conditions ◽  
Navier Stokes Equations ◽  
Slip Boundary ◽  
Parallel Iterative ◽  
Nonlinear Slip Boundary Conditions

AbstractBased on full domain partition, three parallel iterative finite-element algorithms are proposed and analyzed for the Navier–Stokes equations with nonlinear slip boundary conditions. Since the nonlinear slip boundary conditions include the subdifferential property, the variational formulation of these equations is variational inequalities of the second kind. In these parallel algorithms, each subproblem is defined on a global composite mesh that is fine with size h on its subdomain and coarse with size H (H ≫ h) far away from the subdomain, and then we can solve it in parallel with other subproblems by using an existing sequential solver without extensive recoding. All of the subproblems are nonlinear and are independently solved by three kinds of iterative methods. Compared with the corresponding serial iterative finite-element algorithms, the parallel algorithms proposed in this paper can yield an approximate solution with a comparable accuracy and a substantial decrease in computational time. Contributions of this paper are as follows: (1) new parallel algorithms based on full domain partition are proposed for the Navier–Stokes equations with nonlinear slip boundary conditions; (2) nonlinear iterative methods are studied in the parallel algorithms; (3) new theoretical results about the stability, convergence and error estimates of the developed algorithms are obtained; (4) some numerical results are given to illustrate the promise of the developed algorithms.

Sign in / Sign up

Export Citation Format

Share Document