Relationship between accuracy and number of velocity particles of the finite-difference lattice Boltzmann method in velocity slip simulations

2010 ◽  
Vol 132 (10) ◽  
Author(s):  
Minoru Watari
Keyword(s):  
Relaxation Process ◽  
Stokes Equations ◽  
Velocity Slip ◽  
Simulation Models ◽  
Knudsen Layer ◽  
Navier Stokes ◽  
Flow Area ◽  
Navier Stokes Equations ◽  
Conservation Of Momentum ◽  
Boltzmann Method

Relationship between accuracy and number of velocity particles in velocity slip phenomena was investigated by numerical simulations and theoretical considerations. Two types of 2D models were used: the octagon family and the D2Q9 model. Models have to possess the following four prerequisites to accurately simulate the velocity slip phenomena: (a) equivalency to the Navier–Stokes equations in the N-S flow area, (b) conservation of momentum flow Pxy in the whole area, (c) appropriate relaxation process in the Knudsen layer, and (d) capability to properly express the mass and momentum flows on the wall. Both the octagon family and the D2Q9 model satisfy conditions (a) and (b). However, models with fewer velocity particles do not sufficiently satisfy conditions (c) and (d). The D2Q9 model fails to represent a relaxation process in the Knudsen layer and shows a considerable fluctuation in the velocity slip due to the model’s angle to the wall. To perform an accurate velocity slip simulation, models with sufficient velocity particles, such as the triple octagon model with moving particles of 24 directions, are desirable.

scholarly journals Involving the Navier–Stokes equations in the derivation of boundary conditions for the lattice Boltzmann method

2011 ◽  
Vol 369 (1944) ◽  
pp. 2184-2192 ◽  
Author(s):  
Joris C. G. Verschaeve
Keyword(s):  
Boundary Condition ◽  
Lattice Boltzmann Method ◽  
Lattice Boltzmann ◽  
Stokes Equations ◽  
Slip Boundary Condition ◽  
Navier Stokes ◽  
Grid Spacing ◽  
Navier Stokes Equations ◽  
Slip Boundary ◽  
Boltzmann Method

By means of the continuity equation of the incompressible Navier–Stokes equations, additional physical arguments for the derivation of a formulation of the no-slip boundary condition for the lattice Boltzmann method for straight walls at rest are obtained. This leads to a boundary condition that is second-order accurate with respect to the grid spacing and conserves mass. In addition, the boundary condition is stable for relaxation frequencies close to two.

scholarly journals Plane poloidal-toroidal decomposition of doubly periodic vector fields. Part 2. The Stokes equations

2005 ◽  
Vol 47 (1) ◽  
pp. 39-50 ◽  
Author(s):  
G. D. McBain
Keyword(s):  
Vector Fields ◽  
Stokes Equations ◽  
Momentum Equation ◽  
Integral Representations ◽  
Navier Stokes ◽  
Field Equations ◽  
Navier Stokes Equations ◽  
Time Discretisation ◽  
Mass Constraint ◽  
Conservation Of Momentum

AbstractWe continue our study of the adaptation from spherical to doubly periodic slot domains of the poloidal-toroidal representation of vector fields. Building on the successful construction of an orthogonal quinquepartite decomposition of doubly periodic vector fields of arbitrary divergence with integral representations for the projections of known vector fields and equivalent scalar representations for unknown vector fields (Part 1), we now present a decomposition of vector field equations into an equivalent set of scalar field equations. The Stokes equations for slow viscous incompressible fluid flow in an arbitrary force field are treated as an example, and for them the application of the decomposition uncouples the conservation of momentum equation from the conservation of mass constraint. The resulting scalar equations are then solved by elementary methods. The extension to generalised Stokes equations resulting from the application of various time discretisation schemes to the Navier-Stokes equations is also solved.

Influence of Inertia Terms on High Pressure Gap Flow Applications in Hydraulics

2019 ◽  
Author(s):  
Felix Fischer ◽  
Andreas Rhein ◽  
Katharina Schmitz
Keyword(s):  
High Pressure ◽  
Reynolds Equation ◽  
Stokes Equations ◽  
Fluid Phase ◽  
Simulation Models ◽  
Rigid Bodies ◽  
Navier Stokes ◽  
Dependent Manner ◽  
Fluid Inertia ◽  
Navier Stokes Equations

Abstract Hydraulic pumps, which reach pressures up to 3000 bar, are often realized as plunger-piston type pumps. In the case of a common-rail pump for diesel injection systems, the plunger is driven by a cam-tappet construction and the contact during suction stroke is maintained by a helical spring. Many hydraulic piston-based high pressure pumps include gap seals, which are formed by small clearances between the two surfaces of the piston and the bushing. Usually the gap height is in the magnitude of several micrometers. Typical radial gaps are between 0.5 and 1 per mil of the nominal diameter. These gap seals are used to allow and maintain pressure build up in the piston chamber. When the gap is pressurized, a special flow regime is reached. For the description of this particular flow the Reynolds equation, which is a simplification of the Navier-Stokes equations, can be used as done in the state of the art. Furthermore, if the pressure in the gap is high enough — 500 bar and above — fluid-structure interactions must be taken into account. Pressure levels above 1500 or 2000 bar indicate the necessity for solving the energy equation of the fluid phase and the rigid bodies surrounding it. In any case, the fluid properties such as density and viscosity, have to be modelled in a pressure dependent manner. This means, a compressible flow is described in the sealing gap. Viscosity changes in magnitudes while density remains in the same magnitude, but nevertheless changes about 30 %. These facts must be taken into account when solving the Reynolds equation. In this paper the authors work out that the Reynolds equation is not suitable for every piston-bushing gap seal in hydraulic applications. It will be shown that remarkable errors are made, when the inertia terms in the Navier-Stokes equations are neglected, especially in high pressure applications. To work out the influence of the inertia terms in these flows, two simulation models are built up and calculated for the physical problem. One calculates the compressible Reynolds equation neglecting the fluid inertia. The other model, taking the fluid inertia into account, calculates the coupled Navier-Stokes equations on the same geometrical boundaries. Here, the so called SIMPLE (Semi-Implicit Method for Pressure Linked Equations) algorithm is used. The discretization is realized with the Finite Volume Method. Afterwards, the solutions of both models are compared to investigate the influence of the inertia terms on the flow in these specific high pressure applications.

A New Partial Slip Boundary Condition for the Lattice-Boltzmann Method

2013 ◽  
Author(s):  
Marc-Florian Uth ◽  
Alf Crüger ◽  
Heinz Herwig
Keyword(s):  
Boundary Condition ◽  
Lattice Boltzmann ◽  
Stokes Equations ◽  
Slip Boundary Condition ◽  
Navier Stokes ◽  
Slip Model ◽  
Navier Stokes Equations ◽  
Navier Slip ◽  
Slip Boundary ◽  
Boltzmann Method

In micro or nano flows a slip boundary condition is often needed to account for the special flow situation that occurs at this level of refinement. A common model used in the Finite Volume Method (FVM) is the Navier-Slip model which is based on the velocity gradient at the wall. It can be implemented very easily for a Navier-Stokes (NS) Solver. Instead of directly solving the Navier-Stokes equations, the Lattice-Boltzmann method (LBM) models the fluid on a particle basis. It models the streaming and interaction of particles statistically. The pressure and the velocity can be calculated at every time step from the current particle distribution functions. The resulting fields are solutions of the Navier-Stokes equations. Boundary conditions in LBM always not only have to define values for the macroscopic variables but also for the particle distribution function. Therefore a slip model cannot be implemented in the same way as in a FVM-NS solver. An additional problem is the structure of the grid. Curved boundaries or boundaries that are non-parallel to the grid have to be approximated by a stair-like step profile. While this is no problem for no-slip boundaries, any other velocity boundary condition such as a slip condition is difficult to implement. In this paper we will present two different implementations of slip boundary conditions for the Lattice-Boltzmann approach. One will be an implementation that takes advantage of the microscopic nature of the method as it works on a particle basis. The other one is based on the Navier-Slip model. We will compare their applicability for different amounts of slip and different shapes of walls relative to the numerical grid. We will also show what limits the slip rate and give an outlook of how this can be avoided.

scholarly journals Turbulence Modeling a Review for Different Used Methods

2021 ◽  
Vol 39 (1) ◽  
pp. 227-234
Author(s):  
Khelifa Hami
Keyword(s):  
Turbulence Modeling ◽  
Stokes Equations ◽  
Turbulence Models ◽  
Simulation Models ◽  
Navier Stokes ◽  
Future Research ◽  
Detached Eddy Simulation ◽  
Navier Stokes Equations ◽  
Eddy Simulation ◽  
Complicated Geometries

This contribution represents a critical view of the advantages and limits of the set of mathematical models of the physical phenomena of turbulence. Turbulence models can be grouped into two categories, depending on how turbulent quantities are calculated: direct numerical simulations (DNS) and RANS (Reynolds Averaged Navier-Stokes Equations) models. The disadvantage of these models is that they require enormous computing power, inaccessible, especially for large and complicated geometries. For this reason, hybrid models (combinations between DNS and RANS methods) have been developed, for example, the LES (“Large Eddy Simulation”) or DES (“Detached Eddy Simulation”) models. They represent a compromise - are less precise than DNS, but more precise than RANS models. The results presented in this contribution will allow and facilitate future research in the field the choice of the model approach necessary for the case studies whatever their difficulty factor.

Multifidelity modeling similarity conditions for airfoil dynamic stall prediction with manifold mapping

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Vishal Raul ◽  
Leifur Leifsson
Keyword(s):  
Prediction Error ◽  
Stokes Equations ◽  
Trust Region ◽  
Simulation Models ◽  
Time Discretization ◽  
Dynamic Stall ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Content Type ◽  
Manifold Mapping

PurposeThe purpose of this work is to investigate the similarity requirements for the application of multifidelity modeling (MFM) for the prediction of airfoil dynamic stall using computational fluid dynamics (CFD) simulations.Design/methodology/approachDynamic stall is modeled using the unsteady Reynolds-averaged Navier–Stokes equations and Menter's shear stress transport turbulence model. Multifidelity models are created by varying the spatial and temporal discretizations. The effectiveness of the MFM method depends on the similarity between the high- (HF) and low-fidelity (LF) models. Their similarity is tested by computing the prediction error with respect to the HF model evaluations. The proposed approach is demonstrated on three airfoil shapes under deep dynamic stall at a Mach number 0.1 and Reynolds number 135,000.FindingsThe results show that varying the trust-region (TR) radius (λ) significantly affects the prediction accuracy of the MFM. The HF and LF simulation models hold similarity within small (λ ≤ 0.12) to medium (0.12 ≤ λ ≤ 0.23) TR radii producing a prediction error less than 5%, whereas for large TR radii (0.23 ≤ λ ≤ 0.41), the similarity is strongly affected by the time discretization and minimally by the spatial discretization.Originality/valueThe findings of this work present new knowledge for the construction of accurate MFMs for dynamic stall performance prediction using LF model spatial- and temporal discretization setup and the TR radius size. The approach used in this work is general and can be used for other unsteady applications involving CFD-based MFM and optimization.

Lattice Boltzmann Method for the Simulating Extrudate Swell of Viscoelastic Fluid

2015 ◽  
Vol 799-800 ◽  
pp. 784-787
Author(s):  
Wen Qin Liu ◽  
Yong Li
Keyword(s):  
Lattice Boltzmann Method ◽  
Lattice Boltzmann ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Extrudate Swell ◽  
Navier Stokes Equations ◽  
New Approach ◽  
Moving Interface ◽  
Good Accord ◽  
Boltzmann Method

The main objective of this work is to develop a new approach based on the Lattice Boltzmann method (LBM) to simulate the extrudate swell of an Oldroyd B viscoelatic fluid. Two lattice Boltzmann equations are used to solve the Navier-Stokes equations and constitutive equation simultaneously at each time iteration. The single LBM model is used to track the moving interface in this paper. To validate the accuracy and stability of this new scheme, we study the steady 2D Poiseuille flow firstly, finding the numerical results be in good accord with the analytical solution. Then the die-swell phenomenon is solved, we successfully acquire the different swelling state of an Oldroyd B fluid at different time.

scholarly journals Head-on collision of internal waves with trapped cores

2017 ◽  
Vol 24 (4) ◽  
pp. 751-762 ◽  
Author(s):  
Vladimir Maderich ◽  
Kyung Tae Jung ◽  
Kateryna Terletska ◽  
Kyeong Ok Kim
Keyword(s):  
Wave Amplitude ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Class Iii ◽  
Navier Stokes Equations ◽  
Linear Waves ◽  
Non Linear ◽  
Conservation Of Momentum ◽  
Run Up ◽  
The Waves

Abstract. The dynamics and energetics of a head-on collision of internal solitary waves (ISWs) with trapped cores propagating in a thin pycnocline were studied numerically within the framework of the Navier–Stokes equations for a stratified fluid. The peculiarity of this collision is that it involves trapped masses of a fluid. The interaction of ISWs differs for three classes of ISWs: (i) weakly non-linear waves without trapped cores, (ii) stable strongly non-linear waves with trapped cores, and (iii) shear unstable strongly non-linear waves. The wave phase shift of the colliding waves with equal amplitude grows as the amplitudes increase for colliding waves of classes (i) and (ii) and remains almost constant for those of class (iii). The excess of the maximum run-up amplitude, normalized by the amplitude of the waves, over the sum of the amplitudes of the equal colliding waves increases almost linearly with increasing amplitude of the interacting waves belonging to classes (i) and (ii); however, it decreases somewhat for those of class (iii). The colliding waves of class (ii) lose fluid trapped by the wave cores when amplitudes normalized by the thickness of the pycnocline are in the range of approximately between 1 and 1.75. The interacting stable waves of higher amplitude capture cores and carry trapped fluid in opposite directions with little mass loss. The collision of locally shear unstable waves of class (iii) is accompanied by the development of instability. The dependence of loss of energy on the wave amplitude is not monotonic. Initially, the energy loss due to the interaction increases as the wave amplitude increases. Then, the energy losses reach a maximum due to the loss of potential energy of the cores upon collision and then start to decrease. With further amplitude growth, collision is accompanied by the development of instability and an increase in the loss of energy. The collision process is modified for waves of different amplitudes because of the exchange of trapped fluid between colliding waves due to the conservation of momentum.

Orifice Contraction Coefficient for Inviscid Incompressible Flow

1985 ◽  
Vol 107 (1) ◽  
pp. 36-43 ◽  
Author(s):  
R. D. Grose
Keyword(s):  
Stokes Equations ◽  
Inviscid Flow ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Contraction Coefficient ◽  
Two Dimensional Flow ◽  
Conservation Of Momentum ◽  
Orifice Flow ◽  
Contraction Effect ◽  
Control Volumes

The theory for steady flow of an incompressible fluid through an orifice has been semi-empirically established for only certain flow conditions. In this paper, the development of a more rigorous theory for the prediction of the orifice flow contraction effect is presented. This theory is based on the conservation of momentum and mass principles applied to global control volumes for continuum flow. The control volumes are chosen to have a particular geometric construction which is based on certain characteristics of the Navier-Stokes equations for incompressible and, in the limit, inviscid flow. The treatment is restricted to steady incompressible, single phase, single component, inviscid Newtonian flow, but the principles that are developed hold for more general conditions. The resultant equations predict the orifice contraction coefficient as a function of the upstream geometry ratio for both axisymmetric and two-dimensional flow fields. The predicted contraction coefficient values agree with experimental orifice discharge coefficient data without the need for empirical adjustment.

Hydrodynamic Behaviors of the Gas-Lubricated Film in Wedge-Shaped Microchannel

2016 ◽  
Vol 138 (3) ◽  
Author(s):  
Xueqing Zhang ◽  
Qinghua Chen ◽  
Juanfang Liu
Keyword(s):  
Wall Temperature ◽  
Reynolds Equation ◽  
Energy Equation ◽  
Stokes Equations ◽  
Load Capacity ◽  
Velocity Slip ◽  
Navier Stokes ◽  
Hydrodynamic Properties ◽  
Obvious Effect ◽  
Navier Stokes Equations

As for the micro gas bearing operating at a high temperature and speed, one wedge-shaped microchannel is established, and the hydrodynamic properties of the wedge-shaped gas film are comprehensively investigated. The Reynolds equation, modified Reynolds equation, energy equation, and Navier–Stokes equations are employed to describe and analyze the hydrodynamics of the gas film. Furthermore, the comparisons among the hydrodynamic properties predicted by various models were performed for the different wedge factors and the different wall temperatures. The results show that coupling the simplified energy equation with the Reynolds or modified Reynolds equations has an obvious effect on the change of the friction force acting on the horizontal plate and the load capacity of the gas film at the higher wedge factor and the lower wall temperature. The velocity slip weakens the squeeze of the gas film and strengths the gas backflow. A larger wedge factor or a higher wall temperature leads to a higher gas film temperature and thus enhances the rarefaction effect. As the wall temperature is elevated, the load capacity obtained by the Reynolds equation increases, while the results by the Navier–Stokes equations coupled with the full energy equation rapidly decrease. Additionally, the vertical flow across the gas film in the Navier–Stokes equations weakens the squeeze between the gas film and the tilt plate and the gas backflow.

Sign in / Sign up

Export Citation Format

Share Document