A modification to the PANS(partially averaged Navier-Stokes) model is proposed to simulate unsteady cavitating flows. In the model, the parameter fk is modified to vary as a function of the ratios between the water de...A modification to the PANS(partially averaged Navier-Stokes) model is proposed to simulate unsteady cavitating flows. In the model, the parameter fk is modified to vary as a function of the ratios between the water density and the mixture density in the local flows. The objective of this study is to validate the modified model and further understand the interaction between turbulence and cavitation around a Clark-Y hydrofoil. The comparisons between the numerical and experiment results show that the modified model can be improved to predict the cavity evolution, vortex shedding frequency and the lift force fluctuating in time fairly well, as it can effectively modulate the eddy viscosity in the cavitating region and various levels of physical turbulent fluctuations are resolved. In addition, from the computational results, it is proved that cavitation phenomenon physically influences the turbulent level, especially by the vortex shedding behaviors. Also, the mean u-velocity profiles demonstrate that the attached cavity thickness can alter the local turbulent shear layer.展开更多
Effective management of a river reach requires a sound understanding of flow and sediment transport generated by varying natural and artificial runoff conditions. Flow and sediment transport within the Ning-Meng reach...Effective management of a river reach requires a sound understanding of flow and sediment transport generated by varying natural and artificial runoff conditions. Flow and sediment transport within the Ning-Meng reach of the Yellow River（NMRYR）, northern China are controlled by a complex set of factors/processes, mainly including four sets of factors：（1） aeolian sediments from deserts bordering the main stream;（2） inflow of water and sediment from numerous tributaries;（3） impoundment of water by reservoir/hydro-junction; and（4） complex diversion and return of irrigation water. In this study, the 1-D flow ＆ sediment transport model developed by the Yellow River Institute of Hydraulic Research was used to simulate the flow and sediment transport within the NMRYR from 2001 to 2012. All four sets of factors that primarily control the flow and sediment transport mentioned above were considered in this model. Compared to the measured data collected from the hydrological stations along the NMRYR, the simulated flow and sediment transport values were generally acceptable, with relative mean deviation between measured and simulated values of 〈15%. However, simulated sediment concentration and siltation values within two sub-reaches（i.e., Qingtongxia Reservoir to Bayan Gol Hydrological Station and Bayan Gol Hydrological Station to Toudaoguai Hydrological Station） for some periods exhibited relatively large errors（the relative mean deviations between measured and simulated values of 18% and 25%, respectively）. These errors are presumably related to the inability to accurately determine the quantity of aeolian sediment influx to the river reach and the inflow of water from the ten ephemeral tributaries. This study may provide some valuable insights into the numerical simulations of flow and sediment transport in large watersheds and also provide a useful model for the effective management of the NMRYR.展开更多
Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assu...Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assumptions and approximations. In this kind of groundwater systems, studying contamination spreading and analyzing risk are challenging tasks. The main difficulty stems from determining both the travel times and the maximum distances covered by pollutants. In this context, the risk of contamination in the deep carbonate aquifer of the Augusta coastal area is presented. We used a geostatistical approach and numerical codes (MODFLOW-2000, MT3DMS) to reconstruct the complex geological framework of the study area where several contamination scenarios of hypothetical point source in a risk assessment framework were simulated. Results on the contaminant spreading are discussed and the effect of the geological structures characterizing the zone under study, namely horst and graben, are described.展开更多
Cavitating flows inside a diesel injection nozzle hole were simulated using a two-fluid model. Attention was focused on the complex cavitation processes and flow characteristics under constant inlet pressure and fluct...Cavitating flows inside a diesel injection nozzle hole were simulated using a two-fluid model. Attention was focused on the complex cavitation processes and flow characteristics under constant inlet pressure and fluctuant inlet pressure modes. To validate the two-fluid model, model predictions were compared with the experimental data available in the literatures, and good agreement was achieved. The numerical results show that the appearance of supercavitation in the diesel nozzle hole induces obvious changes of flow field structures and exit flow conditions. The distributions of liquid phase turbulent kinetic energy and exit velocity profiles corresponding to the supercavitation regime indicate the potential for promoting the primary breakup of a diesel jet. Furthermore, the upstream pressure fluctuations significantly influence the cavitation processes. Both partial cavitation and supercavitation show unsteady behaviors as the rapid rise or fall of upstream pressure.展开更多
Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster managem...Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.展开更多
In the present paper,the unsteady cavitating turbulent flow over the twisted NACA66 hydrofoil is investigated based on an modified shear stress transfer k-ωpartially averaged Navier-Stokes(MSST PANS)model,i.e.,new MS...In the present paper,the unsteady cavitating turbulent flow over the twisted NACA66 hydrofoil is investigated based on an modified shear stress transfer k-ωpartially averaged Navier-Stokes(MSST PANS)model,i.e.,new MSST PANS(NMSST PANS)model,where the production term of kinetic energy in the turbulence model is modified with helicity.Compared with the experimental data,cavitation evolution and its characteristic frequency are satisfactorily predicted by the proposed NMSST PANS model.It is revealed that the interaction among the main flow,the reentrant jets,and sheet cavitation causes the formation of the primary shedding cavity near the mid-span and the secondary shedding cavity at each side of the twisted hydrofoil,and further induces the remarkable pressure gradient around shedding cavities.Along with the development of the primary and the secondary shedding cavities,the great pressure gradient associated with large cavity volume variation promotes the vortical flow generation and the spatial deformation of vortex structure during cavitation evolution,and results in the primary and the secondary U-type vortices.Further,dynamic mode decomposition(DMD)analysis is utilized to confirm the interaction among the main flow,the main reentrant jet and two side reentrant jets,and cavitation.These results indicate that the proposed NMSST PANS model is suitable to simulate the complicated cavitating turbulent flow for various engineering applications.展开更多
The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the n...The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.展开更多
Mud flows are common phenomena in mountainous areas,which can threaten human safety and cause property losses under certain extreme circumstances.Studying the dynamic characteristics of mud flows,especially in the ver...Mud flows are common phenomena in mountainous areas,which can threaten human safety and cause property losses under certain extreme circumstances.Studying the dynamic characteristics of mud flows,especially in the vertical direction,is helpful for risk reduction and hazard mitigation.In this study,a 2D depth-resolved numerical model based on Herschel-Bulkley rheology was developed to study the vertical structures of unsteady mud flows with a free-surface.The numerical model was solved by the projection method,and the free surface of mud flows was captured through the VOF method.To fully validate this new model,a series of laboratory experiments involving dam break mud flows were conducted,and the mud flow heights,bottom pressures and envelopes of mud residuum were measured.The numerical model proposed in this study was first validated by the steady-state solution for uniform flows of Herschel-Bulkley fluid on an inclined plane.Additionally,the simulated and measured mud flow heights,bottom pressures at different x locations and envelopes with different bed slopes showed good agreement.Furthermore,the numerical results for a Herschel-Bulkley fluid dam break flow were used to validate the proposed model,which further revealed good agreements.After that,the scenarios in which mud flows impact on a structure were numerically studied,and the vertical profiles of the front velocity and impact pressure on the structure were analyzed and discussed.The results show that a plug layer was formed in the mud flow under unsteady and nonuniform flow conditions,and the impact pressure on the structure was dominated by the dynamic pressure.In addition,the vertical position with the maximum impact pressure acting on the structure was not at the bottom or the surface of the mud flows,and the normalized vertical position rose as the yield stress and consistency coefficient increase for Herschel-Bulkley fluids.展开更多
Cavitation typically occurs when the fluid pressure is lower than the vapor pressure in a local thermodynamic state,and the flow is frequently unsteady and turbulent.The Reynolds-Averaged Navier-Stokes（RANS）approach...Cavitation typically occurs when the fluid pressure is lower than the vapor pressure in a local thermodynamic state,and the flow is frequently unsteady and turbulent.The Reynolds-Averaged Navier-Stokes（RANS）approach has been popular for turbulent flow computations.The most widely used ones,such as the standard k-εmodel,have well-recognized deficiencies when treating time dependent flow field.To identify ways to improve the predictive capability of the current RANS-based engineering turbulence closures,conditional averaging is adopted for the Navier-Stokes equation,and one more parameter,based on the filter size,is introduced into the k-εmodel.In the Partially Averaged Navier-Stokes（PANS）model,the filter width is mainly controlled by the ratio of unresolved-to-total kinetic energy1f.This model is assessed in unsteady cavitating flows over a Clark-Y hydrofoil.From the experimental validations regarding the forces,frequencies,cavity visualizations and velocity distributions,the PANS model is shown to improve the predictive capability considerably,in comparison to the standard k-ε model,and also,it is observed the value of1f in the PANS model has substantial influence on the predicting result.As the filter width1f is decreased,the PANS model can effectively reduce the eddy viscosity near the closure region which can significantly influence the capture of the detach cavity,and this model can reproduce the time-averaged velocity quantitatively around the hydrofoil.展开更多
Computational Fluid Dynamics (CFD) simulations of cavitating flow through water hydraulic poppet valves were performed using advanced RNG k-epsilon turbulence model. The flow was turbulent, incompressible and unsteady...Computational Fluid Dynamics (CFD) simulations of cavitating flow through water hydraulic poppet valves were performed using advanced RNG k-epsilon turbulence model. The flow was turbulent, incompressible and unsteady, for Reynolds numbers greater than 43 000. The working fluid was water, and the structure of the valve was simplified as a two dimensional axisymmetric geometrical model. Flow field visualization was numerically achieved. The effects of inlet velocity, outlet pressure, opening size as well as poppet angle on cavitation intensity in the poppet valve were numerically investigated. Experimental flow visualization was conducted to capture cavitation images near the orifice in the poppet valve with 30° poppet angle using high speed video camera. The binary cavitating flow field distribution obtained from digital processing of the original cavitation image showed a good agreement with the numerical result.展开更多
On the basis of the full velocity difference （FVD） model, an improved multiple car-following （MCF） model is proposed by taking into account multiple information inputs from preceding vehicles. The linear stability...On the basis of the full velocity difference （FVD） model, an improved multiple car-following （MCF） model is proposed by taking into account multiple information inputs from preceding vehicles. The linear stability condition of the model is obtained by using the linear stability theory. Through nonlinear analysis, a modified Korteweg-de Vries equation is constructed and solved. The traffic jam can thus be described by the klnk-antikink soliton solution for the mKdV equation. The improvement of this new model over the previous ones lies in the fact that it not only theoretically retains many strong points of the previous ones, but also performs more realistically than others in the dynamical evolution of congestion. Furthermore, numerical simulation of traffic dynamics shows that the proposed model can avoid the disadvantage of negative velocity that occurs at small sensitivity coefficients λ in the FVD model by adjusting the information on the multiple leading vehicles. No collision occurs and no unrealistic deceleration appears in the improved model.展开更多
A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The res...A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.展开更多
Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which ...Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which the wet steam is regarded as mixture comprisingtwo coupled systems: the vapor phase and the liquid phase. These two systems are both described byconservation equations. High resolution TVD scheme is employed to capture condensing phenomena inwet steam flow. This model has been validated by numerical simulations of condensing flows in 1D and2D nozzles. Compared with experimental data, a good agreement is observed. This Eulerian/Eulerianmodel can be extended to 3D calculation of condensing flow.展开更多
The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid ada...The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid adaptation techniques. The investigation focuses on two configurations, related to hypersonic and subsonic experiments. The applicability tests of different turbu- lence models are conducted on the level of two-equation models calculating the steady state solution of the Reynolds-averaged Navier-Stokes（RANS） equations. All used models, the original Wilcox k-co, the Menter shear-stress transport （SST） and the ex- plicit algebraic Reynolds stress model（EARSM） formulation, predict an asymmetric base flow in both cases caused by the support of the models. A comparison with preliminary experimental results indicates a preference for the SST and EARSM results over the results from the older k-co model. Sensitivity studies show no significant influence of the grid topology or the location of the laminar to turbulent transition on the base flow field, but a strong influence of even small angles of attack is reported from the related experiments.展开更多
A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-por...A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-porosity technique was used to represent the solidification,and Darcy law was adopted to simulate the flow in mushy region.The prediction for steel flow and solidification was evaluated by the comparison of two turbulence models.It is found that both Lam-Bremhorst low-Reynolds-number and high-Reynolds-number k-ε models predict the same trend of the steel flow and temperature distribution.However,due to the effect of turbulent flow on heat transfer,the low-Reynolds-number turbulence model predicts longer penetration depth of molten steel in sub-mold region,less shell growth and higher shell surface temperature at the narrow face compared with standard k-ε model.展开更多
ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged u...ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.展开更多
Based on Xue＇s lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is o...Based on Xue＇s lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is obtained by employing the linear stability theory. The density wave is investigated analytically with the perturbation method. The results show that the occurrence of traffic jamming transitions can be described by the kink-antikink solution of the modified Korteweg-de Vries （mKdV） equation. The simulation results are in good agreement with the analytical results, showing that the stability of traffic flow can be enhanced when the relative current of next-nearest-neighbour sites ahead is considered.展开更多
A mathematical modei of two-dimensional turbulent gas-particle twophase flow based on the modified diffusion flux modei (DFM) and a numerical simulation method to analyze the gas-particle flow structures are developed...A mathematical modei of two-dimensional turbulent gas-particle twophase flow based on the modified diffusion flux modei (DFM) and a numerical simulation method to analyze the gas-particle flow structures are developed. The modified diffusion flux modei, in which the acceleration due to various forces is taken into account for the calculation of the diffusion velocity of particles, is applicable to the analysis of multi-dimensional gas-particle two-phase turbulent flow. In order to verify its accuracy and efficiency, the numerical simulation by DFM is compared with experimental studies and the prediction by k-ε-kp two-fluid modei, which shows a reasonable agreement. It is confirmed that the modified diffusion flux modei is suitable for simulating the multi-dimensional gas-particle two-phase flow.展开更多
Some flow-like landslides tend to lose materials while moving on a relatively dry sliding surface.This phenomenon is called bed deposition.In contrast to the bed entrainment phenomenon,bed deposition is relatively poo...Some flow-like landslides tend to lose materials while moving on a relatively dry sliding surface.This phenomenon is called bed deposition.In contrast to the bed entrainment phenomenon,bed deposition is relatively poorly understood.Therefore,an improved depth-averaged model is proposed to quantify this phenomenon.The deposition depth is calculated according to the momentum conservation of the deposited mass,and the rheological property of the sliding mass on the bottom is modified considering an abrupt increase in the depth of the sliding surface after deposition.Utilizing the proposed model,the Gaolou landslide,a typical flowlike landside occurring on October 6 th,2006 with an obvious bed deposition phenomenon in Shaanxi Province of China,is simulated to investigate the influence of bed deposition on its propagation process.The results indicate that the proposed model can effectively depict the bed deposition phenomenon in the Gaolou landslide.Bed deposition dissipates part of the kinetic energy of this landslide;thus,the simulated debris inundation area would likely be overestimated when this effect is neglected.On the other hand,the thin liquefied layer formed by the bed deposition process reduces the friction energy dissipation,contributing to the high mobility of this landslide.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.11172040 and 51239005)the Beijing Municipal Natural Science Foundation(Grant No.3144043)
文摘A modification to the PANS(partially averaged Navier-Stokes) model is proposed to simulate unsteady cavitating flows. In the model, the parameter fk is modified to vary as a function of the ratios between the water density and the mixture density in the local flows. The objective of this study is to validate the modified model and further understand the interaction between turbulence and cavitation around a Clark-Y hydrofoil. The comparisons between the numerical and experiment results show that the modified model can be improved to predict the cavity evolution, vortex shedding frequency and the lift force fluctuating in time fairly well, as it can effectively modulate the eddy viscosity in the cavitating region and various levels of physical turbulent fluctuations are resolved. In addition, from the computational results, it is proved that cavitation phenomenon physically influences the turbulent level, especially by the vortex shedding behaviors. Also, the mean u-velocity profiles demonstrate that the attached cavity thickness can alter the local turbulent shear layer.
基金financially supported by the National Natural Science Foundation of China(51579113,51309111,51309113)
文摘Effective management of a river reach requires a sound understanding of flow and sediment transport generated by varying natural and artificial runoff conditions. Flow and sediment transport within the Ning-Meng reach of the Yellow River（NMRYR）, northern China are controlled by a complex set of factors/processes, mainly including four sets of factors：（1） aeolian sediments from deserts bordering the main stream;（2） inflow of water and sediment from numerous tributaries;（3） impoundment of water by reservoir/hydro-junction; and（4） complex diversion and return of irrigation water. In this study, the 1-D flow ＆ sediment transport model developed by the Yellow River Institute of Hydraulic Research was used to simulate the flow and sediment transport within the NMRYR from 2001 to 2012. All four sets of factors that primarily control the flow and sediment transport mentioned above were considered in this model. Compared to the measured data collected from the hydrological stations along the NMRYR, the simulated flow and sediment transport values were generally acceptable, with relative mean deviation between measured and simulated values of 〈15%. However, simulated sediment concentration and siltation values within two sub-reaches（i.e., Qingtongxia Reservoir to Bayan Gol Hydrological Station and Bayan Gol Hydrological Station to Toudaoguai Hydrological Station） for some periods exhibited relatively large errors（the relative mean deviations between measured and simulated values of 18% and 25%, respectively）. These errors are presumably related to the inability to accurately determine the quantity of aeolian sediment influx to the river reach and the inflow of water from the ten ephemeral tributaries. This study may provide some valuable insights into the numerical simulations of flow and sediment transport in large watersheds and also provide a useful model for the effective management of the NMRYR.
文摘Simulation of flow in fractured aquifers is a complex issue. The problem of mathematical modeling of highly heterogeneous porous media, typical of natural systems, couples with the needs making proper simplifying assumptions and approximations. In this kind of groundwater systems, studying contamination spreading and analyzing risk are challenging tasks. The main difficulty stems from determining both the travel times and the maximum distances covered by pollutants. In this context, the risk of contamination in the deep carbonate aquifer of the Augusta coastal area is presented. We used a geostatistical approach and numerical codes (MODFLOW-2000, MT3DMS) to reconstruct the complex geological framework of the study area where several contamination scenarios of hypothetical point source in a risk assessment framework were simulated. Results on the contaminant spreading are discussed and the effect of the geological structures characterizing the zone under study, namely horst and graben, are described.
基金Supported by the National Natural Science Foundation of China, Key Project Fund-ing (Grant No. 50636040)Major State Basic Research Development Program (Grant No. 2007CB210001)
文摘Cavitating flows inside a diesel injection nozzle hole were simulated using a two-fluid model. Attention was focused on the complex cavitation processes and flow characteristics under constant inlet pressure and fluctuant inlet pressure modes. To validate the two-fluid model, model predictions were compared with the experimental data available in the literatures, and good agreement was achieved. The numerical results show that the appearance of supercavitation in the diesel nozzle hole induces obvious changes of flow field structures and exit flow conditions. The distributions of liquid phase turbulent kinetic energy and exit velocity profiles corresponding to the supercavitation regime indicate the potential for promoting the primary breakup of a diesel jet. Furthermore, the upstream pressure fluctuations significantly influence the cavitation processes. Both partial cavitation and supercavitation show unsteady behaviors as the rapid rise or fall of upstream pressure.
基金financially supported by Department of Space,India(Grant No.ISRO/RES/4/663/18-19)。
文摘Debris flows are rapid mass movements with a mixture of rock,soil and water.High-intensity rainfall events have triggered multiple debris flows around the globe,making it an important concern from the disaster management perspective.This study presents a numerical model called debris flow simulation 2D(DFS 2D)and applicability of the proposed model is investigated through the values of the model parameters used for the reproduction of an occurred debris flow at Yindongzi gully in China on 13 August 2010.The model can be used to simulate debris flows using three different rheologies and has a userfriendly interface for providing the inputs.Using DFS 2D,flow parameters can be estimated with respect to space and time.The values of the flow resistance parameters of model,dry-Coulomb and turbulent friction,were calibrated through the back analysis and the values obtained are 0.1 and 1000 m/s^(2),respectively.Two new methods of calibration are proposed in this study,considering the crosssectional area of flow and topographical changes induced by the debris flow.The proposed methods of calibration provide an effective solution to the cumulative errors induced by coarse-resolution digital elevation models(DEMs)in numerical modelling of debris flows.The statistical indices such as Willmott's index of agreement,mean-absolute-error,and normalized-root-mean-square-error of the calibrated model are 0.5,1.02 and 1.44,respectively.The comparison between simulated and observed values of topographic changes indicates that DFS 2D provides satisfactory results and can be used for dynamic modelling of debris flows.
基金Project supported by the National Natural Science Foundation of China (Grant No.52336001).
文摘In the present paper,the unsteady cavitating turbulent flow over the twisted NACA66 hydrofoil is investigated based on an modified shear stress transfer k-ωpartially averaged Navier-Stokes(MSST PANS)model,i.e.,new MSST PANS(NMSST PANS)model,where the production term of kinetic energy in the turbulence model is modified with helicity.Compared with the experimental data,cavitation evolution and its characteristic frequency are satisfactorily predicted by the proposed NMSST PANS model.It is revealed that the interaction among the main flow,the reentrant jets,and sheet cavitation causes the formation of the primary shedding cavity near the mid-span and the secondary shedding cavity at each side of the twisted hydrofoil,and further induces the remarkable pressure gradient around shedding cavities.Along with the development of the primary and the secondary shedding cavities,the great pressure gradient associated with large cavity volume variation promotes the vortical flow generation and the spatial deformation of vortex structure during cavitation evolution,and results in the primary and the secondary U-type vortices.Further,dynamic mode decomposition(DMD)analysis is utilized to confirm the interaction among the main flow,the main reentrant jet and two side reentrant jets,and cavitation.These results indicate that the proposed NMSST PANS model is suitable to simulate the complicated cavitating turbulent flow for various engineering applications.
基金Supported by the National Natural Science Foundation of China （50606026, 50736006）.
文摘The second-order moment combustion model, proposed by the authors is validated using the direct numerical simulation （DNS） of incompressible turbulent reacting channel flows. The instantaneous DNS results show the near-wall strip structures of concentration and temperature fluctuations. The DNS statistical results give the budget of the terms in the correlation equations, showing that the production and dissipation terms are most important. The DNS statistical data are used to validate the closure model in RANS second-order moment （SOM） combustion model. It is found that the simulated diffusion and production terms are in agreement with the DNS data in most flow regions, except in the near-wall region, where the near-wall modification should be made, and the closure model for the dissipation term needs further improvement. The algebraic second-order moment （ASOM） combustion model is well validated by DNS.
基金Fundings from National Natural Science Foundation of China(Grant No.41941017 and U20A20112)Sichuan Science and Technology Program(Grant No.2021YFH0009)are gratefully acknowledged。
文摘Mud flows are common phenomena in mountainous areas,which can threaten human safety and cause property losses under certain extreme circumstances.Studying the dynamic characteristics of mud flows,especially in the vertical direction,is helpful for risk reduction and hazard mitigation.In this study,a 2D depth-resolved numerical model based on Herschel-Bulkley rheology was developed to study the vertical structures of unsteady mud flows with a free-surface.The numerical model was solved by the projection method,and the free surface of mud flows was captured through the VOF method.To fully validate this new model,a series of laboratory experiments involving dam break mud flows were conducted,and the mud flow heights,bottom pressures and envelopes of mud residuum were measured.The numerical model proposed in this study was first validated by the steady-state solution for uniform flows of Herschel-Bulkley fluid on an inclined plane.Additionally,the simulated and measured mud flow heights,bottom pressures at different x locations and envelopes with different bed slopes showed good agreement.Furthermore,the numerical results for a Herschel-Bulkley fluid dam break flow were used to validate the proposed model,which further revealed good agreements.After that,the scenarios in which mud flows impact on a structure were numerically studied,and the vertical profiles of the front velocity and impact pressure on the structure were analyzed and discussed.The results show that a plug layer was formed in the mud flow under unsteady and nonuniform flow conditions,and the impact pressure on the structure was dominated by the dynamic pressure.In addition,the vertical position with the maximum impact pressure acting on the structure was not at the bottom or the surface of the mud flows,and the normalized vertical position rose as the yield stress and consistency coefficient increase for Herschel-Bulkley fluids.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 50679001, 50979004)the Fundation from State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology
文摘Cavitation typically occurs when the fluid pressure is lower than the vapor pressure in a local thermodynamic state,and the flow is frequently unsteady and turbulent.The Reynolds-Averaged Navier-Stokes（RANS）approach has been popular for turbulent flow computations.The most widely used ones,such as the standard k-εmodel,have well-recognized deficiencies when treating time dependent flow field.To identify ways to improve the predictive capability of the current RANS-based engineering turbulence closures,conditional averaging is adopted for the Navier-Stokes equation,and one more parameter,based on the filter size,is introduced into the k-εmodel.In the Partially Averaged Navier-Stokes（PANS）model,the filter width is mainly controlled by the ratio of unresolved-to-total kinetic energy1f.This model is assessed in unsteady cavitating flows over a Clark-Y hydrofoil.From the experimental validations regarding the forces,frequencies,cavity visualizations and velocity distributions,the PANS model is shown to improve the predictive capability considerably,in comparison to the standard k-ε model,and also,it is observed the value of1f in the PANS model has substantial influence on the predicting result.As the filter width1f is decreased,the PANS model can effectively reduce the eddy viscosity near the closure region which can significantly influence the capture of the detach cavity,and this model can reproduce the time-averaged velocity quantitatively around the hydrofoil.
文摘Computational Fluid Dynamics (CFD) simulations of cavitating flow through water hydraulic poppet valves were performed using advanced RNG k-epsilon turbulence model. The flow was turbulent, incompressible and unsteady, for Reynolds numbers greater than 43 000. The working fluid was water, and the structure of the valve was simplified as a two dimensional axisymmetric geometrical model. Flow field visualization was numerically achieved. The effects of inlet velocity, outlet pressure, opening size as well as poppet angle on cavitation intensity in the poppet valve were numerically investigated. Experimental flow visualization was conducted to capture cavitation images near the orifice in the poppet valve with 30° poppet angle using high speed video camera. The binary cavitating flow field distribution obtained from digital processing of the original cavitation image showed a good agreement with the numerical result.
基金Project supported by the National High Tech Research and Development Program of China (Grant No 511-0910-1031)the National "10th Five-year" Science and Technique Important Program of China (Grant No 2002BA404A07)
文摘On the basis of the full velocity difference （FVD） model, an improved multiple car-following （MCF） model is proposed by taking into account multiple information inputs from preceding vehicles. The linear stability condition of the model is obtained by using the linear stability theory. Through nonlinear analysis, a modified Korteweg-de Vries equation is constructed and solved. The traffic jam can thus be described by the klnk-antikink soliton solution for the mKdV equation. The improvement of this new model over the previous ones lies in the fact that it not only theoretically retains many strong points of the previous ones, but also performs more realistically than others in the dynamical evolution of congestion. Furthermore, numerical simulation of traffic dynamics shows that the proposed model can avoid the disadvantage of negative velocity that occurs at small sensitivity coefficients λ in the FVD model by adjusting the information on the multiple leading vehicles. No collision occurs and no unrealistic deceleration appears in the improved model.
文摘A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.
基金This project is supported by National Natural Science Foundation of China (No.50176035)Doctorate Foundation of Xi'an Jiaotong University,China (No.DFXJU2000-18)
文摘Wet steam condensing flow in low-pressure steam turbine leads to efficiencylosses and blade erosions. In order to investigate this problem by numerical approach, anEulerian/Eulerian model has been developed, in which the wet steam is regarded as mixture comprisingtwo coupled systems: the vapor phase and the liquid phase. These two systems are both described byconservation equations. High resolution TVD scheme is employed to capture condensing phenomena inwet steam flow. This model has been validated by numerical simulations of condensing flows in 1D and2D nozzles. Compared with experimental data, a good agreement is observed. This Eulerian/Eulerianmodel can be extended to 3D calculation of condensing flow.
基金German Research Foundation (Deutsche Forschungs-gemeinschaft-DFG) Sonderforschungsbereich Transregio 40
文摘The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid adaptation techniques. The investigation focuses on two configurations, related to hypersonic and subsonic experiments. The applicability tests of different turbu- lence models are conducted on the level of two-equation models calculating the steady state solution of the Reynolds-averaged Navier-Stokes（RANS） equations. All used models, the original Wilcox k-co, the Menter shear-stress transport （SST） and the ex- plicit algebraic Reynolds stress model（EARSM） formulation, predict an asymmetric base flow in both cases caused by the support of the models. A comparison with preliminary experimental results indicates a preference for the SST and EARSM results over the results from the older k-co model. Sensitivity studies show no significant influence of the grid topology or the location of the laminar to turbulent transition on the base flow field, but a strong influence of even small angles of attack is reported from the related experiments.
文摘A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-porosity technique was used to represent the solidification,and Darcy law was adopted to simulate the flow in mushy region.The prediction for steel flow and solidification was evaluated by the comparison of two turbulence models.It is found that both Lam-Bremhorst low-Reynolds-number and high-Reynolds-number k-ε models predict the same trend of the steel flow and temperature distribution.However,due to the effect of turbulent flow on heat transfer,the low-Reynolds-number turbulence model predicts longer penetration depth of molten steel in sub-mold region,less shell growth and higher shell surface temperature at the narrow face compared with standard k-ε model.
基金Project(2006BAJ04B04)supported by the National Science and Technology Pillar Program in the Eleventh Five-year Plan PeriodProject(2006AA05Z229)supported by the National High Technology Research and Development Program of China+1 种基金Project supportedby the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education MinistryProject(06wk3023)supported by Hunan Science and Technology Office
文摘ABE-KONDOH-NAGANO,ABID,YANG-SHIH and LAUNDER-SHARMA low-Reynolds number turbulence models were applied to simulating unsteady turbulence flow around a square cylinder in different phases flow field and time-averaged unsteady flow field.Meanwhile,drag and lift coefficients of the four different low-Reynolds number turbulence models were analyzed.The simulated results of YANG-SHIH model are close to the large eddy simulation results and experimental results,and they are significantly better than those of ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMR models.The modification of the generation of turbulence kinetic energy is the key factor to a successful simulation for YANG-SHIH model,while the correction of the turbulence near the wall has minor influence on the simulation results.For ABE-KONDOH-NAGANO,ABID and LAUNDER-SHARMA models satisfactory simulation results cannot be obtained due to lack of the modification of the generation of turbulence kinetic energy.With the joint force of wall function and the turbulence models with the adoption of corrected swirl stream,flow around a square cylinder can be fully simulated with less grids by the near-wall.
基金Project supported by the National High Technology Research and Development Program of China(Grant No.511-0910-1031)
文摘Based on Xue＇s lattice model, an extended lattice model is proposed by considering the relative current information about next-nearest-neighbour sites ahead. The linear stability condition of the presented model is obtained by employing the linear stability theory. The density wave is investigated analytically with the perturbation method. The results show that the occurrence of traffic jamming transitions can be described by the kink-antikink solution of the modified Korteweg-de Vries （mKdV） equation. The simulation results are in good agreement with the analytical results, showing that the stability of traffic flow can be enhanced when the relative current of next-nearest-neighbour sites ahead is considered.
基金Special Funds for Major State Basic Research Projects of China(G1999022200)
文摘A mathematical modei of two-dimensional turbulent gas-particle twophase flow based on the modified diffusion flux modei (DFM) and a numerical simulation method to analyze the gas-particle flow structures are developed. The modified diffusion flux modei, in which the acceleration due to various forces is taken into account for the calculation of the diffusion velocity of particles, is applicable to the analysis of multi-dimensional gas-particle two-phase turbulent flow. In order to verify its accuracy and efficiency, the numerical simulation by DFM is compared with experimental studies and the prediction by k-ε-kp two-fluid modei, which shows a reasonable agreement. It is confirmed that the modified diffusion flux modei is suitable for simulating the multi-dimensional gas-particle two-phase flow.
基金fund received from the National Key R&D Program of China(2017YFC1501302)the National Natural Science Foundation of China(No.41877242)the China Scholarship Council(CSC)–University of Bologna Joint Scholarship(File No.201806560011)which supported this study。
文摘Some flow-like landslides tend to lose materials while moving on a relatively dry sliding surface.This phenomenon is called bed deposition.In contrast to the bed entrainment phenomenon,bed deposition is relatively poorly understood.Therefore,an improved depth-averaged model is proposed to quantify this phenomenon.The deposition depth is calculated according to the momentum conservation of the deposited mass,and the rheological property of the sliding mass on the bottom is modified considering an abrupt increase in the depth of the sliding surface after deposition.Utilizing the proposed model,the Gaolou landslide,a typical flowlike landside occurring on October 6 th,2006 with an obvious bed deposition phenomenon in Shaanxi Province of China,is simulated to investigate the influence of bed deposition on its propagation process.The results indicate that the proposed model can effectively depict the bed deposition phenomenon in the Gaolou landslide.Bed deposition dissipates part of the kinetic energy of this landslide;thus,the simulated debris inundation area would likely be overestimated when this effect is neglected.On the other hand,the thin liquefied layer formed by the bed deposition process reduces the friction energy dissipation,contributing to the high mobility of this landslide.
