3,788 research outputs found
Moment-based formulation of Navier–Maxwell slip boundary conditions for lattice Boltzmann simulations of rarefied flows in microchannels
We present an implementation of first-order Navier–Maxwell slip boundary conditions for simulating near-continuum rarefied flows in microchannels with the lattice Boltzmann method. Rather than imposing boundary conditions directly on the particle velocity distribution functions, following the existing discrete analogs of the specular and diffuse reflection conditions from continuous kinetic theory, we use a moment-based method to impose the Navier–Maxwell slip boundary conditions that relate the velocity and the strain rate at the boundary. We use these conditions to solve for the unknown distribution functions that propagate into the\ud
domain across the boundary. We achieve second-order accuracy by reformulating these conditions for the second set of distribution functions that arise in the derivation of the lattice Boltzmann method by an integration along characteristics. The results are in excellent agreement with asymptotic solutions of the compressible Navier-Stokes equations for microchannel flows in the slip regime. Our moment formalism is also valuable for analysing the existing boundary conditions, and explains the origin of numerical slip in the bounce-back and other common boundary conditions that impose explicit conditions on the higher moments instead of on the local tangential velocity
Mesoscopic simulation of diffusive contaminant spreading in gas flows at low pressure
Many modern production and measurement facilities incorporate multiphase
systems at low pressures. In this region of flows at small, non-zero Knudsen-
and low Mach numbers the classical mesoscopic Monte Carlo methods become
increasingly numerically costly. To increase the numerical efficiency of
simulations hybrid models are promising. In this contribution, we propose a
novel efficient simulation approach for the simulation of two phase flows with
a large concentration imbalance in a low pressure environment in the low
intermediate Knudsen regime. Our hybrid model comprises a lattice-Boltzmann
method corrected for the lower intermediate Kn regime proposed by Zhang et al.
for the simulation of an ambient flow field. A coupled event-driven
Monte-Carlo-style Boltzmann solver is employed to describe particles of a
second species of low concentration. In order to evaluate the model, standard
diffusivity and diffusion advection systems are considered.Comment: 9 pages, 8 figure
Inertial Frame Independent Forcing for Discrete Velocity Boltzmann Equation: Implications for Filtered Turbulence Simulation
We present a systematic derivation of a model based on the central moment
lattice Boltzmann equation that rigorously maintains Galilean invariance of
forces to simulate inertial frame independent flow fields. In this regard, the
central moments, i.e. moments shifted by the local fluid velocity, of the
discrete source terms of the lattice Boltzmann equation are obtained by
matching those of the continuous full Boltzmann equation of various orders.
This results in an exact hierarchical identity between the central moments of
the source terms of a given order and the components of the central moments of
the distribution functions and sources of lower orders. The corresponding
source terms in velocity space are then obtained from an exact inverse
transformation due to a suitable choice of orthogonal basis for moments.
Furthermore, such a central moment based kinetic model is further extended by
incorporating reduced compressibility effects to represent incompressible flow.
Moreover, the description and simulation of fluid turbulence for full or any
subset of scales or their averaged behavior should remain independent of any
inertial frame of reference. Thus, based on the above formulation, a new
approach in lattice Boltzmann framework to incorporate turbulence models for
simulation of Galilean invariant statistical averaged or filtered turbulent
fluid motion is discussed.Comment: 37 pages, 1 figur
Detailed analysis of the lattice Boltzmann method on unstructured grids
The lattice Boltzmann method has become a standard for efficiently solving
problems in fluid dynamics. While unstructured grids allow for a more efficient
geometrical representation of complex boundaries, the lattice Boltzmann methods
is often implemented using regular grids. Here we analyze two implementations
of the lattice Boltzmann method on unstructured grids, the standard forward
Euler method and the operator splitting method. We derive the evolution of the
macroscopic variables by means of the Chapman-Enskog expansion, and we prove
that it yields the Navier-Stokes equation and is first order accurate in terms
of the temporal discretization and second order in terms of the spatial
discretization. Relations between the kinetic viscosity and the integration
time step are derived for both the Euler method and the operator splitting
method. Finally we suggest an improved version of the bounce-back boundary
condition. We test our implementations in both standard benchmark geometries
and in the pore network of a real sample of a porous rock.Comment: 42 page
Ludwig: A parallel Lattice-Boltzmann code for complex fluids
This paper describes `Ludwig', a versatile code for the simulation of
Lattice-Boltzmann (LB) models in 3-D on cubic lattices. In fact `Ludwig' is not
a single code, but a set of codes that share certain common routines, such as
I/O and communications. If `Ludwig' is used as intended, a variety of complex
fluid models with different equilibrium free energies are simple to code, so
that the user may concentrate on the physics of the problem, rather than on
parallel computing issues. Thus far, `Ludwig''s main application has been to
symmetric binary fluid mixtures. We first explain the philosophy and structure
of `Ludwig' which is argued to be a very effective way of developing large
codes for academic consortia. Next we elaborate on some parallel implementation
issues such as parallel I/O, and the use of MPI to achieve full portability and
good efficiency on both MPP and SMP systems. Finally, we describe how to
implement generic solid boundaries, and look in detail at the particular case
of a symmetric binary fluid mixture near a solid wall. We present a novel
scheme for the thermodynamically consistent simulation of wetting phenomena, in
the presence of static and moving solid boundaries, and check its performance.Comment: Submitted to Computer Physics Communication
- …