It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development ...It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent bound- ary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation （DNS）. It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart Allmaras （S A） model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S-A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S-A model are obviously improved. Thus it can be concluded that the mod- ification of S-A model with the pressure gradient can improve the predictive accuracy for simulat- ing the shock-wave/turbulent boundary laver interaction.展开更多
The sand dunes are typical bed forms of natural alluvial rivers. In this article, a vertical 2-D Reynolds stress model is established for the simulation of turbulent flows around sand dunes, and water-sand boundary co...The sand dunes are typical bed forms of natural alluvial rivers. In this article, a vertical 2-D Reynolds stress model is established for the simulation of turbulent flows around sand dunes, and water-sand boundary conditions are set with particular attention. By numerical simulations, the following conclusions can be drawn. （1） The flow resistance in rivers with sand dunes could be divided into the sand-grain resistance and the sand dune resistance, and the sand-grain resistance coefficient mainly depends on Reynolds number, relative sand grain roughness and sand dune steepness. This coefficient in rivers with sand dunes would be larger than that calculated in a flat riverbed, and the steeper the sand dunes, the larger the sand-grain resistance coefficient. （2） The sand dune resistance coefficient mainly depends on the relative sand dune height and sand dune steepness, the steeper the sand dunes, the larger the sand dune resistance coefficient. （3） For the flat riverbed, the turbulent eddy viscosity coefficient and the sediment diffusion coefficient are approximately identical, but for the sand dune riverbed, in the vertical position, where the sediment diffusion coefficient reaches its maximum, it would be higher than the turbulent eddy viscosity coefficient.展开更多
Time sequence signals of streamwise and normal velocity components,as well as velocity strain rate,at different vertical locations in the turbulent boundary layer over a smooth flat plate in a wind tunnel have been fi...Time sequence signals of streamwise and normal velocity components,as well as velocity strain rate,at different vertical locations in the turbulent boundary layer over a smooth flat plate in a wind tunnel have been finely examined by the use of double-sensor hot-wire anemometry.The local module maximum for wavelet coefficient of longitudinal velocity component,as a detecting index,is employed to educe the ejection and sweep process of the coherent structure burst in the turbulent boundary layer from the random fluctuating background.The coherent waveforms of Reynolds stress residual contribution term for random fluctuations to coherent structure,as well as the velocity strain rate of coherent structure,are extracted by the conditional phase average technique.Based on the theoretical analysis of eddy viscosity coefficient in complex eddy viscosity model for coherent structure,the macro-relaxation effect between Reynolds stress residual contribution term of random fluctuations to coherent structure and the velocity strain rate of coherent structure is studied and the variations of the phase difference between them across the turbulent boundary layer are investigated experimentally.The rationality of complex eddy viscosity model for coherent structure is confirmed through the investigation.展开更多
Breaking wave induced nearsurface turbulence has important consequences for many physical and biochemical processes including water column and nutrients mixing, heat and gases exchange across air-sea interface. The en...Breaking wave induced nearsurface turbulence has important consequences for many physical and biochemical processes including water column and nutrients mixing, heat and gases exchange across air-sea interface. The energy loss from wave breaking and the bubble plume penetration depth are estimated. As a consequence, the vertical distribution of the turbulent kinetic energy （TKE）, the TKE dissipation rate and the eddy viscosity induced by wave breaking are also provided. It is indicated that model results are found to be consistent with the observational evidence that most TKE generated by wave breaking is lost within a depth of a few meters near the sea surface. High turbulence level with intensities of eddy viscosity induced by breaking is nearly four orders larger than vw1（ = κu ＊wz）, the value predicted for the wall layer scaling close to the surface, where u ＊w is the friction velocity in water, κ with 0. 4 is the yon Kármán constant, and z is the water depth, and the strength of the eddy viscosity depends both on wind speed and sea state, and decays rapidly through the depth. This leads to the conclusion that the breaking wave induced vertical mixing is mainly limited to the near surface layer, well above the classical values expected from the similarity theory. Deeper down, however, the effects of wave breaking on the vertical mixing become less important.展开更多
This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transp...This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transport equations, separately, to compute laminar kinetic energy, turbulent kinetic energy, and dissipation rate in a flow field. It needs neither correlations of intermittency factors nor knowledge of the transition onset. Two transition tests are carried out: flat plate boundary layer under zero ...展开更多
The relationship of lateral eddy viscosity depending on length scale is estimated with the decay rate of mesoscale eddies identified from sea level anomaly of satellite observations. The eddy viscosity is expressed in...The relationship of lateral eddy viscosity depending on length scale is estimated with the decay rate of mesoscale eddies identified from sea level anomaly of satellite observations. The eddy viscosity is expressed in terms of the mesoscale eddy parameters according to vortex dynamics. The census of mesoscale eddies shows, in general, that the eddy numbers obey the e-folding decay laws in terms of their amplitude, area and lifetime. The intrinsic values in the e-folding laws are used to estimate the lateral eddy viscosity. Dislike the previous theory that diffusivities are proportional to the length square, the eddy mixing rates (diffusivity and viscosity) from satellite mesoscale eddy datasets are proportional to rs to power of 1.8 (slightly less than 2), where rs is the radius of eddy with radius larger than the Batchelor scale. Additionally, the extrapolation of the eddy mixing to the molecule scale implies that the above power laws may hold until the value of rs is less than O (1 m). These mixing rates with the new parameterizations are suggested to use in numerical schemes. Finally, the climatological distributions of eddy viscosity are calculated.展开更多
基金supported by the National Natural Science Foundation of China (No.11302012,51376001,51136003)the National Basic Research Program of China (No.2012CB720205)+3 种基金the National Magnetic Confinement Fusion Research Program of China (No.2012GB102006)the Aeronautical Science Foundation of China (No.2012ZB51014)the ‘‘111’’ Project(No.B08009)the Astronautical Technology Innovation Foundation of China
文摘It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent bound- ary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation （DNS）. It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart Allmaras （S A） model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S-A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S-A model are obviously improved. Thus it can be concluded that the mod- ification of S-A model with the pressure gradient can improve the predictive accuracy for simulat- ing the shock-wave/turbulent boundary laver interaction.
基金support by the National Natural Science Foundation of China (Grant No.50539060)
文摘The sand dunes are typical bed forms of natural alluvial rivers. In this article, a vertical 2-D Reynolds stress model is established for the simulation of turbulent flows around sand dunes, and water-sand boundary conditions are set with particular attention. By numerical simulations, the following conclusions can be drawn. （1） The flow resistance in rivers with sand dunes could be divided into the sand-grain resistance and the sand dune resistance, and the sand-grain resistance coefficient mainly depends on Reynolds number, relative sand grain roughness and sand dune steepness. This coefficient in rivers with sand dunes would be larger than that calculated in a flat riverbed, and the steeper the sand dunes, the larger the sand-grain resistance coefficient. （2） The sand dune resistance coefficient mainly depends on the relative sand dune height and sand dune steepness, the steeper the sand dunes, the larger the sand dune resistance coefficient. （3） For the flat riverbed, the turbulent eddy viscosity coefficient and the sediment diffusion coefficient are approximately identical, but for the sand dune riverbed, in the vertical position, where the sediment diffusion coefficient reaches its maximum, it would be higher than the turbulent eddy viscosity coefficient.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10832001 and 10872145)Opening Subject of State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences
文摘Time sequence signals of streamwise and normal velocity components,as well as velocity strain rate,at different vertical locations in the turbulent boundary layer over a smooth flat plate in a wind tunnel have been finely examined by the use of double-sensor hot-wire anemometry.The local module maximum for wavelet coefficient of longitudinal velocity component,as a detecting index,is employed to educe the ejection and sweep process of the coherent structure burst in the turbulent boundary layer from the random fluctuating background.The coherent waveforms of Reynolds stress residual contribution term for random fluctuations to coherent structure,as well as the velocity strain rate of coherent structure,are extracted by the conditional phase average technique.Based on the theoretical analysis of eddy viscosity coefficient in complex eddy viscosity model for coherent structure,the macro-relaxation effect between Reynolds stress residual contribution term of random fluctuations to coherent structure and the velocity strain rate of coherent structure is studied and the variations of the phase difference between them across the turbulent boundary layer are investigated experimentally.The rationality of complex eddy viscosity model for coherent structure is confirmed through the investigation.
基金This research was supported by the National Natural Science Foundation of China under contract Nos 40576021 and 40531005.
文摘Breaking wave induced nearsurface turbulence has important consequences for many physical and biochemical processes including water column and nutrients mixing, heat and gases exchange across air-sea interface. The energy loss from wave breaking and the bubble plume penetration depth are estimated. As a consequence, the vertical distribution of the turbulent kinetic energy （TKE）, the TKE dissipation rate and the eddy viscosity induced by wave breaking are also provided. It is indicated that model results are found to be consistent with the observational evidence that most TKE generated by wave breaking is lost within a depth of a few meters near the sea surface. High turbulence level with intensities of eddy viscosity induced by breaking is nearly four orders larger than vw1（ = κu ＊wz）, the value predicted for the wall layer scaling close to the surface, where u ＊w is the friction velocity in water, κ with 0. 4 is the yon Kármán constant, and z is the water depth, and the strength of the eddy viscosity depends both on wind speed and sea state, and decays rapidly through the depth. This leads to the conclusion that the breaking wave induced vertical mixing is mainly limited to the near surface layer, well above the classical values expected from the similarity theory. Deeper down, however, the effects of wave breaking on the vertical mixing become less important.
文摘This article presents a linear eddy-viscosity turbulence model for predicting bypass and natural transition in boundary layers by using Reynolds-averaged Navier-Stokes (RANS) equations. The model includes three transport equations, separately, to compute laminar kinetic energy, turbulent kinetic energy, and dissipation rate in a flow field. It needs neither correlations of intermittency factors nor knowledge of the transition onset. Two transition tests are carried out: flat plate boundary layer under zero ...
文摘The relationship of lateral eddy viscosity depending on length scale is estimated with the decay rate of mesoscale eddies identified from sea level anomaly of satellite observations. The eddy viscosity is expressed in terms of the mesoscale eddy parameters according to vortex dynamics. The census of mesoscale eddies shows, in general, that the eddy numbers obey the e-folding decay laws in terms of their amplitude, area and lifetime. The intrinsic values in the e-folding laws are used to estimate the lateral eddy viscosity. Dislike the previous theory that diffusivities are proportional to the length square, the eddy mixing rates (diffusivity and viscosity) from satellite mesoscale eddy datasets are proportional to rs to power of 1.8 (slightly less than 2), where rs is the radius of eddy with radius larger than the Batchelor scale. Additionally, the extrapolation of the eddy mixing to the molecule scale implies that the above power laws may hold until the value of rs is less than O (1 m). These mixing rates with the new parameterizations are suggested to use in numerical schemes. Finally, the climatological distributions of eddy viscosity are calculated.
