To study the development of imbibition such as the imbibition front and phase distribution in shale,the Lattice Boltzmann Method(LBM)is used to study the imbibition processes in the pore-throat network of shale.Throug...To study the development of imbibition such as the imbibition front and phase distribution in shale,the Lattice Boltzmann Method(LBM)is used to study the imbibition processes in the pore-throat network of shale.Through dimensional analysis,four dimensionless parameters affecting the imbibition process were determined.A color gradient model of LBM was used in computation based on a real core pore size distribution.The numerical results show that the four factors have great effects on imbibition.The impact of each factor is not monotonous.The imbibition process is the comprehensive effect of all aspects.The imbibition front becomes more and more non-uniform with time in a heterogeneous pore-throat network.Some non-wetting phases(oil here)cannot be displaced out.The displacement efficiency and velocity do not change monotonously with any factor.The development of the average imbibition length with time is not smooth and not linear in a heterogeneous pore-throat network.Two fitting relations between the four dimensionless parameters and the imbibition velocity and efficiency are obtained,respectively.展开更多
The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that o...The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that of conventional reservoirs.At the early stage of flowback,there is no single-phase flow of the liquid phase in shale,but rather a gas-water two-phase flow,such that the single-phase flow model for tight oil and gas reservoirs is not applicable.In this study,pores and microfractures are extracted based on the experimental results of computed tomography(CT)scanning,and a spatial model of microfractures is established.Then,the influence of rough microfracture surfaces on the flow is corrected using the modified cubic law,which was modified by introducing the average deviation of the microfracture height as a roughness factor to consider the influence of microfracture surface roughness.The flow in the fracture network is simulated using the modified cubic law and the lattice Boltzmann method(LBM).The results obtained demonstrate that most of the fracturing fluid is retained in the shale microfractures,which explains the low fracturing fluid flowback rate in shale hydraulic fracturing.展开更多
The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive comp...The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate parti...Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate particle-fluid interaction problems involving heat transfer at the grain level.In this extended technique,an immersed moving boundary(IMB)scheme is used to couple the discrete element method(DEM)and lattice Boltzmann method(LBM),while a recently proposed Dirichlet-type thermal boundary condition is also adapted to account for heat transfer between fluid phase and solid particles.The resulting DEM-IBM-LBM model is robust to simulate moving curved boundaries with constant temperature in thermal flows.To facilitate the understanding and implementation of this coupled model for non-isothermal problems,a complete list is given for the conversion of relevant physical variables to lattice units.Then,benchmark tests,including a single-particle sedimentation and a two-particle drafting-kissing-tumbling(DKT)simulation with heat transfer,are carried out to validate the accuracy of our coupled technique.To further investigate the role of heat transfer in particle-laden flows,two multiple-particle problems with heat transfer are performed.Numerical examples demonstrate that the proposed coupling model is a promising high-resolution approach for simulating the heat-particle-fluid coupling at the grain level.展开更多
The objective of this investigation is to assess the effect of obstacles on numerical heat transfer and fluid flow momentum in a rectangular microchannel(MC).Two distinct configurations were studied:one without obstac...The objective of this investigation is to assess the effect of obstacles on numerical heat transfer and fluid flow momentum in a rectangular microchannel(MC).Two distinct configurations were studied:one without obstacles and the other with alternating obstacles placed on the upper and lower walls.The research utilized the thermal lattice Boltzmann method(LBM),which solves the energy and momentum equations of fluids with the BGK approximation,implemented in a Python coding environment.Temperature jump and slip velocity conditions were utilized in the simulation for the MC and extended to all obstacle boundaries.The study aims to analyze the rarefaction effect,with Knudsen numbers(Kn)of 0.012,0.02,and 0.05.The outcomes indicate that rarefaction has a significant impact on the velocity and temperature distribution.The presence of nine obstacles led to slower fluid movement inside the microchannel MC,resulting in faster cooling at the outlet.In MCs with obstacles,the rarefaction effect plays a crucial role in decreasing the Nusselt number(Nu)and skin friction coefficient(Cf).Furthermore,the study demonstrated that the obstacles played a crucial role in boosting fluid flow and heat transfer in the MC.The findings suggest that the examined configurations could have potential applications as cooling technologies in micro-electro-mechanical systems and microdevice applications.展开更多
A multiphase field model coupled with a lattice Boltzmann(PF-LBM)model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface,the interaction between dendrites and bubb...A multiphase field model coupled with a lattice Boltzmann(PF-LBM)model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface,the interaction between dendrites and bubbles,and the effects of different temperatures,anisotropic strengths and tilting angles on the solidified organization of the SCN-0.24wt.%butanedinitrile alloy during the solidification process.The model adopts a multiphase field model to simulate the growth of dendrites,calculates the growth motions of dendrites based on the interfacial solute equilibrium;and adopts a lattice Boltzmann model(LBM)based on the Shan-Chen multiphase flow to simulate the growth and motions of bubbles in the liquid phase,which includes the interaction between solid-liquid-gas phases.The simulation results show that during the directional growth of columnar dendrites,bubbles first precipitate out slowly at the very bottom of the dendrites,and then rise up due to the different solid-liquid densities and pressure differences.The bubbles will interact with the dendrite in the process of flow migration,such as extrusion,overflow,fusion and disappearance.In the case of wide gaps in the dendrite channels,bubbles will fuse to form larger irregular bubbles,and in the case of dense channels,bubbles will deform due to the extrusion of dendrites.In the simulated region,as the dendrites converge and diverge,the bubbles precipitate out of the dendrites by compression and diffusion,which also causes physical phenomena such as fusion and spillage of the bubbles.These results reveal the physical mechanisms of bubble nucleation,growth and kinematic evolution during solidification and interaction with dendrite growth.展开更多
This work focuses on numerically studying hydrodynamic interaction between a passive particle and a self-propelled particle,termed a squirmer,by using a two-dimensional lattice Boltzmann method(LBM).It is found that t...This work focuses on numerically studying hydrodynamic interaction between a passive particle and a self-propelled particle,termed a squirmer,by using a two-dimensional lattice Boltzmann method(LBM).It is found that the squirmer can capture a passive particle and propel it simultaneously,provided the passive particle is situated within the squirmer's wake.Our research shows that the critical capture distance,which determines whether the particle is captured,primarily depends on the intensity of the squirmer's dipolarity.The stronger dipolarity of squirmer results in an increased critical capture distance.Conversely,the Reynolds number is found to have minimal influence on this interaction.Interestingly,the passive particle,when driven by the squirmer's wake,contributes to a reduction in the squirmer's drag.This results in a mutual acceleration for both particles.Our findings can provide valuable perspectives for formulating the principles of reducing the drag of micro-swimmers and help to achieve the goal of using micro-swimmers to transport goods without physical tethers.展开更多
Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface...Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface. To address this challenge, we propose a novel numerical scheme that integrates the lattice Boltzmann method(LBM) and lattice spring model(LSM). In this scheme, LBM simulates viscoacoustic wave propagation in the fluid area and LSM simulates elastic wave propagation in the solid area. We also introduce three different LBM-LSM coupling strategies, a standard bounce back scheme, a specular reflection scheme, and a hybrid scheme, to describe wave propagation across fluid-solid boundaries. To demonstrate the accuracy of these LBM-LSM coupling schemes, we simulate wave propagation in a two-layer model containing a fluid-solid interface. We place excitation sources in the fluid layer and the solid layer respectively, to observe the wave phenomena when seismic waves propagate to interface from different sides. The simulated results by LBM-LSM are compared with the reference wavefields obtained by the finite difference method(FDM) and the analytical solution(ANA).Our LBM-LSM coupling scheme was verified effective, as the relative errors between the LBM-LSM solutions and reference solutions were within an acceptable range, sometimes around 1.00%. The coupled LBM-LSM scheme is further used to model seismic wavefields across a more realistic rugged seabed,which reveals the potential applications of the coupled LBM-LSM scheme in marine seismic imaging techniques, such as reverse-time migration and full-waveform inversion. The method also has potential applications in simulating wave propagation in complex two-and multi-phase media.展开更多
基金supported by the National Natural Science Foundation of China(12072347)the Excellent Training Plan of the Institute of Mechanics,Chinese Academy of SciencesCNPC New Energy Key Project(2021DJ4902).
文摘To study the development of imbibition such as the imbibition front and phase distribution in shale,the Lattice Boltzmann Method(LBM)is used to study the imbibition processes in the pore-throat network of shale.Through dimensional analysis,four dimensionless parameters affecting the imbibition process were determined.A color gradient model of LBM was used in computation based on a real core pore size distribution.The numerical results show that the four factors have great effects on imbibition.The impact of each factor is not monotonous.The imbibition process is the comprehensive effect of all aspects.The imbibition front becomes more and more non-uniform with time in a heterogeneous pore-throat network.Some non-wetting phases(oil here)cannot be displaced out.The displacement efficiency and velocity do not change monotonously with any factor.The development of the average imbibition length with time is not smooth and not linear in a heterogeneous pore-throat network.Two fitting relations between the four dimensionless parameters and the imbibition velocity and efficiency are obtained,respectively.
基金supported by the National Natural Science Foundation of China(Grant No.52022087).
文摘The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that of conventional reservoirs.At the early stage of flowback,there is no single-phase flow of the liquid phase in shale,but rather a gas-water two-phase flow,such that the single-phase flow model for tight oil and gas reservoirs is not applicable.In this study,pores and microfractures are extracted based on the experimental results of computed tomography(CT)scanning,and a spatial model of microfractures is established.Then,the influence of rough microfracture surfaces on the flow is corrected using the modified cubic law,which was modified by introducing the average deviation of the microfracture height as a roughness factor to consider the influence of microfracture surface roughness.The flow in the fracture network is simulated using the modified cubic law and the lattice Boltzmann method(LBM).The results obtained demonstrate that most of the fracturing fluid is retained in the shale microfractures,which explains the low fracturing fluid flowback rate in shale hydraulic fracturing.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.12072217 and 42077254)the Natural Science Foundation of Hunan Province,China(Grant No.2022JJ30567).
文摘The high-resolution DEM-IMB-LBM model can accurately describe pore-scale fluid-solid interactions,but its potential for use in geotechnical engineering analysis has not been fully unleashed due to its prohibitive computational costs.To overcome this limitation,a message passing interface(MPI)parallel DEM-IMB-LBM framework is proposed aimed at enhancing computation efficiency.This framework utilises a static domain decomposition scheme,with the entire computation domain being decomposed into multiple subdomains according to predefined processors.A detailed parallel strategy is employed for both contact detection and hydrodynamic force calculation.In particular,a particle ID re-numbering scheme is proposed to handle particle transitions across sub-domain interfaces.Two benchmarks are conducted to validate the accuracy and overall performance of the proposed framework.Subsequently,the framework is applied to simulate scenarios involving multi-particle sedimentation and submarine landslides.The numerical examples effectively demonstrate the robustness and applicability of the MPI parallel DEM-IMB-LBM framework.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
基金financially supported by the Natural Science Foundation of Hunan Province,China(Grant No.2022JJ30567)the support of EPSRC Grant(UK):PURIFY(EP/V000756/1)the Scientific Research Foundation of Education Department of Hunan Province,China(Grant No.20B557).
文摘Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate particle-fluid interaction problems involving heat transfer at the grain level.In this extended technique,an immersed moving boundary(IMB)scheme is used to couple the discrete element method(DEM)and lattice Boltzmann method(LBM),while a recently proposed Dirichlet-type thermal boundary condition is also adapted to account for heat transfer between fluid phase and solid particles.The resulting DEM-IBM-LBM model is robust to simulate moving curved boundaries with constant temperature in thermal flows.To facilitate the understanding and implementation of this coupled model for non-isothermal problems,a complete list is given for the conversion of relevant physical variables to lattice units.Then,benchmark tests,including a single-particle sedimentation and a two-particle drafting-kissing-tumbling(DKT)simulation with heat transfer,are carried out to validate the accuracy of our coupled technique.To further investigate the role of heat transfer in particle-laden flows,two multiple-particle problems with heat transfer are performed.Numerical examples demonstrate that the proposed coupling model is a promising high-resolution approach for simulating the heat-particle-fluid coupling at the grain level.
文摘The objective of this investigation is to assess the effect of obstacles on numerical heat transfer and fluid flow momentum in a rectangular microchannel(MC).Two distinct configurations were studied:one without obstacles and the other with alternating obstacles placed on the upper and lower walls.The research utilized the thermal lattice Boltzmann method(LBM),which solves the energy and momentum equations of fluids with the BGK approximation,implemented in a Python coding environment.Temperature jump and slip velocity conditions were utilized in the simulation for the MC and extended to all obstacle boundaries.The study aims to analyze the rarefaction effect,with Knudsen numbers(Kn)of 0.012,0.02,and 0.05.The outcomes indicate that rarefaction has a significant impact on the velocity and temperature distribution.The presence of nine obstacles led to slower fluid movement inside the microchannel MC,resulting in faster cooling at the outlet.In MCs with obstacles,the rarefaction effect plays a crucial role in decreasing the Nusselt number(Nu)and skin friction coefficient(Cf).Furthermore,the study demonstrated that the obstacles played a crucial role in boosting fluid flow and heat transfer in the MC.The findings suggest that the examined configurations could have potential applications as cooling technologies in micro-electro-mechanical systems and microdevice applications.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.52161002,51661020,and 11364024)the Postdoctoral Science Foundation of China(Grant No.2014M560371)the Funds for Distinguished Young Scientists of Lanzhou University of Technology of China(Grant No.J201304).
文摘A multiphase field model coupled with a lattice Boltzmann(PF-LBM)model is proposed to simulate the distribution mechanism of bubbles and solutes at the solid-liquid interface,the interaction between dendrites and bubbles,and the effects of different temperatures,anisotropic strengths and tilting angles on the solidified organization of the SCN-0.24wt.%butanedinitrile alloy during the solidification process.The model adopts a multiphase field model to simulate the growth of dendrites,calculates the growth motions of dendrites based on the interfacial solute equilibrium;and adopts a lattice Boltzmann model(LBM)based on the Shan-Chen multiphase flow to simulate the growth and motions of bubbles in the liquid phase,which includes the interaction between solid-liquid-gas phases.The simulation results show that during the directional growth of columnar dendrites,bubbles first precipitate out slowly at the very bottom of the dendrites,and then rise up due to the different solid-liquid densities and pressure differences.The bubbles will interact with the dendrite in the process of flow migration,such as extrusion,overflow,fusion and disappearance.In the case of wide gaps in the dendrite channels,bubbles will fuse to form larger irregular bubbles,and in the case of dense channels,bubbles will deform due to the extrusion of dendrites.In the simulated region,as the dendrites converge and diverge,the bubbles precipitate out of the dendrites by compression and diffusion,which also causes physical phenomena such as fusion and spillage of the bubbles.These results reveal the physical mechanisms of bubble nucleation,growth and kinematic evolution during solidification and interaction with dendrite growth.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.12132015 and 11972336)。
文摘This work focuses on numerically studying hydrodynamic interaction between a passive particle and a self-propelled particle,termed a squirmer,by using a two-dimensional lattice Boltzmann method(LBM).It is found that the squirmer can capture a passive particle and propel it simultaneously,provided the passive particle is situated within the squirmer's wake.Our research shows that the critical capture distance,which determines whether the particle is captured,primarily depends on the intensity of the squirmer's dipolarity.The stronger dipolarity of squirmer results in an increased critical capture distance.Conversely,the Reynolds number is found to have minimal influence on this interaction.Interestingly,the passive particle,when driven by the squirmer's wake,contributes to a reduction in the squirmer's drag.This results in a mutual acceleration for both particles.Our findings can provide valuable perspectives for formulating the principles of reducing the drag of micro-swimmers and help to achieve the goal of using micro-swimmers to transport goods without physical tethers.
基金supported in part by R & D Department of China National Petroleum Corporation (2022DQ0604-01)National Natural Science Foundation of China (42204132)+3 种基金the China Postdoctoral Science Foundations (2020M680667, 2021T140661)Harvard-CUP Joint Laboratory on Petroleum Science“111” project (B13010)the financial support from the CAS Special Research Assistant Project。
文摘Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface. To address this challenge, we propose a novel numerical scheme that integrates the lattice Boltzmann method(LBM) and lattice spring model(LSM). In this scheme, LBM simulates viscoacoustic wave propagation in the fluid area and LSM simulates elastic wave propagation in the solid area. We also introduce three different LBM-LSM coupling strategies, a standard bounce back scheme, a specular reflection scheme, and a hybrid scheme, to describe wave propagation across fluid-solid boundaries. To demonstrate the accuracy of these LBM-LSM coupling schemes, we simulate wave propagation in a two-layer model containing a fluid-solid interface. We place excitation sources in the fluid layer and the solid layer respectively, to observe the wave phenomena when seismic waves propagate to interface from different sides. The simulated results by LBM-LSM are compared with the reference wavefields obtained by the finite difference method(FDM) and the analytical solution(ANA).Our LBM-LSM coupling scheme was verified effective, as the relative errors between the LBM-LSM solutions and reference solutions were within an acceptable range, sometimes around 1.00%. The coupled LBM-LSM scheme is further used to model seismic wavefields across a more realistic rugged seabed,which reveals the potential applications of the coupled LBM-LSM scheme in marine seismic imaging techniques, such as reverse-time migration and full-waveform inversion. The method also has potential applications in simulating wave propagation in complex two-and multi-phase media.