Do you want to publish a course? Click here

Fluctuating Multicomponent Lattice Boltzmann Model

163   0   0.0 ( 0 )
 Added by Mauro Sbragaglia Dr
 Publication date 2014
  fields Physics
and research's language is English




Ask ChatGPT about the research

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.



rate research

Read More

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.
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.
We investigate the phase diagram and critical behavior of three-dimensional multicomponent Abelian-Higgs models, in which an N-component complex field z_x^a of unit length and charge is coupled to compact quantum electrodynamics in the usual Wilson lattice formulation. We determine the phase diagram and study the nature of the transition line for N=2 and N=4. Two phases are identified, specified by the behavior of the gauge-invariant local composite operator Q_x^{ab} = bar{z}_x^a z_x^b - delta^{ab}/N, which plays the role of order parameter. In one phase, we have langle Q_x^{ab}rangle =0, while in the other Q_x^{ab} condenses. Gauge correlations are never critical: gauge excitations are massive for any finite coupling. The two phases are separated by a transition line. Our numerical data are consistent with the simple scenario in which the nature of the transition is independent of the gauge coupling. Therefore, for any finite positive value of the gauge coupling, we predict a continuous transition in the Heisenberg universality class for N=2 and a first-order transition for N=4. However, notable crossover phenomena emerge for large gauge couplings, when gauge fluctuations are suppressed. Such crossover phenomena are related to the unstable O(2N) fixed point, describing the behavior of the model in the infinite gauge-coupling limit.
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.
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.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا