Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was c...Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.展开更多
The variable salinity in stored reservoirs connected by a long channel attracts the attention of scientists worldwide,having applications in environmental and geophysical engineering.This study explores the impact of ...The variable salinity in stored reservoirs connected by a long channel attracts the attention of scientists worldwide,having applications in environmental and geophysical engineering.This study explores the impact of Navier slip conditions on exchange flows within a long channel connecting two large reservoirs of differing salinity.These horizontal density gradients drive the flow.We modify the recent one-dimensional theory,developed to avoid runaway stratification,to account for the presence of uniform slip walls.By adjusting the parameters of the horizontal density gradient based on the slip factor,we resolve analytically various flow regimes ranging from high diffusion to transitional high advection.These regimes are governed by physical parameters like channel aspect ratio,slip factor,Schmidt number,and gravitational Reynolds number.Our solutions align perfectly with ones in the no-slip limit.More importantly,under the conditions of no net flow across the channel and high Schmidt number(where stratification is concentrated near the channel’s mid-layer),we derive a closed-form solution for the slip parameter,aspect ratio,and gravitational Reynolds number that describes the interface’s behavior as a sharp interface separating two distinct zones.This interface,arising from hydrostatic wall gradients,ultimately detaches the low-and high-density regimes throughout the channel when the gravitational Reynolds number is inversely proportional to the aspect ratio for a fixed slip parameter.This phenomenon,observed previously in 2D numerical simulations with no-slip walls in the literature,is thus confirmed by our theoretical results.Our findings further demonstrate that wall slip leads to distinct and diverse flow regimes.展开更多
It is common for the rotating pin largely plunging into lower sheet to break up lap interface of Friction Stir Lap Welding(FSLW)joint,but the unavoidable up-bending morphology of hook outside Nugget Zone(NZ)largely re...It is common for the rotating pin largely plunging into lower sheet to break up lap interface of Friction Stir Lap Welding(FSLW)joint,but the unavoidable up-bending morphology of hook outside Nugget Zone(NZ)largely reduces the joint bearing ability.Based on the novel Opposite-directions Flowing FSLW(OF-FSLW)by the self-developed rotating tool with an Xshaped right-left thread pin,the 2024 aluminum alloys lap joint was successfully welded in this study.The migration law of lap interface during welding was investigated by the experimental and numerical methods,and then how the rotating pin and its rotating velocity affect the formation and strength of OF-FSLW joint was further analyzed.The results show that the Material Concentrated Zone(MCZ)which formed above the original lap interface made the hook bend downward,the NZ greatly enlarged and the beginning part of cold lap compressed and thickened,thereby heightening the joint bearing ability.For the OF-FSLW joint,its maximum tensile strength was403 MPa,and the corresponding joint efficiency of 90.8%was an incredible and superb value for the 2000 series heat-treatment strengthened aluminum alloys friction stir welded joint.The OF-FSLW technology by the rotating tool with an X-shaped right-left thread pin is proven to be a greatly effective approach for manufacturing the aluminum alloys lap joint with superb strength.展开更多
A hierarchical model is developed to predict the streaming potential （SP） in the canaliculi of a loaded os teon. Canaliculi are assumed to run straight across the os teon annular cylinder wall, while disregarding th...A hierarchical model is developed to predict the streaming potential （SP） in the canaliculi of a loaded os teon. Canaliculi are assumed to run straight across the os teon annular cylinder wall, while disregarding the effect of lacuna. SP is generalized by the canalicular fluid flow. Ana lytical solutions are obtained for the canalicular fluid veloc ity, pressure, and SP. Results demonstrate that SP amplitude （SPA） is proportional to the pressure difference, strain am plitude, frequency, and strain rate amplitude. However, the key loading factor governing SP is the strain rate, which is a representative loading parameter under the specific phys iological state. Moreover, SPA is independent of canalicu lar length. This model links external loads to the canalicu lar fluid pressure, velocity, and SP, which can facilitate fur ther understanding of the mechanotransduction and electro mechanotransduction mechanisms of bones.展开更多
Oscillatory flow facilitates gas exchange in human respiration system. In the present study, both numerical calculation and PIV (Particle Image Velocimetry) measurement indicate that, under the application of HFOV (Hi...Oscillatory flow facilitates gas exchange in human respiration system. In the present study, both numerical calculation and PIV (Particle Image Velocimetry) measurement indicate that, under the application of HFOV (High Frequency Oscillatory Ventilation), apparent steady streaming is caused and augmented in distal airways by the continuous oscillation, i.e., the core air moves downwards and the peripheral air evacuates upwards within bronchioles. The net flow of steady streaming serves to overcome the lack of tidal volume in HFOV and delivers fresh air into deeper lung region. Also, numerical calculations reveal that the intensity of steady streaming is mainly influenced by the geometry of airways with provided oscillatory frequency and tidal volume, and it rises with Re and Wo up to a Re of about 124 and Wo of about 5. Steady streaming is considered as an important factor for the ventilation efficiency of HFOV.展开更多
Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the ...Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.展开更多
This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunzia...This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.展开更多
基金the financial support from the National Natural Science Foundation of China (No.42102127)the Postdoctoral Research Foundation of China (No.2024 M751860)。
文摘Cleat serves as the primary flow pathway for coalbed methane(CBM)and water.However,few studies consider the impact of local contact on two-phase flow within cleats.A visual generalized model of endogenous cleats was constructed based on microfluidics.A microscopic and mesoscopic observation technique was proposed to simultaneously capture gas-liquid interface morphology of pores and throat and the two-phase flow characteristics in entire cleat system.The local contact characteristics of cleats reduced absolute permeability,which resulted in a sharp increase in the starting pressure.The reduced gas flow capacity narrowed the co-infiltration area and decreased water saturation at the isotonic point in a hydrophilic environment.The increased local contact area of cleats weakened gas phase flow capacity and narrowed the co-infiltration area.Jumping events occurred in methane-water flow due to altered porosity caused by local contact in cleats.The distribution of residual phases changed the jumping direction on the micro-scale as well as the dominant channel on the mesoscale.Besides,jumping events caused additional energy dissipation,which was ignored in traditional two-phase flow models.This might contribute to the overestimation of relative permeability.The work provides new methods and insights for investigating unsaturated flow in complex porous media.
文摘The variable salinity in stored reservoirs connected by a long channel attracts the attention of scientists worldwide,having applications in environmental and geophysical engineering.This study explores the impact of Navier slip conditions on exchange flows within a long channel connecting two large reservoirs of differing salinity.These horizontal density gradients drive the flow.We modify the recent one-dimensional theory,developed to avoid runaway stratification,to account for the presence of uniform slip walls.By adjusting the parameters of the horizontal density gradient based on the slip factor,we resolve analytically various flow regimes ranging from high diffusion to transitional high advection.These regimes are governed by physical parameters like channel aspect ratio,slip factor,Schmidt number,and gravitational Reynolds number.Our solutions align perfectly with ones in the no-slip limit.More importantly,under the conditions of no net flow across the channel and high Schmidt number(where stratification is concentrated near the channel’s mid-layer),we derive a closed-form solution for the slip parameter,aspect ratio,and gravitational Reynolds number that describes the interface’s behavior as a sharp interface separating two distinct zones.This interface,arising from hydrostatic wall gradients,ultimately detaches the low-and high-density regimes throughout the channel when the gravitational Reynolds number is inversely proportional to the aspect ratio for a fixed slip parameter.This phenomenon,observed previously in 2D numerical simulations with no-slip walls in the literature,is thus confirmed by our theoretical results.Our findings further demonstrate that wall slip leads to distinct and diverse flow regimes.
基金supported by the National Natural Science Foundation of China(Nos.51874201 and 52074184)Natural Science Foundation project of Liaoning Province(2023-MS238)+1 种基金Liaoning Revitalization talents program(XLYC2203009)the Fundamental Research Funds for the Universities of Liaoning Province。
文摘It is common for the rotating pin largely plunging into lower sheet to break up lap interface of Friction Stir Lap Welding(FSLW)joint,but the unavoidable up-bending morphology of hook outside Nugget Zone(NZ)largely reduces the joint bearing ability.Based on the novel Opposite-directions Flowing FSLW(OF-FSLW)by the self-developed rotating tool with an Xshaped right-left thread pin,the 2024 aluminum alloys lap joint was successfully welded in this study.The migration law of lap interface during welding was investigated by the experimental and numerical methods,and then how the rotating pin and its rotating velocity affect the formation and strength of OF-FSLW joint was further analyzed.The results show that the Material Concentrated Zone(MCZ)which formed above the original lap interface made the hook bend downward,the NZ greatly enlarged and the beginning part of cold lap compressed and thickened,thereby heightening the joint bearing ability.For the OF-FSLW joint,its maximum tensile strength was403 MPa,and the corresponding joint efficiency of 90.8%was an incredible and superb value for the 2000 series heat-treatment strengthened aluminum alloys friction stir welded joint.The OF-FSLW technology by the rotating tool with an X-shaped right-left thread pin is proven to be a greatly effective approach for manufacturing the aluminum alloys lap joint with superb strength.
基金supported by the OIT of Higher Learning Institutions of Shanxithe National Natural Science Foundation of China(11302143,11472185)Natural Science Fund of Shanxi(2014021013)
文摘A hierarchical model is developed to predict the streaming potential （SP） in the canaliculi of a loaded os teon. Canaliculi are assumed to run straight across the os teon annular cylinder wall, while disregarding the effect of lacuna. SP is generalized by the canalicular fluid flow. Ana lytical solutions are obtained for the canalicular fluid veloc ity, pressure, and SP. Results demonstrate that SP amplitude （SPA） is proportional to the pressure difference, strain am plitude, frequency, and strain rate amplitude. However, the key loading factor governing SP is the strain rate, which is a representative loading parameter under the specific phys iological state. Moreover, SPA is independent of canalicu lar length. This model links external loads to the canalicu lar fluid pressure, velocity, and SP, which can facilitate fur ther understanding of the mechanotransduction and electro mechanotransduction mechanisms of bones.
文摘Oscillatory flow facilitates gas exchange in human respiration system. In the present study, both numerical calculation and PIV (Particle Image Velocimetry) measurement indicate that, under the application of HFOV (High Frequency Oscillatory Ventilation), apparent steady streaming is caused and augmented in distal airways by the continuous oscillation, i.e., the core air moves downwards and the peripheral air evacuates upwards within bronchioles. The net flow of steady streaming serves to overcome the lack of tidal volume in HFOV and delivers fresh air into deeper lung region. Also, numerical calculations reveal that the intensity of steady streaming is mainly influenced by the geometry of airways with provided oscillatory frequency and tidal volume, and it rises with Re and Wo up to a Re of about 124 and Wo of about 5. Steady streaming is considered as an important factor for the ventilation efficiency of HFOV.
基金Project supported by the National Natural Science Foundation of China(No.52109068)the Water Conservancy Technology Project of Jiangsu Province of China(No.2022060)。
文摘Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.
基金funded by the Italian Ministry of Education,University and Research(MIUR)in the frame of the Departments of Excellence Initiative 2018-2027 attributed to DICAM of the University of Trento(grant L.232/2016)in the frame of the PRIN 2017 project Innovative numerical methods for evolutionary partial differential equations and applications,the PRIN 2022 project High order structure-preserving semi-implicit schemes for hyperbolic equations.D.is member of INdAM GNCS and was also co-funded by the European Union NextGenerationEU(PNRR,Spoke 7 CN HPC).Views and opinions expressed are however those of the author(s)only and do not necessarily reflect those of the European Union or the European Research Council.Neither the European Union nor the granting authority can be held responsible for them.
文摘This paper presents a mass and momentum conservative semi-implicit finite volume(FV)scheme for complex non-hydrostatic free surface flows,interacting with moving solid obstacles.A simplified incompressible Baer-Nunziato type model is considered for two-phase flows containing a liquid phase,a solid phase,and the surrounding void.According to the so-called diffuse interface approach,the different phases and consequently the void are described by means of a scalar volume fraction function for each phase.In our numerical scheme,the dynamics of the liquid phase and the motion of the solid are decoupled.The solid is assumed to be a moving rigid body,whose motion is prescribed.Only after the advection of the solid volume fraction,the dynamics of the liquid phase is considered.As usual in semi-implicit schemes,we employ staggered Cartesian control volumes and treat the nonlinear convective terms explicitly,while the pressure terms are treated implicitly.The non-conservative products arising in the transport equation for the solid volume fraction are treated by a path-conservative approach.The resulting semi-implicit FV discretization of the mass and momentum equations leads to a mildly nonlinear system for the pressure which can be efficiently solved with a nested Newton-type technique.The time step size is only limited by the velocities of the two phases contained in the domain,and not by the gravity wave speed nor by the stiff algebraic relaxation source term,which requires an implicit discretization.The resulting semi-implicit algorithm is first validated on a set of classical incompressible Navier-Stokes test problems and later also adds a fixed and moving solid phase.
