Boundary condition at a two-phase interface in the lattice Boltzmann method for the convection-diffusion equation

Phase-field-based lattice Boltzmann method for two-phase flows with interfacial mass or heat transfer

In this work, we develop a phase-field-based lattice Boltzmann (LB) method for a two-scalar model of the two-phase flows with interfacial mass or heat transfer. Through the Chapman-Enskog analysis, we show that the present LB method can correctly recover the governing equations for phase field, flow field, and concentration or temperature field. In particular, to derive the two-scalar equations for the mass or heat transfer, we propose a new LB model with an auxiliary source distribution function to describe the extra flux terms, and the discretizations of some derivative terms can be avoided. The accuracy and efficiency of the present LB method are also tested through several benchmark problems, and the influence of mass or heat transfer on the fluid viscosity is further considered by introducing an exponential relation. The numerical results show that the present LB method is suitable for the two-phase flows with interfacial mass or heat transfer.

Off-lattice Boltzmann simulation of conjugate heat transfer for natural convection in two-dimensional cavities

This study addresses the inadequacy of isothermal wall conditions in predicting accurate flow features and thermal effects in multicomponent systems. A finite-difference characteristic-based off-lattice Boltzmann method (OLBM) with a source term-based conjugate heat transfer (CHT) model is utilized to analyze buoyancy-driven flows in two-dimensional enclosures. The source term-based CHT model [Karani and Huber, Phys. Rev. E 91, 023304 (2015)1539-375510.1103/PhysRevE.91.023304] is extended to handle curved conjugate boundaries. The proposed CHT-OLBM solver is verified using analytical solutions and reference data from the literature. The effects of wall conduction on conjugate natural convection (CNC) problems in square and horizontal annular cavities are systematically examined with a solid wall of a nondimensional thickness of 0.2. For the square cavity problem, the governing parameters considered are 10^{5}≤Gr≤10^{9} and χ=1,5,10, while for the horizontal annulus problem, the governing parameters are taken as 10^{5}≤Ra≤10^{7} and χ=1,5,10,100, where Gr is the Grashof number, Ra is the Rayleigh number, and χ is the thermal conductivity ratio. Qualitative analysis of the simulation results using isotherms and streamlines and quantitative analysis of the local fluid-solid interface temperature, solid wall temperature distribution, Nusselt number profiles, overall Nusselt number, and effective Grashof/Rayleigh number is conducted. The overall heat transfer reduction inside the cavities due to a solid wall is quantified, and correlations are obtained to represent the overall heat transfer inside both enclosures. The findings demonstrate the significance of CHT analysis, and the CHT-OLBM solver developed in this study can successfully investigate steady, unsteady, and chaotic CNC flows.

Comparative investigation of a lattice Boltzmann boundary treatment of multiphase mass transport with heterogeneous chemical reactions

Multiphase reactive transport in porous media is an important component of many natural and engineering processes. In the present study, boundary schemes for the continuum species transport-lattice Boltzmann (CST-LB) mass transport model and the multicomponent pseudopotential model are proposed to simulate heterogeneous chemical reactions in a multiphase system. For the CST-LB model, a lattice-interface-tracking scheme for the heterogeneous chemical reaction boundary is provided. Meanwhile, a local-average virtual density boundary scheme for the multicomponent pseudopotential model is formulated based on the work of Li et al. [Li, Yu, and Luo, Phys. Rev. E 100, 053313 (2019)10.1103/PhysRevE.100.053313]. With these boundary treatments, a numerical implementation is put forward that couples the multiphase fluid flow, interfacial species transport, heterogeneous chemical reactions, and porous matrix structural evolution. A series of comparison benchmark cases are investigated to evaluate the numerical performance for different pseudopotential wetting boundary treatments, and an application case of multiphase dissolution in porous media is conducted to validate the present models' ability to solve complex problems. By applying the present LB models with reasonable boundary treatments, multiphase reactive transport in various natural or engineering scenarios can be simulated accurately.

Analysis of the typical unified lattice Boltzmann models and a comprehensive multiphase model for convection-diffusion problems in multiphase systems

The present paper analyzes the typical unified lattice Boltzmann (LB) models for different convection-diffusion (CD) problems in multiphase systems. The CD problems in multiphase systems can be roughly classified into three groups: CD problems with a continuous scalar value and a continuous flux, a discontinuous scalar value and a continuous flux, a continuous scalar value and a discontinuous flux. The characteristics of the corresponding unified LB models for the three kinds of CD problems are analyzed and the equivalence between the LB models based on different perspectives or numerical schemes is revealed. Finally, a comprehensive multiphase LB model (CMLBM) capable of solving different isotropic and anisotropic CD problems in multiphase systems is proposed. Four typical CD problems in multiphase systems are calculated to validate the CMLBM; the results show that it performs well against the typical isotropic and anisotropic CD problems in multiphase systems.

Role of thermal disequilibrium on natural convection in porous media: Insights from pore-scale study

The present study investigates the role of thermal nonequilibrium on natural convection in a fluid-saturated porous medium heated from below. We conduct high-resolution direct numerical simulation at the pore scale in a two-dimensional regular porous structure by means of the thermal lattice-Boltzmann method (LBM). We perform a combination of linear stability analysis of continuum-scale heat transfer models, and pore-scale and continuum-scale simulations to study the role of thermal conductivity contrasts among phases on natural convection. The comparison of pore-scale lattice-Boltzmann simulations with linear stability analysis reveals that traditional continuum-scale models fail to capture the correct onset of convection, convection mode, and heat transfer when the thermal conductivity of the solid obstacles does not match that of the fluid.

Prediction of the moments in advection-diffusion lattice Boltzmann method. I. Truncation dispersion, skewness, and kurtosis

The effect of the heterogeneity in the soil structure or the nonuniformity of the velocity field on the modeled resident time distribution (RTD) and breakthrough curves is quantified by their moments. While the first moment provides the effective velocity, the second moment is related to the longitudinal dispersion coefficient (k_{T}) in the developed Taylor regime; the third and fourth moments are characterized by their normalized values skewness (Sk) and kurtosis (Ku), respectively. The purpose of this investigation is to examine the role of the truncation corrections of the numerical scheme in k_{T}, Sk, and Ku because of their interference with the second moment, in the form of the numerical dispersion, and in the higher-order moments, by their definition. Our symbolic procedure is based on the recently proposed extended method of moments (EMM). Originally, the EMM restores any-order physical moments of the RTD or averaged distributions assuming that the solute concentration obeys the advection-diffusion equation in multidimensional steady-state velocity field, in streamwise-periodic heterogeneous structure. In our work, the EMM is generalized to the fourth-order-accurate apparent mass-conservation equation in two- and three-dimensional duct flows. The method looks for the solution of the transport equation as the product of a long harmonic wave and a spatially periodic oscillating component; the moments of the given numerical scheme are derived from a chain of the steady-state fourth-order equations at a single cell. This mathematical technique is exemplified for the truncation terms of the two-relaxation-time lattice Boltzmann scheme, using plug and parabolic flow in straight channel and cylindrical capillary with the d2Q9 and d3Q15 discrete velocity sets as simple but illustrative examples. The derived symbolic dependencies can be readily extended for advection by another, Newtonian or non-Newtonian, flow profile in any-shape open-tabular conduits. It is established that the truncation errors in the three transport coefficients k_{T}, Sk, and Ku decay with the second-order accuracy. While the physical values of the three transport coefficients are set by Péclet number, their truncation corrections additionally depend on the two adjustable relaxation rates and the two adjustable equilibrium weight families which independently determine the convective and diffusion discretization stencils. We identify flow- and dimension-independent optimal strategies for adjustable parameters and confront them to stability requirements. Through specific choices of two relaxation rates and weights, we expect our results be directly applicable to forward-time central differences and leap-frog central-convective Du Fort-Frankel-diffusion schemes. In straight channel, a quasi-exact validation of the truncation predictions through the numerical moments becomes possible thanks to the specular-forward no-flux boundary rule. In the staircase description of a cylindrical capillary, we account for the spurious boundary-layer diffusion and dispersion because of the tangential constraint of the bounce-back no-flux boundary rule.

Effects of tangential-type boundary condition discontinuities on the accuracy of the lattice Boltzmann method for heat and mass transfer

We present a systematic study on the effects of tangential-type boundary condition discontinuities on the accuracy of the lattice Boltzmann equation (LBE) method for Dirichlet and Neumann problems in heat and mass transfer modeling. The second-order accurate boundary condition treatments for continuous Dirichlet and Neumann problems are directly implemented for the corresponding discontinuous boundary conditions. Results from three numerical tests, including both straight and curved boundaries, are presented to show the accuracy and order of convergence of the LBE computations. Detailed error assessments are conducted for the interior temperature or concentration (denoted as a scalar ϕ) and the interior derivatives of ϕ for both types of boundary conditions, for the boundary flux in the Dirichlet problem and for the boundary ϕ values in the Neumann problem. When the discontinuity point on the straight boundary is placed at the center of the unit lattice in the Dirichlet problem, it yields only first-order accuracy for the interior distribution of ϕ, first-order accuracy for the boundary flux, and zeroth-order accuracy for the interior derivatives compared with the second-order accuracy of all quantities of interest for continuous boundary conditions. On the lattice scale, the LBE solution for the interior derivatives near the singularity is largely independent of the resolution and correspondingly the local distribution of the absolute errors is almost invariant with the changing resolution. For Neumann problems, when the discontinuity is placed at the lattice center, second-order accuracy is preserved for the interior distribution of ϕ; and a "superlinear" convergence order of 1.5 for the boundary ϕ values and first-order accuracy for the interior derivatives are obtained. For straight boundaries with the discontinuity point arbitrarily placed within the lattice and curved boundaries, the boundary flux becomes zeroth-order accurate for Dirichlet problems; and all three quantities, including the interior and boundary ϕ values and the interior derivatives, are only first-order accurate for Neumann problems.

Conjugate heat transfer with the entropic lattice Boltzmann method

A conjugate heat-transfer model is presented based on the two-population entropic lattice Boltzmann method. The present approach relies on the extension of Grad's boundary conditions to the two-population model for thermal flows, as well as on the appropriate exact conjugate heat-transfer condition imposed at the fluid-solid interface. The simplicity and efficiency of the lattice Boltzmann method (LBM), and in particular of the entropic multirelaxation LBM, are retained in the present approach, thus enabling simulations of turbulent high Reynolds number flows and complex wall boundaries. The model is validated by means of two-dimensional parametric studies of various setups, including pure solid conduction, conjugate heat transfer with a backward-facing step flow, and conjugate heat transfer with the flow past a circular heated cylinder. Further validations are performed in three dimensions for the case of a turbulent flow around a heated mounted cube.

Lattice Boltzmann method for convection-diffusion equations with general interfacial conditions

In this work, we propose an interfacial scheme accompanying the lattice Boltzmann method for convection-diffusion equations with general interfacial conditions, including conjugate conditions with or without jumps in heat and mass transfer, continuity of macroscopic variables and normal fluxes in ion diffusion in porous media with different porosity, and the Kapitza resistance in heat transfer. The construction of this scheme is based on our boundary schemes [Huang and Yong, J. Comput. Phys. 300, 70 (2015)JCTPAH0021-999110.1016/j.jcp.2015.07.045] for Robin boundary conditions on straight or curved boundaries. It gives second-order accuracy for straight interfaces and first-order accuracy for curved ones. In addition, the new scheme inherits the advantage of the boundary schemes in which only the current lattice nodes are involved. Such an interfacial scheme is highly desirable for problems with complex geometries or in porous media. The interfacial scheme is numerically validated with several examples. The results show the utility of the constructed scheme and very well support our theoretical predications.

Counter-extrapolation method for conjugate interfaces in computational heat and mass transfer

In this paper a conjugate interface method is developed by performing extrapolations along the normal direction. Compared to other existing conjugate models, our method has several technical advantages, including the simple and straightforward algorithm, accurate representation of the interface geometry, applicability to any interface-lattice relative orientation, and availability of the normal gradient. The model is validated by simulating the steady and unsteady convection-diffusion system with a flat interface and the steady diffusion system with a circular interface, and good agreement is observed when comparing the lattice Boltzmann results with respective analytical solutions. A more general system with unsteady convection-diffusion process and a curved interface, i.e., the cooling process of a hot cylinder in a cold flow, is also simulated as an example to illustrate the practical usefulness of our model, and the effects of the cylinder heat capacity and thermal diffusivity on the cooling process are examined. Results show that the cylinder with a larger heat capacity can release more heat energy into the fluid and the cylinder temperature cools down slower, while the enhanced heat conduction inside the cylinder can facilitate the cooling process of the system. Although these findings appear obvious from physical principles, the confirming results demonstrates the application potential of our method in more complex systems. In addition, the basic idea and algorithm of the counter-extrapolation procedure presented here can be readily extended to other lattice Boltzmann models and even other computational technologies for heat and mass transfer systems.

Analysis and improvement of Brinkman lattice Boltzmann schemes: bulk, boundary, interface. Similarity and distinctness with finite elements in heterogeneous porous media

This work focuses on the numerical solution of the Stokes-Brinkman equation for a voxel-type porous-media grid, resolved by one to eight spacings per permeability contrast of 1 to 10 orders in magnitude. It is first analytically demonstrated that the lattice Boltzmann method (LBM) and the linear-finite-element method (FEM) both suffer from the viscosity correction induced by the linear variation of the resistance with the velocity. This numerical artefact may lead to an apparent negative viscosity in low-permeable blocks, inducing spurious velocity oscillations. The two-relaxation-times (TRT) LBM may control this effect thanks to free-tunable two-rates combination Λ. Moreover, the Brinkman-force-based BF-TRT schemes may maintain the nondimensional Darcy group and produce viscosity-independent permeability provided that the spatial distribution of Λ is fixed independently of the kinematic viscosity. Such a property is lost not only in the BF-BGK scheme but also by "partial bounce-back" TRT gray models, as shown in this work. Further, we propose a consistent and improved IBF-TRT model which vanishes viscosity correction via simple specific adjusting of the viscous-mode relaxation rate to local permeability value. This prevents the model from velocity fluctuations and, in parallel, improves for effective permeability measurements, from porous channel to multidimensions. The framework of our exact analysis employs a symbolic approach developed for both LBM and FEM in single and stratified, unconfined, and bounded channels. It shows that even with similar bulk discretization, BF, IBF, and FEM may manifest quite different velocity profiles on the coarse grids due to their intrinsic contrasts in the setting of interface continuity and no-slip conditions. While FEM enforces them on the grid vertexes, the LBM prescribes them implicitly. We derive effective LBM continuity conditions and show that the heterogeneous viscosity correction impacts them, a property also shared by FEM for shear stress. But, in contrast with FEM, effective velocity conditions in LBM give rise to slip velocity jumps which depend on (i) neighbor permeability values, (ii) resolution, and (iii) control parameter Λ, ranging its reliable values from Poiseuille bounce-back solution in open flow to zero in Darcy's limit. We suggest an "upscaling" algorithm for Λ, from multilayers to multidimensions in random extremely dispersive samples. Finally, on the positive side for LBM besides its overall versatility, the implicit boundary layers allow for smooth accommodation of the flat discontinuous Darcy profiles, quite deficient in FEM.