Multi-block Lattice Boltzmann method: Extension to 3D and validation in turbulence

In this paper, a multi-block method is applied to 3D problems. In this application, interpolation schemes are carefully addressed to avoid the inconsistency when information is transferred from coarse to fine blocks. Two test cases are employed to assess information transfer scheme and accuracy improvement with respect to grid refinement.     

Lattice Boltzmann model for three-dimensional decaying homogeneous isotropic turbulence

We implement a lattice Boltzmann method (LBM) for decaying homogeneous isotropic turbulence based on an analogous Galerkin filter and focus on the fundamental statistical isotropic property. This regularized method is constructed based on orthogonal Hermite polynomial space. For decaying homogeneous isotropic turbulence, this regularized method can simulate the isotropic property very well. Numerical studies demonstrate that the novel regularized LBM is a promising approximation of turbulent fluid flows, which paves the way for coupling various turbulent models with LBM.     

Lattice Boltzmann representation of neutral turbulence in the cold gas blanket divertor regime

In the divertor plasma detachment regime one expects to encounter turbulent collisional neutral gas dynamics. These conditions are ideally simulated by lattice Boltzmann methods — methods ideal for application on multi-parallel processors. As illustrative of this method, we specifically show the effect of conductivity in the Rayleigh-Benard convection cells, using an octagonal velocity lattice. This work was supported by DoE and a joint U.S-Czech grant. These results were presented at the 6th International Workshop on Plasma Edge Theory in Fusion Devices, 15–17 September 1997.     

Multiple-Relaxation-Time Lattice Boltzmann Approach to Richtmyer-Meshkov Instability

The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time (MRT) Lattice Boltzmann (LB) model proposed in [Europhys. Lett. 90 (2010) 54003]. We discuss the reason why the Gram-Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka-Tsutahara model [Phys. Rev. E 69 (2004) 035701(R)] is only valid in subsonic flows.The von Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer-Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang-Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.     

Entropy Function Approach to the Lattice Boltzmann Method

In this paper, we present the construction of the Lattice Boltzmann method equipped with the H-theorem. Based on entropy functions whose local equilibria are suitable to recover the Navier–Stokes equations in the framework of the Lattice Boltzmann method, we derive a collision integral which enables simple identification of transport coefficients, and which circumvents construction of the equilibrium. We discuss performance of this approach as compared to the standard realizations.     

Gyrokinetic simulations of solar wind turbulence from ion to electron scales

A three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms. The resulting energy spectra are quantitatively consistent with a magnetic power spectrum scaling of k(-2.8) as observed in in situ spacecraft measurements of the "dissipation range" of solar wind turbulence. Despite the strongly nonlinear nature of the turbulence, the linear kinetic Alfvén wave mode quantitatively describes the polarization of the turbulent fluctuations. The collisional ion heating is measured at subion-Larmor radius scales, which provides evidence of the ion entropy cascade in an electromagnetic turbulence simulation.     

Mesoscale simulations: Lattice Boltzmann and particle algorithms

I introduce two mesoscale algorithms, lattice Boltzmann and stochastic rotation dynamics, and show how they can be used to investigate the hydrodynamics of complex fluids. For each method I describe the algorithm, show that it solves the Navier–Stokes equations, and then discuss physical problems where it is particularly applicable. For lattice Boltzmann the examples I choose are phase ordering in a binary fluid and drop dynamics on a chemically patterned surface. For stochastic rotation dynamics I consider the hydrodynamics of dilute polymer solutions, concentrating on shear thinning and translocation across a barrier.     

A Lattice Boltzmann Model for Anisotropic Crystal Growth from Melt

We coupled the lattice Boltzmann method with enhanced collisions for hydrodynamics with a model for the anisotropic liquid/solid phase transition. The model is based on a simple reaction model. As a test we have performed calculations for dendritic growth of a crystal into an undercooled melt.     

Discretized Boltzmann equation: Lattice limit and non-Maxwellian gases

We continue the study of a discrete model of the Boltzmann equation, in which the spatial variable is replaced by a finite periodic lattice. Using a weak compactness criterion forL ₁, the existence of a lattice limit as the lattice spacing tends to zero is proved. The case of unbounded collision kernels (non-Maxwellian gases) is also treated.This work was supported in part by National Science Foundation grant no. ENG-75 15882.On leave from Mathematisches Institut der Universität München, Federal Republic of Germany     

Boundary Slip and Surface Interaction： A Lattice Boltzmann Simulation

The factors affecting slip length in Couette geometry flows are analysed by means of a two-phase mesoscopic lattice Boltzmann model including non-ideal fluid-fluid and fluid-wall interactions. The main factors influencing the boundary slip are the strength of interactions between fluid-fluid and fluid-wall particles. Other factors, such as fluid viscosity, bulk pressure may also change the slip length. We find that boundary slip only occurs under a certain density （bulk pressure）. If the density is large enough, the slip length will tend to zero. In our simulations, a low density layer near the wall does not need to be postulated a priori but emerges naturally from the underlying non-ideal mesoscopic dynamics. It is the low density layer that induces the boundary slip. The results may be helpful to understand recent experimental observations on the slippage of micro flows.     

Lattice Boltzmann Solver for Multiphase Flows: Application to High Weber and Reynolds Numbers

The lattice Boltzmann method, now widely used for a variety of applications, has also been extended to model multiphase flows through different formulations. While already applied to many different configurations in low Weber and Reynolds number regimes, applications to higher Weber/Reynolds numbers or larger density/viscosity ratios are still the topic of active research. In this study, through a combination of a decoupled phase-field formulation—the conservative Allen–Cahn equation—and a cumulant-based collision operator for a low-Mach pressure-based flow solver, we present an algorithm that can be used for higher Reynolds/Weber numbers. The algorithm was validated through a variety of test cases, starting with the Rayleigh–Taylor instability in both 2D and 3D, followed by the impact of a droplet on a liquid sheet. In all simulations, the solver correctly captured the flow dynamics andmatched reference results very well. As the final test case, the solver was used to model droplet splashing on a thin liquid sheet in 3D with a density ratio of 1000 and kinematic viscosity ratio of 15, matching the water/air system at We = 8000 and Re = 1000. Results showed that the solver correctly captured the fingering instabilities at the crown rim and their subsequent breakup, in agreement with experimental and numerical observations reported in the literature.     

Application of Lattice Boltzmann Method to sensitivity analysis via complex differentiation

Sensitivity analysis of flows computed by Lattice Boltzmann Method via Complex Differentiation is proposed. The theoretical work is illustrated and the proposed method assessed considering the D2Q9 scheme, along with the differentiation of the Lattice Boltzmann Equation. Boundary condition implementation is also detailed. Some examples illustrate the capability of the proposed method.     

Continuous and Discrete Adjoint Approach Based on Lattice Boltzmann Method in Aerodynamic Optimization Part I: Mathematical Derivation of Adjoint Lattice Boltzmann Equations

The significance of flow optimization utilizing the lattice Boltzmann (LB) method becomes obvious regarding its advantages as a novel flow field solution method compared to the other conventional computational fluid dynamics techniques. These unique characteristics of the LB method form the main idea of its application to optimization problems. In this research, for the first time, both continuous and discrete adjoint equations were extracted based on the LB method using a general procedure with low implementation cost. The proposed approach could be performed similarly for any optimization problem with the corresponding cost function and design variables vector. Moreover, this approach was not limited to flow fields and could be employed for steady as well as unsteady flows. Initially, the continuous and discrete adjoint LB equations and the cost function gradient vector were derived mathematically in detail using the continuous and discrete LB equations in space and time, respectively. Meanwhile, new adjoint concepts in lattice space were introduced. Finally, the analytical evaluation of the adjoint distribution functions and the cost function gradients was carried out.     

Highly Efficient Lattice Boltzmann Model for Compressible Fluids:Two-Dimensional Case

We present a highly efficient lattice Boltzmann model for simulatingcompressible flows. This model is based on the combination of an appropriatefinite difference scheme, a 16-discrete-velocity model [Kataoka andTsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion anddissipation terms. The dispersion term effectively reduces the oscillationat the discontinuity and enhances numerical precision. The dissipation termmakes the new model more easily meet with the von Neumann stabilitycondition. This model works for both high-speed and low-speed flows witharbitrary specific-heat-ratio. With the new model simulation results for thewell-known benchmark problems get a high accuracy compared with the analytic or experimental ones. 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. With a more realistic equation ofstate or free-energy functional, the new model has the potential tostudythe complex procedure of shock wave reaction on porous materials.     

Finite-Difference Lattice Boltzmann Scheme for High-Speed Compressible Flow: Two-Dimensional Case

Lattice Boltzmann (LB) modeling of high-speed compressible flows has long been attempted by various authors. One common weakness of most of previous models is the instability problem when the Mach number of the flow is large. In this paper we present a finite-difference LB model, which works for flows with flexible ratios of specific heats and a wide range of Mach number, from $0$ to 30 or higher. Besides the discrete-velocity-model by Watari [Physica A 382 (2007) 502], a modified Lax--Wendroff finite difference scheme and an artificial viscosity are introduced. The combination of the finite-difference scheme and the adding of artificial viscosity must find a balance of numerical stability versus accuracy. The proposed model is validated by recovering results of some well-known benchmark tests: shock tubes and shock reflections. The new model may be used to track shock waves and/or to study the non-equilibrium procedure in the transition between the regular and Mach reflections of shock waves, etc.     

Lattice Boltzmann simulations of capillary filling: Finite vapour density effects

Numerical simulations of two-dimensional capillary filling using the pseudo-potential lattice Boltzmann model for multiphase fluids are presented. It is shown that whenever the density of the light-phase exceeds about ten percent of the dense phase, the front motion proceeds through a combined effect of capillary advection and condensation.     

Lattice Boltzmann Simulations of Blood Flow: Non-Newtonian Rheology and Clotting Processes

The numerical simulation of thrombosis in stented aneurysms is an important issue to estimate the efficiency of a stent. In this paper, we consider a Lattice Boltzmann (LB) approach to bloodflow modeling and we implement a non-Newtonian correction in order to reproduce more realistic flow profiles. We obtain a good agreement between simulations and Casson’s model of blood rheology in a simple geometry. Finally we discuss how, by using a passive scalar suspension model with aggregation on top of the LB dynamics, we can describe the clotting processes in the aneurysm     

Study of Kinetic Equation for Non-Ideal Gases Using Lattice Boltzmann Method: Interparticle Forces from BBGKY Hierarchy

The first equation of the BBGKY hierarchy is closed with a truncation scheme which conserves mass, momentum and total energy. Resulting kinetic equation is solved numerically using the lattice Boltzmann method over an isothermal \(D2Q9\) lattice. The inter particle forces are derived using a density gradient approximation to the collision term of the BBGKY hierarchy and the resulting non-ideal equation of state is calculated. In this paper, we show analytically that the time rate of change of Boltzmann’s H-function i.e., \(\frac{dH \left( t \right) }{dt}\) is zero for this truncation scheme, implying reversibility in time. Therefore, a BGK term is added to the collision integral to describe the macroscopic irreversibility. Both the local and global \(H\) -functions are calculated numerically for significant variations in the density. We also have simulated a switching experiment over this system and obtained the non-equilibrium work distributions.     

格子Boltzmann方法模拟填充纳米流体三角形腔内的自然对流

采用格子Boltzmann方法研究填充水-氧化铝纳米流体的等腰直角三角形腔体中的自然对流.讨论瑞利数、颗粒体积分数、热源位置等因素对对流换热的影响,以及不同纳米流体模型对模拟结果的影响.结果表明:在低瑞利数下,随着热源在左壁面向上移动,换热效率逐渐增加.而在高瑞利数(Ra=106)时,观察到相反的现象;采用单相纳米流体...