A debris flow descending through an erodible convex colluvial bed,originating from a landslide dam and its upstream deposits,can entrain massive amounts of sediment,dramatically increasing the debris flow volume.Most ...A debris flow descending through an erodible convex colluvial bed,originating from a landslide dam and its upstream deposits,can entrain massive amounts of sediment,dramatically increasing the debris flow volume.Most existing erosion models assume that bed sediments are fully saturated,although this condition is rarely observed in nature.Therefore,a thorough understanding of debris flow overtopping erosion on a convex unsaturated bed is crucial for quantifying disaster risk.In this study,we experimentally investigated the effects of sediment composition,specifically coarse-grain size distribution and fine particle content,on the pore pressure evolution and entrainment of debris flows overriding a convex unsaturated colluvial bed.The average entrainment rate at convex sites for continuously graded bed sediment was higher than its discontinuous counterpart.The measured pore pressures within the unsaturated bed sediments were primarily generated by the passing debris flows.Furthermore,it was found that these pressures decreased as the fine particle content increased and the coarse-grain size of the erodible substrates decreased.When the coarse-grain size of the debris flow was smaller than that of the bed sediment,only a portion of the eroded material was entrained by the moving debris flow.In contrast,when the coarse-grain size of the debris flow was equal to or greater than that of the bed sediment,nearly all of the eroded material was entrained.The findings of this study could contribute to the assessment of hazard amplification and inform the design of mitigation and prevention strategies.展开更多
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.展开更多
During the re-entry of a hypersonic aircraft into the earth’s atmosphere,the surrounding air experiences dissociation,ionization,and other complex chemical phenomena due to extreme temperature by shock wave.To ensure...During the re-entry of a hypersonic aircraft into the earth’s atmosphere,the surrounding air experiences dissociation,ionization,and other complex chemical phenomena due to extreme temperature by shock wave.To ensure thermal safety,the thermochemical non-equilibrium effects resulting from real-gas behavior should be taken into account.In this paper,the characteristics of a double-cone hypersonic laminar flow,including distributions of wall pressure,heat flux,and species dissociation are numerically analyzed with incoming enthalpy of 9.65-21.77 MJ/kg.The thermochemical non-equilibrium flow at different enthalpy and wall temperatures is performed with two-temperature model and Park’s seven chemical reaction model.It is found that the doublecone flow features complex shock-shock interactions to form triple points.The flow topology is further brought out from the analysis of streamlines.At the lowest incoming enthalpy with isothermal wall conditions,two foci points appear.While others highlight only one focal point.As the increment of incoming enthalpy,the heat flux and dissociation of nitrogen and oxygen also increase.An increasing wall temperature leads to a larger separation bubble and a lower value of heat flux and pressure peak,while massive dissociation occurs without obvious ionization under considered cases.展开更多
Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the...Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.展开更多
By use of the two-layer adiabatic globe spectral model and the zonally averaged climatic data of winter season as initial values, 10-day integrations are carried out based on three kinds of model topography (i.e., (1)...By use of the two-layer adiabatic globe spectral model and the zonally averaged climatic data of winter season as initial values, 10-day integrations are carried out based on three kinds of model topography (i.e., (1) the averaged topography; (2) the envelope topography; (3) the modified envelope topography). The results show that the orography of the Northern Hemisphere plays an important role in the simulation of large-scale weather patterns in winter season. The simulation based on the envelope topography developed by Wallace et al. has some improvements in the Rocky Mountains area. But this scheme causes very serious horizontal expansion around the Tibetan Plateau (hereafter referred to as the TV). A modified envelope topography scheme has been worked out that increases the slope of the TP by decreasing the horizontal expansion while keeping the maximum altitude. The results show some improvements of the scheme around the TP. By analysis of the mechanical effects of the large-scale orography on the currents, the different forcings of the air flow over and around the TP and the Rocky Mountain (the RM) are investigated.展开更多
The water flow over a semicircular weir is investigated numerically and experimentally in this paper. The numerical model solves the Reynolds equation for a mean flow field with thek-ε-turbulent model. To trace the m...The water flow over a semicircular weir is investigated numerically and experimentally in this paper. The numerical model solves the Reynolds equation for a mean flow field with thek-ε-turbulent model. To trace the motion of the free surface, the COF method with geometric reconstruction is employed. The velocity of the flow is measured by means of LDV technique. Four types of flow patterns, the position of the separation and reattachment point, the distribution of shear stress on the bed at downstream of the weir are presented and discussed. The numerical results agree well with the experiment data.展开更多
The numerical simulation using the multiple relaxation time lattice Boltzmann method （MRT-LBM） is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these th...The numerical simulation using the multiple relaxation time lattice Boltzmann method （MRT-LBM） is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the moving cylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/ D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D = 0.5 and fv= 0.4 （where A is the forced amplitude, D is the cylinder diameter, and fv stands for the vibrating frequency, respectively）. The calculated results not only indicate that the spacing ratios T/D （T is the center-to-center spacing between the two upstream cylinders） have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D （where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder）.展开更多
The flow over mountain is quite complicated. There are a lot of papers on this problem and a lot of progresses have been made. However, in the most of these papers, just the dynamics contributions of mountain have bee...The flow over mountain is quite complicated. There are a lot of papers on this problem and a lot of progresses have been made. However, in the most of these papers, just the dynamics contributions of mountain have been analysed; the effect of the friction is often neglected. Since the frictional effect is always associated with flow, especially when it flows over the mountain. The study shows that the friction is small in the magnitude but it is not a negligible effect because it changes the features of the flow. In the case of super-or sub-critical flow, there are two extremes: one maximum, one minimum of the fluid surface on the lee-side of the mountain, while in the inviscid fluid, there is just one extreme. The frictional effect should neither be too strong nor too weak to make the situation happened according to the investigation of this paper.展开更多
On the basis of a self-similar solution as well as of the assumption of the 'Transserse Motion',a general linear theory on hypersonic flow over a general slender body is set up in this paper By means of this t...On the basis of a self-similar solution as well as of the assumption of the 'Transserse Motion',a general linear theory on hypersonic flow over a general slender body is set up in this paper By means of this theory, the problem concerned can he put into a universal system of O.D Eqs .which can be integrated manerically in adyance展开更多
Based on the observation data and the reanalysis datasets, the variability and the circulation features influencing precipitation in the Tibetan Plateau (TP) are investigated. Taking into account the effects of topogr...Based on the observation data and the reanalysis datasets, the variability and the circulation features influencing precipitation in the Tibetan Plateau (TP) are investigated. Taking into account the effects of topography, surface winds are deconstructed into flow-around and flow-over components relative to the TP. Climatologically, the flow-around component mainly represents cyclonic circulation in the TP during the summer. The transition zone of total precipitation in the summer parallels the convergence belt between the southerlies and the northerlies of the flow-over component. The leading mode of rainfall anomalies in the TP has a meridional dipole structure, and the first principal component (PC1) mainly depicts the variation of rainfall in the southern TP. The wet southern TP experiences strengthened flow-over, which in turn mechanistically favors intensified ascent forced by the flow-over component. In addition, variations in the Indian summer monsoon (ISM) have an important role in influencing the flow over the southern TP, and the ISM ultimately impacts the precipitation over southern TP.展开更多
The magnetohydrodynamic（MHD） flow induced by a stretching or shrinking sheet under slip conditions is studied.Analytical solutions based on the boundary layer assumption are obtained in a closed form and can be appl...The magnetohydrodynamic（MHD） flow induced by a stretching or shrinking sheet under slip conditions is studied.Analytical solutions based on the boundary layer assumption are obtained in a closed form and can be applied to a flow configuration with any arbitrary velocity distributions. Seven typical sheet velocity profiles are employed as illustrating examples. The solutions to the slip MHD flow are derived from the general solution and discussed in detail. Different from self-similar boundary layer flows, the flows studied in this work have solutions in explicit analytical forms. However, the current flows require special mass transfer at the wall, which is determined by the moving velocity of the sheet. The effects of the slip parameter, the mass transfer at the wall, and the magnetic field on the flow are also demonstrated.展开更多
The studies of asymmetric vortices flow over slender body and its active control at high angles of attack have significant importance for both academic field and engineering area.This paper attempts to provide an upda...The studies of asymmetric vortices flow over slender body and its active control at high angles of attack have significant importance for both academic field and engineering area.This paper attempts to provide an update state of art to the investigations on the fields of forebody asymmetric vortices.This review emphasizes the correlation between micro-perturbation on the model nose and its response and evolution behaviors of the asymmetric vortices.The critical issues are discussed, which include the formation and evolution mechanism of asymmetric multi-vortices;main behaviors of asymmetric vortices flow including its deterministic feature and vortices flow structure;the evolution and development of asymmetric vortices under the perturbation on the model nose;forebody vortex active control especially discussed micro-perturbation active control concept and technique in more detail.However present understanding in this area is still very limited and this paper tries to identify the key unknown problems in the concluding remarks.展开更多
Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical d...Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors.WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction,and its advantage is to reduce dissipation in smooth regions of flow while preserving high-resolution around discontinuities.Here WENO-K scheme is adopted to reconstruct left and right state values within the Roe Riemann solver updating the inviscid fluxes on a structured dynamic overset grid.To minimize the accuracy loss for high-order reconstruction on artificial boundaries of overset grid,a method of multilayer fringes is proposed to carry out interpolation between background grid and blade grid.Massively parallel computing considering automatic load balance on averagely partitioned overset grid is developed to reduce the wall-clock time of an unsteady simulation.Numerical results for Caradonna-Tung(C-T)rotor in hover at the conditions of subsonic and transonic tip Mach numbers show that the thrust coefficient error for the result of WENO-K scheme is no more than 3%.Compared with WENO-JS scheme,WENO-K scheme achieves about 40%improvement on accuracy of predicting rotor thrust with only 4.1%extra computational cost.More importantly,WENO-K scheme can capture more sophisticated unsteady flow structures and resolve tip vortices to a larger wake age with an increment of about 270°compared to WENO-JS scheme.展开更多
The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Prov...The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.展开更多
The summer mean water vapor transport （WVT） and cross-equatorial flow （CEF） over the Asian- Australian monsoon region simulated by 22 coupled atmospheric-oceanic general circulation models （AOGCMs） from the Worl...The summer mean water vapor transport （WVT） and cross-equatorial flow （CEF） over the Asian- Australian monsoon region simulated by 22 coupled atmospheric-oceanic general circulation models （AOGCMs） from the World Climate Research Programme＇s Coupled Model Intercomparison Project Phase 5 （CMIP5） were evaluated. Based on climatology of the twentieth-century simulations, most of models have a reason- ably realistic representation of summer monsoon WVT characterized by southeast water vapor conveyor belt over the South Indian Ocean and southwest belt from the Arabian Sea to the East Asian. The correlation coefficients between NCEP reanalysis and simulations of BCC-CSMI-1, BNU-ESM, CanESM2, FGOALS-s2, MIROC4h and MPI-ESM-LR are up to 0.8. The simulated CEF depicted by the meridional wind along the equator includes the Somali jet and eastern CEF in low atmosphere and the reverse circulation in upper atmosphere, which were generally consistent with NCEP reanalysis. Multi-model ensemble means （MME） can reproduce more reasonable climatological features in spatial distribution both of WVT and CEF. Ten models with more reasonable WVT simulations were selected for future projection studies, including BCC- CSMI-1, BNU-ESM, CanESM2, CCSM4, FGOALS-s2, FIO-ESM, GFDL-ESM2G, MRIOCS, MPI-ESM-LR and NorESM-1M. Analysis based on the future projection experiments in RCP （Representative Concentra- tion Pathway） 2.6, RCP4.5, RCP6 and RCP8.5 show that the global warming forced by different RCP scenarios will results in enhanced WVT over the Indian area and the west Pacific and weaken WVT in the low latitudes of tropical Indian Ocean.展开更多
In this study,magneto-hydrodynamics (MHD) mixed convection effects of Al2O3-water nanofluid flow over a backward-facing step were investigated numerically for various electrical conductivity models of nanofluids.A uni...In this study,magneto-hydrodynamics (MHD) mixed convection effects of Al2O3-water nanofluid flow over a backward-facing step were investigated numerically for various electrical conductivity models of nanofluids.A uniform external magnetic field was applied to the flow and strength of magnetic field was varied with different values of dimensionless parameter Hartmann number (Ha=0,10,20,30,40).Three different electrical conductivity models were used to see the effects of MHD nanofluid flow.Besides,five different inclination angles between 0°-90° is used for the external magnetic field.The problem geometry is a backward-facing step which is used in many engineering applications where flow separation and reattachment phenomenon occurs.Mixed type convective heat transfer of backward-facing step was examined with various values of Richardson number (Ri=0.01,0.1,1,10) and four different nanoparticle volume fractions (Ф=0.01,0.015,0.020,0.025) considering different electrical conductivity models.Finite element method via commercial code COMSOL was used for computations.Results indicate that the addition of nanoparticles enhanced heat transfer significantly.Also increasing magnetic field strength and inclination angle increased heat transfer rate.Effects of different electrical conductivity models were also investigated and it was observed that they have significant effects on the fluid flow and heat transfer characteristics in the presence of magnetic field.展开更多
Experimental and numerical methods were applied to investigating high subsonic and supersonic flows over a 60°swept delta wing in fixed state and pitching oscillation. Static pressure coefficient distributions ov...Experimental and numerical methods were applied to investigating high subsonic and supersonic flows over a 60°swept delta wing in fixed state and pitching oscillation. Static pressure coefficient distributions over the wing leeward surface and the hysteresis loops of pressure coefficient versus angle of attack at the sensor locations were obtained by wind tunnel tests. Similar results were obtained by numerical simulations which agreed well with the experiments. Flow structure around the wing was also demonstrated by the numerical simulation. Effects of Mach number and angle of attack on pressure distribution curves in static tests were investigated. Effects of various oscillation parameters including Mach number, mean angle of attack, pitching amplitude and frequency on hysteresis loops were investigated in dynamic tests and the associated physical mechanisms were discussed. Vortex breakdown phenomenon over the wing was identified at high angles of attack using the pressure coefficient curves and hysteresis loops, and its effects on the flow features were discussed.展开更多
Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panM...Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.展开更多
The phenomenon of debris flow is intermediate between mass movement and solid transport. Flows can be sudden, severe and destructive. Understanding debris flow erosion processes is the key to providing geomorphic expl...The phenomenon of debris flow is intermediate between mass movement and solid transport. Flows can be sudden, severe and destructive. Understanding debris flow erosion processes is the key to providing geomorphic explanations, but progress has been limited because the physical-mechanical properties, movement laws and erosion characteristics are different from those of sediment-laden flow. Using infinite slope theory, this research examines the process and mechanism of downcutting erosion over a moveable bed in a viscous debris flow gully. It focuses specifically on the scour depth and the critical slope for viscous debris flow,and formulas for both calculations are presented.Both scour depth and the critical conditions of downcutting erosion are related to debris flow properties(sand volume concentration and flow depth) and gully properties(longitudinal slope,viscous and internal friction angle of gully materials,and coefficient of kinetic friction). In addition, a series of flume experiments was carried out to characterize the scouring process of debris flows with different properties. The calculated values agreed well with the experimental data. These theoretical formulas are reasonable, and using infinite slope theory to analyze down cutting erosion from viscous debris flow is feasible.展开更多
The present article has been fine-tuned with the investigation of mixed convection Darcy-Forchheimer flow of ZnO-SAE50 oil nanolubricant over an inclined rotating disk under the influence of uniform applied magnetic f...The present article has been fine-tuned with the investigation of mixed convection Darcy-Forchheimer flow of ZnO-SAE50 oil nanolubricant over an inclined rotating disk under the influence of uniform applied magnetic field applied to various industries.The current study has been enriched with additional consideration of slip flow,thermal radiation,viscous dissipation,Joulian dissipation and internal heating.In view of augmentation of thermal conductivity of nanolubricant,a new micro-nano-convection model namely Patel model has been invoked.The specialty of this model involves the effects of specific surface area and nano-convection due to Brownian motion of nanoparticles,kinetic theory based micro-convection,liquid layering and particle concentration.Suitably transformed governing equations have been solved numerically by using Runge-Kutta-Fehlberg scheme.An analysis of the present study has shown that applied magnetic field,porosity of the medium,velocity slip and inertia coefficient account for the slowing down of radial as well as tangential flow of ZnO-SAE50 oil nanolubricant,thereby leading to an improvement in velocity and thermal boundary layers.展开更多
基金supported by the National Key R&D Program of China(Grant No.2018YFC1505205)the Science and Technology Research Program of the Institute of Mountain Hazards and Environment,Chinese Academy of Sciences(Grant No.IMHE-ZDRW-01)Sichuan Science and Technology Program(Grant No.2024NSFSC0781).
文摘A debris flow descending through an erodible convex colluvial bed,originating from a landslide dam and its upstream deposits,can entrain massive amounts of sediment,dramatically increasing the debris flow volume.Most existing erosion models assume that bed sediments are fully saturated,although this condition is rarely observed in nature.Therefore,a thorough understanding of debris flow overtopping erosion on a convex unsaturated bed is crucial for quantifying disaster risk.In this study,we experimentally investigated the effects of sediment composition,specifically coarse-grain size distribution and fine particle content,on the pore pressure evolution and entrainment of debris flows overriding a convex unsaturated colluvial bed.The average entrainment rate at convex sites for continuously graded bed sediment was higher than its discontinuous counterpart.The measured pore pressures within the unsaturated bed sediments were primarily generated by the passing debris flows.Furthermore,it was found that these pressures decreased as the fine particle content increased and the coarse-grain size of the erodible substrates decreased.When the coarse-grain size of the debris flow was smaller than that of the bed sediment,only a portion of the eroded material was entrained by the moving debris flow.In contrast,when the coarse-grain size of the debris flow was equal to or greater than that of the bed sediment,nearly all of the eroded material was entrained.The findings of this study could contribute to the assessment of hazard amplification and inform the design of mitigation and prevention strategies.
基金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.
基金supported by the National Natural Science Foundation of China(Grant Nos.12090030 and 12002261)National Postdoctoral Program for Innovative Talents,China(Grant No.BX20200267)the Fundamental Research Funds for the Central Universities,China(Grant No.xzy012024019)。
文摘During the re-entry of a hypersonic aircraft into the earth’s atmosphere,the surrounding air experiences dissociation,ionization,and other complex chemical phenomena due to extreme temperature by shock wave.To ensure thermal safety,the thermochemical non-equilibrium effects resulting from real-gas behavior should be taken into account.In this paper,the characteristics of a double-cone hypersonic laminar flow,including distributions of wall pressure,heat flux,and species dissociation are numerically analyzed with incoming enthalpy of 9.65-21.77 MJ/kg.The thermochemical non-equilibrium flow at different enthalpy and wall temperatures is performed with two-temperature model and Park’s seven chemical reaction model.It is found that the doublecone flow features complex shock-shock interactions to form triple points.The flow topology is further brought out from the analysis of streamlines.At the lowest incoming enthalpy with isothermal wall conditions,two foci points appear.While others highlight only one focal point.As the increment of incoming enthalpy,the heat flux and dissociation of nitrogen and oxygen also increase.An increasing wall temperature leads to a larger separation bubble and a lower value of heat flux and pressure peak,while massive dissociation occurs without obvious ionization under considered cases.
基金Supported by the National Natural Defense Basic Scientific Research Program of China(A262006-1288)the Key Disciplines Program of Shanghai Municipal Commission of Education(J50501)~~
文摘Supersonic axisymmetric jet flow over a missile afterbody containing exhaust jet is simulated using the second order accurate positive schemes method developed for solving the axisymmetric Euler equations based on the 2-D conservation laws.Comparisons between the numerical results and the experimental measurements show excellent agreements.The computed results are in good agreement with the numerical solutions obtained by using third order accurate RKDG finite element method.The results show larger gradient at discontinuous points compared with those obtained by second order accurate TVD schemes.It indicates that the presented method is efficient and reliable for solving the axisymmetric jet with external freestream flows,and shows that the method captures shocks well without numerical noise.
文摘By use of the two-layer adiabatic globe spectral model and the zonally averaged climatic data of winter season as initial values, 10-day integrations are carried out based on three kinds of model topography (i.e., (1) the averaged topography; (2) the envelope topography; (3) the modified envelope topography). The results show that the orography of the Northern Hemisphere plays an important role in the simulation of large-scale weather patterns in winter season. The simulation based on the envelope topography developed by Wallace et al. has some improvements in the Rocky Mountains area. But this scheme causes very serious horizontal expansion around the Tibetan Plateau (hereafter referred to as the TV). A modified envelope topography scheme has been worked out that increases the slope of the TP by decreasing the horizontal expansion while keeping the maximum altitude. The results show some improvements of the scheme around the TP. By analysis of the mechanical effects of the large-scale orography on the currents, the different forcings of the air flow over and around the TP and the Rocky Mountain (the RM) are investigated.
文摘The water flow over a semicircular weir is investigated numerically and experimentally in this paper. The numerical model solves the Reynolds equation for a mean flow field with thek-ε-turbulent model. To trace the motion of the free surface, the COF method with geometric reconstruction is employed. The velocity of the flow is measured by means of LDV technique. Four types of flow patterns, the position of the separation and reattachment point, the distribution of shear stress on the bed at downstream of the weir are presented and discussed. The numerical results agree well with the experiment data.
基金Support by the National Natural Science Foundation of China under Grant Nos.10932010 and 11072220the Natural Science Foundation of Zhejiang Province under Grant Nos.Y607425,Z6090556the Foundation Project for Youths of Zhijiang Normal University under Grant No.KJ20090102
文摘The numerical simulation using the multiple relaxation time lattice Boltzmann method （MRT-LBM） is carried out for the purpose of investigating the two-dimensional flow around three circular cylinders. Among these three circular cylinders, one of the three cylinders on which a forced in-line vibrating is used to do this research and attempt to find out the effects of the moving cylinder and the other two rigid cylinders on the wake characteristics and vortex formation. As a benchmark problem to discuss the problem of lift coefficient r.m.s for these cylinders with spacing ratios T/ D between other rigid side-by-side cylinders, and the calculation is carried out with two compared cases at Reynolds number of 100, two of the cylinders are rigid and the other one is an in-line vibrated cylinder lying downstream, in addition, forced vibrating amplitude and frequency are A/D = 0.5 and fv= 0.4 （where A is the forced amplitude, D is the cylinder diameter, and fv stands for the vibrating frequency, respectively）. The calculated results not only indicate that the spacing ratios T/D （T is the center-to-center spacing between the two upstream cylinders） have influence on the wake patterns and the formation of vortex shedding, but also analyze the lift coefficient r.m.s for the three cylinders with the spacing ratios S/D （where S is the center-to-center spacing between the center of upstream two side-by-side cylinders and downstream cylinder）.
基金This work was supported by the National Science Foundation of U.S.A.National Natural Science Foundation of China.
文摘The flow over mountain is quite complicated. There are a lot of papers on this problem and a lot of progresses have been made. However, in the most of these papers, just the dynamics contributions of mountain have been analysed; the effect of the friction is often neglected. Since the frictional effect is always associated with flow, especially when it flows over the mountain. The study shows that the friction is small in the magnitude but it is not a negligible effect because it changes the features of the flow. In the case of super-or sub-critical flow, there are two extremes: one maximum, one minimum of the fluid surface on the lee-side of the mountain, while in the inviscid fluid, there is just one extreme. The frictional effect should neither be too strong nor too weak to make the situation happened according to the investigation of this paper.
文摘On the basis of a self-similar solution as well as of the assumption of the 'Transserse Motion',a general linear theory on hypersonic flow over a general slender body is set up in this paper By means of this theory, the problem concerned can he put into a universal system of O.D Eqs .which can be integrated manerically in adyance
基金supported by the National Basic Research Program of China (973 program, Grant No.2010CB950400)the National Natural Science Foundation of China (Grant No. 41030961)
文摘Based on the observation data and the reanalysis datasets, the variability and the circulation features influencing precipitation in the Tibetan Plateau (TP) are investigated. Taking into account the effects of topography, surface winds are deconstructed into flow-around and flow-over components relative to the TP. Climatologically, the flow-around component mainly represents cyclonic circulation in the TP during the summer. The transition zone of total precipitation in the summer parallels the convergence belt between the southerlies and the northerlies of the flow-over component. The leading mode of rainfall anomalies in the TP has a meridional dipole structure, and the first principal component (PC1) mainly depicts the variation of rainfall in the southern TP. The wet southern TP experiences strengthened flow-over, which in turn mechanistically favors intensified ascent forced by the flow-over component. In addition, variations in the Indian summer monsoon (ISM) have an important role in influencing the flow over the southern TP, and the ISM ultimately impacts the precipitation over southern TP.
文摘The magnetohydrodynamic（MHD） flow induced by a stretching or shrinking sheet under slip conditions is studied.Analytical solutions based on the boundary layer assumption are obtained in a closed form and can be applied to a flow configuration with any arbitrary velocity distributions. Seven typical sheet velocity profiles are employed as illustrating examples. The solutions to the slip MHD flow are derived from the general solution and discussed in detail. Different from self-similar boundary layer flows, the flows studied in this work have solutions in explicit analytical forms. However, the current flows require special mass transfer at the wall, which is determined by the moving velocity of the sheet. The effects of the slip parameter, the mass transfer at the wall, and the magnetic field on the flow are also demonstrated.
基金The project supported by the National Natural Science Foundation of China(10172017)Aeronautical Science Foundation of China(02A51048)Foundation of National Key Laboratory of Aerodynamic Design and Research(51462020504HK0101)
文摘The studies of asymmetric vortices flow over slender body and its active control at high angles of attack have significant importance for both academic field and engineering area.This paper attempts to provide an update state of art to the investigations on the fields of forebody asymmetric vortices.This review emphasizes the correlation between micro-perturbation on the model nose and its response and evolution behaviors of the asymmetric vortices.The critical issues are discussed, which include the formation and evolution mechanism of asymmetric multi-vortices;main behaviors of asymmetric vortices flow including its deterministic feature and vortices flow structure;the evolution and development of asymmetric vortices under the perturbation on the model nose;forebody vortex active control especially discussed micro-perturbation active control concept and technique in more detail.However present understanding in this area is still very limited and this paper tries to identify the key unknown problems in the concluding remarks.
基金co-supported by the National Natural Science Foundation of China(No.12072285)Shaanxi Science foundation for Distinguished Young Scholars,China(No.2020JC-13)。
文摘Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors.WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction,and its advantage is to reduce dissipation in smooth regions of flow while preserving high-resolution around discontinuities.Here WENO-K scheme is adopted to reconstruct left and right state values within the Roe Riemann solver updating the inviscid fluxes on a structured dynamic overset grid.To minimize the accuracy loss for high-order reconstruction on artificial boundaries of overset grid,a method of multilayer fringes is proposed to carry out interpolation between background grid and blade grid.Massively parallel computing considering automatic load balance on averagely partitioned overset grid is developed to reduce the wall-clock time of an unsteady simulation.Numerical results for Caradonna-Tung(C-T)rotor in hover at the conditions of subsonic and transonic tip Mach numbers show that the thrust coefficient error for the result of WENO-K scheme is no more than 3%.Compared with WENO-JS scheme,WENO-K scheme achieves about 40%improvement on accuracy of predicting rotor thrust with only 4.1%extra computational cost.More importantly,WENO-K scheme can capture more sophisticated unsteady flow structures and resolve tip vortices to a larger wake age with an increment of about 270°compared to WENO-JS scheme.
基金Under the auspices of National Natural Science Foundation of China (No 40571019)
文摘The structure of wind-sand flow under different total sand transport rates was measured with field vertical anemometer and sand trap on the crest of typical coastal transverse ridge in Changli Gold Coast of Hebei Province, which is one of the most typical coastal aeolian distribution regions in China and famous for the tall and typical coastal transverse ridges. The measurement results show that, on the conditions of approximate wind velocities and same surface materials and environments, some changes happen to the structure of wind-sand flow with the increase of total sand transport rate on the crest of coastal transverse ridge. First, the sand transport rates of layers at different heights in the wind-sand flow increase, with the maximum increase at the height layer of 4-8cm. Second, the ratios of sand trans-port rates of layers at different heights to total sand transport rate decrease at the low height layer (0-4cm), but increase at the high height layer (4-60cm). Third, the distribution of the sand transport rate in the wind-sand flow can be expressed by an exponential function at the height layer of 0-40cm, but it changes from power function model to ex-ponential function model in the whole height layer (0-60cm) and changes into polynomial function model at the height layer of 40-60cm with the increase of total sand transport rate. Those changes have a close relationship with the limit of sand grain size of wind flow transporting and composition of sand grain size in the wind-sand flow.
基金supported by the National Basic Research Program of China(973 Program)under Grant 2010CB950500 and 2010CB950304the Scientific Research Foundation of the First Institute of Oceanography,State Oceanic Administration(Grant No.GY02-2001G26)the National Natural Science Foundation of China under Grant No.41206026
文摘The summer mean water vapor transport （WVT） and cross-equatorial flow （CEF） over the Asian- Australian monsoon region simulated by 22 coupled atmospheric-oceanic general circulation models （AOGCMs） from the World Climate Research Programme＇s Coupled Model Intercomparison Project Phase 5 （CMIP5） were evaluated. Based on climatology of the twentieth-century simulations, most of models have a reason- ably realistic representation of summer monsoon WVT characterized by southeast water vapor conveyor belt over the South Indian Ocean and southwest belt from the Arabian Sea to the East Asian. The correlation coefficients between NCEP reanalysis and simulations of BCC-CSMI-1, BNU-ESM, CanESM2, FGOALS-s2, MIROC4h and MPI-ESM-LR are up to 0.8. The simulated CEF depicted by the meridional wind along the equator includes the Somali jet and eastern CEF in low atmosphere and the reverse circulation in upper atmosphere, which were generally consistent with NCEP reanalysis. Multi-model ensemble means （MME） can reproduce more reasonable climatological features in spatial distribution both of WVT and CEF. Ten models with more reasonable WVT simulations were selected for future projection studies, including BCC- CSMI-1, BNU-ESM, CanESM2, CCSM4, FGOALS-s2, FIO-ESM, GFDL-ESM2G, MRIOCS, MPI-ESM-LR and NorESM-1M. Analysis based on the future projection experiments in RCP （Representative Concentra- tion Pathway） 2.6, RCP4.5, RCP6 and RCP8.5 show that the global warming forced by different RCP scenarios will results in enhanced WVT over the Indian area and the west Pacific and weaken WVT in the low latitudes of tropical Indian Ocean.
文摘In this study,magneto-hydrodynamics (MHD) mixed convection effects of Al2O3-water nanofluid flow over a backward-facing step were investigated numerically for various electrical conductivity models of nanofluids.A uniform external magnetic field was applied to the flow and strength of magnetic field was varied with different values of dimensionless parameter Hartmann number (Ha=0,10,20,30,40).Three different electrical conductivity models were used to see the effects of MHD nanofluid flow.Besides,five different inclination angles between 0°-90° is used for the external magnetic field.The problem geometry is a backward-facing step which is used in many engineering applications where flow separation and reattachment phenomenon occurs.Mixed type convective heat transfer of backward-facing step was examined with various values of Richardson number (Ri=0.01,0.1,1,10) and four different nanoparticle volume fractions (Ф=0.01,0.015,0.020,0.025) considering different electrical conductivity models.Finite element method via commercial code COMSOL was used for computations.Results indicate that the addition of nanoparticles enhanced heat transfer significantly.Also increasing magnetic field strength and inclination angle increased heat transfer rate.Effects of different electrical conductivity models were also investigated and it was observed that they have significant effects on the fluid flow and heat transfer characteristics in the presence of magnetic field.
文摘Experimental and numerical methods were applied to investigating high subsonic and supersonic flows over a 60°swept delta wing in fixed state and pitching oscillation. Static pressure coefficient distributions over the wing leeward surface and the hysteresis loops of pressure coefficient versus angle of attack at the sensor locations were obtained by wind tunnel tests. Similar results were obtained by numerical simulations which agreed well with the experiments. Flow structure around the wing was also demonstrated by the numerical simulation. Effects of Mach number and angle of attack on pressure distribution curves in static tests were investigated. Effects of various oscillation parameters including Mach number, mean angle of attack, pitching amplitude and frequency on hysteresis loops were investigated in dynamic tests and the associated physical mechanisms were discussed. Vortex breakdown phenomenon over the wing was identified at high angles of attack using the pressure coefficient curves and hysteresis loops, and its effects on the flow features were discussed.
文摘Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.
基金supported by the National Natural Science Foundation of China (Grant Nos.50979103 and 40901007)
文摘The phenomenon of debris flow is intermediate between mass movement and solid transport. Flows can be sudden, severe and destructive. Understanding debris flow erosion processes is the key to providing geomorphic explanations, but progress has been limited because the physical-mechanical properties, movement laws and erosion characteristics are different from those of sediment-laden flow. Using infinite slope theory, this research examines the process and mechanism of downcutting erosion over a moveable bed in a viscous debris flow gully. It focuses specifically on the scour depth and the critical slope for viscous debris flow,and formulas for both calculations are presented.Both scour depth and the critical conditions of downcutting erosion are related to debris flow properties(sand volume concentration and flow depth) and gully properties(longitudinal slope,viscous and internal friction angle of gully materials,and coefficient of kinetic friction). In addition, a series of flume experiments was carried out to characterize the scouring process of debris flows with different properties. The calculated values agreed well with the experimental data. These theoretical formulas are reasonable, and using infinite slope theory to analyze down cutting erosion from viscous debris flow is feasible.
文摘The present article has been fine-tuned with the investigation of mixed convection Darcy-Forchheimer flow of ZnO-SAE50 oil nanolubricant over an inclined rotating disk under the influence of uniform applied magnetic field applied to various industries.The current study has been enriched with additional consideration of slip flow,thermal radiation,viscous dissipation,Joulian dissipation and internal heating.In view of augmentation of thermal conductivity of nanolubricant,a new micro-nano-convection model namely Patel model has been invoked.The specialty of this model involves the effects of specific surface area and nano-convection due to Brownian motion of nanoparticles,kinetic theory based micro-convection,liquid layering and particle concentration.Suitably transformed governing equations have been solved numerically by using Runge-Kutta-Fehlberg scheme.An analysis of the present study has shown that applied magnetic field,porosity of the medium,velocity slip and inertia coefficient account for the slowing down of radial as well as tangential flow of ZnO-SAE50 oil nanolubricant,thereby leading to an improvement in velocity and thermal boundary layers.
