A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-por...A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-porosity technique was used to represent the solidification,and Darcy law was adopted to simulate the flow in mushy region.The prediction for steel flow and solidification was evaluated by the comparison of two turbulence models.It is found that both Lam-Bremhorst low-Reynolds-number and high-Reynolds-number k-ε models predict the same trend of the steel flow and temperature distribution.However,due to the effect of turbulent flow on heat transfer,the low-Reynolds-number turbulence model predicts longer penetration depth of molten steel in sub-mold region,less shell growth and higher shell surface temperature at the narrow face compared with standard k-ε model.展开更多
This paper presents the gas distribution analysis by injecting air fountain into the containment and simulations with the HYDRAGON code. Turbulence models of standard k-ε(SKE), re-normalization group k-ε(RNG) and a ...This paper presents the gas distribution analysis by injecting air fountain into the containment and simulations with the HYDRAGON code. Turbulence models of standard k-ε(SKE), re-normalization group k-ε(RNG) and a realizable k-ε(RLZ) are used to assess the effects on the gas distribution analysis during a severe accident in a nuclear power plant. By comparing with experimental data,the simulation results of the RNG and SKE turbulence models agree well with the experimental data on the prediction of dimensionless density distributions. The results illustrate that the turbulence model choice had a small effect on the simulation results, particularly the region near to the air fountain source.展开更多
The present study proposes a sub-grid scale model for the one-dimensional Burgers turbulence based on the neuralnetwork and deep learning method.The filtered data of the direct numerical simulation is used to establis...The present study proposes a sub-grid scale model for the one-dimensional Burgers turbulence based on the neuralnetwork and deep learning method.The filtered data of the direct numerical simulation is used to establish thetraining data set,the validation data set,and the test data set.The artificial neural network(ANN)methodand Back Propagation method are employed to train parameters in the ANN.The developed ANN is applied toconstruct the sub-grid scale model for the large eddy simulation of the Burgers turbulence in the one-dimensionalspace.The proposed model well predicts the time correlation and the space correlation of the Burgers turbulence.展开更多
A model that uses the operator splitting technique based on Eulerian-Lagrangian method and embeds the k-ε turbulence mode is developed. The model is used to simulate a two dimensional flow over a backward-facing step...A model that uses the operator splitting technique based on Eulerian-Lagrangian method and embeds the k-ε turbulence mode is developed. The model is used to simulate a two dimensional flow over a backward-facing step with triangular mesh discrete flow field. The convection operator is solved with the characteristic method and the diffusion operator and pressure-Poisson equation are solved by finite element method. The results from the example showed that numerical results are well agreed with the experimental data,and that the method can be adapted to simulate complex turbulent flow with strong non-linear convection. The treatment of boundary conditions is also studied in the paper.展开更多
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 internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence mod...The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved.展开更多
The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations fo...The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations for the electromagnetic braking were used to calculate the velocity distribution in the continuous casting mold of steel. The results show that the revised k-ε model considering the effect of magnetic field application tends to suppress the production of turbulence and difference between the conventional and revised k-e model is small.展开更多
A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carr...A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.展开更多
Numerical study on turbulent mixed convection in inclined plane channels,from 15° to 90° (vertical),was carried out to examine the effect of inclination on fluid flow and heat transfer distributions.The turb...Numerical study on turbulent mixed convection in inclined plane channels,from 15° to 90° (vertical),was carried out to examine the effect of inclination on fluid flow and heat transfer distributions.The turbulent air flows upward or downward into the duct with one wall heated from bottom.Calculation results with several kinds of k-εtype turbulence models were used to compare the experimental data with those in literatures to determine suitable model.The dependents of Nusselt number on the inclination angle of both the buoyancy-aided and buoyancy-opposed flow are discussed.展开更多
In gas turbine engines, laminar-turbulent transition occurs. However, generally, the turbulence models to describe such transition results in too early and too short transition. Combining a turbulence model with a des...In gas turbine engines, laminar-turbulent transition occurs. However, generally, the turbulence models to describe such transition results in too early and too short transition. Combining a turbulence model with a description of intermittency, i.e. the fraction of time the flow is turbulent during the transition phase, can improve it. By letting grow the intermittency from zero to unity, start and evolution of transition can be imposed. In this paper, a method where a dynamic equation of intermittency combining with a two-equation k-ωturbulence model is described. This intermittency factor is a premultiplicator of the turbulent viscosity computed by the turbulence model. Following a suggestion by Menter et al.[1], the start of transition is computed based on local variables.展开更多
An unstructured Reynolds-averaged Navier-Stokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equation turbulence model (k-...An unstructured Reynolds-averaged Navier-Stokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equation turbulence model (k-ω model) is also developed for hybrid grids to compute the turbulence flow. The turbulence flow past NACA0012 airfoil and the double ellipsolids are computed, and the numerical results show that the above methods are very efficient.展开更多
LRN （low-Reynolds number） modifications to the NR （Norris-Reynolds） k-equation turbulence model are proposed and evaluated. The k and e that render the hybrid time scale are determined using the k-transport equati...LRN （low-Reynolds number） modifications to the NR （Norris-Reynolds） k-equation turbulence model are proposed and evaluated. The k and e that render the hybrid time scale are determined using the k-transport equation together with the Bradshaw and other algebraic relations. The eddy-viscosity coefficient Cμ and the empirical damping function are constructed such as to preserve the anisotropic characteristics of turbulence for application to non-equilibrium turbulent flows. The MNR （modified NR） model is applied to calculate two well-documented flows, yielding predictions in good agreement with the DNS （direct numerical simulation） and experimental data. Comparisons demonstrate that the MNR model offers a significant improvement over the original NR model and competitiveness with the Spalart-Allmaras one-equation turbulence model. The performance evaluation dictates that unlike the original NR model, the MNR model can be employed as a single-equation model instead of associating it with the two-layer model of turbulence.展开更多
The numerical simulation based on Reynolds time-averaged equation is one of the approved methods to evaluate the aerodynamic performance of trains in crosswind.However,there are several turbulence models,trains may pr...The numerical simulation based on Reynolds time-averaged equation is one of the approved methods to evaluate the aerodynamic performance of trains in crosswind.However,there are several turbulence models,trains may present different aerodynamic performances in crosswind using different turbulence models.In order to select the most suitable turbulence model,the inter-city express 2(ICE2)model is chosen as a research object,6 different turbulence models are used to simulate the flow characteristics,surface pressure and aerodynamic forces of the train in crosswind,respectively.6 turbulence models are the standard k-ε,Renormalization Group(RNG)k-ε,Realizable k-ε,Shear Stress Transport(SST)k-ω,standard k-ωand Spalart-Allmaras(SPA),respectively.The numerical results and the wind tunnel experimental data are compared.The results show that the most accurate model for predicting the surface pressure of the train is SST k-ω,followed by Realizable k-ε.Compared with the experimental result,the error of the side force coefficient obtained by SST k-ωand Realizable k-εturbulence model is less than 1%.The most accurate prediction for the lift force coefficient is achieved by SST k-ω,followed by RNG k-ε.By comparing 6 different turbulence models,the SST k-ωmodel is most suitable for the numerical simulation of the aerodynamic behavior of trains in crosswind.展开更多
Three-dimensional corner separation is a common phenomenon that significantly affects compressor performance. Turbulence model is still a weakness for RANS method on predicting corner separation flow accurately. In th...Three-dimensional corner separation is a common phenomenon that significantly affects compressor performance. Turbulence model is still a weakness for RANS method on predicting corner separation flow accurately. In the present study, numerical study of corner separation in a linear highly loaded prescribed velocity distribution （PVD） compressor cascade has been investigated using seven frequently used turbulence models. The seven turbulence models include Spalart Allmaras model, standard k-e model, realizable k-e model, standard k-to model, shear stress transport k co model, v2-fmodel and Reynolds stress model. The results of these turbulence models have been compared and analyzed in detail with available experimental data. It is found the standard k-1： model, realizable k-e model, v2-f model and Reynolds stress model can provide reasonable results for predicting three dimensional corner separation in the compressor cascade. The Spalart-Allmaras model, standard k-to model and shear stress transport k-w model overesti- mate corner separation region at incidence of 0°. The turbulence characteristics are discussed and turbulence anisotropy is observed to be stronger in the corner separating region.展开更多
Despite dedicated effort for many decades,statistical description of highly technologically important wall turbulence remains a great challenge.Current models are unfortunately incomplete,or empirical,or qualitative.A...Despite dedicated effort for many decades,statistical description of highly technologically important wall turbulence remains a great challenge.Current models are unfortunately incomplete,or empirical,or qualitative.After a review of the existing theories of wall turbulence,we present a new framework,called the structure ensemble dynamics （SED）,which aims at integrating the turbulence dynamics into a quantitative description of the mean flow.The SED theory naturally evolves from a statistical physics understanding of non-equilibrium open systems,such as fluid turbulence, for which mean quantities are intimately coupled with the fluctuation dynamics.Starting from the ensemble-averaged Navier-Stokes（EANS） equations,the theory postulates the existence of a finite number of statistical states yielding a multi-layer picture for wall turbulence.Then,it uses order functions（ratios of terms in the mean momentum as well as energy equations） to characterize the states and transitions between states.Application of the SED analysis to an incompressible channel flow and a compressible turbulent boundary layer shows that the order functions successfully reveal the multi-layer structure for wall-bounded turbulence, which arises as a quantitative extension of the traditional view in terms of sub-layer,buffer layer,log layer and wake. Furthermore,an idea of using a set of hyperbolic functions for modeling transitions between layers is proposed for a quantitative model of order functions across the entire flow domain.We conclude that the SED provides a theoretical framework for expressing the yet-unknown effects of fluctuation structures on the mean quantities,and offers new methods to analyze experimental and simulation data.Combined with asymptotic analysis,it also offers a way to evaluate convergence of simulations.The SED approach successfully describes the dynamics at both momentum and energy levels, in contrast with all prevalent approaches describing the mean velocity profile only.Moreover,the SED theoretical framework is general,independent of the flow system to study, while the actual functional form of the order functions may vary from flow to flow.We assert that as the knowledge of order functions is accumulated and as more flows are analyzed, new principles（such as hierarchy,symmetry,group invariance,etc.） governing the role of turbulent structures in the mean flow properties will be clarified and a viable theory of turbulence might emerge.展开更多
A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibrati...A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R ＆ D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.展开更多
A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kin...A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.展开更多
Presently developed two-phase turbulence models under-predict the gas turbulent fluctuation, because their turbulence modification models cannot fully reflect the effect of particles. In this paper, a two-time-scale d...Presently developed two-phase turbulence models under-predict the gas turbulent fluctuation, because their turbulence modification models cannot fully reflect the effect of particles. In this paper, a two-time-scale dis- sipation model of turbulence modification, developed for the two-phase velocity correlation and for the dissipation rate of gas turbulent kinetic energy, is proposed and used to simulate sudden-expansion and swirling gas-particle flows. The proposed two-time scale model gives better results than the single-time scale model. Besides, a gas tur- bulence augmentation model accounting for the finite-size particle wake effect in the gas Reynolds stress equation is proposed. The proposed turbulence modification models are used to simulate two-phase pipe flows. It can prop- erly predict both turbulence reduction and turbulence enhancement for a certain size of particles observed in ex- periments.展开更多
In previous attempts of rational subgrid-scale (SGS) modeling by employing the Kolmogorov equation of filtered (KEF) quantities, it was necessary to assume that the resolved-scale second-order structure function is st...In previous attempts of rational subgrid-scale (SGS) modeling by employing the Kolmogorov equation of filtered (KEF) quantities, it was necessary to assume that the resolved-scale second-order structure function is stationary. Forced isotropic turbulence is often used as a framework for establishing and validating such SGS models based on stationary restrictions, for it generates statistical stationary samples. However, traditional forcing method at low wavenumbers cannot provide an analytic form of forcing term for a complete KEF in physical space, which has been illustrated to be essential in the modeling of such SGS models. Thus, an alternative forcing method giving an analytic forcing term in physical space is needed for rational SGS modeling. Giving an analytic linear driving term in physical space, linearly forced isotropic turbulence should be considered an ideal theoretical framework for rational SGS modeling. In this paper, we demonstrate the feasibility of establishing a rational SGS model with stationary restriction based on linearly forced isotropic turbulence. The performance of this rational SGS model is validated. We, therefore, propose the use of linearly forced isotropic turbulence as a complement to free-decaying isotropic turbulence and low-wavenumber forced isotropic turbulence for SGS model validations.展开更多
The mechanism of local scour around submarine pipelines is studied numerically based on a renormalized group （RNG） turbulence model. To validate the numerical model, the equilibrium profiles of local scour for two c...The mechanism of local scour around submarine pipelines is studied numerically based on a renormalized group （RNG） turbulence model. To validate the numerical model, the equilibrium profiles of local scour for two cases are simulated and compared with the experimental data. It shows that the RNG turbulence model can give an appropriate prediction for the configuration of equilibrium scour hole, and it is applicable to this situation. The local scour mechanism around submarine pipelines including the flow structure, shear stress distribution and pressure field is then analyzed and compared with experiments. For further comparison and validation, especially for the flow structure, a numerical calculation employing the large eddy simulation （LES） is also conducted. The numerical results of RNG demonstrate that the critical factor governing the equilibrium profile is the seabed shear stress distribution in the case of bed load sediment transport, and the two-equation RNG turbulence model coupled with the law of wall is capable of giving a satisfying estimation for the bed shear stress. Moreover, the piping phenomena due to the great difference of pressure between the upstream and downstream parts of pipelines and the vortex structure around submarine pipelines are also simulated successfully, which are believed to be the important factor that lead to the onset of local scour.展开更多
文摘A steady three-dimensional fluid flow and solidification model was built based on CFD software by high-Reynolds-number and Lam-Bremhorst low-Reynolds-number k-ε model.During the simulation,the fixed-grid enthalpy-porosity technique was used to represent the solidification,and Darcy law was adopted to simulate the flow in mushy region.The prediction for steel flow and solidification was evaluated by the comparison of two turbulence models.It is found that both Lam-Bremhorst low-Reynolds-number and high-Reynolds-number k-ε models predict the same trend of the steel flow and temperature distribution.However,due to the effect of turbulent flow on heat transfer,the low-Reynolds-number turbulence model predicts longer penetration depth of molten steel in sub-mold region,less shell growth and higher shell surface temperature at the narrow face compared with standard k-ε model.
基金support of the National key Lab of Reactor System Design Technology Chengdu,Chinathe Chinese Scholarship Council for the award of Doctoral study
文摘This paper presents the gas distribution analysis by injecting air fountain into the containment and simulations with the HYDRAGON code. Turbulence models of standard k-ε(SKE), re-normalization group k-ε(RNG) and a realizable k-ε(RLZ) are used to assess the effects on the gas distribution analysis during a severe accident in a nuclear power plant. By comparing with experimental data,the simulation results of the RNG and SKE turbulence models agree well with the experimental data on the prediction of dimensionless density distributions. The results illustrate that the turbulence model choice had a small effect on the simulation results, particularly the region near to the air fountain source.
基金supported by the National Key R&D Program of China(Grant No.2022YFB3303500).
文摘The present study proposes a sub-grid scale model for the one-dimensional Burgers turbulence based on the neuralnetwork and deep learning method.The filtered data of the direct numerical simulation is used to establish thetraining data set,the validation data set,and the test data set.The artificial neural network(ANN)methodand Back Propagation method are employed to train parameters in the ANN.The developed ANN is applied toconstruct the sub-grid scale model for the large eddy simulation of the Burgers turbulence in the one-dimensionalspace.The proposed model well predicts the time correlation and the space correlation of the Burgers turbulence.
基金supported by the Key Technologies Research&Development Program of China in the11th Five-Year Plan(Grant No.2008BAB29B04)
文摘A model that uses the operator splitting technique based on Eulerian-Lagrangian method and embeds the k-ε turbulence mode is developed. The model is used to simulate a two dimensional flow over a backward-facing step with triangular mesh discrete flow field. The convection operator is solved with the characteristic method and the diffusion operator and pressure-Poisson equation are solved by finite element method. The results from the example showed that numerical results are well agreed with the experimental data,and that the method can be adapted to simulate complex turbulent flow with strong non-linear convection. The treatment of boundary conditions is also studied in the paper.
基金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.
基金Projects(59375211,10771178,10676031) supported by the National Natural Science Foundation of ChinaProject(07A068) supported by the Key Project of Hunan Education CommissionProject(2005CB321702) supported by the National Key Basic Research Program of China
文摘The internal turbulent flow in conical diffuser is a very complicated adverse pressure gradient flow.DLR k-ε turbulence model was adopted to study it.The every terms of the Laplace operator in DLR k-ε turbulence model and pressure Poisson equation were discretized by upwind difference scheme.A new full implicit difference scheme of 5-point was constructed by using finite volume method and finite difference method.A large sparse matrix with five diagonals was formed and was stored by three arrays of one dimension in a compressed mode.General iterative methods do not work wel1 with large sparse matrix.With algebraic multigrid method(AMG),linear algebraic system of equations was solved and the precision was set at 10-6.The computation results were compared with the experimental results.The results show that the computation results have a good agreement with the experiment data.The precision of computational results and numerical simulation efficiency are greatly improved.
文摘The research is motivated by the ongoing the electromagnetic continuous casting of molten metal. The revised k-ε model considering the effect of magnetic field application was derived. The specific model equations for the electromagnetic braking were used to calculate the velocity distribution in the continuous casting mold of steel. The results show that the revised k-ε model considering the effect of magnetic field application tends to suppress the production of turbulence and difference between the conventional and revised k-e model is small.
文摘A two-equation turbulence model has been dereloped for predicting two-phase flow the two equations describe the conserration of turbulence kinetic energy and dissipation rate of that energy for the incompressible carrier fluid in a two-phase flow The continuity, the momentum, K and εequations are modeled. In this model,the solid-liquid slip veloeites, the particle-particte interactions and the interactions between two phases are considered,The sandy water pipe turbulent flows are sueeessfuly predicted by this turbulince model.
文摘Numerical study on turbulent mixed convection in inclined plane channels,from 15° to 90° (vertical),was carried out to examine the effect of inclination on fluid flow and heat transfer distributions.The turbulent air flows upward or downward into the duct with one wall heated from bottom.Calculation results with several kinds of k-εtype turbulence models were used to compare the experimental data with those in literatures to determine suitable model.The dependents of Nusselt number on the inclination angle of both the buoyancy-aided and buoyancy-opposed flow are discussed.
文摘In gas turbine engines, laminar-turbulent transition occurs. However, generally, the turbulence models to describe such transition results in too early and too short transition. Combining a turbulence model with a description of intermittency, i.e. the fraction of time the flow is turbulent during the transition phase, can improve it. By letting grow the intermittency from zero to unity, start and evolution of transition can be imposed. In this paper, a method where a dynamic equation of intermittency combining with a two-equation k-ωturbulence model is described. This intermittency factor is a premultiplicator of the turbulent viscosity computed by the turbulence model. Following a suggestion by Menter et al.[1], the start of transition is computed based on local variables.
文摘An unstructured Reynolds-averaged Navier-Stokes flow solver using the finite volume method is studied. The spatial discretisation is based on the Osher approximate Riemann solvers. A two-equation turbulence model (k-ω model) is also developed for hybrid grids to compute the turbulence flow. The turbulence flow past NACA0012 airfoil and the double ellipsolids are computed, and the numerical results show that the above methods are very efficient.
文摘LRN （low-Reynolds number） modifications to the NR （Norris-Reynolds） k-equation turbulence model are proposed and evaluated. The k and e that render the hybrid time scale are determined using the k-transport equation together with the Bradshaw and other algebraic relations. The eddy-viscosity coefficient Cμ and the empirical damping function are constructed such as to preserve the anisotropic characteristics of turbulence for application to non-equilibrium turbulent flows. The MNR （modified NR） model is applied to calculate two well-documented flows, yielding predictions in good agreement with the DNS （direct numerical simulation） and experimental data. Comparisons demonstrate that the MNR model offers a significant improvement over the original NR model and competitiveness with the Spalart-Allmaras one-equation turbulence model. The performance evaluation dictates that unlike the original NR model, the MNR model can be employed as a single-equation model instead of associating it with the two-layer model of turbulence.
基金Supported by National Natural Science Foundation of China(Grant No.51605397)Sichuan Provincial Science and Technology Program of China(Grant No.2019YJ0227)Self-determined Project of State Key Laboratory of Traction Power(Grant No.2019TPL_T02)
文摘The numerical simulation based on Reynolds time-averaged equation is one of the approved methods to evaluate the aerodynamic performance of trains in crosswind.However,there are several turbulence models,trains may present different aerodynamic performances in crosswind using different turbulence models.In order to select the most suitable turbulence model,the inter-city express 2(ICE2)model is chosen as a research object,6 different turbulence models are used to simulate the flow characteristics,surface pressure and aerodynamic forces of the train in crosswind,respectively.6 turbulence models are the standard k-ε,Renormalization Group(RNG)k-ε,Realizable k-ε,Shear Stress Transport(SST)k-ω,standard k-ωand Spalart-Allmaras(SPA),respectively.The numerical results and the wind tunnel experimental data are compared.The results show that the most accurate model for predicting the surface pressure of the train is SST k-ω,followed by Realizable k-ε.Compared with the experimental result,the error of the side force coefficient obtained by SST k-ωand Realizable k-εturbulence model is less than 1%.The most accurate prediction for the lift force coefficient is achieved by SST k-ω,followed by RNG k-ε.By comparing 6 different turbulence models,the SST k-ωmodel is most suitable for the numerical simulation of the aerodynamic behavior of trains in crosswind.
基金supported by the National Natural Science Foundation of China(No.51376001,No.51420105008,No.51306013,No.51136003)the National Basic Research Program of China(2012CB720205,2014CB046405)+2 种基金the Beijing Higher Education Young Elite Teacher Projectthe Fundamental Research Funds for the Central Universitiessupported by the Innovation Foundation of BUAA for Ph.D.Graduates
文摘Three-dimensional corner separation is a common phenomenon that significantly affects compressor performance. Turbulence model is still a weakness for RANS method on predicting corner separation flow accurately. In the present study, numerical study of corner separation in a linear highly loaded prescribed velocity distribution （PVD） compressor cascade has been investigated using seven frequently used turbulence models. The seven turbulence models include Spalart Allmaras model, standard k-e model, realizable k-e model, standard k-to model, shear stress transport k co model, v2-fmodel and Reynolds stress model. The results of these turbulence models have been compared and analyzed in detail with available experimental data. It is found the standard k-1： model, realizable k-e model, v2-f model and Reynolds stress model can provide reasonable results for predicting three dimensional corner separation in the compressor cascade. The Spalart-Allmaras model, standard k-to model and shear stress transport k-w model overesti- mate corner separation region at incidence of 0°. The turbulence characteristics are discussed and turbulence anisotropy is observed to be stronger in the corner separating region.
基金supported by the National Natural Science Foundation of China(90716008)the National Basic Research Program of China(2009CB724100).
文摘Despite dedicated effort for many decades,statistical description of highly technologically important wall turbulence remains a great challenge.Current models are unfortunately incomplete,or empirical,or qualitative.After a review of the existing theories of wall turbulence,we present a new framework,called the structure ensemble dynamics （SED）,which aims at integrating the turbulence dynamics into a quantitative description of the mean flow.The SED theory naturally evolves from a statistical physics understanding of non-equilibrium open systems,such as fluid turbulence, for which mean quantities are intimately coupled with the fluctuation dynamics.Starting from the ensemble-averaged Navier-Stokes（EANS） equations,the theory postulates the existence of a finite number of statistical states yielding a multi-layer picture for wall turbulence.Then,it uses order functions（ratios of terms in the mean momentum as well as energy equations） to characterize the states and transitions between states.Application of the SED analysis to an incompressible channel flow and a compressible turbulent boundary layer shows that the order functions successfully reveal the multi-layer structure for wall-bounded turbulence, which arises as a quantitative extension of the traditional view in terms of sub-layer,buffer layer,log layer and wake. Furthermore,an idea of using a set of hyperbolic functions for modeling transitions between layers is proposed for a quantitative model of order functions across the entire flow domain.We conclude that the SED provides a theoretical framework for expressing the yet-unknown effects of fluctuation structures on the mean quantities,and offers new methods to analyze experimental and simulation data.Combined with asymptotic analysis,it also offers a way to evaluate convergence of simulations.The SED approach successfully describes the dynamics at both momentum and energy levels, in contrast with all prevalent approaches describing the mean velocity profile only.Moreover,the SED theoretical framework is general,independent of the flow system to study, while the actual functional form of the order functions may vary from flow to flow.We assert that as the knowledge of order functions is accumulated and as more flows are analyzed, new principles（such as hierarchy,symmetry,group invariance,etc.） governing the role of turbulent structures in the mean flow properties will be clarified and a viable theory of turbulence might emerge.
基金supported by the National Natural Science Foundation of China(Grant No.51079157)
文摘A growing interest has been devoted to the contra-rotating propellers （CRPs） due to their high propulsive efficiency, torque balance, low fuel consumption, low cavitations, low noise performance and low hull vibration. Compared with the single-screw system, it is more difficult for the open water performance prediction because forward and aft propellers interact with each other and generate a more complicated flow field around the CRPs system. The current work focuses on the open water performance prediction of contra-rotating propellers by RANS and sliding mesh method considering the effect of computational time step size and turbulence model. The validation study has been performed on two sets of contra-rotating propellers developed by David W Taylor Naval Ship R ＆ D center. Compared with the experimental data, it shows that RANS with sliding mesh method and SST k-ω turbulence model has a good precision in the open water performance prediction of contra-rotating propellers, and small time step size can improve the level of accuracy for CRPs with the same blade number of forward and aft propellers, while a relatively large time step size is a better choice for CRPs with different blade numbers.
基金the Special Funds for Major State Basic Research of China(G-1999-0222-08)the National Natural Science Foundation of China(50376004)Ph.D.Program Foundation,Ministry of Education of China(20030007028)
文摘A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.
基金State Key Development Program for Basic Research of China (No.2006CB200305), the National Natural Sci-ence Foundation of China (No.50376004), and Ph.D. Program Foundation of Ministry of Education of China (No.20030007028).
文摘Presently developed two-phase turbulence models under-predict the gas turbulent fluctuation, because their turbulence modification models cannot fully reflect the effect of particles. In this paper, a two-time-scale dis- sipation model of turbulence modification, developed for the two-phase velocity correlation and for the dissipation rate of gas turbulent kinetic energy, is proposed and used to simulate sudden-expansion and swirling gas-particle flows. The proposed two-time scale model gives better results than the single-time scale model. Besides, a gas tur- bulence augmentation model accounting for the finite-size particle wake effect in the gas Reynolds stress equation is proposed. The proposed turbulence modification models are used to simulate two-phase pipe flows. It can prop- erly predict both turbulence reduction and turbulence enhancement for a certain size of particles observed in ex- periments.
基金the National Natural Science Foundation of China (Grant 11772128)the Fundamental Research Funds for the Central Universities (Grants 2017MS022 and 2018ZD09).
文摘In previous attempts of rational subgrid-scale (SGS) modeling by employing the Kolmogorov equation of filtered (KEF) quantities, it was necessary to assume that the resolved-scale second-order structure function is stationary. Forced isotropic turbulence is often used as a framework for establishing and validating such SGS models based on stationary restrictions, for it generates statistical stationary samples. However, traditional forcing method at low wavenumbers cannot provide an analytic form of forcing term for a complete KEF in physical space, which has been illustrated to be essential in the modeling of such SGS models. Thus, an alternative forcing method giving an analytic forcing term in physical space is needed for rational SGS modeling. Giving an analytic linear driving term in physical space, linearly forced isotropic turbulence should be considered an ideal theoretical framework for rational SGS modeling. In this paper, we demonstrate the feasibility of establishing a rational SGS model with stationary restriction based on linearly forced isotropic turbulence. The performance of this rational SGS model is validated. We, therefore, propose the use of linearly forced isotropic turbulence as a complement to free-decaying isotropic turbulence and low-wavenumber forced isotropic turbulence for SGS model validations.
基金supported by the Program for Changjiang Scholars and Innovative Research Team in University of China under contract No,IRT0420the National Natural Science Foundation of China under contract No.50409015.
文摘The mechanism of local scour around submarine pipelines is studied numerically based on a renormalized group （RNG） turbulence model. To validate the numerical model, the equilibrium profiles of local scour for two cases are simulated and compared with the experimental data. It shows that the RNG turbulence model can give an appropriate prediction for the configuration of equilibrium scour hole, and it is applicable to this situation. The local scour mechanism around submarine pipelines including the flow structure, shear stress distribution and pressure field is then analyzed and compared with experiments. For further comparison and validation, especially for the flow structure, a numerical calculation employing the large eddy simulation （LES） is also conducted. The numerical results of RNG demonstrate that the critical factor governing the equilibrium profile is the seabed shear stress distribution in the case of bed load sediment transport, and the two-equation RNG turbulence model coupled with the law of wall is capable of giving a satisfying estimation for the bed shear stress. Moreover, the piping phenomena due to the great difference of pressure between the upstream and downstream parts of pipelines and the vortex structure around submarine pipelines are also simulated successfully, which are believed to be the important factor that lead to the onset of local scour.
