Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuse...Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuses on the L block to investigate the mechanisms influencing efficient fracture propagation and enhanced stimulated reservoir volume(SRV)in fracturing.To explore the mechanisms influencing effective fracture propagation and enhanced SRV,the L block was selected as the research object,with a comprehensive consideration of geological background,reservoir properties,and dynamic production data.By combining the discrete lattice method with numer-ical analysis and true triaxial experimental simulation,the fracture morphology of a single cluster and the propagation patterns of multiple clusters of complex fractures were obtained.Additionally,the optimization of temporary plugging timing and the fracture map under multiple factors were innovatively proposed.Results indicate that greater flow rate and viscosity can effectively overcome the stress shadow effect of the outermost fractures(1st and 6th clusters),increasing the fracture pressure of the single cluster and the equilibrium degree of multiple fracture propagation,thus forming a more complex fracture network.Moreover,when viscosity exceeds 45 pressure concentrates at fracture mPa⋅s,tips,promoting discontinuous propagation and reducing flow resistance.Conversely,increased gangue thickness and spacing between horizontal wells increase the vertical propagation pressure,suppressing fracture growth and reducing central flow velocity.This study provides a multi-cluster fracture propagation map for optimizing volumetric fracturing in coal seams and suggests that the optimal temporary plugging time significantly enhances the SRV.展开更多
A rapid and efficient method for static aeroelastic analysis of a flexible slender wing when considering the structural geometric nonlinearity has been developed in this paper. A non-planar vortex lattice method herei...A rapid and efficient method for static aeroelastic analysis of a flexible slender wing when considering the structural geometric nonlinearity has been developed in this paper. A non-planar vortex lattice method herein is used to compute the non-planar aerodynamics of flexible wings with large deformation. The finite element method is introduced for structural nonlinear statics analysis. The surface spline method is used for structure/aerodynamics coupling. The static aeroelastic characteristics of the wind tunnel model of a flexible wing are studied by the nonlinear method presented, and the nonlinear method is also evaluated by comparing the results with those obtained from two other methods and the wind tunnel test. The results indicate that the traditional linear method of static aeroelastic analysis is not applicable for cases with large deformation because it produces results that are not realistic. However, the nonlinear methodology, which involves combining the structure finite element method with the non-planar vortex lattice method, could be used to solve the aeroelastic deformation with considerable accuracy, which is in fair agreement with the test results. Moreover, the nonlinear finite element method could consider complex structures. The non-planar vortex lattice method has advantages in both the computational accuracy and efficiency. Consequently, the nonlinear method presented is suitable for the rapid and efficient analysis requirements of engineering practice. It could be used in the preliminary stage and also in the detailed stage of aircraft design.展开更多
The Unsteady Vortex Lattice Method(UVLM) is a medium-fidelity aerodynamic tool that has been widely used in aeroelasticity and flight dynamics simulations. The most timeconsuming step is the evaluation of the induced ...The Unsteady Vortex Lattice Method(UVLM) is a medium-fidelity aerodynamic tool that has been widely used in aeroelasticity and flight dynamics simulations. The most timeconsuming step is the evaluation of the induced velocity. Supposing that the number of bound and wake lattices is N and the computational cost is O (N2), we present an OeNT Dipole Panel Fast Multipole Method(DPFMM) for the rapid evaluation of the induced velocity in UVLM. The multipole expansion coefficients of a quadrilateral dipole panel have been derived in spherical coordinates, whose accuracy is the same as that of the Biot-Savart kernel at the same truncation degree P.Two methods(the loosening method and the shrinking method) are proposed and tested for space partitioning volumetric panels. Compared with FMM for vortex filaments(with three harmonics),DPFMM is approximately two times faster for N2 [103,106]. The simulation time of a multirotor(N~104) is reduced from 100 min(with unaccelerated direct solver) to 2 min(with DPFMM).展开更多
In this paper,the liquid–vapor phase separation under viscous shear is investigated by using a pseudopotential central moment lattice Boltzmann method.Physically,the multiphase shear flow is governed by two competing...In this paper,the liquid–vapor phase separation under viscous shear is investigated by using a pseudopotential central moment lattice Boltzmann method.Physically,the multiphase shear flow is governed by two competing mechanisms:surface tension and shear force.It is interesting to find that the liquid tends to form a droplet when the surface tension dominates under conditions of low temperature,shear velocity,and viscosity,and in larger domain size.Otherwise,the liquid tends to form a band if shear force dominates.Moreover,the average density gradient is used as a physical criterion to distinguish the spinodal decomposition and domain growth.Both spatial and temporal changes of density are studied during the phase separation under shear.展开更多
Flow and heat transfer characteristics during droplet impact on hot walls are pivotal for elucidating the mechanisms of spray cooling and exploring pathways for heat transfer enhancement.When the wall temperature exce...Flow and heat transfer characteristics during droplet impact on hot walls are pivotal for elucidating the mechanisms of spray cooling and exploring pathways for heat transfer enhancement.When the wall temperature exceeds the Leidenfrost point,a vapor film forms between the droplet and the wall,rendering the heat transfer process highly complex.Furthermore,for droplet impact on curved walls,the presence of curvature introduces additional factors that modify the spreading behavior of the droplet and necessitate in-depth analysis.Therefore,this work investigates the flow and heat transfer dynamics of droplet impact on hot planes and curved surfaces numerically via a pseudopotential multiple-relaxation-time Lattice Boltzmann model.The results reveal that the maximum spreading factor increases with the Weber number,diameter ratio,and Bond number,and marginally with the contact angle.Moreover,the time required to achieve the maximumspreading factor increaseswith the contact angle.This relationship exhibits a V-shaped trend due to gravitational effects.Furthermore,the total surface heat flux increases with theWeber number but decreases with the contact angle.The results advance the fundamental understanding of droplet impact dynamics on hot curved surfaces,providing practical insights for optimizing spray cooling performance and thermal management systems.展开更多
Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss fu...Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss functions.However,when addressing complex flow problems,PINNs still face some challenges such as activation saturation and vanishing gradients in deep network training,leading to slow convergence and insufficient prediction accuracy.We present physics-informed neural networks incorporating lattice Boltzmann method optimized by tanh robust weight initialization(T-PINN-LBM)to address these challenges.This approach fuses the mesoscopic lattice Boltzmann model with the automatic differentiation framework of PINNs.It also implements a tanh robust weight initialization method derived from fixed point analysis.This model effectively mitigates activation and gradient decay in deep networks,improving convergence speed and data efficiency in multiscale flow simulations.We validate the effectiveness of the model on the classical arithmetic example of lid-driven cavity flow.Compared to the traditional Xavier initialized PINN and PINN-LBM,T-PINNLBM reduces the mean absolute error(MAE)by one order of magnitude at the same network depth and maintains stable convergence in deeper networks.The results demonstrate that this model can accurately capture complex flow structures without prior data,providing a new feasible pathway for data-free driven fluid simulation.展开更多
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.展开更多
This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating som...This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating some indices including the number of failed edges, dynamic edge tolerance capacity and the perturbation of edge. In general, it assumes that the perturbation is attributed to the normal distribution in adopted simulations. By investigating the effectiveness of edge tolerance in scale-free and random networks, it finds that the larger tolerance parameter λ can more efficiently delay the cascading failure process for scale-free networks than random networks. These results indicate that the cascading failure process can be effectively controlled by increasing the tolerance parameter λ. Moreover, the simulations also show that, larger variance of perturbation can easily trigger the cascading failures than the smaller one. This study may be useful for evaluating efficiency of whole traffic systems, and for alleviating cascading failure in such systems.展开更多
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.展开更多
One of the main issues concerning the standard Vortex Lattice Method is its application to partially or fully detached flow conditions,where non-linear aerodynamic characteristics appear as the angle of attack increas...One of the main issues concerning the standard Vortex Lattice Method is its application to partially or fully detached flow conditions,where non-linear aerodynamic characteristics appear as the angle of attack increases and/or the aspect ratio decreases.In order to solve such limitations,a pure numerical approach based entirely on the Vortex Lattice Method concepts has been developed.The so-called steady“Full Multi-wake Vortex Lattice Method”comes from the main hypothesis that each discretized element on the body’s surface detaches their own wakes downstream.The obtained results match for lift,drag and moment coefficients for the entire aspect ratio range configurations(under straight wakes and inviscid assumptions).Future unsteady versions of such a multi-wake approach could improve the current results obtained through Vortex Element Methods(as vortons or isolated vortex filaments).展开更多
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.展开更多
Polymer electrolyte fuel cells(PEFCs)being employed in fuel cell electric vehicles(FCEVs)are promising power generators producing electric power from fuel stream via porous electrodes.Structure of carbon paper gas dif...Polymer electrolyte fuel cells(PEFCs)being employed in fuel cell electric vehicles(FCEVs)are promising power generators producing electric power from fuel stream via porous electrodes.Structure of carbon paper gas diffusion layers(GDLs)applying in the porous electrodes can greatly affect the PEFC performance,especially at the cathode side where electrochemical reaction is more sluggish.To discover the role of carbon paper GDL structure on the mass transfer properties,different cathode electrodes with dissimilar structural parameters are simulated via lattice Boltzmann method(LBM).3D contours of oxygen and water vapor concentration through the GDL as well as the 2D contours of current density on the catalyst layer are illustrated and examined.The results indicate that the carbon fiber diameter has a negligible impact on the current density while the impact of carbon paper thickness and porosity is significant.In fact,increasing of carbon paper thickness or porosity leads to lack of cell performance.展开更多
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...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 practical fluid dynamic simulations,the bou n dary condition should be treated carefully because it always has crucial influen ce on the numerical accuracy,stability and efficiency.Two types of boundary tr eatment ...In practical fluid dynamic simulations,the bou n dary condition should be treated carefully because it always has crucial influen ce on the numerical accuracy,stability and efficiency.Two types of boundary tr eatment methods for lattice Boltzmann method(LBM)are proposed.One is for the treatment of boundaries situated at lattice nodes,and the other is for the appr oximation of boundaries that are not located at the regular lattice nodes.The f irst type of boundary treatment method can deal with various dynamic boundaries on complex geometries by using a general set of formulas,which can maintain sec ond-order accuracy.Based on the fact that the fluid flows simulated by LBM are not far from equilibrium,the unknown distributions at a boundary node are expr essed as the analogous forms of their corresponding equilibrium distributions.T herefore,the number of unknowns can be reduced and an always-closed set of equ ations can be obtained for the solutions to pressure,velocity and special bound ary conditions on various geometries.The second type of boundary treatment is a complete interpolation scheme to treat curved boundaries.It comes from careful analysis of the relations between distribution functions at boundary nodes and their neighboring lattice nodes.It is stable for all situations and of second-order accuracy.Basic ideas,implementation procedures and verifications with ty pical examples for the both treatments are presented.Numerical simulations and analyses show that they are accurate,stable,general and efficient for practica l simulations.展开更多
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 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.展开更多
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.展开更多
This work presents a numerical study on the dynamic high velocity compaction of the metal powder. The analysis of the process is based on a mesoscopic approach using multi-speed lattice Boltzmann method. The boundary ...This work presents a numerical study on the dynamic high velocity compaction of the metal powder. The analysis of the process is based on a mesoscopic approach using multi-speed lattice Boltzmann method. The boundary condition and the relaxation time are tailored to the situation. The dynamic compaction process is vividly presented and the shock wave can be easily found in the simulation. The density is analyzed in order to explore the mechanism of the high velocity compaction.展开更多
基金the project of the State Key Laboratory of Petroleum Resources and Engineering(No.PRE/open-2307)the CNOOC Research Institute(No.2020PFS-03).
文摘Hydraulic fracturing,an effective method for enhancing coal seam productivity,largely determines coalbed methane(CBM)production,which is significantly influenced by geological and engineering factors.This study focuses on the L block to investigate the mechanisms influencing efficient fracture propagation and enhanced stimulated reservoir volume(SRV)in fracturing.To explore the mechanisms influencing effective fracture propagation and enhanced SRV,the L block was selected as the research object,with a comprehensive consideration of geological background,reservoir properties,and dynamic production data.By combining the discrete lattice method with numer-ical analysis and true triaxial experimental simulation,the fracture morphology of a single cluster and the propagation patterns of multiple clusters of complex fractures were obtained.Additionally,the optimization of temporary plugging timing and the fracture map under multiple factors were innovatively proposed.Results indicate that greater flow rate and viscosity can effectively overcome the stress shadow effect of the outermost fractures(1st and 6th clusters),increasing the fracture pressure of the single cluster and the equilibrium degree of multiple fracture propagation,thus forming a more complex fracture network.Moreover,when viscosity exceeds 45 pressure concentrates at fracture mPa⋅s,tips,promoting discontinuous propagation and reducing flow resistance.Conversely,increased gangue thickness and spacing between horizontal wells increase the vertical propagation pressure,suppressing fracture growth and reducing central flow velocity.This study provides a multi-cluster fracture propagation map for optimizing volumetric fracturing in coal seams and suggests that the optimal temporary plugging time significantly enhances the SRV.
基金National Natural Science Foundation of China(Nos.11172025,91116005)Research Fund for the Doctoral Program of Higher Education of China(No.20091102110015)
文摘A rapid and efficient method for static aeroelastic analysis of a flexible slender wing when considering the structural geometric nonlinearity has been developed in this paper. A non-planar vortex lattice method herein is used to compute the non-planar aerodynamics of flexible wings with large deformation. The finite element method is introduced for structural nonlinear statics analysis. The surface spline method is used for structure/aerodynamics coupling. The static aeroelastic characteristics of the wind tunnel model of a flexible wing are studied by the nonlinear method presented, and the nonlinear method is also evaluated by comparing the results with those obtained from two other methods and the wind tunnel test. The results indicate that the traditional linear method of static aeroelastic analysis is not applicable for cases with large deformation because it produces results that are not realistic. However, the nonlinear methodology, which involves combining the structure finite element method with the non-planar vortex lattice method, could be used to solve the aeroelastic deformation with considerable accuracy, which is in fair agreement with the test results. Moreover, the nonlinear finite element method could consider complex structures. The non-planar vortex lattice method has advantages in both the computational accuracy and efficiency. Consequently, the nonlinear method presented is suitable for the rapid and efficient analysis requirements of engineering practice. It could be used in the preliminary stage and also in the detailed stage of aircraft design.
文摘The Unsteady Vortex Lattice Method(UVLM) is a medium-fidelity aerodynamic tool that has been widely used in aeroelasticity and flight dynamics simulations. The most timeconsuming step is the evaluation of the induced velocity. Supposing that the number of bound and wake lattices is N and the computational cost is O (N2), we present an OeNT Dipole Panel Fast Multipole Method(DPFMM) for the rapid evaluation of the induced velocity in UVLM. The multipole expansion coefficients of a quadrilateral dipole panel have been derived in spherical coordinates, whose accuracy is the same as that of the Biot-Savart kernel at the same truncation degree P.Two methods(the loosening method and the shrinking method) are proposed and tested for space partitioning volumetric panels. Compared with FMM for vortex filaments(with three harmonics),DPFMM is approximately two times faster for N2 [103,106]. The simulation time of a multirotor(N~104) is reduced from 100 min(with unaccelerated direct solver) to 2 min(with DPFMM).
基金supported by National Natural Science Foundation of China under Grant No.51806116Guangdong Basic and Applied Basic Research Foundation under Grant No.2024A1515010927+2 种基金China Scholarship Council under Grant No.202306380288Humanities and Social Science Foundation of the Ministry of Education in China under Grant No.24YJCZH163Fundamental Research Funds for the Central Universities,Sun Yat-sen University under Grant No.24qnpy044。
文摘In this paper,the liquid–vapor phase separation under viscous shear is investigated by using a pseudopotential central moment lattice Boltzmann method.Physically,the multiphase shear flow is governed by two competing mechanisms:surface tension and shear force.It is interesting to find that the liquid tends to form a droplet when the surface tension dominates under conditions of low temperature,shear velocity,and viscosity,and in larger domain size.Otherwise,the liquid tends to form a band if shear force dominates.Moreover,the average density gradient is used as a physical criterion to distinguish the spinodal decomposition and domain growth.Both spatial and temporal changes of density are studied during the phase separation under shear.
基金funded by the National Natural Science Foundation of China,grant numbers 52276055 and 52476065.
文摘Flow and heat transfer characteristics during droplet impact on hot walls are pivotal for elucidating the mechanisms of spray cooling and exploring pathways for heat transfer enhancement.When the wall temperature exceeds the Leidenfrost point,a vapor film forms between the droplet and the wall,rendering the heat transfer process highly complex.Furthermore,for droplet impact on curved walls,the presence of curvature introduces additional factors that modify the spreading behavior of the droplet and necessitate in-depth analysis.Therefore,this work investigates the flow and heat transfer dynamics of droplet impact on hot planes and curved surfaces numerically via a pseudopotential multiple-relaxation-time Lattice Boltzmann model.The results reveal that the maximum spreading factor increases with the Weber number,diameter ratio,and Bond number,and marginally with the contact angle.Moreover,the time required to achieve the maximumspreading factor increaseswith the contact angle.This relationship exhibits a V-shaped trend due to gravitational effects.Furthermore,the total surface heat flux increases with theWeber number but decreases with the contact angle.The results advance the fundamental understanding of droplet impact dynamics on hot curved surfaces,providing practical insights for optimizing spray cooling performance and thermal management systems.
基金supported by the National Natural Science Foundation of China(Grant Nos.12474453,12174085,and 12404530).
文摘Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss functions.However,when addressing complex flow problems,PINNs still face some challenges such as activation saturation and vanishing gradients in deep network training,leading to slow convergence and insufficient prediction accuracy.We present physics-informed neural networks incorporating lattice Boltzmann method optimized by tanh robust weight initialization(T-PINN-LBM)to address these challenges.This approach fuses the mesoscopic lattice Boltzmann model with the automatic differentiation framework of PINNs.It also implements a tanh robust weight initialization method derived from fixed point analysis.This model effectively mitigates activation and gradient decay in deep networks,improving convergence speed and data efficiency in multiscale flow simulations.We validate the effectiveness of the model on the classical arithmetic example of lid-driven cavity flow.Compared to the traditional Xavier initialized PINN and PINN-LBM,T-PINNLBM reduces the mean absolute error(MAE)by one order of magnitude at the same network depth and maintains stable convergence in deeper networks.The results demonstrate that this model can accurately capture complex flow structures without prior data,providing a new feasible pathway for data-free driven fluid simulation.
基金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 National Basic Research Program of China (Grant No 2006CB705500)Chang-Jiang Scholars and Innovative Research Team in University of China (Grant No IRT0605)the National Natural Science Foundation of China (Grant No70631001)
文摘This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating some indices including the number of failed edges, dynamic edge tolerance capacity and the perturbation of edge. In general, it assumes that the perturbation is attributed to the normal distribution in adopted simulations. By investigating the effectiveness of edge tolerance in scale-free and random networks, it finds that the larger tolerance parameter λ can more efficiently delay the cascading failure process for scale-free networks than random networks. These results indicate that the cascading failure process can be effectively controlled by increasing the tolerance parameter λ. Moreover, the simulations also show that, larger variance of perturbation can easily trigger the cascading failures than the smaller one. This study may be useful for evaluating efficiency of whole traffic systems, and for alleviating cascading failure in such systems.
基金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.
基金the National Council for Science and Technology of Mexico(Consejo Nacional de Ciencia y Tecnología-CONACyT).
文摘One of the main issues concerning the standard Vortex Lattice Method is its application to partially or fully detached flow conditions,where non-linear aerodynamic characteristics appear as the angle of attack increases and/or the aspect ratio decreases.In order to solve such limitations,a pure numerical approach based entirely on the Vortex Lattice Method concepts has been developed.The so-called steady“Full Multi-wake Vortex Lattice Method”comes from the main hypothesis that each discretized element on the body’s surface detaches their own wakes downstream.The obtained results match for lift,drag and moment coefficients for the entire aspect ratio range configurations(under straight wakes and inviscid assumptions).Future unsteady versions of such a multi-wake approach could improve the current results obtained through Vortex Element Methods(as vortons or isolated vortex filaments).
文摘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.
文摘Polymer electrolyte fuel cells(PEFCs)being employed in fuel cell electric vehicles(FCEVs)are promising power generators producing electric power from fuel stream via porous electrodes.Structure of carbon paper gas diffusion layers(GDLs)applying in the porous electrodes can greatly affect the PEFC performance,especially at the cathode side where electrochemical reaction is more sluggish.To discover the role of carbon paper GDL structure on the mass transfer properties,different cathode electrodes with dissimilar structural parameters are simulated via lattice Boltzmann method(LBM).3D contours of oxygen and water vapor concentration through the GDL as well as the 2D contours of current density on the catalyst layer are illustrated and examined.The results indicate that the carbon fiber diameter has a negligible impact on the current density while the impact of carbon paper thickness and porosity is significant.In fact,increasing of carbon paper thickness or porosity leads to lack of cell performance.
基金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.
基金Supported by the National Natural Science Foundation of China(50009007)
文摘In practical fluid dynamic simulations,the bou n dary condition should be treated carefully because it always has crucial influen ce on the numerical accuracy,stability and efficiency.Two types of boundary tr eatment methods for lattice Boltzmann method(LBM)are proposed.One is for the treatment of boundaries situated at lattice nodes,and the other is for the appr oximation of boundaries that are not located at the regular lattice nodes.The f irst type of boundary treatment method can deal with various dynamic boundaries on complex geometries by using a general set of formulas,which can maintain sec ond-order accuracy.Based on the fact that the fluid flows simulated by LBM are not far from equilibrium,the unknown distributions at a boundary node are expr essed as the analogous forms of their corresponding equilibrium distributions.T herefore,the number of unknowns can be reduced and an always-closed set of equ ations can be obtained for the solutions to pressure,velocity and special bound ary conditions on various geometries.The second type of boundary treatment is a complete interpolation scheme to treat curved boundaries.It comes from careful analysis of the relations between distribution functions at boundary nodes and their neighboring lattice nodes.It is stable for all situations and of second-order accuracy.Basic ideas,implementation procedures and verifications with ty pical examples for the both treatments are presented.Numerical simulations and analyses show that they are accurate,stable,general and efficient for practica l simulations.
基金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.
基金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.
基金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.
基金supported by the National Natural Science Foundation of China(Nos. 50874123 and 51174236)National Basic Research Program of China(No. 2011CB606306)
文摘This work presents a numerical study on the dynamic high velocity compaction of the metal powder. The analysis of the process is based on a mesoscopic approach using multi-speed lattice Boltzmann method. The boundary condition and the relaxation time are tailored to the situation. The dynamic compaction process is vividly presented and the shock wave can be easily found in the simulation. The density is analyzed in order to explore the mechanism of the high velocity compaction.
