The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation(LES)method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in...The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation(LES)method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation(DNS),which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.展开更多
A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energ...A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.展开更多
Gas solid two-phase turbulent plane jet is applied to many natural s it uations and in engineering systems. To predict the particle dispersion in the ga s jet is of great importance in industrial applications and in ...Gas solid two-phase turbulent plane jet is applied to many natural s it uations and in engineering systems. To predict the particle dispersion in the ga s jet is of great importance in industrial applications and in the designing of engineering systems. A large eddy simulation of the two-phase plane jet was con d ucted to investigate the particle dispersion patterns. The particles with Stokes numbers equal to 0 0028, 0 3, 2 5, 28 (corresponding to particle diameter 1 μm , 10 μm, 30 μm, 100 μm, respectively) in \%Re\%=11 300 gas flow were studied. The simulation results of gas phase motion agreed well with previous experimental re sults. And the simulation results of the solid particles motion showed that part icles with different Stokes number have different spatial dispersion; and that p articles with intermediate Stokes number have the largest dispersion ratio.展开更多
This paper describes the design of a new kind of miniature abrading sphere, which is magnetically mounted inside a spherical gap and set in rotation pneumatically with air. Large eddy simulation is performed in conjun...This paper describes the design of a new kind of miniature abrading sphere, which is magnetically mounted inside a spherical gap and set in rotation pneumatically with air. Large eddy simulation is performed in conjunction with the compressible Smagorinsky model. Minimal temperature variation allows for the assumption of adiabatic walls. Fluid-solid interaction is modeled using the law of the wall for compressible turbulent flow. A parametric study is done to determine optimal geometric layout while taking physical restrictions into account. The resulting optimal configuration is then examined in detail in order to determine demands to be met by the computerized control of the magnetic bearing as well as to quantify the force available to the abrasion process. Finally, a mathematical relation is given that determines available abrasion force depending on standard volumetric flow rate and rotation frequency. The findings presented here provide a basis for further development of smaller versions of the tool.展开更多
Large eddy simulations(LES) were performed to study the non-reacting flow fields of a Cambridge swirl burner. The dynamic Smagorinsky eddy viscosity model is used as the sub-grid scale turbulence model. Comparisons of...Large eddy simulations(LES) were performed to study the non-reacting flow fields of a Cambridge swirl burner. The dynamic Smagorinsky eddy viscosity model is used as the sub-grid scale turbulence model. Comparisons of experimental data show that the LES results are capable of predicting mean and root-mean-square velocity profiles. The LES results show that the annular swirling flow has a minor impact on the formation of the bluff-body recirculation zone. The vortex structures near the shear layers, visualized by the iso-surface of Q-criterion, display ring structures in non-swirling flow and helical structures in swirling flow near the burner exit. Spectral analysis was employed to predict the occurrence of flow oscillations induced by vortex shedding and precessing vortex core(PVC). In order to extract accurately the unsteady large-scale structures in swirling flow, a three-dimensional proper orthogonal decomposition(POD) method was developed to reconstruct turbulent fluctuating velocity fields. POD analysis reveals that flow fields contain co-existing helical and toroidal shaped coherent structures. The helical structure associated with the PVC is the most energetic dynamic flow structure. The latter toroidal structure associated with vortex shedding has lower energy content which indicates that it is a secondary structure.展开更多
A three-dimensional large eddy simulation (LES) of a spatially developing round jet is carried out in cylindrical coordinates using a dynamic subgrid model with strong inflow instability. Evolutions of large-scale v...A three-dimensional large eddy simulation (LES) of a spatially developing round jet is carried out in cylindrical coordinates using a dynamic subgrid model with strong inflow instability. Evolutions of large-scale vortex structures represented by tangential vortices are obtained and compared with flow visualization. Also presented are three-dimensional spatial evolutions of coherent structure, which are of quasi two-dimensional Kelvin-Helmholtz instability and vortex rings as well as breaking up of the vortex rings with fully three-dimensional characteristics. Predicted results of mean velocity and turbulent intensity agree well with experiments. They are also compared with the results predicted by LES using standard Smagorinsky model and show good self-similarity. Turbulence spectrum of the predicted velocity shows the -5/3 decay for higher wave number, as expected for turbulent round jet flows. In addition, fl-test and y-test are carded out for the instantaneous velocity, showing that the present LES method can successfully predict the hierarchical structure of round jet.展开更多
Large eddy simulation(LES)is used to calculate the in-cylinder turbulent flow field in a direct injection spark ignition(DISI)engine.The computations are carried out for three different maximum valve lifts(MVL)and thr...Large eddy simulation(LES)is used to calculate the in-cylinder turbulent flow field in a direct injection spark ignition(DISI)engine.The computations are carried out for three different maximum valve lifts(MVL)and throughout 100 consecutive engine cycles.The simulated results as well as corresponding particle image velocimetry(PIV)measurement database are analyzed by the proper orthogonal decomposition(POD)method.Through a new developed POD quadruple decomposition the instantaneous in-cylinder flow fields are decomposed into four parts,named mean field,coherent field,transition field and turbulent field,respectively.Then the in-cylinder turbulent flow characteristics and cycle-to-cycle variations(CCV)are studied separately upon the four part flow fields.Results indicate that each part exhibits its specific characteristics and has close connection with others.The mean part contains more than 50%of the total kinetic energy and the energy cascade phenomenon occurs among the four part fields;the coherent field part possesses the highest CCV level which dominates CCV of the bulk flow.In addition,it is observed that a change in MVL affects significantly the in-cylinder flow behavior including CCV,especially for the coherent part.Furthermore,the POD analysis demonstrates that at least 25 sample cycles for the mean velocity and 50 sample cycles for the RMS velocity are necessary for obtaining converged and correct results in CCV.展开更多
The large eddy simulation (LES) of compressible turbulent channel flows at three different Mach numbers is performed in the present work,by extending the dynamic mixed subgrid-scale (SGS) model to compressible flows.T...The large eddy simulation (LES) of compressible turbulent channel flows at three different Mach numbers is performed in the present work,by extending the dynamic mixed subgrid-scale (SGS) model to compressible flows.The turbulent statistics agree well with those from the existing direct numerical simulation (DNS) results,indicating that the LES method established in the present work is reliable.The analysis of the turbulent fluctuations computed by the present LES reveals that the flows considered in this work follow the Morkovin's hypothesis.Thus,the compressibility effects are dominated by the mean field properties,and the relevant statistical ratios are invariant to the variation of Mach number.The near-wall streamwise streaks are more coherent and the spacing between streaks is wider as the Mach number increases.This can be regarded as a direct feature characterized by the compressibility effects.The restrained influences of compressibility effects on the production and dissipation of the turbulence kinetic energy are also identified based on the present LES results.展开更多
The motion of micro-particles with different mass flow rate in the planer turbulent jet flow has been simulated, using LES method to obtain the flow vorticity evolution and Lagrangian method to track micro-particles. ...The motion of micro-particles with different mass flow rate in the planer turbulent jet flow has been simulated, using LES method to obtain the flow vorticity evolution and Lagrangian method to track micro-particles. The re- suits showed that when the flow rate is small, the particles more likely to present in the vortex periphery, the dis- tribution pattern is similar to the flow pattern. When the flow rate is high, some particles will escape from the mo- tion region to the original static region, so that in the jet region, particles are relatively evenly distributed. When the flow field is full developed, the particles average concentration in the y direction affected by the mass flow rate relative slightly, the normalized mean particles concentrations at different flow rate were similar to Gaussian shape.展开更多
基金Supported by the National Natural Science Foundation of China(10772082)~~
文摘The low-Reynolds-number full developed turbulent flow in channels is simulated using large eddy simulation(LES)method with the preconditioned algorithm and the dynamic subgrid-scale model,with a given disturbance in inlet boundary,after a short development section.The inlet Reynolds number based on momentum thickness is 670.The computed results show good agreement with direct numerical simulation(DNS),which include root mean square fluctuated velocity distribution and average velocity distribution.It is also found that the staggered phenomenon of the coherent structures is caused by sub-harmonic.The results clearly show the formation and evolution of horseshoe vortex in the turbulent boundary layer,including horseshoe vortex structure with a pair of streamwise vortexes and one-side leg of horseshoe vortex.Based on the results,the development of the horseshoe-shaped coherent structures is analyzed in turbulent boundary layer.
基金Supported by the National Natural Science Foundation of China(Nos.41206015,41106019)the National Basic Research Program of China(973 Program)(Nos.2011CB403501,2012CB417402)the Fund for Creative Research Groups by NSFC(No.41121064)
文摘A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.
文摘Gas solid two-phase turbulent plane jet is applied to many natural s it uations and in engineering systems. To predict the particle dispersion in the ga s jet is of great importance in industrial applications and in the designing of engineering systems. A large eddy simulation of the two-phase plane jet was con d ucted to investigate the particle dispersion patterns. The particles with Stokes numbers equal to 0 0028, 0 3, 2 5, 28 (corresponding to particle diameter 1 μm , 10 μm, 30 μm, 100 μm, respectively) in \%Re\%=11 300 gas flow were studied. The simulation results of gas phase motion agreed well with previous experimental re sults. And the simulation results of the solid particles motion showed that part icles with different Stokes number have different spatial dispersion; and that p articles with intermediate Stokes number have the largest dispersion ratio.
文摘This paper describes the design of a new kind of miniature abrading sphere, which is magnetically mounted inside a spherical gap and set in rotation pneumatically with air. Large eddy simulation is performed in conjunction with the compressible Smagorinsky model. Minimal temperature variation allows for the assumption of adiabatic walls. Fluid-solid interaction is modeled using the law of the wall for compressible turbulent flow. A parametric study is done to determine optimal geometric layout while taking physical restrictions into account. The resulting optimal configuration is then examined in detail in order to determine demands to be met by the computerized control of the magnetic bearing as well as to quantify the force available to the abrasion process. Finally, a mathematical relation is given that determines available abrasion force depending on standard volumetric flow rate and rotation frequency. The findings presented here provide a basis for further development of smaller versions of the tool.
基金supported by the National Natural Science Foundation of China(Grant Nos.51176178&91441117)the Key Program of National Natural Science Foundation of China(Grant No.50936005)
文摘Large eddy simulations(LES) were performed to study the non-reacting flow fields of a Cambridge swirl burner. The dynamic Smagorinsky eddy viscosity model is used as the sub-grid scale turbulence model. Comparisons of experimental data show that the LES results are capable of predicting mean and root-mean-square velocity profiles. The LES results show that the annular swirling flow has a minor impact on the formation of the bluff-body recirculation zone. The vortex structures near the shear layers, visualized by the iso-surface of Q-criterion, display ring structures in non-swirling flow and helical structures in swirling flow near the burner exit. Spectral analysis was employed to predict the occurrence of flow oscillations induced by vortex shedding and precessing vortex core(PVC). In order to extract accurately the unsteady large-scale structures in swirling flow, a three-dimensional proper orthogonal decomposition(POD) method was developed to reconstruct turbulent fluctuating velocity fields. POD analysis reveals that flow fields contain co-existing helical and toroidal shaped coherent structures. The helical structure associated with the PVC is the most energetic dynamic flow structure. The latter toroidal structure associated with vortex shedding has lower energy content which indicates that it is a secondary structure.
基金supported by the National Natural Science Foundation of China (Grant Nos. 50176027 and 50706021)a grant from the Research Committee of The Hong Kong Polytechnic University (Grant No.G-U294)
文摘A three-dimensional large eddy simulation (LES) of a spatially developing round jet is carried out in cylindrical coordinates using a dynamic subgrid model with strong inflow instability. Evolutions of large-scale vortex structures represented by tangential vortices are obtained and compared with flow visualization. Also presented are three-dimensional spatial evolutions of coherent structure, which are of quasi two-dimensional Kelvin-Helmholtz instability and vortex rings as well as breaking up of the vortex rings with fully three-dimensional characteristics. Predicted results of mean velocity and turbulent intensity agree well with experiments. They are also compared with the results predicted by LES using standard Smagorinsky model and show good self-similarity. Turbulence spectrum of the predicted velocity shows the -5/3 decay for higher wave number, as expected for turbulent round jet flows. In addition, fl-test and y-test are carded out for the instantaneous velocity, showing that the present LES method can successfully predict the hierarchical structure of round jet.
基金supported by the National Natural Science Foundation of China(Grant Nos.51176020 and 51376029)
文摘Large eddy simulation(LES)is used to calculate the in-cylinder turbulent flow field in a direct injection spark ignition(DISI)engine.The computations are carried out for three different maximum valve lifts(MVL)and throughout 100 consecutive engine cycles.The simulated results as well as corresponding particle image velocimetry(PIV)measurement database are analyzed by the proper orthogonal decomposition(POD)method.Through a new developed POD quadruple decomposition the instantaneous in-cylinder flow fields are decomposed into four parts,named mean field,coherent field,transition field and turbulent field,respectively.Then the in-cylinder turbulent flow characteristics and cycle-to-cycle variations(CCV)are studied separately upon the four part flow fields.Results indicate that each part exhibits its specific characteristics and has close connection with others.The mean part contains more than 50%of the total kinetic energy and the energy cascade phenomenon occurs among the four part fields;the coherent field part possesses the highest CCV level which dominates CCV of the bulk flow.In addition,it is observed that a change in MVL affects significantly the in-cylinder flow behavior including CCV,especially for the coherent part.Furthermore,the POD analysis demonstrates that at least 25 sample cycles for the mean velocity and 50 sample cycles for the RMS velocity are necessary for obtaining converged and correct results in CCV.
文摘The large eddy simulation (LES) of compressible turbulent channel flows at three different Mach numbers is performed in the present work,by extending the dynamic mixed subgrid-scale (SGS) model to compressible flows.The turbulent statistics agree well with those from the existing direct numerical simulation (DNS) results,indicating that the LES method established in the present work is reliable.The analysis of the turbulent fluctuations computed by the present LES reveals that the flows considered in this work follow the Morkovin's hypothesis.Thus,the compressibility effects are dominated by the mean field properties,and the relevant statistical ratios are invariant to the variation of Mach number.The near-wall streamwise streaks are more coherent and the spacing between streaks is wider as the Mach number increases.This can be regarded as a direct feature characterized by the compressibility effects.The restrained influences of compressibility effects on the production and dissipation of the turbulence kinetic energy are also identified based on the present LES results.
基金supported by National Natural Science Foundation of China (Grant No. 50976107 and 51206149)National Key Technology R&D Program of China (Grant No. 2009BAF39B01)Zhejiang Science Technology Project (Grant No. 2011C11073)
文摘The motion of micro-particles with different mass flow rate in the planer turbulent jet flow has been simulated, using LES method to obtain the flow vorticity evolution and Lagrangian method to track micro-particles. The re- suits showed that when the flow rate is small, the particles more likely to present in the vortex periphery, the dis- tribution pattern is similar to the flow pattern. When the flow rate is high, some particles will escape from the mo- tion region to the original static region, so that in the jet region, particles are relatively evenly distributed. When the flow field is full developed, the particles average concentration in the y direction affected by the mass flow rate relative slightly, the normalized mean particles concentrations at different flow rate were similar to Gaussian shape.