This paper presents a new kind of everage for the loeally-generated large vortexes so that the physieal quantities of the locally-generated large vortexes and the external large vortexes canberigorously separated from...This paper presents a new kind of everage for the loeally-generated large vortexes so that the physieal quantities of the locally-generated large vortexes and the external large vortexes canberigorously separated from the equal ions for the large vortexes proposed in a previous paper[1] To the equations for the two kinds of large vortexes, some auxiliary relations are introduced, and the value, of the length-scale lN of energy dissipation of the external large vortexes may he determined according to the actual circumstances of the disturbance of external sources. Thus the resulting equations of the second moments of turbulent velocity fluctuations for the two kinds of large vortexes can be made closed. Meanwhile, the corresponding coefficients of diffusion in the previous paper[1] are improved,Finally, a closed set of numerically-solvable equations of turbulence model are obtained.展开更多
In view of the fact that large scale vortices play the substantial role of momentum transport in turbulent flows, large eddy simulation （LES） is considered as a better simulation model. However, the sub-grid scale ...In view of the fact that large scale vortices play the substantial role of momentum transport in turbulent flows, large eddy simulation （LES） is considered as a better simulation model. However, the sub-grid scale （SGS） models reported so far have not ascertained under what flow conditions the LES can lapse into the direct nu-merical simulation. To overcome this discrepancy, this paper develops a swirling strength based the SGS model to properly model the turbulence intermittency, with the primary characteristics that when the local swirling strength is zero, the local sub-grid viscosity will be vanished. In this paper, the model is used to investigate the flow characteris-tics of zero-incident incompressible turbulent flows around a single square cylinder （SC） at a low Reynolds number range Re ∈ [103, 104]. The flow characteristics investigated include the Reynolds number dependence of lift and drag coefficients, the distributions of time-spanwise averaged variables such as the sub-grid viscosity and the logarithm of Kolmogorov micro-scale to the base of 10 at Re=2 500 and 104, the contours of spanwise and streamwise vorticity components at t = 170. It is revealed that the peak value of sub-grid viscosity ratio and its root mean square （RMS） values grow with the Reynolds number. The dissipation rate of turbulent kinetic energy is larger near the SC solid walls. The instantaneous factor of swirling strength intermittency （FSI） exhibits some laminated structure involved with vortex shedding.展开更多
A new vortex sheet model was proposed for simulating aircraft wake vortex evolution.Rather than beginning with a pair of counter-rotating cylindrical vortices as in the traditional models, a lift-drag method is used t...A new vortex sheet model was proposed for simulating aircraft wake vortex evolution.Rather than beginning with a pair of counter-rotating cylindrical vortices as in the traditional models, a lift-drag method is used to initialize a vortex sheet so that the roll-up phase is taken into account. The results of this model report a better approximation to a real situation when compared to the measurement data. The roll-up induced structures are proved to influence the far-field decay.On one hand, they lead to an early decay in the diffusion phase. On the other hand, the growth of linear instability such as elliptical instability is suppressed, resulting in a slower decay in the rapid decay phase. This work provides a simple and practicable model for simulating wake vortex evolution, which combines the roll-up process and the far-field phase in simulation. It is also proved that the roll-up phase should not be ignored when simulating the far-field evolution of an aircraft wake vortex pair, which indicates the necessity of this new model.展开更多
In this paper, the three_dimensional vo rtex filament method was used to simulate the evolution of vortex structures in the axisymmetric round jet. The results agree well with the ones given by Chung and Troutt. The...In this paper, the three_dimensional vo rtex filament method was used to simulate the evolution of vortex structures in the axisymmetric round jet. The results agree well with the ones given by Chung and Troutt. Then one_coupling model was employed to calculate the particle motio n based on the computed flows. The results show that the particle motion is affe cted by flows obviously at the case of particle number St1 and negligibly at St1 ,particles distribute around the vortex structures uniformly at St ～1 . When perturbations with wavenumber 5 are introduced to vortex rings, part icles disperse wider along radial direction, which conforms to the experimental results. The degree of particle dispersion is in the direct ratio to the amplitu de of perturbation. The conclusions given in the paper are useful to the practic e usage.展开更多
Large vortices with scales ranging from hundreds meters to tens of kilometers are generally found in the atmospheric convective boundary layer(CBL). These vortices play important roles in the vertical transport of mom...Large vortices with scales ranging from hundreds meters to tens of kilometers are generally found in the atmospheric convective boundary layer(CBL). These vortices play important roles in the vertical transport of momentum,heat,water vapor and other tracers in the boundary layer.On the basis of the view of interaction between the convection in CBL and the gravity waves in the upper stable layer the authors developed a convection-wave theory on the formation of large vortices.According to the theory the wavenumber spectrum of the large vortices mainly depends on the atmospheric conditions in both of the upper and lower layers,such as wind speed,wind direction shear,stratification as well as temperature jump. In the present paper satellite image and weather data in a case of cold air outbreak over warm ocean are analyzed to study every stage of the convective processes,such as cloud street, convective cell as well as their transformation.According to the theory the wavenumber compositions for cloud street and convective cell are calculated,respectively,on the basis of the atmospheric conditions at every stage.The distributions of vertical motions,convergent band and disturbed interface are obtained and compared with the cloud patterns in the convective processes. Thus the study seems to offer a likely explanation for the origin of large vortices in CBL.展开更多
The flow structure of one isothermal swirling case in the Sydney swirl flame database was studied using two numerical methods. Results from the Reynolds-averaged Navier-Stokes (RANS) approach and large eddy simulation...The flow structure of one isothermal swirling case in the Sydney swirl flame database was studied using two numerical methods. Results from the Reynolds-averaged Navier-Stokes (RANS) approach and large eddy simulation (LES) were compared with experimental measurements. The simulations were applied in two different Cartesian grids which were investigated by a grid independence study for RANS and a post-estimator for LES. The RNG k-ε turbulence model was used in RANS and dynamic Smagorinsky-Lilly model was used as the sub-grid scale model in LES. A validation study and cross comparison of ensemble average and root mean square (RMS) results showed LES outperforms RANS statistic results. Flow field results indicated that both approaches could capture dominant flow structures, like vortex breakdown (VB), and precessing vortex core (PVC). Streamlines indicate that the formation mechanisms of VB deducted from the two methods were different. The vorticity field was also studied using a velocity gradient based method. This research gained in-depth understanding of isothermal swirling flow.展开更多
The large eddy simulation（LES） of the flow characteristics in an annular jet pump（AJP） is conducted, and the flow characteristics are systematically analyzed from both time-averaged and instantaneous aspects. The ...The large eddy simulation（LES） of the flow characteristics in an annular jet pump（AJP） is conducted, and the flow characteristics are systematically analyzed from both time-averaged and instantaneous aspects. The jet expansion, the velocity distribution and the energy are considered to analyze the time-averaged evolution of the flow field in the AJP. The transient flow characteristics can also be acquired from the analysis of the turbulence intensity and the Reynolds stress. The simulation demonstrates that in the time-averaged characteristics, the potential cores increase linearly with the increase of the flow ratio. With the flow development, the jet half-width gradually increases and the residual energy coefficient decreases. Compared with the distribution of the time-averaged axial velocity, that of the instantaneous velocity is more complex and disorderly. The high intensity of the axial turbulence mainly occurs in the mixing layer and the near-wall regions of the diffuser. The annular distribution of the Reynolds stress is mainly in the mixing layer and the recirculation region. There is a low-stress zone between the mixing layer and the high-stress region in the wall-boundary layer. The intensity of the spanwise vortexes is larger than that of the streamwise vortexes, and therefore, the former make greater contribution to the total vorticity. This research provides a better understanding of the flow characteristics in the AJP.展开更多
By using the linearized barotropic vorticity equation in polar coordinates the stability of pertur- bations on a large scale circular basic flow is transformed into a generalized eigenvalue problem, yielding the relat...By using the linearized barotropic vorticity equation in polar coordinates the stability of pertur- bations on a large scale circular basic flow is transformed into a generalized eigenvalue problem, yielding the relationship between the growth rate of the amplitude of perturbations and the az- imuthal wave number. Then, numerical experiments whose integration time is 60 model hours are performed in terms of a quasi-geostrophic barotropic model in Cartesian coordinates using the per- turbation stream function field of unstable mode superimposed on a strong and weak circular basic flows as the initial fields. The experimental results reveal that the amplitudes of the initial pertur- bations in the model atmosphere grow with time. The amplitude of the perturbations superimposed on the strong circular basic flow grows quicker and forms a spiral-band-like structure.展开更多
In this paper,the 3-D turbulent flow around two side-by-side circular cylinders of different diameters,at sub-critical Reynolds number（Re=3 900）,is numerically simulated by the large eddy simulation（LES）.The spaci...In this paper,the 3-D turbulent flow around two side-by-side circular cylinders of different diameters,at sub-critical Reynolds number（Re=3 900）,is numerically simulated by the large eddy simulation（LES）.The spacing ratios（T/D）between the two cylinders are considered in four cases（T/D=1.2,1.5,1.8 and 2.7）to study the vortex shedding and turbulent properties in the flow field.The main results are focused on the drag and lift coefficients,the vortex shedding frequency,the coherent structure,and the scale properties.It is shown that when T/D is equal to 1.2,the vortex shedding of the main cylinder is strongly suppressed by the small cylinder,the drag and lift coefficients of the main cylinder are smaller than those in other three cases.While T/D is equal to 1.5,the vortex shedding of the main cylinder can be improved,the drag and lift coefficients of the main cylinder are larger than those in other three cases.The empirical mode decomposition（EMD）method is applied to decompose the velocity signals traced by the LES.It is shown that there is a linear relationship between the mean period and the mode in the semi-log coordinates.The vortex shedding period of the main cylinder is consistent with the period of the restructured coherent structures quantitatively.展开更多
基金The Project Supported by National Natural Science Foundation of China
文摘This paper presents a new kind of everage for the loeally-generated large vortexes so that the physieal quantities of the locally-generated large vortexes and the external large vortexes canberigorously separated from the equal ions for the large vortexes proposed in a previous paper[1] To the equations for the two kinds of large vortexes, some auxiliary relations are introduced, and the value, of the length-scale lN of energy dissipation of the external large vortexes may he determined according to the actual circumstances of the disturbance of external sources. Thus the resulting equations of the second moments of turbulent velocity fluctuations for the two kinds of large vortexes can be made closed. Meanwhile, the corresponding coefficients of diffusion in the previous paper[1] are improved,Finally, a closed set of numerically-solvable equations of turbulence model are obtained.
基金Project supported by the National Natural Science Foundation of China(No.11372303)
文摘In view of the fact that large scale vortices play the substantial role of momentum transport in turbulent flows, large eddy simulation （LES） is considered as a better simulation model. However, the sub-grid scale （SGS） models reported so far have not ascertained under what flow conditions the LES can lapse into the direct nu-merical simulation. To overcome this discrepancy, this paper develops a swirling strength based the SGS model to properly model the turbulence intermittency, with the primary characteristics that when the local swirling strength is zero, the local sub-grid viscosity will be vanished. In this paper, the model is used to investigate the flow characteris-tics of zero-incident incompressible turbulent flows around a single square cylinder （SC） at a low Reynolds number range Re ∈ [103, 104]. The flow characteristics investigated include the Reynolds number dependence of lift and drag coefficients, the distributions of time-spanwise averaged variables such as the sub-grid viscosity and the logarithm of Kolmogorov micro-scale to the base of 10 at Re=2 500 and 104, the contours of spanwise and streamwise vorticity components at t = 170. It is revealed that the peak value of sub-grid viscosity ratio and its root mean square （RMS） values grow with the Reynolds number. The dissipation rate of turbulent kinetic energy is larger near the SC solid walls. The instantaneous factor of swirling strength intermittency （FSI） exhibits some laminated structure involved with vortex shedding.
基金supported by the Boeing-COMAC Aviation Energy Conservation and Emissions Reduction Technology Center (AECER)
文摘A new vortex sheet model was proposed for simulating aircraft wake vortex evolution.Rather than beginning with a pair of counter-rotating cylindrical vortices as in the traditional models, a lift-drag method is used to initialize a vortex sheet so that the roll-up phase is taken into account. The results of this model report a better approximation to a real situation when compared to the measurement data. The roll-up induced structures are proved to influence the far-field decay.On one hand, they lead to an early decay in the diffusion phase. On the other hand, the growth of linear instability such as elliptical instability is suppressed, resulting in a slower decay in the rapid decay phase. This work provides a simple and practicable model for simulating wake vortex evolution, which combines the roll-up process and the far-field phase in simulation. It is also proved that the roll-up phase should not be ignored when simulating the far-field evolution of an aircraft wake vortex pair, which indicates the necessity of this new model.
文摘In this paper, the three_dimensional vo rtex filament method was used to simulate the evolution of vortex structures in the axisymmetric round jet. The results agree well with the ones given by Chung and Troutt. Then one_coupling model was employed to calculate the particle motio n based on the computed flows. The results show that the particle motion is affe cted by flows obviously at the case of particle number St1 and negligibly at St1 ,particles distribute around the vortex structures uniformly at St ～1 . When perturbations with wavenumber 5 are introduced to vortex rings, part icles disperse wider along radial direction, which conforms to the experimental results. The degree of particle dispersion is in the direct ratio to the amplitu de of perturbation. The conclusions given in the paper are useful to the practic e usage.
基金This research was supported by the National Natural Science Foundation of China
文摘Large vortices with scales ranging from hundreds meters to tens of kilometers are generally found in the atmospheric convective boundary layer(CBL). These vortices play important roles in the vertical transport of momentum,heat,water vapor and other tracers in the boundary layer.On the basis of the view of interaction between the convection in CBL and the gravity waves in the upper stable layer the authors developed a convection-wave theory on the formation of large vortices.According to the theory the wavenumber spectrum of the large vortices mainly depends on the atmospheric conditions in both of the upper and lower layers,such as wind speed,wind direction shear,stratification as well as temperature jump. In the present paper satellite image and weather data in a case of cold air outbreak over warm ocean are analyzed to study every stage of the convective processes,such as cloud street, convective cell as well as their transformation.According to the theory the wavenumber compositions for cloud street and convective cell are calculated,respectively,on the basis of the atmospheric conditions at every stage.The distributions of vertical motions,convergent band and disturbed interface are obtained and compared with the cloud patterns in the convective processes. Thus the study seems to offer a likely explanation for the origin of large vortices in CBL.
基金the financial support of DONG Energy [PSO2007-7333]
文摘The flow structure of one isothermal swirling case in the Sydney swirl flame database was studied using two numerical methods. Results from the Reynolds-averaged Navier-Stokes (RANS) approach and large eddy simulation (LES) were compared with experimental measurements. The simulations were applied in two different Cartesian grids which were investigated by a grid independence study for RANS and a post-estimator for LES. The RNG k-ε turbulence model was used in RANS and dynamic Smagorinsky-Lilly model was used as the sub-grid scale model in LES. A validation study and cross comparison of ensemble average and root mean square (RMS) results showed LES outperforms RANS statistic results. Flow field results indicated that both approaches could capture dominant flow structures, like vortex breakdown (VB), and precessing vortex core (PVC). Streamlines indicate that the formation mechanisms of VB deducted from the two methods were different. The vorticity field was also studied using a velocity gradient based method. This research gained in-depth understanding of isothermal swirling flow.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.51179134,11472197)
文摘The large eddy simulation（LES） of the flow characteristics in an annular jet pump（AJP） is conducted, and the flow characteristics are systematically analyzed from both time-averaged and instantaneous aspects. The jet expansion, the velocity distribution and the energy are considered to analyze the time-averaged evolution of the flow field in the AJP. The transient flow characteristics can also be acquired from the analysis of the turbulence intensity and the Reynolds stress. The simulation demonstrates that in the time-averaged characteristics, the potential cores increase linearly with the increase of the flow ratio. With the flow development, the jet half-width gradually increases and the residual energy coefficient decreases. Compared with the distribution of the time-averaged axial velocity, that of the instantaneous velocity is more complex and disorderly. The high intensity of the axial turbulence mainly occurs in the mixing layer and the near-wall regions of the diffuser. The annular distribution of the Reynolds stress is mainly in the mixing layer and the recirculation region. There is a low-stress zone between the mixing layer and the high-stress region in the wall-boundary layer. The intensity of the spanwise vortexes is larger than that of the streamwise vortexes, and therefore, the former make greater contribution to the total vorticity. This research provides a better understanding of the flow characteristics in the AJP.
基金This work was supported by the National Natural Science Foundation of China
文摘By using the linearized barotropic vorticity equation in polar coordinates the stability of pertur- bations on a large scale circular basic flow is transformed into a generalized eigenvalue problem, yielding the relationship between the growth rate of the amplitude of perturbations and the az- imuthal wave number. Then, numerical experiments whose integration time is 60 model hours are performed in terms of a quasi-geostrophic barotropic model in Cartesian coordinates using the per- turbation stream function field of unstable mode superimposed on a strong and weak circular basic flows as the initial fields. The experimental results reveal that the amplitudes of the initial pertur- bations in the model atmosphere grow with time. The amplitude of the perturbations superimposed on the strong circular basic flow grows quicker and forms a spiral-band-like structure.
基金supported by the National Natural Science Foundation of China(Grant Nos.11572203,11332006)the Innovation Program of Shanghai Municipal Education Commission(Grant No.13YZ124)
文摘In this paper,the 3-D turbulent flow around two side-by-side circular cylinders of different diameters,at sub-critical Reynolds number（Re=3 900）,is numerically simulated by the large eddy simulation（LES）.The spacing ratios（T/D）between the two cylinders are considered in four cases（T/D=1.2,1.5,1.8 and 2.7）to study the vortex shedding and turbulent properties in the flow field.The main results are focused on the drag and lift coefficients,the vortex shedding frequency,the coherent structure,and the scale properties.It is shown that when T/D is equal to 1.2,the vortex shedding of the main cylinder is strongly suppressed by the small cylinder,the drag and lift coefficients of the main cylinder are smaller than those in other three cases.While T/D is equal to 1.5,the vortex shedding of the main cylinder can be improved,the drag and lift coefficients of the main cylinder are larger than those in other three cases.The empirical mode decomposition（EMD）method is applied to decompose the velocity signals traced by the LES.It is shown that there is a linear relationship between the mean period and the mode in the semi-log coordinates.The vortex shedding period of the main cylinder is consistent with the period of the restructured coherent structures quantitatively.
