With the continuous increase of aeroengine flight ceiling(>20 km),the thin atmosphere at high altitudes and the size effect all cause the compressor component inlet Reynolds number to decrease rapidly to a critical...With the continuous increase of aeroengine flight ceiling(>20 km),the thin atmosphere at high altitudes and the size effect all cause the compressor component inlet Reynolds number to decrease rapidly to a critical value(approximately 2.0×10^(5)),and the significant transition process on the blade/endwall surface leads to the sharp degradation of compressor performance,which seriously affects the engine fuel consumption and working stability at high altitudes.In this paper,the research progress on the internal flow mechanism and flow control methods of axial compressors at low Reynolds numbers is reviewed from the aspects of quantification and prediction of performance variation,flow loss mechanism related to separation and transition,efficient transition control and flow field organization.The development trend of the low-Reynolds-number effect of axial flow compressors is noted,and the difficulties and application prospects of aerodynamic design and efficient flow control methods for compressors under low Reynolds numbers at high altitudes are discussed.展开更多
Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of un...Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.展开更多
The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the p...The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the part most related to flow ripple,and its structural design is an important topic.In this study,an analytical model for the axial piston pump flow ripple was established and verified using a numerical analysis with computational fluid dynamics(CFD)calculations.Moreover,a parametric analysis of the valve plate was performed to investigate the critical parameters and their ranges.A fast optimization method,the rotation vector optimization method(RVOM),was proposed for the valve plate design and compared with the currently used optimization methods to prove its efficiency.As a constant-pressure pump works in different states of swashplate angle,outlet pressure,and pump speed,an optimization principle for the entire working status was proposed to achieve the overall reduction performance.A test rig for an aircraft hydraulic pump was established,and validation experiments were conducted.It was determined that the optimized pump could achieve reduction at multiple working statuses,and the largest pressure pulsation reduction ratios for the typical speed and speed sweep tests reached 64.7%and 71.7%,respectively.The model and method proposed in this study are proven to be effective and accurate.展开更多
Vortex-induced vibration(VIV)of an underwater manipulator in pulsating flow presents a notable engineering problem in precise control due to the velocity variation in the flow.This study investigates the VIV response ...Vortex-induced vibration(VIV)of an underwater manipulator in pulsating flow presents a notable engineering problem in precise control due to the velocity variation in the flow.This study investigates the VIV response of an underwater manipulator subjected to pulsating flow,focusing on how different postures affect the behavior of the system.The effects of pulsating parameters and manipulator arrangement on the hydrodynamic coefficient,vibration response,motion trajectory,and vortex shedding behaviors were analyzed.Results indicated that the cross flow vibration displacement in pulsating flow increased by 32.14%compared to uniform flow,inducing a shift in the motion trajectory from a crescent shape to a sideward vase shape.In the absence of interference between the upper and lower arms,the lift coefficient of the manipulator substantially increased with rising pulsating frequency,reaching a maximum increment of 67.0%.This increase in the lift coefficient led to a 67.05%rise in the vibration frequency of the manipulator in the in-line direction.As the pulsating amplitude increased,the drag coefficient of the underwater manipulator rose by 36.79%,but the vibration frequency in the cross-flow direction decreased by 56.26%.Additionally,when the upper and lower arms remained in a state of mutual interference,the cross-flow vibration amplitudes of the upper and lower arms were approximately 1.84 and 4.82 times higher in a circular-elliptical arrangement compared to an elliptical-circular arrangement,respectively.Consequently,the flow field shifted from a P+S pattern to a disordered pattern,disrupting the regularity of the motion trajectory.展开更多
Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture...Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.展开更多
An unsteady numerical simulation is conducted to examine the dynamic runback characteristics of a water film flow driven by a boundary layer airflow over a solid surface pertinent to the dynamic glaze ice accretion pr...An unsteady numerical simulation is conducted to examine the dynamic runback characteristics of a water film flow driven by a boundary layer airflow over a solid surface pertinent to the dynamic glaze ice accretion process over aircraft wing surfaces.The multiphase flow simulation results of the wind-driven water runback(WDWR)flow are compared quantitatively with the experimental results in terms of the time-dependent variations of the water film thickness profiles and evolution of the front contact point of the runback water film flow.The underlying mechanism of the intermittent water runback behavior is elucidated by analyzing the time evolution of the airflow velocity and vorticity fields above the runback water film flow over the solid surface.To the best knowledge of the authors,the work presented here is the first successful attempt to numerically examine the transient runback characteristics of WDWR flows.It serves as an excellent benchmark case for the development of best practices to model the important micro-physical processes responsible for the transient water transport over aircraft wing surfaces.展开更多
The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled wi...The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.展开更多
A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three ...The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.展开更多
The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, ...The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.展开更多
Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous f...Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.展开更多
Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks....Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.展开更多
From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship betw...From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship between material flow and the energy intensity is useful to save energy in steel industry. Based on the concept of standard material flow diagram, all possible situations of ferric material flow in steel manufacturing process are analyzed. The expressions of the influence of material flow deviated from standard material flow diagram on energy consumption are put forward.展开更多
In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes eq...In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes equation of fluid mechanics and Maxwell＇s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different a, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann ：number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.展开更多
The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is use...The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field ...The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field and thermal radiation.Initially,a similarity transformation is used to convert momentum and energy conservation equations in partial differential forms into non-linear ordinary differential equations (ODE) applying meaningful boundary conditions.In order to obtain the non-linear ODEs analytically,Galerkin method (GM) is employed.Subsequently,the ODEs are also solved by a reliable numerical solution.In order to test the accuracy,precision and reliability of the analytical method,results of the analytical analysis are compared with the numerical results.With respect to the comparisons,fairly good compatibilities with insignificant errors are observed.Eventually,the impacts of effective parameters including magnetic and radiation parameters and nanofluid volume fraction on the velocity,skin friction coefficient and Nusselt number distributions are comprehensively described.Based on the results,it is revealed that with increasing the role of magnetic force,velocity profile,skin friction coefficient and thermal performance descend.Radiation parameter has insignificant influence on velocity profile while it obviously has augmentative and decreasing effects on skin friction and Nusselt number,respectively.展开更多
The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte ...The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.展开更多
On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through ...On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through the homogeneous reservoir. It is indicate that the pressure distribution of dimensionless reservoir and bottom hole in Laplace space, which take on the radial flow, also shows similar structure, and the internal relationship between the above solutions were illustrated in detail.展开更多
基金co-supported by the National Natural Science Foundation of China(No.52306053)the Science Center for Gas Turbine Project,China(No.P2022-B-Ⅱ-005-001)the National Science and Technology Major Project of China(No.2017-Ⅱ-0010-0024)。
文摘With the continuous increase of aeroengine flight ceiling(>20 km),the thin atmosphere at high altitudes and the size effect all cause the compressor component inlet Reynolds number to decrease rapidly to a critical value(approximately 2.0×10^(5)),and the significant transition process on the blade/endwall surface leads to the sharp degradation of compressor performance,which seriously affects the engine fuel consumption and working stability at high altitudes.In this paper,the research progress on the internal flow mechanism and flow control methods of axial compressors at low Reynolds numbers is reviewed from the aspects of quantification and prediction of performance variation,flow loss mechanism related to separation and transition,efficient transition control and flow field organization.The development trend of the low-Reynolds-number effect of axial flow compressors is noted,and the difficulties and application prospects of aerodynamic design and efficient flow control methods for compressors under low Reynolds numbers at high altitudes are discussed.
基金supported by the National Natural Science Foundation of China(No.92252201)the Fundamental Research Funds for the Central Universitiesthe Academic Excellence Foundation of Beihang University(BUAA)for PhD Students。
文摘Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.
基金Supported by National Natural Science Foundation of China(Grant No.51975025)National Key Research and Development Program of China(Grant No.2019YFB2004500)。
文摘The flow ripple caused by an axial piston pump may lead to pipe vibrations and lower hydraulic component reliability,which are of particular concern in hydraulic systems.The valve plate of the pump is considered the part most related to flow ripple,and its structural design is an important topic.In this study,an analytical model for the axial piston pump flow ripple was established and verified using a numerical analysis with computational fluid dynamics(CFD)calculations.Moreover,a parametric analysis of the valve plate was performed to investigate the critical parameters and their ranges.A fast optimization method,the rotation vector optimization method(RVOM),was proposed for the valve plate design and compared with the currently used optimization methods to prove its efficiency.As a constant-pressure pump works in different states of swashplate angle,outlet pressure,and pump speed,an optimization principle for the entire working status was proposed to achieve the overall reduction performance.A test rig for an aircraft hydraulic pump was established,and validation experiments were conducted.It was determined that the optimized pump could achieve reduction at multiple working statuses,and the largest pressure pulsation reduction ratios for the typical speed and speed sweep tests reached 64.7%and 71.7%,respectively.The model and method proposed in this study are proven to be effective and accurate.
基金Supported by the National Natural Science Foundation of China(No.51905211)A Project of the“20 Regulations for New Universities”Funding Program of Jinan(No.202228116).
文摘Vortex-induced vibration(VIV)of an underwater manipulator in pulsating flow presents a notable engineering problem in precise control due to the velocity variation in the flow.This study investigates the VIV response of an underwater manipulator subjected to pulsating flow,focusing on how different postures affect the behavior of the system.The effects of pulsating parameters and manipulator arrangement on the hydrodynamic coefficient,vibration response,motion trajectory,and vortex shedding behaviors were analyzed.Results indicated that the cross flow vibration displacement in pulsating flow increased by 32.14%compared to uniform flow,inducing a shift in the motion trajectory from a crescent shape to a sideward vase shape.In the absence of interference between the upper and lower arms,the lift coefficient of the manipulator substantially increased with rising pulsating frequency,reaching a maximum increment of 67.0%.This increase in the lift coefficient led to a 67.05%rise in the vibration frequency of the manipulator in the in-line direction.As the pulsating amplitude increased,the drag coefficient of the underwater manipulator rose by 36.79%,but the vibration frequency in the cross-flow direction decreased by 56.26%.Additionally,when the upper and lower arms remained in a state of mutual interference,the cross-flow vibration amplitudes of the upper and lower arms were approximately 1.84 and 4.82 times higher in a circular-elliptical arrangement compared to an elliptical-circular arrangement,respectively.Consequently,the flow field shifted from a P+S pattern to a disordered pattern,disrupting the regularity of the motion trajectory.
基金supported by the Japan Society for the Promotion of Science(JSPS)KAKENHI(Grant Nos.JP23KK0182,JP23K26356,and JP24K00971).
文摘Stony debris flows,characterized by coarse boulders embedded in a sediment-laden matrix,greatly amplify destructive potential by altering flow dynamics and impact forces.Conventional single-phase particle-fluidmixture models often struggle to capture the complexities introduced by coarse boulders and multi-phase interactions,while strong-coupling methods can be computationally prohibitive for practical hazard assessments.In this study,we propose a semi-hybrid,fully resolved coupling numerical framework for modeling boulder-laden debris flows.This framework conceptualizes debris flows as a composite system comprising a continuous viscous fluidphase(including finesediments)and a discrete phase of arbitrarily shaped coarse particles.The continuous phase is treated as a generalized nonlinear Coulomb-viscoplastic fluidusing the smoothed particle hydrodynamics(SPH)method,while coarse particles are modeled via the distributed contact discrete element method(DCDEM).These two phases are coupled through an efficienttwo-way resolved scheme,ensuring accurate simulation of flow-boulder interactions within a unifiedtimeframe.We validate the proposed method against two physical experiments:(1)gravity-driven concrete flows and(2)debris flowinteracting with slit-type barriers.Results confirmthe method's robustness in accurately capturing fluid-solid-structureinteractions and deposition processes.Its capabilities are further showcased through the simulation of a stony debris-flowevent inWenchuan County,China,highlighting its promise for real-world engineering applications and validating the effectiveness of the existing cascade dam system in mitigating debrisflowimpact and energy dissipation.
基金supported by the National Science Foundation(NSF)of the USA(Grant Nos.TIP-2140489,CBET-2313310,and CBET-2415347).
文摘An unsteady numerical simulation is conducted to examine the dynamic runback characteristics of a water film flow driven by a boundary layer airflow over a solid surface pertinent to the dynamic glaze ice accretion process over aircraft wing surfaces.The multiphase flow simulation results of the wind-driven water runback(WDWR)flow are compared quantitatively with the experimental results in terms of the time-dependent variations of the water film thickness profiles and evolution of the front contact point of the runback water film flow.The underlying mechanism of the intermittent water runback behavior is elucidated by analyzing the time evolution of the airflow velocity and vorticity fields above the runback water film flow over the solid surface.To the best knowledge of the authors,the work presented here is the first successful attempt to numerically examine the transient runback characteristics of WDWR flows.It serves as an excellent benchmark case for the development of best practices to model the important micro-physical processes responsible for the transient water transport over aircraft wing surfaces.
基金Project (10964004) supported by the National Natural Science Foundation of ChinaProject (20070731001) supported by Research Fund for the Doctoral Program of ChinaProject (096RJZA104) supported by the Natural Science Foundation of Gansu Province,China
文摘The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
文摘The hydrodynamic characteristics of a rigid, single, circular cylinder in a three dimensional, incompressible, uniform cross flow were calculated using the large-eddy simulation method of CFX5. Solutions to the three dimensional N-S equations were obtained by the finite volume method. The focus of this numerical simulation was to research the characteristics of pressure distribution （drag and litt forces） and vortex tubes at high Reynolds numbers. The results of the calculations showed that the forces at every section in the spanwise direction of the cylinder were symmetrical about the middle section and smaller than the forces calculated in two dimensional cases. Moreover, the flow around the cylinder obviously presents three dimensional characteristics.
基金Supported by the Aeronautical Science Foundation of China(20061453020)Foundation for Basic Research of Northwestern Polytechnical University(03)~~
文摘The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.
基金Supported by the National Natural Science Foundation of China(10372043,11172134)the Fundingof Jiangsu Innovation Program for Graduate Education(CXZZ11-0192)~~
文摘Gridless method is developed for unsteady viscous flows involving moving boundaries. The point distri- bution of gridless method is implemented in an isotropic or anisotropic way according to the features of viscous flows. In the area far away from the body, the traditional cloud of isotropic points is used, while in the adjacent area, the cloud of anisotropic points is distributed. In this way, the point spacing normal to the wall can be small enough for simulating the boundary layer, and meanwhile, the total number of points in the computational do- main can be controlled due to large spacing in other tangential direction through the anisotropic way. A fast mov- ing technique of clouds of points at each time-step is presented based on the attenuation law of disturbed motion for unsteady flows involving moving boundaries. In the mentioned cloud of points, a uniform weighted least- square curve fit method is utilized to discretize the spatial derivatives of the Navier-Stokes equations. The pro- posed gridless method, coupled with a dual time-stepping method and the Spalart-Allmaras turbulence model, is implemented for the Navier-Stokes equations. The computational results of unsteady viscous flows around a NLR7301 airfoil with an oscillating flap and a pitching NACA0012 airfoil are presented in a good agreement with the available experimental data.
文摘Two methods based on a slight modification of the regular traffic assignmentalgorithms are proposed to directly compute turn flows instead of estimating them from link flows orobtaining them by expanding the networks. The first one is designed on the path-turn incidencerelationship, and it is similar to the computational procedure of link flows. It applies to thetraffic assignment algorithms that can provide detailed path structures. The second utilizes thelink-turn incidence relationship and the conservation of flow on links, a law deriving from thisrelationship. It is actually an improved version of Dial's logit assignment algorithm. The proposedapproaches can avoid the shortcomings both of the estimation methods, e. g. Furness's model andFrator's model, and of the network-expanding method in precision, stability and computation scale.Finally, they are validated by numerical examples.
基金Item Sponsored by National Basic Research Programof China (200002600)
文摘From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship between material flow and the energy intensity is useful to save energy in steel industry. Based on the concept of standard material flow diagram, all possible situations of ferric material flow in steel manufacturing process are analyzed. The expressions of the influence of material flow deviated from standard material flow diagram on energy consumption are put forward.
文摘In this study, the effects of magnetic field and nanoparticle on the Jeffery- Hamel flow are studied using a powerful analytical method called the Adomian decomposition method （ADM）. The traditional Navier-Stokes equation of fluid mechanics and Maxwell＇s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different a, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann ：number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.
基金The project supported by the National Natural Science Foundation of China (50506028)Action Scheme for Invigorating Education Towards the twenty-first century
文摘The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
文摘The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field and thermal radiation.Initially,a similarity transformation is used to convert momentum and energy conservation equations in partial differential forms into non-linear ordinary differential equations (ODE) applying meaningful boundary conditions.In order to obtain the non-linear ODEs analytically,Galerkin method (GM) is employed.Subsequently,the ODEs are also solved by a reliable numerical solution.In order to test the accuracy,precision and reliability of the analytical method,results of the analytical analysis are compared with the numerical results.With respect to the comparisons,fairly good compatibilities with insignificant errors are observed.Eventually,the impacts of effective parameters including magnetic and radiation parameters and nanofluid volume fraction on the velocity,skin friction coefficient and Nusselt number distributions are comprehensively described.Based on the results,it is revealed that with increasing the role of magnetic force,velocity profile,skin friction coefficient and thermal performance descend.Radiation parameter has insignificant influence on velocity profile while it obviously has augmentative and decreasing effects on skin friction and Nusselt number,respectively.
文摘The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.
文摘On the basis of similar structure of solutions of ordinary differential equation (ODE) boundary value problem, the similar construction method was put forward by solving problems of fluid flow in porous media through the homogeneous reservoir. It is indicate that the pressure distribution of dimensionless reservoir and bottom hole in Laplace space, which take on the radial flow, also shows similar structure, and the internal relationship between the above solutions were illustrated in detail.
