The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is use...The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.展开更多
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展开更多
A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently so...A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently solved by the modified incomplete LU decomposition conjugate gradient scheme ( MILU-CG). The effects of surface suction or blowing' s position and strength on the vortex structures in the cylinder wake, as well as on the drag and lift forces at Reynoldes number Re = 100 were investigated numerically. The results show that the suction on the shoulder of the cylinder or the blowing on the rear of the cylinder can effeciently suppress the asymmetry of the vortex wake in the transverse direction and greatly reduce the lift force; the suction on the shoulder of the cylinder, when its strength is properly chosen, can reduce the drag force significantly, too.展开更多
For deposit body medium, the internal structural properties may be the controlling factors for the strength of the material and the mechanical response. Based on the results of soil-rock meso-statistics using digital ...For deposit body medium, the internal structural properties may be the controlling factors for the strength of the material and the mechanical response. Based on the results of soil-rock meso-statistics using digital imaging, a simulated annealing algorithm is adopted to expand the meso-structural features of deposit bodies in 3D. The construction of the 3D meso-structure of a deposit body is achieved, and then the particle flow analysis program PFC3 D is used to simulate the mechanical properties of the deposit body. It is shown that with a combination of the simulated annealing algorithm and the statistical feature functions, the randomness and heterogeneity of the rock distribution in the 3D inner structure of deposit body medium can be realized, and the reconstructed structural features of the deposit medium can match the features of the digital images well. The spatial utilizations and the compacting effects of the body-centered cubic, hexagonal close and face-centered packing models are high, so these structures can be applied in the simulations of the deposit structures. However, the shear features of the deposit medium vary depending on the different model constructive modes. Rocks, which are the backbone of the deposit, are the factors that determine the shear strength and deformation modulus of the deposit body. The modeling method proposed is useful for the construction of 3D meso-scope models from 2D meso-scope statistics and can be used for studying the mechanical properties of mixed media, such as deposit bodies.展开更多
The coupled motion of two flexible bodies with different lengths immersed in moving fluid is studied numerically. The flapping frequency, flapping amplitude and average drag coefficient of each body are calculated and...The coupled motion of two flexible bodies with different lengths immersed in moving fluid is studied numerically. The flapping frequency, flapping amplitude and average drag coefficient of each body are calculated and the influences of the arranging manner and separation distance are analyzed. In our simulation, when placed in the flow individually, the flexible body with a longer length will flap in period and the shorter one will maintain still straightly in the flow direction. The numerical results show that, two different flexible structures near placed in moving flow would strongly interact. When they are placed side by side, the existence of the stable shorter flexible body will restrain the flapping of the longer one while the existence of the longer flexible body may also induce the shorter one to flap synchronously. When placed in tandem with the shorter flexible body in upstream, the flapping of the longer one in downstream will be obviously enhanced. In the situation for the longer flexible body placed in upstream of the shorter one, the coupled flapping amplitude and average drag coefficients increase and decrease periodically with increasing the arranging space, and peak values appear as a result of the mediate of the tail wakes.展开更多
Large eddy simulation cooperated with the second order full extension ETG(Euler-Taylor-Galerkin) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ra...Large eddy simulation cooperated with the second order full extension ETG(Euler-Taylor-Galerkin) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ratio of (1.5.) The second order full extension ETG finite element method was developed by Wang and He. By means of Taylor expansion of terms containing time derivative, time derivative is replaced by space derivative. The function of it is equal to introducing an artificial viscosity term. The streamlines of the flow at different moments were obtained. The time history of drag coefficient, lift coefficient and the streamwise velocity on the symmetrical points were presented. Furthermore, the symmetrical problem of the frequency spectrum of flow around two square cylinders arranged side by side were studied by using the spectral analysis technology. The data obtained at the initial stage are excluded in order to avoid the influence of initial condition on the results. The power spectrums of drag coefficient, lift coefficient, the streamwise velocity on the symmetrical points were analyzed respectively. The results show that although the time domain process of dynamic parameters is non-symmetrical, the frequency domain process of them is symmetrical under the symmetrical boundary conditions.展开更多
The aerodynamic characteristics are vital for short cylindrical Terminal Sensitive Bullets(TSB)with lowaspect ratio,especially in terminal trajectory.Currently,there is little research in terms of the TSB andshort cyl...The aerodynamic characteristics are vital for short cylindrical Terminal Sensitive Bullets(TSB)with lowaspect ratio,especially in terminal trajectory.Currently,there is little research in terms of the TSB andshort cylinder with two free ends,and particularly in this trajectory,where the scanning angleβand rollangle a vary over a broad range between 0°and 180.In this work wind tunnel experiments are firstconducted to learn the effects of Reynolds number and scanning angle on aerodynamic parameters forshort cylinder with aspect ratio L/D=1.Similar to infinite cylinder,for the short cylinder with two freeends,the drag crisis phenomenon still exists in the critical regime 1.7×10^(5)≤Re≤6.8×10^(5).Then 3Dsimulations are performed to demonstrate the aerodynamic characteristics of short cylinder and TSBover a broad range of Re,L/D,a and 6.The sensitivity analysis of time step and grid are presented as well.whenβ3=0,for short cylinder,the drag crisis phenomenon was also observed in the simulation,but notas obvious as in the wind tunnel test.In some attitudes,there is an obvious Karman vortexin the wake ofshort cylinder and TSB.The correlation between time-averaged aerodynamic coefficients and L/D,Re,a&l is discussed.The vortex shedding frequency and shear layer behavior are obtained for quasi-steadyand unsteady flow.Finally,the effect of end's shape on drag reduction and vortex shedding frequency isanalyzed.展开更多
Based on integral equation method, the study of unsteady supercavitating flow past cones is presented. The shape and length of supercavity are calculated respectively using the finite difference time discretization me...Based on integral equation method, the study of unsteady supercavitating flow past cones is presented. The shape and length of supercavity are calculated respectively using the finite difference time discretization method. The characteristics of the shape and length of supercavities, which vary with the cone＇s angle and cavitation number, are investigated respectively, the varied features of some supercavity scales are analyzed when the flow field is perturbed periodically. The curves relationship between cavity length and cavitation number, which are based respectively on present method and other theories, are discussed and compared. It is obviously shown that the supercavity changes have two characteristics： retardance and wave. These results obtained would be useful in the case of design and analysis of cavitator under water.展开更多
When there are several bodies with relative motion in a flow field,such as the flow in the compression system of modern aero-engine,the flow field will have certain special features,one of which is that the time-space...When there are several bodies with relative motion in a flow field,such as the flow in the compression system of modern aero-engine,the flow field will have certain special features,one of which is that the time-space structure of such multi-bodies unsteady vorticity flow field would be either of unsteady natural flow(UNF)pattern or of unsteady cooperative flow(UCF)pattern.If we further examine the aerodynamic design system of aero-engine,there is no mechanism for the unsteady cooperative flow to occur,in other words the flow field must be of the unsteady natural flow type.If certain technical measures can be adopted to transform UNF into UCF,the aerodynamic performances will surely be improved.This is the main task the author and their colleague have been devoted to and the results are reviewed in the present paper with emphases laid on basic ideas,technical approaches and experimental verifications.展开更多
Flows containing steady or nearly steady strong shocks on parts of the flow field,and unsteady turbulence with shocklets on other parts of the flow field are difficult to capture accurately and efficiently employing t...Flows containing steady or nearly steady strong shocks on parts of the flow field,and unsteady turbulence with shocklets on other parts of the flow field are difficult to capture accurately and efficiently employing the same numerical scheme,even under the multiblock grid or adaptive grid refinement framework.While sixthorder or higher-order shock-capturing methods are appropriate for unsteady turbulence with shocklets,third-order or lower shock-capturing methods are more effective for strong steady or nearly steady shocks in terms of convergence.In order to minimize the short comings of low order and high order shock-capturing schemes for the subject flows,a multiblock overlapping grid with different types of spatial schemes and orders of accuracy on different blocks is proposed.The recently developed single block high order filter scheme in generalized geometries for Navier Stokes and magnetohydrodynamics systems is extended to multiblock overlapping grid geometries.The first stage in validating the high order overlapping approach with several test cases is included.展开更多
基金The project supported by the National Natural Science Foundation of China (50506028)Action Scheme for Invigorating Education Towards the twenty-first century
文摘The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.
文摘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
基金Foundation item:the Natural Science Foundation of Jiangsu Province(BK97056109)
文摘A vorticity-velocity method was used to study the incompressible viscous fluid flow around a circular cylinder with surface suction or blowing. The resulted high order implicit difference equations were effeciently solved by the modified incomplete LU decomposition conjugate gradient scheme ( MILU-CG). The effects of surface suction or blowing' s position and strength on the vortex structures in the cylinder wake, as well as on the drag and lift forces at Reynoldes number Re = 100 were investigated numerically. The results show that the suction on the shoulder of the cylinder or the blowing on the rear of the cylinder can effeciently suppress the asymmetry of the vortex wake in the transverse direction and greatly reduce the lift force; the suction on the shoulder of the cylinder, when its strength is properly chosen, can reduce the drag force significantly, too.
基金Projects(51309089,11202063)supported by the National Natural Science Foundation of ChinaProject(2013BAB06B01)supported by the National High Technology Research and Development Program of China+1 种基金Project(2015CB057903)supported by the National Basic Research Program of ChinaProject(BK20130846)supported by Natural Science Foundation of Jiangsu Province,China
文摘For deposit body medium, the internal structural properties may be the controlling factors for the strength of the material and the mechanical response. Based on the results of soil-rock meso-statistics using digital imaging, a simulated annealing algorithm is adopted to expand the meso-structural features of deposit bodies in 3D. The construction of the 3D meso-structure of a deposit body is achieved, and then the particle flow analysis program PFC3 D is used to simulate the mechanical properties of the deposit body. It is shown that with a combination of the simulated annealing algorithm and the statistical feature functions, the randomness and heterogeneity of the rock distribution in the 3D inner structure of deposit body medium can be realized, and the reconstructed structural features of the deposit medium can match the features of the digital images well. The spatial utilizations and the compacting effects of the body-centered cubic, hexagonal close and face-centered packing models are high, so these structures can be applied in the simulations of the deposit structures. However, the shear features of the deposit medium vary depending on the different model constructive modes. Rocks, which are the backbone of the deposit, are the factors that determine the shear strength and deformation modulus of the deposit body. The modeling method proposed is useful for the construction of 3D meso-scope models from 2D meso-scope statistics and can be used for studying the mechanical properties of mixed media, such as deposit bodies.
基金Supported by the National Natural Science Foundation of China under Grant Nos 51479007,51309017,and 11102027the Natural Science Foundation of Hubei Province under Grant No 2015CFA026the Fundamental Research Fund for State Public-Benefic Scientific Institutes of CRSRI under Grant No CKSF2015026/SL
文摘The coupled motion of two flexible bodies with different lengths immersed in moving fluid is studied numerically. The flapping frequency, flapping amplitude and average drag coefficient of each body are calculated and the influences of the arranging manner and separation distance are analyzed. In our simulation, when placed in the flow individually, the flexible body with a longer length will flap in period and the shorter one will maintain still straightly in the flow direction. The numerical results show that, two different flexible structures near placed in moving flow would strongly interact. When they are placed side by side, the existence of the stable shorter flexible body will restrain the flapping of the longer one while the existence of the longer flexible body may also induce the shorter one to flap synchronously. When placed in tandem with the shorter flexible body in upstream, the flapping of the longer one in downstream will be obviously enhanced. In the situation for the longer flexible body placed in upstream of the shorter one, the coupled flapping amplitude and average drag coefficients increase and decrease periodically with increasing the arranging space, and peak values appear as a result of the mediate of the tail wakes.
文摘Large eddy simulation cooperated with the second order full extension ETG(Euler-Taylor-Galerkin) finite element method was applied to simulate the flow around two square cylinders arranged side by side at a spacing ratio of (1.5.) The second order full extension ETG finite element method was developed by Wang and He. By means of Taylor expansion of terms containing time derivative, time derivative is replaced by space derivative. The function of it is equal to introducing an artificial viscosity term. The streamlines of the flow at different moments were obtained. The time history of drag coefficient, lift coefficient and the streamwise velocity on the symmetrical points were presented. Furthermore, the symmetrical problem of the frequency spectrum of flow around two square cylinders arranged side by side were studied by using the spectral analysis technology. The data obtained at the initial stage are excluded in order to avoid the influence of initial condition on the results. The power spectrums of drag coefficient, lift coefficient, the streamwise velocity on the symmetrical points were analyzed respectively. The results show that although the time domain process of dynamic parameters is non-symmetrical, the frequency domain process of them is symmetrical under the symmetrical boundary conditions.
基金The authors would like to acknowledge the financial supports from the National Natural Science Foundation of China(Grant No.11372136)the Special fund for basic scientific research of Central University(Grant No.30916011306)the Postgraduate Research&Practice Innovation Program of Jiangsu Province(Grant No.KYCX17_0386).
文摘The aerodynamic characteristics are vital for short cylindrical Terminal Sensitive Bullets(TSB)with lowaspect ratio,especially in terminal trajectory.Currently,there is little research in terms of the TSB andshort cylinder with two free ends,and particularly in this trajectory,where the scanning angleβand rollangle a vary over a broad range between 0°and 180.In this work wind tunnel experiments are firstconducted to learn the effects of Reynolds number and scanning angle on aerodynamic parameters forshort cylinder with aspect ratio L/D=1.Similar to infinite cylinder,for the short cylinder with two freeends,the drag crisis phenomenon still exists in the critical regime 1.7×10^(5)≤Re≤6.8×10^(5).Then 3Dsimulations are performed to demonstrate the aerodynamic characteristics of short cylinder and TSBover a broad range of Re,L/D,a and 6.The sensitivity analysis of time step and grid are presented as well.whenβ3=0,for short cylinder,the drag crisis phenomenon was also observed in the simulation,but notas obvious as in the wind tunnel test.In some attitudes,there is an obvious Karman vortexin the wake ofshort cylinder and TSB.The correlation between time-averaged aerodynamic coefficients and L/D,Re,a&l is discussed.The vortex shedding frequency and shear layer behavior are obtained for quasi-steadyand unsteady flow.Finally,the effect of end's shape on drag reduction and vortex shedding frequency isanalyzed.
基金Project supported by the National Nature Science Foundation of China (Grant No:10572045) and Distinguished Young Scholar Science Foundation of Heilongjiang Province of China(Grant No: JC-9)
文摘Based on integral equation method, the study of unsteady supercavitating flow past cones is presented. The shape and length of supercavity are calculated respectively using the finite difference time discretization method. The characteristics of the shape and length of supercavities, which vary with the cone＇s angle and cavitation number, are investigated respectively, the varied features of some supercavity scales are analyzed when the flow field is perturbed periodically. The curves relationship between cavity length and cavitation number, which are based respectively on present method and other theories, are discussed and compared. It is obviously shown that the supercavity changes have two characteristics： retardance and wave. These results obtained would be useful in the case of design and analysis of cavitator under water.
基金the financial support received from the National Natural Science Foundation of China(Grant no.10902008)for this project.
文摘When there are several bodies with relative motion in a flow field,such as the flow in the compression system of modern aero-engine,the flow field will have certain special features,one of which is that the time-space structure of such multi-bodies unsteady vorticity flow field would be either of unsteady natural flow(UNF)pattern or of unsteady cooperative flow(UCF)pattern.If we further examine the aerodynamic design system of aero-engine,there is no mechanism for the unsteady cooperative flow to occur,in other words the flow field must be of the unsteady natural flow type.If certain technical measures can be adopted to transform UNF into UCF,the aerodynamic performances will surely be improved.This is the main task the author and their colleague have been devoted to and the results are reviewed in the present paper with emphases laid on basic ideas,technical approaches and experimental verifications.
基金This work performed under the auspices of the U.S.Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344。
文摘Flows containing steady or nearly steady strong shocks on parts of the flow field,and unsteady turbulence with shocklets on other parts of the flow field are difficult to capture accurately and efficiently employing the same numerical scheme,even under the multiblock grid or adaptive grid refinement framework.While sixthorder or higher-order shock-capturing methods are appropriate for unsteady turbulence with shocklets,third-order or lower shock-capturing methods are more effective for strong steady or nearly steady shocks in terms of convergence.In order to minimize the short comings of low order and high order shock-capturing schemes for the subject flows,a multiblock overlapping grid with different types of spatial schemes and orders of accuracy on different blocks is proposed.The recently developed single block high order filter scheme in generalized geometries for Navier Stokes and magnetohydrodynamics systems is extended to multiblock overlapping grid geometries.The first stage in validating the high order overlapping approach with several test cases is included.
