To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced v...To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.展开更多
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.展开更多
A two-dimensional（2-D） incompressible plane jet is investigated using the lattice Boltzmann method（LBM） for low Reynolds numbers of 42 and 65 based on the jet-exit-width and the maximum jet-exit-velocity. The resu...A two-dimensional（2-D） incompressible plane jet is investigated using the lattice Boltzmann method（LBM） for low Reynolds numbers of 42 and 65 based on the jet-exit-width and the maximum jet-exit-velocity. The results show that the mean centerline velocity decays as x-1/3 and the jet spreads as x2/3 in the self-similar region, which are consistent with the theoretical predictions and the experimental data. The time histories and PSD analyses of the instantaneous centerline velocities indicate the periodic behavior and the interaction between periodic components of velocities should not be neglected in the far field region, although it is invisible in the near field region.展开更多
The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S）...The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.展开更多
The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, ...The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.展开更多
In this paper, a novel model is proposed to investigate the neutron transport in scattering and absorbing medium. This solution to the linear Boltzmann equation is expanded from the idea of lattice Boltzmann method(LB...In this paper, a novel model is proposed to investigate the neutron transport in scattering and absorbing medium. This solution to the linear Boltzmann equation is expanded from the idea of lattice Boltzmann method(LBM) with the collision and streaming process. The theoretical derivation of lattice Boltzmann model for transient neutron transport problem is proposed for the first time.The fully implicit backward difference scheme is used to ensure the numerical stability, and relaxation time and equilibrium particle distribution function are obtained. To validate the new lattice Boltzmann model, the LBM formulation is tested for a homogenous media with different sources, and both transient and steady-state LBM results get a good agreement with the benchmark solutions.展开更多
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.展开更多
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 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 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.展开更多
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.展开更多
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.展开更多
This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential...This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential momentum accommodation coefficient on the rough boundary condition, the lattice Boltzmann simulations of nitrogen and helium flows are performed in a two-dimensional microchannel with rough boundaries. The surface roughness effects in the microchannel on the velocity field, the mass flow rate and the friction coefficient are studied and analysed. Numerical results for the two gases in micro scale show different characteristics from macroscopic flows and demonstrate the feasibility of the lattice Boltzmann model in rarefied gas dynamics.展开更多
Active boundary layer flow control and boundary layer manipulation in the channel flow that was based on low temperature plasma were studied by means of a lattice Boltzmann method. Two plasma actuators were placed in ...Active boundary layer flow control and boundary layer manipulation in the channel flow that was based on low temperature plasma were studied by means of a lattice Boltzmann method. Two plasma actuators were placed in a row to obtain the influence rule of their separation distance on the velocity profile at three locations and maximum velocity in the flow field. Two plasma actuators were placed symmetrically inside a channel to examine the effect of channel height and voltage on the velocity profile and flow rate. It was found that the channel height controls the distribution of flow velocity, which affected the flow rate and its direction. Increasing plasma voltage had a negative effect on the flow rate due to the generation of a larger and stronger flow vortex.展开更多
For comprehensive characteristics of flow in a gas bearing,lattice Boltzmann method(LBM)is applied for study of the two-dimensional flow between two eccentric cylinders with the inner one rotating at a high speed.The ...For comprehensive characteristics of flow in a gas bearing,lattice Boltzmann method(LBM)is applied for study of the two-dimensional flow between two eccentric cylinders with the inner one rotating at a high speed.The flow pattern and circumferential pressure distribution are discussed based on critical issues such as eccentricity ranging from 0.2 to 0.9,clearance ratio varying from 0.005 to 0.01 and rotating speed in the range of 3×104—1.8×105 r/min.The analysis and discussion on the circumferential pressure distribution affirmed the quasilinear relation between the extremum pressure and rotating speed.Furthermore,a high eccentricity and small clearance ratio contributes most to the fluctuation of the circumferential pressure distribution.The flow pattern inside the channel exhibits separation vortex under a large eccentricity.The conclusions drawn in this work give rise to prediction of the flow pattern in the gas bearing which is beneficial for evaluating the performance of as well as instructing the design and development.展开更多
A new method for selecting dimensionless relaxation time in the lattice Boltzmann model was proposed based on similarity criterion and gas true physical parameters.At the same time,the dimensionless relaxation time wa...A new method for selecting dimensionless relaxation time in the lattice Boltzmann model was proposed based on similarity criterion and gas true physical parameters.At the same time,the dimensionless relaxation time was modified by considering the influence of the boundary Knudsen layer.On this basis,the second-order slip boundary condition of the wall was considered,and the key parameters in the corresponding combined bounce-back/specular-reflection boundary condition were deduced to build a new model of unconventional gas microscale flow simulation based on the lattice Boltzmann method suitable for high temperatures and high pressures.The simulation results of methane gas flow driven by body force in infinite micro-channels and flow driven by inlet-outlet pressure differential in long straight channels were compared with the numerical and analytical solutions in the literature to verify the accuracy of the model,and the dimensionless relaxation time modification was formally optimized.The results show that the new model can effectively characterize the slippage effect,compression effect,gas density and the effect of boundary Knudsen layer in the micro-scale flow of unconventional natural gas.The new model can achieve a more comprehensive characterization of the real gas flow conditions and can be used as a basic model for the simulation of unconventional gas flow on the micro-nano scale.展开更多
The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass cons...The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.展开更多
The thermal lattice Boltzmann method （TLBM）, which was proposed by J. G. M. Eggels and J. A. Somers previously, has been improved in this paper. The improved method has introduced a new equilibrium solution for the ...The thermal lattice Boltzmann method （TLBM）, which was proposed by J. G. M. Eggels and J. A. Somers previously, has been improved in this paper. The improved method has introduced a new equilibrium solution for the temperature distribution function on the assumption that flow is incompressible, and it can correct the effect of compressibility on the macroscopic temperature computed. Compared to the previous method, where the half- way bounce back boundary condition was used for non-slip velocity and temperature, a non-equilibrium extrapolation scheme has been adopted for both velocity and temperature boundary conditions in this paper. Its second-order accuracy coincides with the ensemble accuracy of lattice Boltzmann method. In order to validate the improved thermal scheme, the natural convection of air in a square cavity is simulated by using this method. The results obtained in the simulation agree very well with the data of other numerical methods and benchmark data. It is indicated that the improved TLBM is also successful for the simulations of non-isothermal flows. Moreover, this thermal scheme can be applied to simulate the natural convection in a non-uniform high magnetic field. The simulation has been completed in a square cavity filled with the aqueous solutions of KC1 （llwt%）, which is considered as a diamagnetic fluid with electrically low-conducting, with Grashof number Gr=4.64~104 and Prandtl number Pr----7.0. And three cases, with different cavity locations in the magnetic field, have been studied. In the presence of a high magnetic field, the natural convection is quenched by the body forces exerted on the electrically low-conducting fluids, such as the magnetization force and the Lorentz force. From the results obtained, it can be seen that the quenching efficiencies decrease with the variation of location from left, symmetrical line, to the right. These phenomena originate from the different distributions of the magnetic field strengths in the zones of the symmetrical central line of the magnetic fields. The results are also compared with those without a magnetic field. Finally, we can conclude that the improved TLBM will enable effective simulation of the natural convection under a high magnetic field.展开更多
We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eu...We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.展开更多
The central problem of the lattice Boltzmann method （LBM） is to construct a discrete equilibrium. In this paper, a multi-speed 1D cell-model of Boltzmann equation is proposed, in which the cell-population equilibriu...The central problem of the lattice Boltzmann method （LBM） is to construct a discrete equilibrium. In this paper, a multi-speed 1D cell-model of Boltzmann equation is proposed, in which the cell-population equilibrium, a direct non- negative approximation to the continuous Maxwellian distribution, plays an important part. By applying the explicit one-order Chapman-Enskog distribution, the model reduces the transportation and collision, two basic evolution steps in LBM, to the transportation of the non-equilibrium distribution. Furthermore, 1D dam-break problem is performed and the numerical results agree well with the analytic solutions.展开更多
基金Supported by the National Natural Science Foundation of China(52201350,52201394,and 52271301)the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai)(Grant No.SML2022008).
文摘To explore the relationship between dynamic characteristics and wake patterns,numerical simulations were conducted on three equal-diameter cylinders arranged in an equilateral triangle.The simulations varied reduced velocities and gap spacing to observe flow-induced vibrations(FIVs).The immersed boundary–lattice Boltzmann flux solver(IB–LBFS)was applied as a numerical solution method,allowing for straightforward application on a simple Cartesian mesh.The accuracy and rationality of this method have been verified through comparisons with previous numerical results,including studies on flow past three stationary circular cylinders arranged in a similar pattern and vortex-induced vibrations of a single cylinder across different reduced velocities.When examining the FIVs of three cylinders,numerical simulations were carried out across a range of reduced velocities(3.0≤Ur≤13.0)and gap spacing(L=3D,4D,and 5D).The observed vibration response included several regimes:the desynchronization regime,the initial branch,and the lower branch.Notably,the transverse amplitude peaked,and a double vortex street formed in the wake when the reduced velocity reached the lower branch.This arrangement of three cylinders proved advantageous for energy capture as the upstream cylinder’s vibration response mirrored that of an isolated cylinder,while the response of each downstream cylinder was significantly enhanced.Compared to a single cylinder,the vibration and flow characteristics of this system are markedly more complex.The maximum transverse amplitudes of the downstream cylinders are nearly identical and exceed those observed in a single-cylinder set-up.Depending on the gap spacing,the flow pattern varied:it was in-phase for L=3D,antiphase for L=4D,and exhibited vortex shedding for L=5D.The wake configuration mainly featured double vortex streets for L=3D and evolved into two pairs of double vortex streets for L=5D.Consequently,it well illustrates the coupling mechanism that dynamics characteristics and wake vortex change with gap spacing and reduced velocities.
文摘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 National Nature Science Foundation of China(10472046)the Scientific Innova-tion Research of College Graduate in Jiangsu Province(CX08B-035Z)the Innovation and Excellence Foundation of Doctoral Dissertation of Nanjing University of Aeronautics and Astronautics(BCXJ08-01)~~
文摘A two-dimensional（2-D） incompressible plane jet is investigated using the lattice Boltzmann method（LBM） for low Reynolds numbers of 42 and 65 based on the jet-exit-width and the maximum jet-exit-velocity. The results show that the mean centerline velocity decays as x-1/3 and the jet spreads as x2/3 in the self-similar region, which are consistent with the theoretical predictions and the experimental data. The time histories and PSD analyses of the instantaneous centerline velocities indicate the periodic behavior and the interaction between periodic components of velocities should not be neglected in the far field region, although it is invisible in the near field region.
基金Supported by the Aeronautical Science Foundation of China(20111453012)the National Defense Pre-Research Foundation of China(9140A13040111HK0329)~~
文摘The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.
基金Supported by the Aeronautical Science Foundation of China(20061453020)Foundation for Basic Research of Northwestern Polytechnical University(03)~~
文摘The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.
基金supported by the Foundation of National Key Laboratory of Reactor System Design Technology(No.HT-LW-02-2014003)the State Key Program of National Natural Science of China(No.51436009)
文摘In this paper, a novel model is proposed to investigate the neutron transport in scattering and absorbing medium. This solution to the linear Boltzmann equation is expanded from the idea of lattice Boltzmann method(LBM) with the collision and streaming process. The theoretical derivation of lattice Boltzmann model for transient neutron transport problem is proposed for the first time.The fully implicit backward difference scheme is used to ensure the numerical stability, and relaxation time and equilibrium particle distribution function are obtained. To validate the new lattice Boltzmann model, the LBM formulation is tested for a homogenous media with different sources, and both transient and steady-state LBM results get a good agreement with the benchmark solutions.
基金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.
基金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.
基金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.
基金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.
基金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.
基金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.
基金Project supported by the National Natural Science Foundation of China and NSAF (Grant No 10576010)the Creation Foundation of Fudan University (Grant No 2126003)
文摘This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential momentum accommodation coefficient on the rough boundary condition, the lattice Boltzmann simulations of nitrogen and helium flows are performed in a two-dimensional microchannel with rough boundaries. The surface roughness effects in the microchannel on the velocity field, the mass flow rate and the friction coefficient are studied and analysed. Numerical results for the two gases in micro scale show different characteristics from macroscopic flows and demonstrate the feasibility of the lattice Boltzmann model in rarefied gas dynamics.
文摘Active boundary layer flow control and boundary layer manipulation in the channel flow that was based on low temperature plasma were studied by means of a lattice Boltzmann method. Two plasma actuators were placed in a row to obtain the influence rule of their separation distance on the velocity profile at three locations and maximum velocity in the flow field. Two plasma actuators were placed symmetrically inside a channel to examine the effect of channel height and voltage on the velocity profile and flow rate. It was found that the channel height controls the distribution of flow velocity, which affected the flow rate and its direction. Increasing plasma voltage had a negative effect on the flow rate due to the generation of a larger and stronger flow vortex.
基金partially supported by the Aeronautical Science Foundation of China (No.201928052008)
文摘For comprehensive characteristics of flow in a gas bearing,lattice Boltzmann method(LBM)is applied for study of the two-dimensional flow between two eccentric cylinders with the inner one rotating at a high speed.The flow pattern and circumferential pressure distribution are discussed based on critical issues such as eccentricity ranging from 0.2 to 0.9,clearance ratio varying from 0.005 to 0.01 and rotating speed in the range of 3×104—1.8×105 r/min.The analysis and discussion on the circumferential pressure distribution affirmed the quasilinear relation between the extremum pressure and rotating speed.Furthermore,a high eccentricity and small clearance ratio contributes most to the fluctuation of the circumferential pressure distribution.The flow pattern inside the channel exhibits separation vortex under a large eccentricity.The conclusions drawn in this work give rise to prediction of the flow pattern in the gas bearing which is beneficial for evaluating the performance of as well as instructing the design and development.
基金Supported by National Natural Science Foundation of China(Key Program)(51534006)National Natural Science Foundation of China(51874251)。
文摘A new method for selecting dimensionless relaxation time in the lattice Boltzmann model was proposed based on similarity criterion and gas true physical parameters.At the same time,the dimensionless relaxation time was modified by considering the influence of the boundary Knudsen layer.On this basis,the second-order slip boundary condition of the wall was considered,and the key parameters in the corresponding combined bounce-back/specular-reflection boundary condition were deduced to build a new model of unconventional gas microscale flow simulation based on the lattice Boltzmann method suitable for high temperatures and high pressures.The simulation results of methane gas flow driven by body force in infinite micro-channels and flow driven by inlet-outlet pressure differential in long straight channels were compared with the numerical and analytical solutions in the literature to verify the accuracy of the model,and the dimensionless relaxation time modification was formally optimized.The results show that the new model can effectively characterize the slippage effect,compression effect,gas density and the effect of boundary Knudsen layer in the micro-scale flow of unconventional natural gas.The new model can achieve a more comprehensive characterization of the real gas flow conditions and can be used as a basic model for the simulation of unconventional gas flow on the micro-nano scale.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11862003 and 81860635)the Key Project of the Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2017GXNSFDA198038)+3 种基金the Project of Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2018GXNSFAA281302)the Project for Promotion of Young and Middle-aged Teachers’Basic Scientific Research Ability in Guangxi Universities,China(Grant No.2019KY0084)the“Bagui Scholar”Teams for Innovation and Research Project of Guangxi Zhuang Autonomous Region,Chinathe Graduate Innovation Program of Guangxi Normal University,China(Grant No.JXYJSKT-2019-007)。
文摘The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.
基金Project supported by the National Natural Science Foundation of China (Grant No 10772150)the Aeronautical Science Fund of China (Grant No 20061453020)Foundation for Basic Research of Northwestern Polytechnical University
文摘The thermal lattice Boltzmann method （TLBM）, which was proposed by J. G. M. Eggels and J. A. Somers previously, has been improved in this paper. The improved method has introduced a new equilibrium solution for the temperature distribution function on the assumption that flow is incompressible, and it can correct the effect of compressibility on the macroscopic temperature computed. Compared to the previous method, where the half- way bounce back boundary condition was used for non-slip velocity and temperature, a non-equilibrium extrapolation scheme has been adopted for both velocity and temperature boundary conditions in this paper. Its second-order accuracy coincides with the ensemble accuracy of lattice Boltzmann method. In order to validate the improved thermal scheme, the natural convection of air in a square cavity is simulated by using this method. The results obtained in the simulation agree very well with the data of other numerical methods and benchmark data. It is indicated that the improved TLBM is also successful for the simulations of non-isothermal flows. Moreover, this thermal scheme can be applied to simulate the natural convection in a non-uniform high magnetic field. The simulation has been completed in a square cavity filled with the aqueous solutions of KC1 （llwt%）, which is considered as a diamagnetic fluid with electrically low-conducting, with Grashof number Gr=4.64~104 and Prandtl number Pr----7.0. And three cases, with different cavity locations in the magnetic field, have been studied. In the presence of a high magnetic field, the natural convection is quenched by the body forces exerted on the electrically low-conducting fluids, such as the magnetization force and the Lorentz force. From the results obtained, it can be seen that the quenching efficiencies decrease with the variation of location from left, symmetrical line, to the right. These phenomena originate from the different distributions of the magnetic field strengths in the zones of the symmetrical central line of the magnetic fields. The results are also compared with those without a magnetic field. Finally, we can conclude that the improved TLBM will enable effective simulation of the natural convection under a high magnetic field.
基金supported by the National Natural Science Foundation of China(Grant No.11572062)the Fundamental Research Funds for the Central Universities,China(Grant No.CDJZR13248801)+2 种基金the Program for Changjiang Scholars and Innovative Research Team in University,China(Grant No.IRT13043)Key Laboratory of Functional Crystals and Laser Technology,TIPCChinese Academy of Sciences
文摘We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop, and the drop deformation in a shear flow, were simulated to validate the present model.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 70271069 and 60773195)
文摘The central problem of the lattice Boltzmann method （LBM） is to construct a discrete equilibrium. In this paper, a multi-speed 1D cell-model of Boltzmann equation is proposed, in which the cell-population equilibrium, a direct non- negative approximation to the continuous Maxwellian distribution, plays an important part. By applying the explicit one-order Chapman-Enskog distribution, the model reduces the transportation and collision, two basic evolution steps in LBM, to the transportation of the non-equilibrium distribution. Furthermore, 1D dam-break problem is performed and the numerical results agree well with the analytic solutions.
