In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the phys...In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the physical domain to resolvethe relevant scales in multiscale physical systems while minimizing computationalcosts. The algorithm is a generalization of the moving mesh methods basedon harmonic mappings developed by Li et al. [J. Comput. Phys., 170 (2001), pp. 562-588, and 177 (2002), pp. 365-393]. To make 3D moving mesh simulations possible,the key is to develop an efficient mesh redistribution procedure so that this part willcost as little as possible comparing with the solution evolution part. Since the meshredistribution procedure normally requires to solve large size matrix equations, wewill describe a procedure to decouple the matrix equation to a much simpler blocktridiagonaltype which can be efficiently solved by a particularly designed multi-gridmethod. To demonstrate the performance of the proposed 3D moving mesh strategy,the algorithm is implemented in finite element simulations of fluid-fluid interface interactionsin multiphase flows. To demonstrate the main ideas, we consider the formationof drops by using an energetic variational phase field model which describesthe motion of mixtures of two incompressible fluids. Numerical results on two- andthree-dimensional simulations will be presented.展开更多
A thermodynamically complete multi-phase equation of state(EOS)applicable to both dense and porous metals at wide ranges of temperature and pressure is constructed.A standard three-term decomposition of the Helmholtz ...A thermodynamically complete multi-phase equation of state(EOS)applicable to both dense and porous metals at wide ranges of temperature and pressure is constructed.A standard three-term decomposition of the Helmholtz free energy as a function of specific volume and temperature is presented,where the cold component models both compression and expansion states,the thermal ion component introduces the Debye approximation and melting entropy,and the thermal electron component employs the Thomas-Fermi-Kirzhnits(TFK)model.The porosity of materials is considered by introducing the dynamic porosity coefficientαand the constitutive P-αrelation,connecting the thermodynamic properties between dense and porous systems,allowing for an accurate description of the volume decrease caused by void collapse while maintaining the quasi-static thermodynamic properties of porous systems identical to the dense ones.These models enable the EOS applicable and robust at wide ranges of temperature,pressure and porosity.A systematic evaluation of the new EOS is conducted with aluminum(Al)as an example.300 K isotherm,shock Hugoniot,as well as melting curves of both dense and porous Al are calculated,which shows great agreements with experimental data and validates the effectiveness of the models and the accuracy of parameterizations.Notably,it is for the first time Hugoniot P-σcurves up to 10~6 GPa and shock melting behaviors of porous Al are derived from analytical EOS models,which predict much lower compression limit and shock melting temperatures than those of dense Al.展开更多
This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with ...This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with a filter-based density correction turbulence model and a modified Zwart cavitation model.The study investigates the dynamic cavitation features of the thermal fluid around the hydrofoil at various incoming flow velocities.It systematically elucidates the evolution of cavitation and vortex dynamics corresponding to each velocity condition.The results indicate that with increasing incoming flow velocity,distinct cavitation processes take place in the flow field.展开更多
The impacts of natural boulders carried by debris flows pose serious risks to the safety and reliability of structures and buildings.Natural boulders can be highly random and unpredictable.Consequently,boulder control...The impacts of natural boulders carried by debris flows pose serious risks to the safety and reliability of structures and buildings.Natural boulders can be highly random and unpredictable.Consequently,boulder control during debris flows is crucial but difficult.Herein,an eco-friendly control system featuring anchoring natural boulders(NBs)with(negative Poisson's ratio)NPR anchor cables is proposed to form an NB-NPR baffle.A series of flume experiments are conducted to verify the effect of NB-NPR baffles on controlling debris flow impact.The deployment of NB-NPR baffles substantially influences the kinematic behavior of a debris flow,primarily in the form of changes in the depositional properties and impact intensities.The results show that the NB-NPR baffle matrix successfully controls boulder mobility and exhibits positive feedback on solid particle deposition.The NB-NPR baffle group exhibits a reduction in peak impact force ranging from 29%to 79%compared to that of the control group in the basic experiment.The NPR anchor cables play a significant role in the NB-NPR baffle by demonstrating particular characteristics,including consistent resistance,large deformation,and substantial energy absorption.The NB-NPR baffle innovatively utilizes the natural boulders in a debris flow gully by converting destructive boulders into constructive boulders.Overall,this research serves as a basis for future field experiments and applications.展开更多
In the petroleum industry,detection of multi-phase fluid flow is very important in both surface and down-hole measurements.Accurate measurement of high rate of water or gas multi-phase flow has always been an academic...In the petroleum industry,detection of multi-phase fluid flow is very important in both surface and down-hole measurements.Accurate measurement of high rate of water or gas multi-phase flow has always been an academic and industrial focus.NMR is an efficient and accurate technique for the detection of fluids;it is widely used in the determination of fluid compositions and properties.This paper is aimed to quantitatively detect multi-phase flow in oil and gas wells and pipelines and to propose an innovative method for online nuclear magnetic resonance(NMR)detection.The online NMR data acquisition,processing and interpretation methods are proposed to fill the blank of traditional methods.A full-bore straight tube design without pressure drop,a Halbach magnet structure design with zero magnetic leakage outside the probe,a separate antenna structure design without flowing effects on NMR measurement and automatic control technology will achieve unattended operation.Through the innovation of this work,the application of NMR for the real-time and quantitative detection of multi-phase flow in oil and gas wells and pipelines can be implemented.展开更多
Most current lattice Boltzmann （LBM） models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular ...Most current lattice Boltzmann （LBM） models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular dynamics （MD） simulation and LBM to avoid such defect. The basic idea is to first construct a molecular model based on the actual components of the rock-fluid system, then to compute the interaction force between the rock and the fluid of different densities through the MD simulation. This calculated rock-fluid interaction force, combined with the fluid-fluid force determined from the equation of state, is then used in LBM modeling. Without parameter fitting, this study presents a new systematic approach for pore-scale modeling of multi-phase flow. We have validated this ap- proach by simulating a two-phase separation process and gas-liquid-solid three-phase contact angle. Based on an actual X-ray CT image of a reservoir core, we applied our workflow to calculate the absolute permeability of the core, vapor-liquid H20 relative permeability, and capillary pressure curves.展开更多
The transient friction in channel mean flows is the sum of two contributions,i.e.,the underlying laminar flow(ULF)and the purely turbulent component(PTC),and the contributions are analyzed separately by theoretical ex...The transient friction in channel mean flows is the sum of two contributions,i.e.,the underlying laminar flow(ULF)and the purely turbulent component(PTC),and the contributions are analyzed separately by theoretical experiments.It is found that,the transient friction may be higher or remarkably lower than that in equal-Reynolds number steady-state flows.The universal time constant for plane-parallel laminar flows is reported,and the role of the time constant in a turbulent mean flow is examined.It is shown that the time constant is related to the turbulence's frozen time.Finally,a study of the logarithmic layer during the transient flow is accomplished,which shows that the logarithmic layer is destroyed.展开更多
Two-phase pipe flow occurs frequently in oil&gas industry,nuclear power plants,and CCUS.Reliable calculations of gas void fraction(or liquid holdup)play a central role in two-phase pipe flow models.In this paper w...Two-phase pipe flow occurs frequently in oil&gas industry,nuclear power plants,and CCUS.Reliable calculations of gas void fraction(or liquid holdup)play a central role in two-phase pipe flow models.In this paper we apply the fractional flow theory to multiphase flow in pipes and present a unified modeling framework for predicting the fluid phase volume fractions over a broad range of pipe flow conditions.Compared to existing methods and correlations,this new framework provides a simple,approximate,and efficient way to estimate the phase volume fraction in two-phase pipe flow without invoking flow patterns.Notably,existing correlations for estimating phase volume fraction can be transformed and expressed under this modeling framework.Different fractional flow models are applicable to different flow conditions,and they demonstrate good agreement against experimental data within 5%errors when compared with an experimental database comprising of 2754 data groups from 14literature sources,covering various pipe geometries,flow patterns,fluid properties and flow inclinations.The gas void fraction predicted by the framework developed in this work can be used as inputs to reliably model the hydraulic and thermal behaviors of two-phase pipe flows.展开更多
Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 t...Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 to 44 km from the source volcano. The blockage impounded a body of water which then released producing four stratigraphic units in the downstream debris fan. The four stratigraphic units are a boulder core comprised of locally sourced bedrock boulders and three sand-rich units including a fine-grained sand unit, a sandy pumice gravel (±basalt/hydrovolcanic tuff) unit, and a pumice pebble-bearing, crystal-rich sand unit. Hand-drilled auger holes up to ~1.6 m deep were used to obtain samples of the sand-rich units. Units were delimited using surface and down-hole observations, composition and texture, estimated density, statistical parameters of grain size, and vertical and lateral distribution of properties. Overtopping followed by rapid incision into the ash-rich pyroclastic flows progressively cleared the canyon, but a bedrock knickpoint near the head of the canyon limited the volume of debris available for transport to about 0.04 km<sup>3</sup> to 0.08 km<sup>3</sup>. Co-deposition of bedrock boulders and lithic-rich sand was followed by rapid deposition with minimal reworking of remobilized pyroclastics. Continued draining of the impounded lake sent hyperconcentrated flows onto the debris fan depositing pumice-rich gravels that graded upward to crystal-rich sands.展开更多
Based on the theory of superimposed deformation and the regional tectonic background,the multi-phase non-coaxial superimposed structures in Junggar Basin were systematically analyzed using seismic interpretation,field...Based on the theory of superimposed deformation and the regional tectonic background,the multi-phase non-coaxial superimposed structures in Junggar Basin were systematically analyzed using seismic interpretation,field outcrop observation,and paleo-stress field recovery methods according to the characteristics of the current tectonic framework.Moreover,the tectonic evolution process of the basin was reconstructed using sandbox analogue modelling technology.The results showed that the study area has experienced five phases of non-coaxial deformation with superimposition:The first phase of deformation(D1)is characterized by NNE-SSW extension during late Carboniferous to early Permian,which formed large graben,half graben and other extensional structure style around the basin.The second phase of deformation(D2)is represented by NE-SW compression during the middle to late Permian,and it comprised numerous contraction structures that developed based on D1.The basic form of the entire basin is alternating uplift and depression.The third phase of deformation(D3)is the NW-SE transpressional strike-slip in the Triassic-Jurassic,which produced numerous strike-slip structural styles in the middle part of the basin.The fourth phase of deformation(D4)is the uniform sedimentation during Cretaceous,and the fifth phase(D5)is the compression along NNE-SSW due to the North Tianshan northward thrust,which produced three rows of fold thrust belts and tear faults in the front of the mountain in the southern margin of the basin.The newly established three-dimensional tectonic evolution model shows that,based on the large number of NW-trending grabens and half grabens in the Carboniferous basement of Junggar Basin,multiple level NE trending uplifts have formed with the joint superposition of the late structural inversion and multiple stress fields.This has resulted in the current tectonic units of alternating uplifts and depressions in different directions in the study area.展开更多
Conventional conductivity methods for measuring the void fraction in gas-liquid multiphase systems are typically affected by accuracy problems due to the presence of fluid flow and salinity.This study presents a novel...Conventional conductivity methods for measuring the void fraction in gas-liquid multiphase systems are typically affected by accuracy problems due to the presence of fluid flow and salinity.This study presents a novel approach for determining the void fraction based on a reciprocating dynamic conductivity probe used to measure the liquid film thickness under forced annular-flow conditions.The measurement system comprises a cyclone,a conductivity probe,a probe reciprocating device,and a data acquisition and processing system.This method ensures that the flow pattern is adjusted to a forced annular flow,thereby minimizing the influence of complex and variable gas-liquid flow patterns on the measurement results;Moreover,it determines the liquid film thickness solely according to circuit connectivity rather than specific conductivity values,thereby mitigating the impact of salinity.The reliability of the measurement system is demonstrated through laboratory experiments.The experimental results indicate that,in a range of gas phase superficial velocities 5–20 m/s and liquid phase superficial velocities 0.079–0.48 m/s,the maximum measurement deviation for the void fraction is 4.23%.展开更多
Capturing elaborated flow structures and phenomena is required for well-solved numerical flows.The finite difference methods allow simple discretization of mesh and model equations.However,they need simpler meshes,e.g...Capturing elaborated flow structures and phenomena is required for well-solved numerical flows.The finite difference methods allow simple discretization of mesh and model equations.However,they need simpler meshes,e.g.,rectangular.The inverse Lax-Wendroff(ILW)procedure can handle complex geometries for rectangular meshes.High-resolution and high-order methods can capture elaborated flow structures and phenomena.They also have strong mathematical and physical backgrounds,such as positivity-preserving,jump conditions,and wave propagation concepts.We perceive an effort toward direct numerical simulation,for instance,regarding weighted essentially non-oscillatory(WENO)schemes.Thus,we propose to solve a challenging engineering application without turbulence models.We aim to verify and validate recent high-resolution and high-order methods.To check the solver accuracy,we solved vortex and Couette flows.Then,we solved inviscid and viscous nozzle flows for a conical profile.We employed the finite difference method,positivity-preserving Lax-Friedrichs splitting,high-resolution viscous terms discretization,fifth-order multi-resolution WENO,ILW,and third-order strong stability preserving Runge-Kutta.We showed the solver is high-order and captured elaborated flow structures and phenomena.One can see oblique shocks in both nozzle flows.In the viscous flow,we also captured a free-shock separation,recirculation,entrainment region,Mach disk,and the diamond-shaped pattern of nozzle flows.展开更多
This study investigates the impact of pore network characteristics on fluid flow through complex and heterogeneous porous media,providing insights into the factors affecting fluid propagation in such systems.Specifica...This study investigates the impact of pore network characteristics on fluid flow through complex and heterogeneous porous media,providing insights into the factors affecting fluid propagation in such systems.Specifically,high-resolution or micro X-ray computed tomography(CT)imaging techniques were utilized to examine outcrop stromatolite samples of the Lagoa Salgada,considered flow analogous to the Brazilian Pre-salt carbonate reservoirs.The petrophysical results comprised two distinct stromatolite depositional facies,the columnar and the fine-grained facies.By generating pore network model(PNM),the study quantified the relationship between key features of the porous system,including pore and throat radius,throat length,coordination number,shape factor,and pore volume.The study found that the less dense pore network of the columnar sample is typically characterized by larger pores and wider and longer throats but with a weaker connection of throats to pores.Both facies exhibited less variability in the radius of the pores and throats in comparison to throat length.Additionally,a series of core flooding experiments coupled with medical CT scanning was designed and conducted in the plug samples to assess flow propagation and saturation fields.The study revealed that the heterogeneity and presence of disconnected or dead-end pores significantly impacted the flow patterns and saturation.Two-phase flow patterns and oil saturation distribution reveal a preferential and heterogeneous displacement that mainly swept displaced fluid in some regions of plugs and bypassed it in others.The relation between saturation profiles,porosity profiles,and the number of fluid flow patterns for the samples was evident.Only for the columnar plug sample was the enhancement in recovery factor after shifting to lower salinity water injection(SB)observed.展开更多
Monocular 3D object detection is challenging due to the lack of accurate depth information.Some methods estimate the pixel-wise depth maps from off-the-shelf depth estimators and then use them as an additional input t...Monocular 3D object detection is challenging due to the lack of accurate depth information.Some methods estimate the pixel-wise depth maps from off-the-shelf depth estimators and then use them as an additional input to augment the RGB images.Depth-based methods attempt to convert estimated depth maps to pseudo-LiDAR and then use LiDAR-based object detectors or focus on the perspective of image and depth fusion learning.However,they demonstrate limited performance and efficiency as a result of depth inaccuracy and complex fusion mode with convolutions.Different from these approaches,our proposed depth-guided vision transformer with a normalizing flows(NF-DVT)network uses normalizing flows to build priors in depth maps to achieve more accurate depth information.Then we develop a novel Swin-Transformer-based backbone with a fusion module to process RGB image patches and depth map patches with two separate branches and fuse them using cross-attention to exchange information with each other.Furthermore,with the help of pixel-wise relative depth values in depth maps,we develop new relative position embeddings in the cross-attention mechanism to capture more accurate sequence ordering of input tokens.Our method is the first Swin-Transformer-based backbone architecture for monocular 3D object detection.The experimental results on the KITTI and the challenging Waymo Open datasets show the effectiveness of our proposed method and superior performance over previous counterparts.展开更多
Forests play an important role in controlling the formation and movement processes of debris flows.They contribute to soil stabilization,regulation of soil water content,and act as robust structures impeding the downs...Forests play an important role in controlling the formation and movement processes of debris flows.They contribute to soil stabilization,regulation of soil water content,and act as robust structures impeding the downstream progression of debris flows.On the positive side,trees,to some extent,can intercept debris flows and effectively mitigate their velocity by increasing flow resistance.On the negative side,trees may suffer damage from debris-flow hazards,characterized by the generation of substantial quantities of wood fragments and consequential ramifications such as river channel blockage,resulting in backwater rise.In extreme cases,this blockage collapse can lead to instantaneous discharge amplification,thereby adversely impacting urban safety and impeding sustainable development.Therefore,in order to grasp the effects of tree characteristics on tree failure modes,the tree failure modes and corresponding parameters,diameters at breast height(DBH)and root-soil plate size,were identified and recorded through the post-event field investigation in Keze Gully,a region prone to debrisflow events in Sichuan,China,respectively.To investigate the impact of spatial variability in tree root distribution on tree failure modes,the root crosssectional area ratio(RAR),root density(RD),root length density(RLD)and soil detachment rate(SDR)were obtained.The findings indicated that:(1)Tree characteristics reflect the interactions of debris flows and trees,and influence the tree failure modes ultimately.The root distribution characteristics influence the size and shape of the root-soil plate to affect the resistance of trees.(2)Compared to burial and abrasion,stem breakage and overturning are the predominant modes of tree failure in debris-flow hazards.Trees with a smaller DBH primarily experience stem breakage and bending,and trees with a larger DBH mostly experience overturning.(3)The root-soil plate shapes of overturned trees,affected by the root architecture and root growth range,are generally semielliptical or semicircular,and the horizontal and vertical radii increase with DBH,but the correlation between the root-soil plate’s breadth-depth ratio and DBH is low.(4)The biomass and RAR decrease with distance.The RAR distribution exhibit the order of upslope direction>downslope direction>lateral direction.The coarse root biomass significantly increases with DBH,but no clear trend in fine root biomass.(5)The roots can significantly enhance the soil erosion resistance,but the erosion resistance of coarse roots is not as significant as that of fine roots.The erosion resistance increases with DBH,and follows the order of upslope direction>downslope direction>lateral direction.The results could provide new insights into the influences of tree and root distribution characteristics on tree failure modes during debris flows.展开更多
A model is proposed for liquid film profile prediction in gas-liquid two-phase flow,which is able to provide the film thickness along the circumferential direction and the pressure gradient in the flow direction.A two...A model is proposed for liquid film profile prediction in gas-liquid two-phase flow,which is able to provide the film thickness along the circumferential direction and the pressure gradient in the flow direction.A two-fluid model is used to calculate both gas and liquid phases’flow characteristics.The secondary flow occurring in the gas phase is taken into account and a sailing boat mechanism is introduced.Moreover,energy conservation is applied for obtaining the liquid film thickness distribution along the circumference.Liquid film thickness distribution is calculated accordingly for different cases;its values are compared with other models and available experimental data.As a result,the newly proposed model is tested and good performances are demonstrated.The liquid film thickness distribution in small pipes and inclined pipes is also studied,and regime transition is revealed by liquid film profile evolution.The observed inflection point demonstrates that the liquid film thickness decreases steeply along the circumference,when the circle angle ranges between 30°and 50°for gas-liquid stratified flow with small superficial velocities.展开更多
The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direc...The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.展开更多
The accelerated method in solving optimization problems has always been an absorbing topic.Based on the fixedtime(FxT)stability of nonlinear dynamical systems,we provide a unified approach for designing FxT gradient f...The accelerated method in solving optimization problems has always been an absorbing topic.Based on the fixedtime(FxT)stability of nonlinear dynamical systems,we provide a unified approach for designing FxT gradient flows(FxTGFs).First,a general class of nonlinear functions in designing FxTGFs is provided.A unified method for designing first-order FxTGFs is shown under Polyak-Łjasiewicz inequality assumption,a weaker condition than strong convexity.When there exist both bounded and vanishing disturbances in the gradient flow,a specific class of nonsmooth robust FxTGFs with disturbance rejection is presented.Under the strict convexity assumption,Newton-based FxTGFs is given and further extended to solve time-varying optimization.Besides,the proposed FxTGFs are further used for solving equation-constrained optimization.Moreover,an FxT proximal gradient flow with a wide range of parameters is provided for solving nonsmooth composite optimization.To show the effectiveness of various FxTGFs,the static regret analyses for several typical FxTGFs are also provided in detail.Finally,the proposed FxTGFs are applied to solve two network problems,i.e.,the network consensus problem and solving a system linear equations,respectively,from the perspective of optimization.Particularly,by choosing component-wisely sign-preserving functions,these problems can be solved in a distributed way,which extends the existing results.The accelerated convergence and robustness of the proposed FxTGFs are validated in several numerical examples stemming from practical applications.展开更多
This paper presents a topology optimization approach for the surface flows on variable design domains.Via this approach,the matching between the pattern of a surface flow and the 2-manifold used to define the pattern ...This paper presents a topology optimization approach for the surface flows on variable design domains.Via this approach,the matching between the pattern of a surface flow and the 2-manifold used to define the pattern can be optimized,where the 2-manifold is implicitly defined on another fixed 2-manifold named as the base manifold.The fiber bundle topology optimization approach is developed based on the description of the topological structure of the surface flow by using the differential geometry concept of the fiber bundle.The material distribution method is used to achieve the evolution of the pattern of the surface flow.The evolution of the implicit 2-manifold is realized via a homeomorphous map.The design variable of the pattern of the surface flow and that of the implicit 2-manifold are regularized by two sequentially implemented surface-PDE filters.The two surface-PDE filters are coupled,because they are defined on the implicit 2-manifold and base manifold,respectively.The surface Navier-Stokes equations,defined on the implicit 2-manifold,are used to describe the surface flow.The fiber bundle topology optimization problem is analyzed using the continuous adjoint method implemented on the first-order Sobolev space.Several numerical examples have been provided to demonstrate this approach,where the combination of the viscous dissipation and pressure drop is used as the design objective.展开更多
基金the Joint Applied Mathematics Research Institute of Peking University and Hong Kong Baptist University.Li was also partially supported by the National Basic Research Program of China under the grant 2005CB321701The research of Tang was supported by CERG Grants of Hong Kong Research Grant Council,FRG grants of Hong Kong Baptist University,and NSAF Grant#10476032 of National Science Foundation of China.He was supported in part by the Chinese Academy of Sciences while visiting its Institute of Computational Mathematics.
文摘In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the physical domain to resolvethe relevant scales in multiscale physical systems while minimizing computationalcosts. The algorithm is a generalization of the moving mesh methods basedon harmonic mappings developed by Li et al. [J. Comput. Phys., 170 (2001), pp. 562-588, and 177 (2002), pp. 365-393]. To make 3D moving mesh simulations possible,the key is to develop an efficient mesh redistribution procedure so that this part willcost as little as possible comparing with the solution evolution part. Since the meshredistribution procedure normally requires to solve large size matrix equations, wewill describe a procedure to decouple the matrix equation to a much simpler blocktridiagonaltype which can be efficiently solved by a particularly designed multi-gridmethod. To demonstrate the performance of the proposed 3D moving mesh strategy,the algorithm is implemented in finite element simulations of fluid-fluid interface interactionsin multiphase flows. To demonstrate the main ideas, we consider the formationof drops by using an energetic variational phase field model which describesthe motion of mixtures of two incompressible fluids. Numerical results on two- andthree-dimensional simulations will be presented.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12205023,U2230401,12374056,U23A20537,11904027)。
文摘A thermodynamically complete multi-phase equation of state(EOS)applicable to both dense and porous metals at wide ranges of temperature and pressure is constructed.A standard three-term decomposition of the Helmholtz free energy as a function of specific volume and temperature is presented,where the cold component models both compression and expansion states,the thermal ion component introduces the Debye approximation and melting entropy,and the thermal electron component employs the Thomas-Fermi-Kirzhnits(TFK)model.The porosity of materials is considered by introducing the dynamic porosity coefficientαand the constitutive P-αrelation,connecting the thermodynamic properties between dense and porous systems,allowing for an accurate description of the volume decrease caused by void collapse while maintaining the quasi-static thermodynamic properties of porous systems identical to the dense ones.These models enable the EOS applicable and robust at wide ranges of temperature,pressure and porosity.A systematic evaluation of the new EOS is conducted with aluminum(Al)as an example.300 K isotherm,shock Hugoniot,as well as melting curves of both dense and porous Al are calculated,which shows great agreements with experimental data and validates the effectiveness of the models and the accuracy of parameterizations.Notably,it is for the first time Hugoniot P-σcurves up to 10~6 GPa and shock melting behaviors of porous Al are derived from analytical EOS models,which predict much lower compression limit and shock melting temperatures than those of dense Al.
文摘This paper aims to numerically explore the characteristics of unsteady cavitating flow around a NACA0015 hydrofoil,with a focus on vorticity attributes.The simulation utilizes a homogeneous mixture model coupled with a filter-based density correction turbulence model and a modified Zwart cavitation model.The study investigates the dynamic cavitation features of the thermal fluid around the hydrofoil at various incoming flow velocities.It systematically elucidates the evolution of cavitation and vortex dynamics corresponding to each velocity condition.The results indicate that with increasing incoming flow velocity,distinct cavitation processes take place in the flow field.
基金financial support from the National Natural Science Foundation of China(Grant No.41941018).
文摘The impacts of natural boulders carried by debris flows pose serious risks to the safety and reliability of structures and buildings.Natural boulders can be highly random and unpredictable.Consequently,boulder control during debris flows is crucial but difficult.Herein,an eco-friendly control system featuring anchoring natural boulders(NBs)with(negative Poisson's ratio)NPR anchor cables is proposed to form an NB-NPR baffle.A series of flume experiments are conducted to verify the effect of NB-NPR baffles on controlling debris flow impact.The deployment of NB-NPR baffles substantially influences the kinematic behavior of a debris flow,primarily in the form of changes in the depositional properties and impact intensities.The results show that the NB-NPR baffle matrix successfully controls boulder mobility and exhibits positive feedback on solid particle deposition.The NB-NPR baffle group exhibits a reduction in peak impact force ranging from 29%to 79%compared to that of the control group in the basic experiment.The NPR anchor cables play a significant role in the NB-NPR baffle by demonstrating particular characteristics,including consistent resistance,large deformation,and substantial energy absorption.The NB-NPR baffle innovatively utilizes the natural boulders in a debris flow gully by converting destructive boulders into constructive boulders.Overall,this research serves as a basis for future field experiments and applications.
基金supported by the National Natural Science Foundation of China(Grant No.51704327)
文摘In the petroleum industry,detection of multi-phase fluid flow is very important in both surface and down-hole measurements.Accurate measurement of high rate of water or gas multi-phase flow has always been an academic and industrial focus.NMR is an efficient and accurate technique for the detection of fluids;it is widely used in the determination of fluid compositions and properties.This paper is aimed to quantitatively detect multi-phase flow in oil and gas wells and pipelines and to propose an innovative method for online nuclear magnetic resonance(NMR)detection.The online NMR data acquisition,processing and interpretation methods are proposed to fill the blank of traditional methods.A full-bore straight tube design without pressure drop,a Halbach magnet structure design with zero magnetic leakage outside the probe,a separate antenna structure design without flowing effects on NMR measurement and automatic control technology will achieve unattended operation.Through the innovation of this work,the application of NMR for the real-time and quantitative detection of multi-phase flow in oil and gas wells and pipelines can be implemented.
文摘Most current lattice Boltzmann （LBM） models suffer from the deficiency that their parameters have to be obtained by fitting experimental results. In this paper, we propose a new method that integrates the molecular dynamics （MD） simulation and LBM to avoid such defect. The basic idea is to first construct a molecular model based on the actual components of the rock-fluid system, then to compute the interaction force between the rock and the fluid of different densities through the MD simulation. This calculated rock-fluid interaction force, combined with the fluid-fluid force determined from the equation of state, is then used in LBM modeling. Without parameter fitting, this study presents a new systematic approach for pore-scale modeling of multi-phase flow. We have validated this ap- proach by simulating a two-phase separation process and gas-liquid-solid three-phase contact angle. Based on an actual X-ray CT image of a reservoir core, we applied our workflow to calculate the absolute permeability of the core, vapor-liquid H20 relative permeability, and capillary pressure curves.
文摘The transient friction in channel mean flows is the sum of two contributions,i.e.,the underlying laminar flow(ULF)and the purely turbulent component(PTC),and the contributions are analyzed separately by theoretical experiments.It is found that,the transient friction may be higher or remarkably lower than that in equal-Reynolds number steady-state flows.The universal time constant for plane-parallel laminar flows is reported,and the role of the time constant in a turbulent mean flow is examined.It is shown that the time constant is related to the turbulence's frozen time.Finally,a study of the logarithmic layer during the transient flow is accomplished,which shows that the logarithmic layer is destroyed.
基金financial support from the Energize Program between the University of Texas at Austin and Southwest Research InstituteHydraulic Fracturing and Sand Control Industrial Affiliates Program at the University of Texas at Austin for financially supporting this research。
文摘Two-phase pipe flow occurs frequently in oil&gas industry,nuclear power plants,and CCUS.Reliable calculations of gas void fraction(or liquid holdup)play a central role in two-phase pipe flow models.In this paper we apply the fractional flow theory to multiphase flow in pipes and present a unified modeling framework for predicting the fluid phase volume fractions over a broad range of pipe flow conditions.Compared to existing methods and correlations,this new framework provides a simple,approximate,and efficient way to estimate the phase volume fraction in two-phase pipe flow without invoking flow patterns.Notably,existing correlations for estimating phase volume fraction can be transformed and expressed under this modeling framework.Different fractional flow models are applicable to different flow conditions,and they demonstrate good agreement against experimental data within 5%errors when compared with an experimental database comprising of 2754 data groups from 14literature sources,covering various pipe geometries,flow patterns,fluid properties and flow inclinations.The gas void fraction predicted by the framework developed in this work can be used as inputs to reliably model the hydraulic and thermal behaviors of two-phase pipe flows.
文摘Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 to 44 km from the source volcano. The blockage impounded a body of water which then released producing four stratigraphic units in the downstream debris fan. The four stratigraphic units are a boulder core comprised of locally sourced bedrock boulders and three sand-rich units including a fine-grained sand unit, a sandy pumice gravel (±basalt/hydrovolcanic tuff) unit, and a pumice pebble-bearing, crystal-rich sand unit. Hand-drilled auger holes up to ~1.6 m deep were used to obtain samples of the sand-rich units. Units were delimited using surface and down-hole observations, composition and texture, estimated density, statistical parameters of grain size, and vertical and lateral distribution of properties. Overtopping followed by rapid incision into the ash-rich pyroclastic flows progressively cleared the canyon, but a bedrock knickpoint near the head of the canyon limited the volume of debris available for transport to about 0.04 km<sup>3</sup> to 0.08 km<sup>3</sup>. Co-deposition of bedrock boulders and lithic-rich sand was followed by rapid deposition with minimal reworking of remobilized pyroclastics. Continued draining of the impounded lake sent hyperconcentrated flows onto the debris fan depositing pumice-rich gravels that graded upward to crystal-rich sands.
基金supported by the National Natural Science Foundation of China,(Grant No.42072144)Shengli Oilfield,SINOPEC,China(Nos.30200018-21-ZC0613-0030 and 30200018-20-ZC0613-0116)。
文摘Based on the theory of superimposed deformation and the regional tectonic background,the multi-phase non-coaxial superimposed structures in Junggar Basin were systematically analyzed using seismic interpretation,field outcrop observation,and paleo-stress field recovery methods according to the characteristics of the current tectonic framework.Moreover,the tectonic evolution process of the basin was reconstructed using sandbox analogue modelling technology.The results showed that the study area has experienced five phases of non-coaxial deformation with superimposition:The first phase of deformation(D1)is characterized by NNE-SSW extension during late Carboniferous to early Permian,which formed large graben,half graben and other extensional structure style around the basin.The second phase of deformation(D2)is represented by NE-SW compression during the middle to late Permian,and it comprised numerous contraction structures that developed based on D1.The basic form of the entire basin is alternating uplift and depression.The third phase of deformation(D3)is the NW-SE transpressional strike-slip in the Triassic-Jurassic,which produced numerous strike-slip structural styles in the middle part of the basin.The fourth phase of deformation(D4)is the uniform sedimentation during Cretaceous,and the fifth phase(D5)is the compression along NNE-SSW due to the North Tianshan northward thrust,which produced three rows of fold thrust belts and tear faults in the front of the mountain in the southern margin of the basin.The newly established three-dimensional tectonic evolution model shows that,based on the large number of NW-trending grabens and half grabens in the Carboniferous basement of Junggar Basin,multiple level NE trending uplifts have formed with the joint superposition of the late structural inversion and multiple stress fields.This has resulted in the current tectonic units of alternating uplifts and depressions in different directions in the study area.
基金the National Natural Science Foundation of China(No.62173049)the Open Fund of the Hubei Key Laboratory of Oil and Gas Drilling and Production Engineering(Yangtze University),YQZC202309.
文摘Conventional conductivity methods for measuring the void fraction in gas-liquid multiphase systems are typically affected by accuracy problems due to the presence of fluid flow and salinity.This study presents a novel approach for determining the void fraction based on a reciprocating dynamic conductivity probe used to measure the liquid film thickness under forced annular-flow conditions.The measurement system comprises a cyclone,a conductivity probe,a probe reciprocating device,and a data acquisition and processing system.This method ensures that the flow pattern is adjusted to a forced annular flow,thereby minimizing the influence of complex and variable gas-liquid flow patterns on the measurement results;Moreover,it determines the liquid film thickness solely according to circuit connectivity rather than specific conductivity values,thereby mitigating the impact of salinity.The reliability of the measurement system is demonstrated through laboratory experiments.The experimental results indicate that,in a range of gas phase superficial velocities 5–20 m/s and liquid phase superficial velocities 0.079–0.48 m/s,the maximum measurement deviation for the void fraction is 4.23%.
基金supported by the AFOSR grant FA9550-20-1-0055 and the NSF grant DMS-2010107.
文摘Capturing elaborated flow structures and phenomena is required for well-solved numerical flows.The finite difference methods allow simple discretization of mesh and model equations.However,they need simpler meshes,e.g.,rectangular.The inverse Lax-Wendroff(ILW)procedure can handle complex geometries for rectangular meshes.High-resolution and high-order methods can capture elaborated flow structures and phenomena.They also have strong mathematical and physical backgrounds,such as positivity-preserving,jump conditions,and wave propagation concepts.We perceive an effort toward direct numerical simulation,for instance,regarding weighted essentially non-oscillatory(WENO)schemes.Thus,we propose to solve a challenging engineering application without turbulence models.We aim to verify and validate recent high-resolution and high-order methods.To check the solver accuracy,we solved vortex and Couette flows.Then,we solved inviscid and viscous nozzle flows for a conical profile.We employed the finite difference method,positivity-preserving Lax-Friedrichs splitting,high-resolution viscous terms discretization,fifth-order multi-resolution WENO,ILW,and third-order strong stability preserving Runge-Kutta.We showed the solver is high-order and captured elaborated flow structures and phenomena.One can see oblique shocks in both nozzle flows.In the viscous flow,we also captured a free-shock separation,recirculation,entrainment region,Mach disk,and the diamond-shaped pattern of nozzle flows.
基金the support of EPIC—Energy Production Innovation Center,hosted by the University of Campinas(UNICAMP)sponsored by FAPESP—Sao Paulo Research Foundation(2017/15736—3 process)+2 种基金the support and funding from Equinor Brazil and the support of ANP(Brazil's National Oil,Natural Gas and Biofuels Agency)through the R&D levy regulationthe Center of Energy and Petroleum Studies(CEPETRO)the School of Mechanical Engineering(FEM)。
文摘This study investigates the impact of pore network characteristics on fluid flow through complex and heterogeneous porous media,providing insights into the factors affecting fluid propagation in such systems.Specifically,high-resolution or micro X-ray computed tomography(CT)imaging techniques were utilized to examine outcrop stromatolite samples of the Lagoa Salgada,considered flow analogous to the Brazilian Pre-salt carbonate reservoirs.The petrophysical results comprised two distinct stromatolite depositional facies,the columnar and the fine-grained facies.By generating pore network model(PNM),the study quantified the relationship between key features of the porous system,including pore and throat radius,throat length,coordination number,shape factor,and pore volume.The study found that the less dense pore network of the columnar sample is typically characterized by larger pores and wider and longer throats but with a weaker connection of throats to pores.Both facies exhibited less variability in the radius of the pores and throats in comparison to throat length.Additionally,a series of core flooding experiments coupled with medical CT scanning was designed and conducted in the plug samples to assess flow propagation and saturation fields.The study revealed that the heterogeneity and presence of disconnected or dead-end pores significantly impacted the flow patterns and saturation.Two-phase flow patterns and oil saturation distribution reveal a preferential and heterogeneous displacement that mainly swept displaced fluid in some regions of plugs and bypassed it in others.The relation between saturation profiles,porosity profiles,and the number of fluid flow patterns for the samples was evident.Only for the columnar plug sample was the enhancement in recovery factor after shifting to lower salinity water injection(SB)observed.
基金supported in part by the Major Project for New Generation of AI (2018AAA0100400)the National Natural Science Foundation of China (61836014,U21B2042,62072457,62006231)the InnoHK Program。
文摘Monocular 3D object detection is challenging due to the lack of accurate depth information.Some methods estimate the pixel-wise depth maps from off-the-shelf depth estimators and then use them as an additional input to augment the RGB images.Depth-based methods attempt to convert estimated depth maps to pseudo-LiDAR and then use LiDAR-based object detectors or focus on the perspective of image and depth fusion learning.However,they demonstrate limited performance and efficiency as a result of depth inaccuracy and complex fusion mode with convolutions.Different from these approaches,our proposed depth-guided vision transformer with a normalizing flows(NF-DVT)network uses normalizing flows to build priors in depth maps to achieve more accurate depth information.Then we develop a novel Swin-Transformer-based backbone with a fusion module to process RGB image patches and depth map patches with two separate branches and fuse them using cross-attention to exchange information with each other.Furthermore,with the help of pixel-wise relative depth values in depth maps,we develop new relative position embeddings in the cross-attention mechanism to capture more accurate sequence ordering of input tokens.Our method is the first Swin-Transformer-based backbone architecture for monocular 3D object detection.The experimental results on the KITTI and the challenging Waymo Open datasets show the effectiveness of our proposed method and superior performance over previous counterparts.
基金supported by the National Natural Science Foundation of China (Grant No.41925030)the Nyingchi National Sustainable Development Experimental Zone Project (2023-SYQ-007)+1 种基金the Science and Technology Research Program of Institute of Mountain Hazards and Environment, Chinese Academy of Sciences (Grant No. IMHE-ZDRW-02)the Chinese Academy of Sciences Light of West China Program
文摘Forests play an important role in controlling the formation and movement processes of debris flows.They contribute to soil stabilization,regulation of soil water content,and act as robust structures impeding the downstream progression of debris flows.On the positive side,trees,to some extent,can intercept debris flows and effectively mitigate their velocity by increasing flow resistance.On the negative side,trees may suffer damage from debris-flow hazards,characterized by the generation of substantial quantities of wood fragments and consequential ramifications such as river channel blockage,resulting in backwater rise.In extreme cases,this blockage collapse can lead to instantaneous discharge amplification,thereby adversely impacting urban safety and impeding sustainable development.Therefore,in order to grasp the effects of tree characteristics on tree failure modes,the tree failure modes and corresponding parameters,diameters at breast height(DBH)and root-soil plate size,were identified and recorded through the post-event field investigation in Keze Gully,a region prone to debrisflow events in Sichuan,China,respectively.To investigate the impact of spatial variability in tree root distribution on tree failure modes,the root crosssectional area ratio(RAR),root density(RD),root length density(RLD)and soil detachment rate(SDR)were obtained.The findings indicated that:(1)Tree characteristics reflect the interactions of debris flows and trees,and influence the tree failure modes ultimately.The root distribution characteristics influence the size and shape of the root-soil plate to affect the resistance of trees.(2)Compared to burial and abrasion,stem breakage and overturning are the predominant modes of tree failure in debris-flow hazards.Trees with a smaller DBH primarily experience stem breakage and bending,and trees with a larger DBH mostly experience overturning.(3)The root-soil plate shapes of overturned trees,affected by the root architecture and root growth range,are generally semielliptical or semicircular,and the horizontal and vertical radii increase with DBH,but the correlation between the root-soil plate’s breadth-depth ratio and DBH is low.(4)The biomass and RAR decrease with distance.The RAR distribution exhibit the order of upslope direction>downslope direction>lateral direction.The coarse root biomass significantly increases with DBH,but no clear trend in fine root biomass.(5)The roots can significantly enhance the soil erosion resistance,but the erosion resistance of coarse roots is not as significant as that of fine roots.The erosion resistance increases with DBH,and follows the order of upslope direction>downslope direction>lateral direction.The results could provide new insights into the influences of tree and root distribution characteristics on tree failure modes during debris flows.
基金support provided by Shandong Provincial Science and Technology Plan Project(No.2023TSGC0625)Natural Resources Defense Council(NRDC,K94).
文摘A model is proposed for liquid film profile prediction in gas-liquid two-phase flow,which is able to provide the film thickness along the circumferential direction and the pressure gradient in the flow direction.A two-fluid model is used to calculate both gas and liquid phases’flow characteristics.The secondary flow occurring in the gas phase is taken into account and a sailing boat mechanism is introduced.Moreover,energy conservation is applied for obtaining the liquid film thickness distribution along the circumference.Liquid film thickness distribution is calculated accordingly for different cases;its values are compared with other models and available experimental data.As a result,the newly proposed model is tested and good performances are demonstrated.The liquid film thickness distribution in small pipes and inclined pipes is also studied,and regime transition is revealed by liquid film profile evolution.The observed inflection point demonstrates that the liquid film thickness decreases steeply along the circumference,when the circle angle ranges between 30°and 50°for gas-liquid stratified flow with small superficial velocities.
基金supported by the National Natural Science Foundation of China(21627813)。
文摘The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.
基金supported by the National Key Research and Development Program of China(2020YFA0714300)the National Natural Science Foundation of China(62003084,62203108,62073079)+3 种基金the Natural Science Foundation of Jiangsu Province of China(BK20200355)the General Joint Fund of the Equipment Advance Research Program of Ministry of Education(8091B022114)Jiangsu Province Excellent Postdoctoral Program(2022ZB131)China Postdoctoral Science Foundation(2022M720720,2023T160105).
文摘The accelerated method in solving optimization problems has always been an absorbing topic.Based on the fixedtime(FxT)stability of nonlinear dynamical systems,we provide a unified approach for designing FxT gradient flows(FxTGFs).First,a general class of nonlinear functions in designing FxTGFs is provided.A unified method for designing first-order FxTGFs is shown under Polyak-Łjasiewicz inequality assumption,a weaker condition than strong convexity.When there exist both bounded and vanishing disturbances in the gradient flow,a specific class of nonsmooth robust FxTGFs with disturbance rejection is presented.Under the strict convexity assumption,Newton-based FxTGFs is given and further extended to solve time-varying optimization.Besides,the proposed FxTGFs are further used for solving equation-constrained optimization.Moreover,an FxT proximal gradient flow with a wide range of parameters is provided for solving nonsmooth composite optimization.To show the effectiveness of various FxTGFs,the static regret analyses for several typical FxTGFs are also provided in detail.Finally,the proposed FxTGFs are applied to solve two network problems,i.e.,the network consensus problem and solving a system linear equations,respectively,from the perspective of optimization.Particularly,by choosing component-wisely sign-preserving functions,these problems can be solved in a distributed way,which extends the existing results.The accelerated convergence and robustness of the proposed FxTGFs are validated in several numerical examples stemming from practical applications.
基金Supported by National Natural Science Foundation of China (Grant No.51875545)Innovation Grant of Changchun Institute of Optics+2 种基金Fine Mechanics and Physics (CIOMP)CAS Project for Young Scientists in Basic Research of China (Grant No.YSBR-066)Science and Technology Development Program of Jilin Province of China (Grant No.SKL202302020)。
文摘This paper presents a topology optimization approach for the surface flows on variable design domains.Via this approach,the matching between the pattern of a surface flow and the 2-manifold used to define the pattern can be optimized,where the 2-manifold is implicitly defined on another fixed 2-manifold named as the base manifold.The fiber bundle topology optimization approach is developed based on the description of the topological structure of the surface flow by using the differential geometry concept of the fiber bundle.The material distribution method is used to achieve the evolution of the pattern of the surface flow.The evolution of the implicit 2-manifold is realized via a homeomorphous map.The design variable of the pattern of the surface flow and that of the implicit 2-manifold are regularized by two sequentially implemented surface-PDE filters.The two surface-PDE filters are coupled,because they are defined on the implicit 2-manifold and base manifold,respectively.The surface Navier-Stokes equations,defined on the implicit 2-manifold,are used to describe the surface flow.The fiber bundle topology optimization problem is analyzed using the continuous adjoint method implemented on the first-order Sobolev space.Several numerical examples have been provided to demonstrate this approach,where the combination of the viscous dissipation and pressure drop is used as the design objective.
