Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStable multispeed lattice Boltzmann methods
110
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We demonstrate how to produce a stable multispeed lattice Boltzmann method (LBM) for a wide range of velocity sets, many of which were previously thought to be intrinsically unstable. We use non-Gauss--Hermitian cubatures. The method operates stably for almost zero viscosity, has second-order accuracy, suppresses typical spurious oscillation (only a modest Gibbs effect is present) and introduces no artificial viscosity. There is almost no computational cost for this innovation. DISCLAIMER: Additional tests and wide discussion of this preprint show that the claimed property of coupled steps: no artificial dissipation and the second-order accuracy of the method are valid only on sufficiently fine grids. For coarse grids the higher-order terms destroy coupling of steps and additional dissipation appears. The equations are true.
rate research
Read More
We construct a system of nonequilibrium entropy limiters for the lattice Boltzmann methods (LBM). These limiters erase spurious oscillations without blurring of shocks, and do not affect smooth solutions. In general, they do the same work for LBM as
flux limiters do for finite differences, finite volumes and finite elements methods, but for LBM the main idea behind the construction of nonequilibrium entropy limiter schemes is to transform a field of a scalar quantity - nonequilibrium entropy. There are two families of limiters: (i) based on restriction of nonequilibrium entropy (entropy trimming) and (ii) based on filtering of nonequilibrium entropy (entropy filtering). The physical properties of LBM provide some additional benefits: the control of entropy production and accurate estimate of introduced artificial dissipation are possible. The constructed limiters are tested on classical numerical examples: 1D athermal shock tubes with an initial density ratio 1:2 and the 2D lid-driven cavity for Reynolds numbers Re between 2000 and 7500 on a coarse 100*100 grid. All limiter constructions are applicable for both entropic and non-entropic quasiequilibria.
Current implementations of fluctuating lattice Boltzmann equations (FLBE) describe single component fluids. In this paper, a model based on the continuum kinetic Boltzmann equation for describing multicomponent fluids is extended to incorporate the effects of thermal fluctuations. The thus obtained fluctuating Boltzmann equation is first linearized to apply the theory of linear fluctuations, and expressions for the noise covariances are determined by invoking the fluctuation-dissipation theorem (FDT) directly at the kinetic level. Crucial for our analysis is the projection of the Boltzmann equation onto the ortho-normal Hermite basis. By integrating in space and time the fluctuating Boltzmann equation with a discrete number of velocities, the FLBE is obtained for both ideal and non-ideal multicomponent fluids. Numerical simulations are specialized to the case where mean-field interactions are introduced on the lattice, indicating a proper thermalization of the system.
We present an energy-conserving multiple-relaxation-time finite difference lattice Boltzmann model for compressible flows. This model is based on a 16-discrete-velocity model. The collision step is first calculated in the moment space and then mapped back to the velocity space. The moment space and corresponding transformation matrix are constructed according to the group representation theory. Equilibria of the nonconserved moments are chosen according to the need of recovering compressible Navier-Stokes equations through the Chapman-Enskog expansion. Numerical experiments showed that compressible flows with strong shocks can be well simulated by the present model. The used benchmark tests include (i) shock tubes, such as the Sod, Lax, Sjogreen, Colella explosion wave and collision of two strong shocks, (ii) regular and Mach shock reflections, and (iii) shock wave reaction on cylindrical bubble problems. The new model works for both low and high speeds compressible flows. It contains more physical information and has better numerical stability and accuracy than its single-relaxation-time version.
The lattice-Boltzmann method (LBM) and its variants have emerged as promising, computationally efficient and increasingly popular numerical methods for modelling complex fluid flow. However, it is acknowledged that the method can demonstrate numerical instabilities, e.g., in the vicinity of shocks. We propose a simple and novel technique to stabilise the lattice-Boltzmann method by monitoring the difference between microscopic and macroscopic entropy. Populations are returned to their equilibrium states if a threshold value is exceeded. We coin the name Ehrenfests steps for this procedure in homage to the vehicle that we use to introduce the procedure, namely, the Ehrenfests idea of coarse-graining. The one-dimensional shock tube for a compressible isothermal fluid is a standard benchmark test for hydrodynamic codes. We observe that, of all the LBMs considered in the numerical experiment with the one-dimensional shock tube, only the method which includes Ehrenfests steps is capable of suppressing spurious post-shock oscillations.
We present a systematic account of recent developments of the relativistic Lattice Boltzmann method (RLBM) for dissipative hydrodynamics. We describe in full detail a unified, compact and dimension-independent procedure to design relativistic LB schemes capable of bridging the gap between the ultra-relativistic regime, $k_{rm B} T gg mc^2$, and the non-relativistic one, $k_{rm B} T ll mc^2$. We further develop a systematic derivation of the transport coefficients as a function of the kinetic relaxation time in $d=1,2,3$ spatial dimensions. The latter step allows to establish a quantitative bridge between the parameters of the kinetic model and the macroscopic transport coefficients. This leads to accurate calibrations of simulation parameters and is also relevant at the theoretical level, as it provides neat numerical evidence of the correctness of the Chapman-Enskog procedure. We present an extended set of validation tests, in which simulation results based on the RLBMs are compared with existing analytic or semi-analytic results in the mildly-relativistic ($k_{rm B} T sim mc^2$) regime for the case of shock propagations in quark-gluon plasmas and laminar electronic flows in ultra-clean graphene samples. It is hoped and expected that the material collected in this paper may allow the interested readers to reproduce the present results and generate new applications of the RLBM scheme.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
