Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCentral Moment Lattice Boltzmann Method on a Rectangular Lattice
148
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Simulating inhomogeneous flows with different characteristic scales in different coordinate directions using the collide-and-stream based lattice Boltzmann methods (LBM) can be accomplished efficiently using rectangular lattice grids. We develop and investigate a new rectangular central moment LBM based on non-orthogonal moment basis (referred to as RC-LBM). The equilibria to which the central moments relax under collision in this approach are obtained from matching with those corresponding to the continuous Maxwell distribution. A Chapman-Enskog analysis is performed to derive the correction terms to the second order moment equilibria involving the grid aspect ratio and velocity gradients that restores the isotropy of the viscous stress tensor and eliminates the non-Galilean invariant cubic velocity terms of the resulting hydrodynamical equations. A special case of this rectangular formulation involving the raw moments (referred to as the RNR-LBM) is also constructed. The resulting schemes represent a considerable simplification, especially for the transformation matrices and isotropy corrections, and improvement over the existing MRT-LB schemes on rectangular lattice grids that use orthogonal moment basis. Numerical validation study of both the RC-LBM and RNR-LBM for a variety of benchmark flow problems are performed that show good accuracy at various grid aspect ratios. The ability of our proposed schemes to simulate flows using relatively lower grid aspect ratios than considered in prior rectangular LB approaches is demonstrated. Furthermore, simulations reveal the superior stability characteristics of the RC-LBM over RNR-LBM in handling shear flows at lower viscosities and/or higher characteristic velocities. In addition, computational advantages of using our rectangular LB formulation in lieu of that based on the square lattice is shown.
rate research
Read More
Lattice Boltzmann (LB) models used for the computation of fluid flows represented by the Navier-Stokes (NS) equations on standard lattices can lead to non-Galilean invariant (GI) viscous stress involving cubic velocity errors. This arises from the dependence of their third order diagonal moments on the first order moments for standard lattices, and strategies have recently been introduced to restore GI without such errors using a modified collision operator involving either corrections to the relaxation times or to the moment equilibria. Convergence acceleration in the simulation of steady flows can be achieved by solving the preconditioned NS equations, which contain a preconditioning parameter that alleviates the numerical stiffness. In the present study, we present a GI formulation of the preconditioned cascaded central moment LB method used to solve the preconditioned NS equations, which is free of cubic velocity errors on a standard lattice. A Chapman-Enskog analysis reveals the structure of the spurious non-GI defect terms and it is demonstrated that the anisotropy of the resulting viscous stress is dependent on the preconditioning parameter, in addition to the fluid velocity. It is shown that partial correction to eliminate the cubic velocity defects is achieved by scaling the cubic velocity terms in the off-diagonal third-order moment equilibria with the square of the preconditioning parameter. Furthermore, we develop additional corrections based on the extended moment equilibria involving gradient terms with coefficients dependent locally on the fluid velocity and the preconditioning parameter. Several conclusions are drawn from the analysis of the structure of the non-GI errors and the associated corrections, with particular emphasis on their dependence on the preconditioning parameter. Improvements in accuracy and convergence acceleration are demonstrated.
We present a multi-scale lattice Boltzmann scheme, which adaptively refines particles velocity space. Different velocity sets, i.e., higher- and lower-order lattices, are consistently and efficiently coupled, allowing us to use the higher-order lattice only when and where needed. This includes regions of either high Mach number or high Knudsen number. The coupling procedure of different lattices consists of either projection of the moments of the higher-order lattice onto the lower-order lattice or lifting of the lower-order lattice to the higher-order velocity space. Both lifting and projection are local operations, which enable a flexible adaptive velocity set. The proposed scheme can be formulated both in a static and an optimal, co-moving reference frame, in the spirit of the recently introduced Particles on Demand method. The multi-scale scheme is first validated through a convected athermal vortex and also studied in a jet flow setup. The performance of the proposed scheme is further investigated through the shock structure problem and a high Knudsen Couette flow, typical examples of highly non-equilibrium flows in which the order of the velocity set plays a decisive role. The results demonstrate that the proposed multi-scale scheme can operate accurately, with flexibility in terms of the underlying models and with reduced computational requirements.
It is well-known that the original lattice Boltzmann (LB) equation deviates from the Navier-Stokes equations due to an unphysical velocity dependent viscosity. This unphysical dependency violates the Galilean invariance and limits the validation domain of the LB method to near incompressible flows. As previously shown, recovery of correct transport phenomena in kinetic equations depends on the higher hydrodynamic moments. In this Letter, we give specific criteria for recovery of various transport coefficients. The Galilean invariance of a general class of LB models is demonstrated via numerical experiments.
When boiling occurs in a liquid flow field, the phenomenon is known as forced-convection boiling. We numerically investigate such a boiling system on a cylinder in a flow at a saturated condition. To deal with the complicated liquid-vapor phase-change phenomenon, we develop a numerical scheme based on the pseudopotential lattice Boltzmann method (LBM). The collision stage is performed in the space of central moments (CMs) to enhance numerical stability for high Reynolds numbers. The adopted forcing scheme, consistent with the CMs-based LBM, leads to a concise yet robust algorithm. Furthermore, additional terms required to ensure thermodynamic consistency are derived in a CMs framework. The effectiveness of the present scheme is successfully tested against a series of boiling processes, including nucleation, growth, and departure of a vapor bubble for Reynolds numbers varying between 30 and 30000. Our CMs-based LBM can reproduce all the boiling regimes, i.e., nucleate boiling, transition boiling, and film boiling, without any artificial input such as initial vapor phase. We find that the typical boiling curve, also known as the Nukiyama curve, appears even though the focused system is not the pool boiling but the forced-convection system. Also, our simulations support experimental observations of intermittent direct solid-liquid contact even in the film-boiling regime. Finally, we provide quantitative comparison with the semi-empirical correlations for the forced-convection film boiling on a cylinder on the Nu-Ja diagram.
A new lattice Boltzmann model for multicomponent ideal gas mixtures is presented. The model development consists of two parts. First, a new kinetic model for Stefan- Maxwell diffusion amongst the species is proposed and realized as a lattice Boltzmann equation on the standard discrete velocity set. Second, a compressible lattice Boltzmann model for the momentum and energy of the mixture is established. Both parts are consistently coupled through mixture composition, momentum, pressure, energy and enthalpy whereby a passive scalar advection-diffusion coupling is obviated, unlike in previous approaches. The proposed model is realized on the standard three-dimensional lattices and is validated with a set of benchmarks highlighting various physical aspects of compressible mixtures. Stefan-Maxwell diffusion is tested against experiment and theory of uphill diffusion of argon and methane in a ternary mixture with hydrogen. The speed of sound is measured in various binary and ternary compositions. We further validate the Stefan-Maxwell diffusion coupling with hydrodynamics by simulating diffusion in opposed jets and the three-dimensional Kelvin-Helmholtz instability of shear layers in a two-component mixture. Apart from the multicomponent compressible mixture, the proposed lattice Boltzmann model also provides an extension of the lattice Boltzmann equation to the compressible flow regime on the standard three-dimensional lattice.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
