Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advecti...Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state...The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state. The calculation results indicate that the motion of the droplet on the cylinder can be divided into three stages: spreading, sliding, and aggregating.The contact length and contact time of a droplet on a cylindrical surface can be affected by factors such as the wettability gradient of the cylindrical wall, the Bond number, and droplet size. Furthermore, phase diagrams showing the relationship between Bond number, cylinder wall wettability gradient, and contact time as well as maximum contact length for three different droplet sizes are given. A theoretical foundation for additional research into the heat and mass transfer process between the droplet and the cylinder can be established by comprehending the variable rules of maximum contact length and contact time.展开更多
With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium...With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.展开更多
Conventional seismic wave forward simulation generally uses mathematical means to solve the macroscopic wave equation,and then obtains the corresponding seismic wavefield.Usually,when the subsurface structure is finel...Conventional seismic wave forward simulation generally uses mathematical means to solve the macroscopic wave equation,and then obtains the corresponding seismic wavefield.Usually,when the subsurface structure is finely constructed and the continuity of media is poor,this strategy is difficult to meet the requirements of accurate wavefield calculation.This paper uses the multiple-relaxation-time lattice Boltzmann method(MRT-LBM)to conduct the seismic acoustic wavefield simulation and verify its computational accuracy.To cope with the problem of severe reflections at the truncated boundaries,we analogize the viscous absorbing boundary and perfectly matched layer(PML)absorbing boundary based on the single-relaxation-time lattice Boltzmann(SRT-LB)equation to the MRT-LB equation,and further,propose a joint absorbing boundary through comparative analysis.We give the specific forms of the modified MRT-LB equation loaded with the joint absorbing boundary in the two-dimensional(2D)and three-dimensional(3D)cases,respectively.Then,we verify the effects of this absorbing boundary scheme on a 2D homogeneous model,2D modified British Petroleum(BP)gas-cloud model,and 3D homogeneous model,respectively.The results reveal that by comparing with the viscous absorbing boundary and PML absorbing boundary,the joint absorbing boundary has the best absorption performance,although it is a little bit complicated.Therefore,this joint absorbing boundary better solves the problem of truncated boundary reflections of MRT-LBM in simulating seismic acoustic wavefields,which is pivotal to its wide application in the field of exploration seismology.展开更多
Sharp phase interfaces and accurate temperature distributions are important criteria in the simulation of solid-liquid phase changes.The multi-relaxation-time lattice Boltzmann method(MRT-LBM)shows great numerical per...Sharp phase interfaces and accurate temperature distributions are important criteria in the simulation of solid-liquid phase changes.The multi-relaxation-time lattice Boltzmann method(MRT-LBM)shows great numerical performance during simulation;however,the value method of the relaxation parameters needs to be specified.Therefore,in this study,a random forest(RF)model is used to discriminate the importance of different relaxation parameters to the convergence,and a support vector machine(SVM)is used to explore the decision boundary of the convergent samples in each dimensional model.The results show that the convergence of the samples is consistent with the sign of the decision number,and two types of the numerical deviations appear,i.e.,the phase mushy zone and the non-physical heat transfer.The relaxation parameters chosen on the decision boundary can further suppress the numerical bias and improve numerical accuracy.展开更多
A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, su...A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.展开更多
An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D...An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.展开更多
The standard lattice Boltzmann method utilizes uniform grids to maintain a compact computational procedure. However, it is often less efficient to perform hydrodynamic and aerodynamic flow simulations when there is a ...The standard lattice Boltzmann method utilizes uniform grids to maintain a compact computational procedure. However, it is often less efficient to perform hydrodynamic and aerodynamic flow simulations when there is a need for high resolution. To resolve these difficulties, a multiple nested lattice Boltzmann method(MNLBM) was developed, which contains several overlapped layers with different resolutions in the computational domain. The data transference of flow field on two layers is accomplished by a Filippova procedure which is proved to satisfy the continuity of mass, momentum, and stresses across the interface. The proposed method is based on the standard lattice Boltzmann method, so it is easily performed.By numerical investigation, the result of present method has been agreed with that of literature, but the computation efficiency is higher than the standard lattice Boltzmann method.展开更多
The lattice Boltzmann method (LBM) is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanopa...The lattice Boltzmann method (LBM) is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanoparticles: TiO2, Ag, Cu, and A1203. The Brinkman and Maxwell-Garnetts models are used to simulate the viscosity and the effective thermal conductivity of nanofluids, respectively. Results from the performed numerical analysis show good agreement with those obtained from other numerical meth- ods. A variety of the Rayleigh number, the nanoparticle volume fraction, and the aspect ratio are examined. According to the results, choosing copper as the nanoparticle leads to obtaining the highest enhancement for this problem. The results also indicate that the maximum value of enhancement occurs at λ =2.5 when Ra = 106 while at A = 1.5 for other Rayleigh numbers.展开更多
By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improv...By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.展开更多
The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method (LBM) was used to study the electrophysical property of rock...The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method (LBM) was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation (I-Sw). The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.展开更多
The upper limit of the solid scatters density ns (x), a key parameter for the simulation of flows in porous media with a gray Lattice Boltzmann Method, is studied by an analytical way for the infiltration Poiseuille...The upper limit of the solid scatters density ns (x), a key parameter for the simulation of flows in porous media with a gray Lattice Boltzmann Method, is studied by an analytical way for the infiltration Poiseuille flow between two infinite parallel plates. Analyses of three different gray Lattice Boltzmann schemes, separately proposed by Gao and Sharma et al., Dardis and McCloskey, and Thorne and Sukop, indicate that the effective domain of Gao and Sharma's scheme is restricted to ns 〈 1/2√3≈0.289, Dardis and McCloskey's scheme is restricted to ns 〈 (√57-1)/28≈0.234, and that there is no extra restriction on ns(x) with Thorne and Sukop's scheme. These results are obtained for the dimensionless relaxation time τ= 1. The above analytical results are verified by our numerical simulations. The use of a gray LBM is further illustrated by simulating the flow at the interface of a porous medium. Simulation results yield velocity profiles which agree very well with Brinkman's prediction.展开更多
This paper describes the application of a three-dimensional lattice Boltzmann method (LBM) to Newtonian and non-Newtonian (Bingham fluid in this work) flows with free surfaces. A mass tracking algorithm was incorp...This paper describes the application of a three-dimensional lattice Boltzmann method (LBM) to Newtonian and non-Newtonian (Bingham fluid in this work) flows with free surfaces. A mass tracking algorithm was incorporated to capture the free surface, whereas Papanastasiou's modified model was used for Bingham fluids. The lattice Boltzmann method was first validated using two benchmarks: Newtonian flow through a square cross-section tube and Bingham flow through a circular cross-section tube. Afterward, the dam-break problem for the Newtonian fluid and the slump test for Bingham fluid were simulated to validate the free-surface-capturing algorithm. The numerical results were in good agreement with analytical results, as well as other simulations, thereby proving the validity and correctness of the current method. The proposed method is a promising substitute for time-consuming and costly physical experiments to solve problems encountered in geotechnical and geological engineering, such as the surge and debris flow induced by a landslide or earthquake.展开更多
In the present study,the process of droplet condensation on superhydrophobic nanoarrays is simulated using a multicomponent multi-phase lattice Boltzmann model.The results indicate that three typical nucleation modes ...In the present study,the process of droplet condensation on superhydrophobic nanoarrays is simulated using a multicomponent multi-phase lattice Boltzmann model.The results indicate that three typical nucleation modes of condensate droplets are produced by changing the geometrical parameters of nanoarrays.Droplets nucleated at the top(top-nucleation mode),or in the upside interpillar space of nanoarrays(side-nucleation mode),generate the non-wetting Cassie state,whereas the ones nucleated at the bottom corners between the nanoarrays(bottom-nucleation mode) present the wetting Wenzel state.Time evolutions of droplet pressures at the upside and downside of the liquid phase are analyzed to understand the wetting behaviors of the droplets condensed from different nucleation modes.The phenomena of droplet condensation on nanoarrays patterned with different hydrophilic and hydrophobic regions are simulated,indicating that the nucleation mode of condensate droplets can also be manipulated by modifying the local intrinsic wettability of nanoarray surface.The simulation results are compared well with the experimental observations reported in the literature.展开更多
The three-dimensional (3D) lattice Boltzmann models, 3DQ15, 3DQ19 and 3DQ27, under different wall boundary conditions and lattice resolutions have been investigated by simulating Poiseuille flow in a circular cylind...The three-dimensional (3D) lattice Boltzmann models, 3DQ15, 3DQ19 and 3DQ27, under different wall boundary conditions and lattice resolutions have been investigated by simulating Poiseuille flow in a circular cylinder for a wide range of Reynolds numbers. The 3DQ19 model with improved Fillippova and Hanel (FH) curved boundary condition represents a good compromise between computational efficiency and reliability. Blood flow in an aortic arch is then simulated as a typical haemodynamic application. Axial and secondary fluid velocity and effective wall shear stress profiles in a 180° bend are obtained, and the results also demonstrate that the lattice Boltzmann method is suitable for simulating the flow in 3D large-curved vessels.展开更多
The numerical investigation of the two-dimensional laminar flow past two ro- tating circular cylinders in the tandem arrangement is conducted by the lattice Boltzmann method. The numerical strategy is used for dealing...The numerical investigation of the two-dimensional laminar flow past two ro- tating circular cylinders in the tandem arrangement is conducted by the lattice Boltzmann method. The numerical strategy is used for dealing with curved and moving boundaries of the second-order accuracy for velocity and temperature fields. The effects of various rotational speed ratios and gap spacing are studied with the Reynolds number of 100 and the Prandtl number of 0.71. A varied range of rotational speed ratios are investigated for four different gap spacing, i.e., 3.0, 1.5, 0.7, and 0.2. The results show that, for the first cylinder, the lift and drag coefficients for large gap spacing are similar to those for a single cylinder; for the second cylinder, the lift coefficient descends with the increase in the angular velocity for all gap spacing, while the drag coefficient ascends except for the gap spacing of 3.0. The results of the averaged periodic Nusselt number on the surface of the cylinders show that, for small distances between the cylinders and low angular velocities, conduction is a dominant mechanism of heat transfer, but for large distances and high angular velocities, convection is the main mechanism of heat transfer.展开更多
In this study,a two-dimensional approach is elaborated to study with the lattice Boltzmann method(LBM)the seepage of water in the pores of a soil.Firstly,the D2Q9 model is selected to account for the discrete velocity...In this study,a two-dimensional approach is elaborated to study with the lattice Boltzmann method(LBM)the seepage of water in the pores of a soil.Firstly,the D2Q9 model is selected to account for the discrete velocity distribution of water flow.In particular,impermeability is considered as macroscopic boundary condition for the left and right domain sides,while the upper and lower boundaries are assumed to behave as pressure boundaries controlled by different densities.The micro-boundary conditions are implemented through the standard rebound strategy and a non-equilibrium extrapolation scheme.Matlab is used for the development of the related algorithm.Finally,the influence of porosity,permeability,osmotic pressure and other factors is assessed with regard to seepage characteristics and the ensuing results are compared with Darcy’s law.The computations show that,for fixed initial conditions,the pore structure has a certain influence on the local velocity of seepage,but the overall state is stable,and the average velocity of each layer is the same.The larger the pore passage is,the faster the flow velocity is,and vice versa.For low permeability,the numerical results are consistent with the Darcy's law.The greater the pressure difference between the inlet and outlet of seepage,the greater the seepage rate.The relationship between them is linear(yet in good agreement with Darcy’s law).展开更多
One dimensional traffic model was studied with the consideration of stochastic deceleration by using the lattice Boltzmann method. The evolution equation for vehicle density was derived, which has the form of the Bu...One dimensional traffic model was studied with the consideration of stochastic deceleration by using the lattice Boltzmann method. The evolution equation for vehicle density was derived, which has the form of the Burgers equation exhibiting the effect of viscosity corresponding to the deceleration. The simulation with the model shows that the variation of vehicle density in space tends to a periodic one, which implies the existence of kinematic waves in the 1D traffic flow and coincides with theoretical prediction.展开更多
This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, where...This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time (SRT) lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.展开更多
基金supported as part of the Center for Hierarchical Waste Form Materials,an Energy Frontier Research Center funded by the U.S.Department of Energy,Office of Science,Basic Energy Sciences under Award No.DE-SC0016574.
文摘Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
文摘The study of a droplet spreading on a circular cylinder under gravity was carried out using the pseudo-potential lattice Boltzmann high-density ratios multiphase model with a non-ideal Peng–Robinson equation of state. The calculation results indicate that the motion of the droplet on the cylinder can be divided into three stages: spreading, sliding, and aggregating.The contact length and contact time of a droplet on a cylindrical surface can be affected by factors such as the wettability gradient of the cylindrical wall, the Bond number, and droplet size. Furthermore, phase diagrams showing the relationship between Bond number, cylinder wall wettability gradient, and contact time as well as maximum contact length for three different droplet sizes are given. A theoretical foundation for additional research into the heat and mass transfer process between the droplet and the cylinder can be established by comprehending the variable rules of maximum contact length and contact time.
基金supported by the Hebei Province Graduate Innovation Funding Project(CXZZBS2022029).
文摘With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.
基金This work is supported in part by the National Natural Science Foundation of China(U19B6003-04-01,42204132,41874130)R&D Department of CNPC(2022DQ0604-01)China Postdoctoral Science Foundation(2020M680667,2021T140661).
文摘Conventional seismic wave forward simulation generally uses mathematical means to solve the macroscopic wave equation,and then obtains the corresponding seismic wavefield.Usually,when the subsurface structure is finely constructed and the continuity of media is poor,this strategy is difficult to meet the requirements of accurate wavefield calculation.This paper uses the multiple-relaxation-time lattice Boltzmann method(MRT-LBM)to conduct the seismic acoustic wavefield simulation and verify its computational accuracy.To cope with the problem of severe reflections at the truncated boundaries,we analogize the viscous absorbing boundary and perfectly matched layer(PML)absorbing boundary based on the single-relaxation-time lattice Boltzmann(SRT-LB)equation to the MRT-LB equation,and further,propose a joint absorbing boundary through comparative analysis.We give the specific forms of the modified MRT-LB equation loaded with the joint absorbing boundary in the two-dimensional(2D)and three-dimensional(3D)cases,respectively.Then,we verify the effects of this absorbing boundary scheme on a 2D homogeneous model,2D modified British Petroleum(BP)gas-cloud model,and 3D homogeneous model,respectively.The results reveal that by comparing with the viscous absorbing boundary and PML absorbing boundary,the joint absorbing boundary has the best absorption performance,although it is a little bit complicated.Therefore,this joint absorbing boundary better solves the problem of truncated boundary reflections of MRT-LBM in simulating seismic acoustic wavefields,which is pivotal to its wide application in the field of exploration seismology.
基金the National Natural Science Foundation of China(Nos.12172017 and 12202021)。
文摘Sharp phase interfaces and accurate temperature distributions are important criteria in the simulation of solid-liquid phase changes.The multi-relaxation-time lattice Boltzmann method(MRT-LBM)shows great numerical performance during simulation;however,the value method of the relaxation parameters needs to be specified.Therefore,in this study,a random forest(RF)model is used to discriminate the importance of different relaxation parameters to the convergence,and a support vector machine(SVM)is used to explore the decision boundary of the convergent samples in each dimensional model.The results show that the convergence of the samples is consistent with the sign of the decision number,and two types of the numerical deviations appear,i.e.,the phase mushy zone and the non-physical heat transfer.The relaxation parameters chosen on the decision boundary can further suppress the numerical bias and improve numerical accuracy.
基金Project supported by the National Natural Science Foundation of China (Grant No. 51775077)。
文摘A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.
基金the Youth Program of the National Natural Science Foundation of China(Grant Nos.11972272,12072246,and 12202331)the National Key Project,China(Grant No.GJXM92579)the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2022JQ-028)。
文摘An additional potential energy distribution function is introduced on the basis of previous D3Q25 model,and the equilibrium distribution function of D3Q25 is obtained by spherical function.A novel three-dimensional(3D)shifted lattice model is proposed,therefore a shifted lattice model is introduced into D3Q25.Under the finite volume scheme,several typical compressible calculation examples are used to verify whether the numerical stability of the D3Q25 model can be improved by adding the shifted lattice model.The simulation results show that the numerical stability is indeed improved after adding the shifted lattice model.
文摘The standard lattice Boltzmann method utilizes uniform grids to maintain a compact computational procedure. However, it is often less efficient to perform hydrodynamic and aerodynamic flow simulations when there is a need for high resolution. To resolve these difficulties, a multiple nested lattice Boltzmann method(MNLBM) was developed, which contains several overlapped layers with different resolutions in the computational domain. The data transference of flow field on two layers is accomplished by a Filippova procedure which is proved to satisfy the continuity of mass, momentum, and stresses across the interface. The proposed method is based on the standard lattice Boltzmann method, so it is easily performed.By numerical investigation, the result of present method has been agreed with that of literature, but the computation efficiency is higher than the standard lattice Boltzmann method.
文摘The lattice Boltzmann method (LBM) is used to examine free convection of nanofluids. The space between the cold outer square and heated inner circular cylinders is filled with water including various kinds of nanoparticles: TiO2, Ag, Cu, and A1203. The Brinkman and Maxwell-Garnetts models are used to simulate the viscosity and the effective thermal conductivity of nanofluids, respectively. Results from the performed numerical analysis show good agreement with those obtained from other numerical meth- ods. A variety of the Rayleigh number, the nanoparticle volume fraction, and the aspect ratio are examined. According to the results, choosing copper as the nanoparticle leads to obtaining the highest enhancement for this problem. The results also indicate that the maximum value of enhancement occurs at λ =2.5 when Ra = 106 while at A = 1.5 for other Rayleigh numbers.
基金Project supported by the National Natural Science Foundation of China (Grant No 70271069).
文摘By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.
基金sponsored by the project No.50404001 from the National Natural Science Foundation of Chinathe National Key Fundamental Research & Development Project(Grant No.2007CB209601)+1 种基金the China National PetroleumCorporation Fundamental Research Program (Grant No.06A30102)the China Postdoctoral Science Foundation(Project No.2004035350)
文摘The electrophysical property of saturated rocks is very important for reservoir identification and evaluation. In this paper, the lattice Boltzmann method (LBM) was used to study the electrophysical property of rock saturated with fluid because of its advantages over conventional numerical approaches in handling complex pore geometry and boundary conditions. The digital core model was constructed through the accumulation of matrix grains based on their radius distribution obtained by the measurements of core samples. The flow of electrical current through the core model saturated with oil and water was simulated on the mesoscopic scale to reveal the non-Archie relationship between resistivity index and water saturation (I-Sw). The results from LBM simulation and laboratory measurements demonstrated that the I-Sw relation in the range of low water saturation was generally not a straight line in the log-log coordinates as described by the Archie equation. We thus developed a new equation based on numerical simulation and physical experiments. This new equation was used to fit the data from laboratory core measurements and previously published data. Determination of fluid saturation and reservoir evaluation could be significantly improved by using the new equation.
基金the National Natural Science Foundation of China(10772097)
文摘The upper limit of the solid scatters density ns (x), a key parameter for the simulation of flows in porous media with a gray Lattice Boltzmann Method, is studied by an analytical way for the infiltration Poiseuille flow between two infinite parallel plates. Analyses of three different gray Lattice Boltzmann schemes, separately proposed by Gao and Sharma et al., Dardis and McCloskey, and Thorne and Sukop, indicate that the effective domain of Gao and Sharma's scheme is restricted to ns 〈 1/2√3≈0.289, Dardis and McCloskey's scheme is restricted to ns 〈 (√57-1)/28≈0.234, and that there is no extra restriction on ns(x) with Thorne and Sukop's scheme. These results are obtained for the dimensionless relaxation time τ= 1. The above analytical results are verified by our numerical simulations. The use of a gray LBM is further illustrated by simulating the flow at the interface of a porous medium. Simulation results yield velocity profiles which agree very well with Brinkman's prediction.
基金support from the Natural Science Foundation of China(Grant Nos.11272048,51239006 and 11572178)the Tsinghua University Initiative Scientific Research Program
文摘This paper describes the application of a three-dimensional lattice Boltzmann method (LBM) to Newtonian and non-Newtonian (Bingham fluid in this work) flows with free surfaces. A mass tracking algorithm was incorporated to capture the free surface, whereas Papanastasiou's modified model was used for Bingham fluids. The lattice Boltzmann method was first validated using two benchmarks: Newtonian flow through a square cross-section tube and Bingham flow through a circular cross-section tube. Afterward, the dam-break problem for the Newtonian fluid and the slump test for Bingham fluid were simulated to validate the free-surface-capturing algorithm. The numerical results were in good agreement with analytical results, as well as other simulations, thereby proving the validity and correctness of the current method. The proposed method is a promising substitute for time-consuming and costly physical experiments to solve problems encountered in geotechnical and geological engineering, such as the surge and debris flow induced by a landslide or earthquake.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51101035,51371051,and 51306037)
文摘In the present study,the process of droplet condensation on superhydrophobic nanoarrays is simulated using a multicomponent multi-phase lattice Boltzmann model.The results indicate that three typical nucleation modes of condensate droplets are produced by changing the geometrical parameters of nanoarrays.Droplets nucleated at the top(top-nucleation mode),or in the upside interpillar space of nanoarrays(side-nucleation mode),generate the non-wetting Cassie state,whereas the ones nucleated at the bottom corners between the nanoarrays(bottom-nucleation mode) present the wetting Wenzel state.Time evolutions of droplet pressures at the upside and downside of the liquid phase are analyzed to understand the wetting behaviors of the droplets condensed from different nucleation modes.The phenomena of droplet condensation on nanoarrays patterned with different hydrophilic and hydrophobic regions are simulated,indicating that the nucleation mode of condensate droplets can also be manipulated by modifying the local intrinsic wettability of nanoarray surface.The simulation results are compared well with the experimental observations reported in the literature.
基金Project supported by the National Natural Science Foundation of China(Grant No10274006)Education Ministry of China(Grant No03011)
文摘The three-dimensional (3D) lattice Boltzmann models, 3DQ15, 3DQ19 and 3DQ27, under different wall boundary conditions and lattice resolutions have been investigated by simulating Poiseuille flow in a circular cylinder for a wide range of Reynolds numbers. The 3DQ19 model with improved Fillippova and Hanel (FH) curved boundary condition represents a good compromise between computational efficiency and reliability. Blood flow in an aortic arch is then simulated as a typical haemodynamic application. Axial and secondary fluid velocity and effective wall shear stress profiles in a 180° bend are obtained, and the results also demonstrate that the lattice Boltzmann method is suitable for simulating the flow in 3D large-curved vessels.
文摘The numerical investigation of the two-dimensional laminar flow past two ro- tating circular cylinders in the tandem arrangement is conducted by the lattice Boltzmann method. The numerical strategy is used for dealing with curved and moving boundaries of the second-order accuracy for velocity and temperature fields. The effects of various rotational speed ratios and gap spacing are studied with the Reynolds number of 100 and the Prandtl number of 0.71. A varied range of rotational speed ratios are investigated for four different gap spacing, i.e., 3.0, 1.5, 0.7, and 0.2. The results show that, for the first cylinder, the lift and drag coefficients for large gap spacing are similar to those for a single cylinder; for the second cylinder, the lift coefficient descends with the increase in the angular velocity for all gap spacing, while the drag coefficient ascends except for the gap spacing of 3.0. The results of the averaged periodic Nusselt number on the surface of the cylinders show that, for small distances between the cylinders and low angular velocities, conduction is a dominant mechanism of heat transfer, but for large distances and high angular velocities, convection is the main mechanism of heat transfer.
文摘In this study,a two-dimensional approach is elaborated to study with the lattice Boltzmann method(LBM)the seepage of water in the pores of a soil.Firstly,the D2Q9 model is selected to account for the discrete velocity distribution of water flow.In particular,impermeability is considered as macroscopic boundary condition for the left and right domain sides,while the upper and lower boundaries are assumed to behave as pressure boundaries controlled by different densities.The micro-boundary conditions are implemented through the standard rebound strategy and a non-equilibrium extrapolation scheme.Matlab is used for the development of the related algorithm.Finally,the influence of porosity,permeability,osmotic pressure and other factors is assessed with regard to seepage characteristics and the ensuing results are compared with Darcy’s law.The computations show that,for fixed initial conditions,the pore structure has a certain influence on the local velocity of seepage,but the overall state is stable,and the average velocity of each layer is the same.The larger the pore passage is,the faster the flow velocity is,and vice versa.For low permeability,the numerical results are consistent with the Darcy's law.The greater the pressure difference between the inlet and outlet of seepage,the greater the seepage rate.The relationship between them is linear(yet in good agreement with Darcy’s law).
文摘One dimensional traffic model was studied with the consideration of stochastic deceleration by using the lattice Boltzmann method. The evolution equation for vehicle density was derived, which has the form of the Burgers equation exhibiting the effect of viscosity corresponding to the deceleration. The simulation with the model shows that the variation of vehicle density in space tends to a periodic one, which implies the existence of kinematic waves in the 1D traffic flow and coincides with theoretical prediction.
基金supported by the Innovation Fund for Aerospace Science and Technology of China(Grant No.2009200066)the Aeronautical Science Fund of China(Grant No.20111453012)
文摘This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time (SRT) lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.