We investigate the role of extended intrinsic mean spin tensor introduced in this work for turbulence modelling in a non-inertial frame of reference. It is described by the Euclidean group of transformations and, in p...We investigate the role of extended intrinsic mean spin tensor introduced in this work for turbulence modelling in a non-inertial frame of reference. It is described by the Euclidean group of transformations and, in particular, its significance and importance in the approach of the algebraic Reynolds stress modelling, such as in a nonlinear K-ε model. To this end and for illustration of the effect of extended intrinsic spin tensor on turbulence modelling, we examine several recently developed nonlinear K-ε models and compare their performance in predicting the homogeneous turbulent shear flow in a rotating frame of reference with LES data. Our results and analysis indicate that, only if the deficiencies of these models and the like be well understood and properly corrected, may in the near future, more sophisticated nonlinear K-ε models be developed to better predict complex turbulent flows in a non-inertial frame of reference.展开更多
A stratified wake has multiple flow regimes,and exhibits different behaviors in these regimes due to the competing physical effects of momentum and buoyancy.This work aims at automated classification of the weakly and...A stratified wake has multiple flow regimes,and exhibits different behaviors in these regimes due to the competing physical effects of momentum and buoyancy.This work aims at automated classification of the weakly and the strongly stratified turbulence regimes based on information available in a full Reynolds stress model.First,we generate a direct numerical simulation database with Reynolds numbers from 10,000 to 50,000 and Froude numbers from 2 to 50.Order(100)independent realizations of temporally evolving wakes are computed to get converged statistics.Second,we train a linear logistic regression classifier with weight thresholding for automated flow regime classification.The classifier is designed to identify the physics critical to classification.Trained against data at one flow condition,the classifier is found to generalize well to other Reynolds and Froude numbers.The results show that the physics governing wake evolution is universal,and that the classifier captures that physics.展开更多
The Speziale,Sarkar and Gatski Reynolds Stress Model(SSG RSM)is utilized to simulate the fluid dynamics in a full baffled stirred tank with a Rushton turbine impeller.Four levels of grid resolutions are chosen to dete...The Speziale,Sarkar and Gatski Reynolds Stress Model(SSG RSM)is utilized to simulate the fluid dynamics in a full baffled stirred tank with a Rushton turbine impeller.Four levels of grid resolutions are chosen to determine an optimised number of grids for further simulations.CFD model data in terms of the flow field,trailing vortex,and the power number are compared with published experimental results.The comparison shows that the global fluid dynamics throughout the stirred tank and the local characteristics of trailing vortices near the blade tips can be captured by the SSG RSM.The predicted mean velocity components in axial,radial and tangential direction are also in good agreement with experiment data.The power number predicted is quite close to the designed value,which demonstrates that this model can accurately calculate the power number in the stirred tank.Therefore,the simulation by using a combination of SSG RSM and MRF impeller rotational model can accurately model turbulent fluid flow in the stirred tank,and it offers an alternative method for design and optimisation of stirred tanks.展开更多
Propeller cavitation is a problematic issue because of its negative effects, such as performances losses, noise,vibration and erosion. Numerical methodology is an effective and efficient technical tool for the study o...Propeller cavitation is a problematic issue because of its negative effects, such as performances losses, noise,vibration and erosion. Numerical methodology is an effective and efficient technical tool for the study of propeller cavitation, however, it is hard to capture tip-vortex cavitation in the previous work by using common turbulence models based on turbulent-viscosity hypothesis. In this work, the Reynolds-Averaged Naiver-Stokes(RANS)approach, adopting the Reynolds stress turbulence model(RSM), is taken to study the unsteady characteristics of the cavitation on the four-bladed INSEAN E779 A model propeller. The numerical simulation was carried out using the commercial CFD software ANSYS Fluent 14.0. One kind of uniform wake flow and two kinds of nonuniform wake flows are considered here. The results in the uniform flow show a good agreement with previous experimental results on both the sheet cavitation and the tip vortex cavitation and prove the ability of the RSM on capturing the tip vortex cavitation. Two kinds of nonuniform wake flows are designed based on the previous experimental researches and the unsteady characteristics of the propeller cavitation are analyzed by comparing the results in the uniform and two nonuniform wake flows together.展开更多
Reynolds-Averaged Navier-Stokes(RANS) Computational Fluid Dynamics(CFD) has been widely used in compressor design and analysis. However, reasonable prediction of compressor flow and its impact on compressor performanc...Reynolds-Averaged Navier-Stokes(RANS) Computational Fluid Dynamics(CFD) has been widely used in compressor design and analysis. However, reasonable prediction of compressor flow and its impact on compressor performance remains challenging. In this study, Menter’s Shear Stress Transport(SST) model and its variants, as well as the ω-based Reynolds stress model(Stress-BSL) are assessed. For a single rotor(Rotor 67), under the peak efficiency operating condition, all studied turbulence models predict its performance with reasonable accuracy;under the off-design conditions, SST with Helicity correction(SST-Helicity) shows superiority in predicting the effect of flow on the spanwise distribution of aerodynamic parameters. For Darmstadt’s 1.5-stage transonic axial compressor, SST-Helicity outperforms SST, SST with the Quadratic Constitutive Relation(SST-QCR) and Stress-BSL in predicting the performance as well as the spanwise distribution of aerodynamic parameters. At the design rotating speed, the stall margin given by SST-Helicity(20.90%) is the closest to the experimental measurement(24.81%), which is more than twice that by SST(8.71%) and 1.5 times that by SST-QCR(14.14%). This paper demonstrates that SSTHelicity model, together with a good quality and sufficiently refined grid, can capture the compressor flow features with reasonable accuracy, which results in a credible prediction of compressor performance and stage matching.展开更多
Floating vegetation island(FVI)provides an effective way to remove excessive nutrition and pollutants in rivers.The Reynolds stress model(RSM)is employed to investigate the hydrodynamic characteristics induced by vari...Floating vegetation island(FVI)provides an effective way to remove excessive nutrition and pollutants in rivers.The Reynolds stress model(RSM)is employed to investigate the hydrodynamic characteristics induced by varied canopy densities of FVI in an open channel.In longitudinal direction,four regions are subdivided according to the flow development process:upstream adjustment region(LUD),diverging flow region(LDF),shear layer growth region(LSD),and flilly developed region.The increasing canopy density accelerates the flow adjustment in the diverging flow region and shear layer growth region,signaling a shorter distance to reach an equilibrium stage,while LUD keeps a constant.The vertical profiles of the normalized velocity are found to be self-similar downstream of the diverging flow region.In the vertical direction,the streamwise velocity profiles in the mixing layer collapse for all densities and obey the hyperbolic tangent law.Normalized penetration depth into the canopy was fitted as a function of dimensionless canopy density given by δc/hc=0.404(CDahc)^-0.316.This finding indicates a large space for rapid water renewal between the canopy region and the underlying water driven by the shear-scale vortices.In the lateral direction,the intensification of secondary current and the increasing number of secondary current cells with increasing canopy density reveal that dense floating canopies contribute to strong momentum exchange.The centers of vortices move as canopy density increases,while the vortices in canopy region do not merge with those in the gap region,as limited by the height and width of the canopy region.The distribution of longitudinal velocity in the transects is significantly influenced by secondary current.展开更多
The oxidation reactor plays a key role in producing rutile titanium dioxide (TiO2) from vapor-phase titanium tetrachloride (TiCl4) by employing a swirling flow operation for enhanced gas mixing. This work aims to ...The oxidation reactor plays a key role in producing rutile titanium dioxide (TiO2) from vapor-phase titanium tetrachloride (TiCl4) by employing a swirling flow operation for enhanced gas mixing. This work aims to understand the effect of reactor configuration on the 3-D swirling flow field using computational fluid dynamics (CFD) simulation. Considering the anisotropic turbulence involved, the Reynolds stress model is applied to describe the complex swirling flow together with the cross-flow mixing of gases. The results show significant effect of the flow angle between the wall jet of air stream (representing TiCl4 in practice) and the axial direction on the initial flow field of cross-flow mixing, where 60° gives smooth profiles of axial velocity development while 90° may provide the fastest mixing between the jet and the axial bulk flow. The pipe shape for the reaction and developing zone, i.e., straight, expanding and shrinking, shows slight influence on the hydrodynamics.展开更多
This paper presents numerical and experimental results of the flow in the tipclearance region of the guide vane row in hydraulic turbine. The 3-D Navier-Stokes equations wereemployed to model the flow in end clearance...This paper presents numerical and experimental results of the flow in the tipclearance region of the guide vane row in hydraulic turbine. The 3-D Navier-Stokes equations wereemployed to model the flow in end clearance region of guide vane cascade, the Reynolds stressdifferential model was used for turbulence closure, and the body-fitted curvilinear coordinates andthe SIMPLE! algorithm were adopted. The governing equations were discretized with the non-staggeredgrids by means of the finite volume method. Detailed comparison of hydrodynamical characteristics ofguide vane in hydraulic turbine with or without tip gap was made. Special attention was paid to theinfluence of leakage flow on the main flow and to the movement of tip leakage flow in the end guidevane. The position and strength of the roll-vortex on the sides of guide vane in were determined.The numerical solutions agree with the experimental results obtained by particle image ve-locimetry.The results help to clarify the loss, wear and cavi-tation erosion between the guide vane andrings, especially for those used in the Yellow River which has a high sediment content.展开更多
By using the Reynolds Stress Closure Model (RSM), turbulentCounter-Gradient-Transport (CGT) phenomenon was numerically investigated in asymmetric flow with ajet, and the computational results were compared with experi...By using the Reynolds Stress Closure Model (RSM), turbulentCounter-Gradient-Transport (CGT) phenomenon was numerically investigated in asymmetric flow with ajet, and the computational results were compared with experimental data. The computational resultsshow that the negative turbulent energy production only appears at some certain stations in CGTregion, this fact indicates that the CGT phenomenon exists more widely than the negative turbulentenergy production; while the CGT region exists all along, it gradually shrinks in the favorablepressure gradient zone until the position of the wing central part is reached, where it vanishes,but it appears in the adverse pressure gradient region; in addition, the location in the flow whereuv = 0 switched sides, relative to where partial deriv U/partial deriv y = 0, from favorablepressure gradient to adverse pressure gradient. The pressure gradient takes an important effect onthe region of negative turbulent energy production and CGT.展开更多
The Reynolds analogy concept has been used in almost all turbulent reacting flow RANS(Reynoldsaveraged Navier–Stokes)simulations,where the turbulence scalar transfers in flow fields are calculated based on the modele...The Reynolds analogy concept has been used in almost all turbulent reacting flow RANS(Reynoldsaveraged Navier–Stokes)simulations,where the turbulence scalar transfers in flow fields are calculated based on the modeled turbulence momentum transfer.This concept,applied to a lean premixed combustion system,was assessed in this paper in terms of exit temperature distribution.Because of the isotropic assumption involved in this analogy,the prediction in some flow condition,such as jet cross flow mixing,would be inaccurate.In this study,using Flamelet Generated Manifold as reaction model,some of the numerical results,obtained from an annular combustor configuration with the turbulent Schmidt number varying from 0.85 to 0.2,were presented and compared with a benchmark atmospheric test results.It was found that the Schmidt numberσt in mean mass fraction f transport equation had significant effect on dilution air mixing process.The mixing between dilution air and reaction products from the primary zone obviously improved asσt decreased on the combustor exit surface.Meanwhile,the sensitivity ofσt in three turbulence models including Realizable k-ε,SST(Shear Stress Transport)and RSM(Reynolds Stress Model)has been compared as well.Since the calculation method of eddy viscosity was different within these three models,RSM was proved to be less sensitive than another two models and can guarantee the best prediction of mixing process condition.On the other hand,the results of dilution air mixing were almost independent of Schmidt number Sct in progress variable c transport equation.This study suggested that for accurate prediction of combustor exit temperature distribution in steady state reacting flow simulation,the turbulent Schmidt number in steady state simulation should be modified to cater to dilution air mixing process.展开更多
文摘We investigate the role of extended intrinsic mean spin tensor introduced in this work for turbulence modelling in a non-inertial frame of reference. It is described by the Euclidean group of transformations and, in particular, its significance and importance in the approach of the algebraic Reynolds stress modelling, such as in a nonlinear K-ε model. To this end and for illustration of the effect of extended intrinsic spin tensor on turbulence modelling, we examine several recently developed nonlinear K-ε models and compare their performance in predicting the homogeneous turbulent shear flow in a rotating frame of reference with LES data. Our results and analysis indicate that, only if the deficiencies of these models and the like be well understood and properly corrected, may in the near future, more sophisticated nonlinear K-ε models be developed to better predict complex turbulent flows in a non-inertial frame of reference.
基金This work is supported by ONR,contract N000142012315.
文摘A stratified wake has multiple flow regimes,and exhibits different behaviors in these regimes due to the competing physical effects of momentum and buoyancy.This work aims at automated classification of the weakly and the strongly stratified turbulence regimes based on information available in a full Reynolds stress model.First,we generate a direct numerical simulation database with Reynolds numbers from 10,000 to 50,000 and Froude numbers from 2 to 50.Order(100)independent realizations of temporally evolving wakes are computed to get converged statistics.Second,we train a linear logistic regression classifier with weight thresholding for automated flow regime classification.The classifier is designed to identify the physics critical to classification.Trained against data at one flow condition,the classifier is found to generalize well to other Reynolds and Froude numbers.The results show that the physics governing wake evolution is universal,and that the classifier captures that physics.
基金support from the Major State Basic Research Development Program of China(973 Program,Grant No.2005CB221205)。
文摘The Speziale,Sarkar and Gatski Reynolds Stress Model(SSG RSM)is utilized to simulate the fluid dynamics in a full baffled stirred tank with a Rushton turbine impeller.Four levels of grid resolutions are chosen to determine an optimised number of grids for further simulations.CFD model data in terms of the flow field,trailing vortex,and the power number are compared with published experimental results.The comparison shows that the global fluid dynamics throughout the stirred tank and the local characteristics of trailing vortices near the blade tips can be captured by the SSG RSM.The predicted mean velocity components in axial,radial and tangential direction are also in good agreement with experiment data.The power number predicted is quite close to the designed value,which demonstrates that this model can accurately calculate the power number in the stirred tank.Therefore,the simulation by using a combination of SSG RSM and MRF impeller rotational model can accurately model turbulent fluid flow in the stirred tank,and it offers an alternative method for design and optimisation of stirred tanks.
基金financially supported by the National Natural Science Foundation of China (Grant Nos. 11872247,11572194,and11472174)。
文摘Propeller cavitation is a problematic issue because of its negative effects, such as performances losses, noise,vibration and erosion. Numerical methodology is an effective and efficient technical tool for the study of propeller cavitation, however, it is hard to capture tip-vortex cavitation in the previous work by using common turbulence models based on turbulent-viscosity hypothesis. In this work, the Reynolds-Averaged Naiver-Stokes(RANS)approach, adopting the Reynolds stress turbulence model(RSM), is taken to study the unsteady characteristics of the cavitation on the four-bladed INSEAN E779 A model propeller. The numerical simulation was carried out using the commercial CFD software ANSYS Fluent 14.0. One kind of uniform wake flow and two kinds of nonuniform wake flows are considered here. The results in the uniform flow show a good agreement with previous experimental results on both the sheet cavitation and the tip vortex cavitation and prove the ability of the RSM on capturing the tip vortex cavitation. Two kinds of nonuniform wake flows are designed based on the previous experimental researches and the unsteady characteristics of the propeller cavitation are analyzed by comparing the results in the uniform and two nonuniform wake flows together.
文摘Reynolds-Averaged Navier-Stokes(RANS) Computational Fluid Dynamics(CFD) has been widely used in compressor design and analysis. However, reasonable prediction of compressor flow and its impact on compressor performance remains challenging. In this study, Menter’s Shear Stress Transport(SST) model and its variants, as well as the ω-based Reynolds stress model(Stress-BSL) are assessed. For a single rotor(Rotor 67), under the peak efficiency operating condition, all studied turbulence models predict its performance with reasonable accuracy;under the off-design conditions, SST with Helicity correction(SST-Helicity) shows superiority in predicting the effect of flow on the spanwise distribution of aerodynamic parameters. For Darmstadt’s 1.5-stage transonic axial compressor, SST-Helicity outperforms SST, SST with the Quadratic Constitutive Relation(SST-QCR) and Stress-BSL in predicting the performance as well as the spanwise distribution of aerodynamic parameters. At the design rotating speed, the stall margin given by SST-Helicity(20.90%) is the closest to the experimental measurement(24.81%), which is more than twice that by SST(8.71%) and 1.5 times that by SST-QCR(14.14%). This paper demonstrates that SSTHelicity model, together with a good quality and sufficiently refined grid, can capture the compressor flow features with reasonable accuracy, which results in a credible prediction of compressor performance and stage matching.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11872285,11672213).
文摘Floating vegetation island(FVI)provides an effective way to remove excessive nutrition and pollutants in rivers.The Reynolds stress model(RSM)is employed to investigate the hydrodynamic characteristics induced by varied canopy densities of FVI in an open channel.In longitudinal direction,four regions are subdivided according to the flow development process:upstream adjustment region(LUD),diverging flow region(LDF),shear layer growth region(LSD),and flilly developed region.The increasing canopy density accelerates the flow adjustment in the diverging flow region and shear layer growth region,signaling a shorter distance to reach an equilibrium stage,while LUD keeps a constant.The vertical profiles of the normalized velocity are found to be self-similar downstream of the diverging flow region.In the vertical direction,the streamwise velocity profiles in the mixing layer collapse for all densities and obey the hyperbolic tangent law.Normalized penetration depth into the canopy was fitted as a function of dimensionless canopy density given by δc/hc=0.404(CDahc)^-0.316.This finding indicates a large space for rapid water renewal between the canopy region and the underlying water driven by the shear-scale vortices.In the lateral direction,the intensification of secondary current and the increasing number of secondary current cells with increasing canopy density reveal that dense floating canopies contribute to strong momentum exchange.The centers of vortices move as canopy density increases,while the vortices in canopy region do not merge with those in the gap region,as limited by the height and width of the canopy region.The distribution of longitudinal velocity in the transects is significantly influenced by secondary current.
文摘The oxidation reactor plays a key role in producing rutile titanium dioxide (TiO2) from vapor-phase titanium tetrachloride (TiCl4) by employing a swirling flow operation for enhanced gas mixing. This work aims to understand the effect of reactor configuration on the 3-D swirling flow field using computational fluid dynamics (CFD) simulation. Considering the anisotropic turbulence involved, the Reynolds stress model is applied to describe the complex swirling flow together with the cross-flow mixing of gases. The results show significant effect of the flow angle between the wall jet of air stream (representing TiCl4 in practice) and the axial direction on the initial flow field of cross-flow mixing, where 60° gives smooth profiles of axial velocity development while 90° may provide the fastest mixing between the jet and the axial bulk flow. The pipe shape for the reaction and developing zone, i.e., straight, expanding and shrinking, shows slight influence on the hydrodynamics.
文摘This paper presents numerical and experimental results of the flow in the tipclearance region of the guide vane row in hydraulic turbine. The 3-D Navier-Stokes equations wereemployed to model the flow in end clearance region of guide vane cascade, the Reynolds stressdifferential model was used for turbulence closure, and the body-fitted curvilinear coordinates andthe SIMPLE! algorithm were adopted. The governing equations were discretized with the non-staggeredgrids by means of the finite volume method. Detailed comparison of hydrodynamical characteristics ofguide vane in hydraulic turbine with or without tip gap was made. Special attention was paid to theinfluence of leakage flow on the main flow and to the movement of tip leakage flow in the end guidevane. The position and strength of the roll-vortex on the sides of guide vane in were determined.The numerical solutions agree with the experimental results obtained by particle image ve-locimetry.The results help to clarify the loss, wear and cavi-tation erosion between the guide vane andrings, especially for those used in the Yellow River which has a high sediment content.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.12002379 and 11972370)the National Key Project(Grant No.GJXM92579).
文摘By using the Reynolds Stress Closure Model (RSM), turbulentCounter-Gradient-Transport (CGT) phenomenon was numerically investigated in asymmetric flow with ajet, and the computational results were compared with experimental data. The computational resultsshow that the negative turbulent energy production only appears at some certain stations in CGTregion, this fact indicates that the CGT phenomenon exists more widely than the negative turbulentenergy production; while the CGT region exists all along, it gradually shrinks in the favorablepressure gradient zone until the position of the wing central part is reached, where it vanishes,but it appears in the adverse pressure gradient region; in addition, the location in the flow whereuv = 0 switched sides, relative to where partial deriv U/partial deriv y = 0, from favorablepressure gradient to adverse pressure gradient. The pressure gradient takes an important effect onthe region of negative turbulent energy production and CGT.
基金funded by the Key Programs of Chinese Academy of Sciences (Project No. ZDRW-CN-2017-2)
文摘The Reynolds analogy concept has been used in almost all turbulent reacting flow RANS(Reynoldsaveraged Navier–Stokes)simulations,where the turbulence scalar transfers in flow fields are calculated based on the modeled turbulence momentum transfer.This concept,applied to a lean premixed combustion system,was assessed in this paper in terms of exit temperature distribution.Because of the isotropic assumption involved in this analogy,the prediction in some flow condition,such as jet cross flow mixing,would be inaccurate.In this study,using Flamelet Generated Manifold as reaction model,some of the numerical results,obtained from an annular combustor configuration with the turbulent Schmidt number varying from 0.85 to 0.2,were presented and compared with a benchmark atmospheric test results.It was found that the Schmidt numberσt in mean mass fraction f transport equation had significant effect on dilution air mixing process.The mixing between dilution air and reaction products from the primary zone obviously improved asσt decreased on the combustor exit surface.Meanwhile,the sensitivity ofσt in three turbulence models including Realizable k-ε,SST(Shear Stress Transport)and RSM(Reynolds Stress Model)has been compared as well.Since the calculation method of eddy viscosity was different within these three models,RSM was proved to be less sensitive than another two models and can guarantee the best prediction of mixing process condition.On the other hand,the results of dilution air mixing were almost independent of Schmidt number Sct in progress variable c transport equation.This study suggested that for accurate prediction of combustor exit temperature distribution in steady state reacting flow simulation,the turbulent Schmidt number in steady state simulation should be modified to cater to dilution air mixing process.