We describe how regularization of lattice Boltzmann methods can be achieved by modifying dissipation. Classes of techniques used to try to improve regularization of LBMs include flux limiters, enforcing the exact correct production of entropy and man
ipulating non-hydrodynamic modes of the system in relaxation. Each of these techniques corresponds to an additional modification of dissipation compared with the standard LBGK model. Using some standard 1D and 2D benchmarks including the shock tube and lid driven cavity, we explore the effectiveness of these classes of methods.
In the absence of fractures, methane bubbles in deep-water sediments can be immovably trapped within a porous matrix by surface tension. The dominant mechanism of transfer of gas mass therefore becomes the diffusion of gas molecules through porewater
. The accurate description of this process requires non-Fickian diffusion to be accounted for, including both thermodiffusion and gravitational action. We evaluate the diffusive flux of aqueous methane considering non-Fickian diffusion and predict the existence of extensive bubble mass accumulation zones within deep-water sediments. The limitation on the hydrate deposit capacity is revealed; too weak deposits cannot reach the base of the hydrate stability zone and form any bubbly horizon.
The Navier--Stokes equations arise naturally as a result of Ehrenfests coarse-graining in phase space after a period of free-flight dynamics. This point of view allows for a very flexible approach to the simulation of fluid flow for high-Reynolds num
ber. We construct regularisers for lattice Boltzmann computational models. These regularisers are based on Ehrenfests coarse-graining idea and could be applied to schemes with either entropic or non-entropic quasiequilibria. We give a numerical scheme which gives good results for the standard test cases of the shock tube and the flow past a square cylinder.
We adapt Schabacks error doubling trick [R. Schaback. Improved error bounds for scattered data interpolation by radial basis functions. Math. Comp., 68(225):201--216, 1999.] to give error estimates for radial interpolation of functions with smoothnes
s lying (in some sense) between that of the usual native space and the subspace with double the smoothness. We do this for both bounded subsets of R^d and spheres. As a step on the way to our ultimate goal we also show convergence of pseudoderivatives of the interpolation error.
In this paper we consider two sources of enhancement for the meshfree Lagrangian particle method smoothed particle hydrodynamics (SPH) by improving the accuracy of the particle approximation. Namely, we will consider shape functions constructed using
: moving least-squares approximation (MLS); radial basis functions (RBF). Using MLS approximation is appealing because polynomial consistency of the particle approximation can be enforced. RBFs further appeal as they allow one to dispense with the smoothing-length -- the parameter in the SPH method which governs the number of particles within the support of the shape function. Currently, only ad hoc methods for choosing the smoothing-length exist. We ensure that any enhancement retains the conservative and meshfree nature of SPH. In doing so, we derive a new set of variationally-consistent hydrodynamic equations. Finally, we demonstrate the performance of the new equations on the Sod shock tube problem.
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.