We present the approaches to implementing the k-√k L turbulence model within the framework of the high-order discontinuous Galerkin(DG)method.We use the DG discretization to solve the full Reynolds-averaged Navier-St...We present the approaches to implementing the k-√k L turbulence model within the framework of the high-order discontinuous Galerkin(DG)method.We use the DG discretization to solve the full Reynolds-averaged Navier-Stokes equations.In order to enhance the robustness of approaches,some effective techniques are designed.The HWENO(Hermite weighted essentially non-oscillatory)limiting strategy is adopted for stabilizing the turbulence model variable k.Modifications have been made to the model equation itself by using the auxiliary variable that is always positive.The 2nd-order derivatives of velocities required in computing the von Karman length scale are evaluated in a way to maintain the compactness of DG methods.Numerical results demonstrate that the approaches have achieved the desirable accuracy for both steady and unsteady turbulent simulations.展开更多
We introduce a framework for statistical inference of the closure coefficients using machine learning methods.The objective of this framework is to quantify the epistemic uncertainty associated with the closure model ...We introduce a framework for statistical inference of the closure coefficients using machine learning methods.The objective of this framework is to quantify the epistemic uncertainty associated with the closure model by using experimental data via Bayesian statistics.The framework is tailored towards cases for which a limited amount of experimental data is available.It consists of two components.First,by treating all latent variables(non-observed variables)in the model as stochastic variables,all sources of uncertainty of the probabilistic closure model are quantified by a fully Bayesian approach.The probabilistic model is defined to consist of the closure coefficients as parameters and other parameters incorporating noise.Then,the uncertainty associated with the closure coefficients is extracted from the overall uncertainty by considering the noise being zero.The overall uncertainty is rigorously evaluated by using Markov-Chain Monte Carlo sampling assisted by surrogate models.We apply the framework to the Spalart-Allmars one-equation turbulence model.Two test cases are considered,including an industrially relevant full aircraft model at transonic flow conditions,the Airbus XRF1.Eventually,we demonstrate that epistemic uncertainties in the closure coefficients result into uncertainties in flow quantities of interest which are prominent around,and downstream,of the shock occurring over the XRF1 wing.This data-driven approach could help to enhance the predictive capabilities of CFD in terms of reliable turbulence modeling at extremes of the flight envelope if measured data is available,which is important in the context of robust design and towards virtual aircraft certification.The plentiful amount of information about the uncertainties could also assist when it comes to estimating the influence of the measured data on the inferred model coefficients.Finally,the developed framework is flexible and can be applied to different test cases and to various turbulence models.展开更多
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.展开更多
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 variety of turbulence models were used to perform numerical simulations of heat transfer for hydrocarbon fuel flowing upward and downward through uniformly heated vertical pipes at supercritical pressure. Inlet temp...A variety of turbulence models were used to perform numerical simulations of heat transfer for hydrocarbon fuel flowing upward and downward through uniformly heated vertical pipes at supercritical pressure. Inlet temperatures varied from 373 K to 663 K, with heat flux rang- ing from 300 kW/m2 to 550 kW/m2. Comparative analyses between predicted and experimental results were used to evaluate the ability of turbulence models to respond to variable thermophysical properties of hydrocarbon fuel at supercritical pressure. It was found that the prediction performance of turbulence models is mainly determined by the damping function, which enables them to respond differently to local flow conditions. Although prediction accuracy for experimental results varied from condition to condition, the shear stress transport （SST） and launder and sharma models performed better than all other models used in the study. For very small buoyancy-influenced runs, the thermal-induced acceleration due to variations in density lead to the impairment of heat transfer occurring in the vicinity of pseudo-critical points, and heat transfer was enhanced at higher temperatures through the combined action of four thermophysical properties： density, viscosity, thermal conductivity and specific heat. For very large buoyancy- influenced runs, the thermal-induced acceleration effect was over predicted by the LS and AB models.展开更多
The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order differ...The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order difference schemes. The compressibility effect of density gradient, pressure dilatation and turbulent Mach number is accounted. In order to reduce confusions between model uncertainties and discretization errors, the formally fifth-order explicit weighted compact nonlinear scheme （WCNS-E-5） is adopted for convection terms, and a fourth-order staggered central difference scheme is applied for viscous terms. The 15° and 34° compression comers at Mach number 9.22 are investigated. Numerical results show that the original SST model is superior to the original S-A model in the resolution of separated regions and predictions of wall pressures and wall heat-flux rates. The capability of the S-A model can be largely improved by blending Catris＇ and Shur＇s compressibility corrections. Among the three corrections of the SST model listed in the present paper, Catris＇ modification brings the best results. However, the dissipation and pressure dilatation corrections result in much larger separated regions than that of the experiment, and are much worse than the original SST model as well as the other two corrections. The correction of turbulent Mach number makes the separated region slightly smaller than that of the original SST model. Some results of low-order schemes are also presented. When compared to the results of the high-order schemes, the separated regions are smaller, and the peak wall pressures and peak heat-flux rates are lower in the region of the reattachment points.展开更多
A new and much simpler rotation and curvature effects factor, which takes the form of Richardson number suggested by Hellsten originally for SST k-ω model, is presented for Spalart and Shur＇s rotation and curvature ...A new and much simpler rotation and curvature effects factor, which takes the form of Richardson number suggested by Hellsten originally for SST k-ω model, is presented for Spalart and Shur＇s rotation and curvature correction in the context of Spalart-Allmaras （SA） turbulence model. The new factor excludes the Lagrangian derivative of the strain rate tensor that exists in the SARC model, resulting in a simple, efficient and easy-to-implement approach for SA turbulence model （denoted as SARCM） to account for the effects of system rotation and curvature, techniquely. And then the SARCM is tested through turbulent curved wall flows： one is the flow over a zeropressure-gradient curved wall and the other is the channel flow in a duct with a U-turn. Predictions of the SARCM model are compared with experimental data and with the results obtained using original SA and SARC models. The numerical results show that SARCM can predict the rotation-curvature effects as accurately as SARC, but considerably more efficiently. Additionally, the accuracy of SARCM might strongly depend on the rotation-curvature model constants. Suggesting values for those constants are given, after some trials and errors.展开更多
Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the ...Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the waveforms and peak-peak values of pressure fluctuations between numerical simulation and moving model test,the structured grid and the SST k-ωturbulence model are selected for numerical simulating the process of high-speed train passing through the tunnel.The largest value of pressure wave amplitudes of numerical simulation and moving model test meet each other.And the locations of the largest value of the initial compression and expansion wave amplitude of numerical simulation are in agreement with that of moving model test.The calculated pressure at the measurement point fully conforms to the propagation law of compression and expansion waves in the tunnel.展开更多
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.展开更多
The Reynolds Averaged Navier-Stokes(RANS) models are still the workhorse in current engineering applications due to its high efficiency and robustness. However, the closure coefficients of RANS turbulence models are d...The Reynolds Averaged Navier-Stokes(RANS) models are still the workhorse in current engineering applications due to its high efficiency and robustness. However, the closure coefficients of RANS turbulence models are determined by model builders according to some simple fundamental flows, and the suggested values may not be applicable to complex flows, especially supersonic jet interaction flow. In this work, the Bayesian method is employed to recalibrate the closure coefficients of Spalart-Allmaras(SA) turbulence model to improve its performance in supersonic jet interaction problem and quantify the uncertainty of wall pressure and separation length. The embedded model error approach is applied to the Bayesian uncertainty analysis. Firstly, the total Sobol index is calculated by non-intrusive polynomial chaos method to represent the sensitivity of wall pressure and separation length to model parameters. Then, the pressure data and the separation length are respectively served as calibration data to get the posterior uncertainty of model parameters and Quantities of Interests(Qo Is). The results show that the relative error of the wall pressure predicted by the SA turbulence model can be reduced from 14.99% to 2.95% through effective Bayesian parameter estimation. Besides, the calibration effects of four likelihood functions are systematically evaluated. The posterior uncertainties of wall pressure and separation length estimated by different likelihood functions are significantly discrepant, and the Maximum a Posteriori(MAP) values of parameters inferred by all functions show better performance than the nominal values. Finally, the closure coefficients are also estimated at different jet total pressures. The similar posterior distributions of model parameters are obtained in different cases, and the MAP values of parameters calibrated in one case are also applicable to other cases.展开更多
Entropy represents the dissipation rate of energy. Through direct numerical simulation （DNS） of supersonic compression ramp flow, we find the value of entropy is monotonously decreasing along the wall-normal directi...Entropy represents the dissipation rate of energy. Through direct numerical simulation （DNS） of supersonic compression ramp flow, we find the value of entropy is monotonously decreasing along the wall-normal direction no matter in the attached or the separated region. Based on this feature, a new version of Baldwin-Lomax turbulence model （BL-entropy） is proposed in this paper. The supersonic compression ramp and cavity-ramp flows in which the original Baldwin-Lomax model fails to get convergent solutions are chosen to evaluate the performance of this model. Results from one-equation Spalart-Allmaras model （SA） and two-equation Wilcox k-x model are also included to compare with available experimental and DNS data. It is shown that BLentropy could conquer the essential deficiency of the original version by providing a more physically meaningful length scale in the complex flows. Moreover, this method is simple, computationally efficient and general, making it applicable to other models related with the supersonic boundary layer.展开更多
The flow field of liquid phase (water) of agitated extraction columns is simulated with the help of computational fluid dynamics (CFD). Four kinds of Reynolds-averaged turbulence models, i.e. the standard k-ε model, ...The flow field of liquid phase (water) of agitated extraction columns is simulated with the help of computational fluid dynamics (CFD). Four kinds of Reynolds-averaged turbulence models, i.e. the standard k-ε model, the RNG (renormalization group) k-s model, the realizable k-ε model and the Reynolds stress model, are compared in detail in order to judge which is the best model in terms of the accuracy, less CPU time and memory required. The performance of the realizable k-s model is obviously improved by reducing the model constant from C2 = 1.90 to C2 = 1.61. It is concluded that the improved realizable k-e model is the optimal model.展开更多
From the controlling equations of atmosphere motion, Prandtl's mixing length theory is used to derive the atmospheric turbulence models, such as Burgers equation model and Burgers-KdV equation model. And then the ...From the controlling equations of atmosphere motion, Prandtl's mixing length theory is used to derive the atmospheric turbulence models, such as Burgers equation model and Burgers-KdV equation model. And then the projective Riccati equations are applied to solve these atmospheric turbulence models, where much more patterns are obtained, including solitary wave pattern, singular pattern, and so on.展开更多
Data-driven turbulence modeling studies have reached such a stage that the basic framework is settled,but several essential issues remain that strongly affect the performance.Two problems are studied in the current re...Data-driven turbulence modeling studies have reached such a stage that the basic framework is settled,but several essential issues remain that strongly affect the performance.Two problems are studied in the current research:(1)the processing of the Reynolds stress tensor and(2)the coupling method between the machine learning model and flow solver.For the Reynolds stress processing issue,we perform the theoretical derivation to extend the relevant tensor arguments of Reynolds stress.Then,the tensor representation theorem is employed to give the complete irreducible invariants and integrity basis.An adaptive regularization term is employed to enhance the representation performance.For the coupling issue,an iterative coupling framework with consistent convergence is proposed and then applied to a canonical separated flow.The results have high consistency with the direct numerical simulation true values,which proves the validity of the current approach.展开更多
The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid ada...The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid adaptation techniques. The investigation focuses on two configurations, related to hypersonic and subsonic experiments. The applicability tests of different turbu- lence models are conducted on the level of two-equation models calculating the steady state solution of the Reynolds-averaged Navier-Stokes（RANS） equations. All used models, the original Wilcox k-co, the Menter shear-stress transport （SST） and the ex- plicit algebraic Reynolds stress model（EARSM） formulation, predict an asymmetric base flow in both cases caused by the support of the models. A comparison with preliminary experimental results indicates a preference for the SST and EARSM results over the results from the older k-co model. Sensitivity studies show no significant influence of the grid topology or the location of the laminar to turbulent transition on the base flow field, but a strong influence of even small angles of attack is reported from the related experiments.展开更多
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.展开更多
Five turbulence models of Reynolds average Navier-Stokes(RANS),including the standard k-ω model,the RNG k-e model taking into account the low Reynolds number effect,the realizable k-ω model,the SST k-ω model,and th...Five turbulence models of Reynolds average Navier-Stokes(RANS),including the standard k-ω model,the RNG k-e model taking into account the low Reynolds number effect,the realizable k-ω model,the SST k-ω model,and the Reynolds stress model(RSM),are employed in the numerical simulations of direct current(DC)arc plasma torches in the range of arc current from 80 A to 240 A and air gas flow rate from 10 m^3 h^-1 to 50 m^3 h^-1.The calculated voltage,electric field intensity,and the heat loss in the arc chamber are compared with the experiments.The results indicate that the arc voltage,the electric field,and the heat loss in the arc chamber calculated by using the standard k-ω model,the RNG k-ωmodel taking into account the low Reynolds number effect,and the realizable k-ω model are much larger than those in the experiments.The RSM predicts relatively close results to the experiments,but fails in the trend of heat loss varying with the gas flow rate.The calculated results of the SST k-ω model are in the best agreement with the experiments,which may be attributed to the reasonable predictions of the turbulence as well as its distribution.展开更多
With the rapid development of artificial intelligence techniques such as neural networks,data-driven machine learning methods are popular in improving and constructing turbulence models.For high Reynolds number turbul...With the rapid development of artificial intelligence techniques such as neural networks,data-driven machine learning methods are popular in improving and constructing turbulence models.For high Reynolds number turbulence in aerodynamics,our previous work built a data-driven model applicable to subsonic airfoil flows with different free stream conditions.The results calculated by the proposed model are encouraging.In this work,we aim to model the turbulence of transonic wing flows with fully connected deep neural networks,where there is less research at present.The proposed model is driven by two flow cases of the ONERA(Office National d'Etudes et de Recherches Aerospatiales)wing and coupled with the Navier-Stokes equation solver.Four subcritical and transonic benchmark cases of different wings are used to evaluate the model performance.The iteration process is stable,and final convergence is achieved.The proposed model can be used to surrogate the traditional Reynolds averaged Navier-Stokes turbulence model.Compared with the data calculated by the Spallart-Allmaras model,the results show that the proposed model can be well generalized to the test cases.The mean relative error of the drag coefficient at different sections is below 4%for each case.This work demonstrates that modeling turbulence by data-driven methods is feasible and that our modeling pattern is effective.展开更多
A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion t...A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion terms were calculated using the second-order AUSM(advection upstream splitting method) scheme and the second-order central difference scheme,respectively.The Runge-Kutta time marching method was employed to solve the governing equations for steady state solutions.Significant flow separation-region which indicates highly non-isotropic turbulence structure has been found in the present work due to intensity interaction under the 24° compression corner.Comparisons between the calculated results and experimental data have been carried out,including surface pressure distribution,boundary-layer static pressure profiles and mean velocity profiles.The numerical results agree well with the experimental values,which indicate Gao-Yong compressible turbulence model is suitable for the prediction of shock wave turbulent boundary layer interaction in two-dimensional compression corner flows.展开更多
The Rotation and Curvature(RC)correction is an important turbulence model modifi-cation approach,and the Spalart-Allmaras model with the RC correction(SA-RC)has been exten-sively studied and used.As a multiplier of th...The Rotation and Curvature(RC)correction is an important turbulence model modifi-cation approach,and the Spalart-Allmaras model with the RC correction(SA-RC)has been exten-sively studied and used.As a multiplier of the modelling equation’s production term,the rotation function f_(r1)should have a cautiously designed value range,but its limit varies in different models and flow solvers.Therefore,the need of restriction is discussed theoretically,and the common range of f_(r1)is explored in Burgers vortexes.Afterwards,the SA-RC model with different limits is tested numerically.Negative f_(r1)always appears in the SA-RC model,and the difference between simula-tion results brought by the limits is not negligible.A lower limit of 0 enhances turbulence produc-tion,and therefore the vortex structures are dissipated faster and shrink in size,while an upper limit plays an opposite role.Considering that the lower limit of 0 usually promotes the simulation accu-racy and fixes the numerical defect,whereas the upper limit worsens the predictive performance in most cases,it is recommended to limit f_(r1)non-negative while utilizing the SA-RC model.In addi-tion,the RC-corrected model has a better prediction of the attached flow near curved walls,while the SA-Helicity model largely improves the simulation accuracy of three-dimensional large-scale vortices.The model combining both corrections has the potential to become more adaptive and more accurate.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.92252201 and 11721202)the Fundamental Research Funds for the Central Universities.
文摘We present the approaches to implementing the k-√k L turbulence model within the framework of the high-order discontinuous Galerkin(DG)method.We use the DG discretization to solve the full Reynolds-averaged Navier-Stokes equations.In order to enhance the robustness of approaches,some effective techniques are designed.The HWENO(Hermite weighted essentially non-oscillatory)limiting strategy is adopted for stabilizing the turbulence model variable k.Modifications have been made to the model equation itself by using the auxiliary variable that is always positive.The 2nd-order derivatives of velocities required in computing the von Karman length scale are evaluated in a way to maintain the compactness of DG methods.Numerical results demonstrate that the approaches have achieved the desirable accuracy for both steady and unsteady turbulent simulations.
文摘We introduce a framework for statistical inference of the closure coefficients using machine learning methods.The objective of this framework is to quantify the epistemic uncertainty associated with the closure model by using experimental data via Bayesian statistics.The framework is tailored towards cases for which a limited amount of experimental data is available.It consists of two components.First,by treating all latent variables(non-observed variables)in the model as stochastic variables,all sources of uncertainty of the probabilistic closure model are quantified by a fully Bayesian approach.The probabilistic model is defined to consist of the closure coefficients as parameters and other parameters incorporating noise.Then,the uncertainty associated with the closure coefficients is extracted from the overall uncertainty by considering the noise being zero.The overall uncertainty is rigorously evaluated by using Markov-Chain Monte Carlo sampling assisted by surrogate models.We apply the framework to the Spalart-Allmars one-equation turbulence model.Two test cases are considered,including an industrially relevant full aircraft model at transonic flow conditions,the Airbus XRF1.Eventually,we demonstrate that epistemic uncertainties in the closure coefficients result into uncertainties in flow quantities of interest which are prominent around,and downstream,of the shock occurring over the XRF1 wing.This data-driven approach could help to enhance the predictive capabilities of CFD in terms of reliable turbulence modeling at extremes of the flight envelope if measured data is available,which is important in the context of robust design and towards virtual aircraft certification.The plentiful amount of information about the uncertainties could also assist when it comes to estimating the influence of the measured data on the inferred model coefficients.Finally,the developed framework is flexible and can be applied to different test cases and to various turbulence models.
基金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(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.
基金funding support from National Natural Science Foundation of China (No.51406005)Defense Industrial Technology Development Program of China (No.B2120132006)
文摘A variety of turbulence models were used to perform numerical simulations of heat transfer for hydrocarbon fuel flowing upward and downward through uniformly heated vertical pipes at supercritical pressure. Inlet temperatures varied from 373 K to 663 K, with heat flux rang- ing from 300 kW/m2 to 550 kW/m2. Comparative analyses between predicted and experimental results were used to evaluate the ability of turbulence models to respond to variable thermophysical properties of hydrocarbon fuel at supercritical pressure. It was found that the prediction performance of turbulence models is mainly determined by the damping function, which enables them to respond differently to local flow conditions. Although prediction accuracy for experimental results varied from condition to condition, the shear stress transport （SST） and launder and sharma models performed better than all other models used in the study. For very small buoyancy-influenced runs, the thermal-induced acceleration due to variations in density lead to the impairment of heat transfer occurring in the vicinity of pseudo-critical points, and heat transfer was enhanced at higher temperatures through the combined action of four thermophysical properties： density, viscosity, thermal conductivity and specific heat. For very large buoyancy- influenced runs, the thermal-induced acceleration effect was over predicted by the LS and AB models.
基金Foundation items： National Basic Research Program of China （2009CB723801） National Natural Science Foundation of China （11072259）
文摘The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order difference schemes. The compressibility effect of density gradient, pressure dilatation and turbulent Mach number is accounted. In order to reduce confusions between model uncertainties and discretization errors, the formally fifth-order explicit weighted compact nonlinear scheme （WCNS-E-5） is adopted for convection terms, and a fourth-order staggered central difference scheme is applied for viscous terms. The 15° and 34° compression comers at Mach number 9.22 are investigated. Numerical results show that the original SST model is superior to the original S-A model in the resolution of separated regions and predictions of wall pressures and wall heat-flux rates. The capability of the S-A model can be largely improved by blending Catris＇ and Shur＇s compressibility corrections. Among the three corrections of the SST model listed in the present paper, Catris＇ modification brings the best results. However, the dissipation and pressure dilatation corrections result in much larger separated regions than that of the experiment, and are much worse than the original SST model as well as the other two corrections. The correction of turbulent Mach number makes the separated region slightly smaller than that of the original SST model. Some results of low-order schemes are also presented. When compared to the results of the high-order schemes, the separated regions are smaller, and the peak wall pressures and peak heat-flux rates are lower in the region of the reattachment points.
基金supported by National Natural Science Foundation of China (No. 11002117)
文摘A new and much simpler rotation and curvature effects factor, which takes the form of Richardson number suggested by Hellsten originally for SST k-ω model, is presented for Spalart and Shur＇s rotation and curvature correction in the context of Spalart-Allmaras （SA） turbulence model. The new factor excludes the Lagrangian derivative of the strain rate tensor that exists in the SARC model, resulting in a simple, efficient and easy-to-implement approach for SA turbulence model （denoted as SARCM） to account for the effects of system rotation and curvature, techniquely. And then the SARCM is tested through turbulent curved wall flows： one is the flow over a zeropressure-gradient curved wall and the other is the channel flow in a duct with a U-turn. Predictions of the SARCM model are compared with experimental data and with the results obtained using original SA and SARC models. The numerical results show that SARCM can predict the rotation-curvature effects as accurately as SARC, but considerably more efficiently. Additionally, the accuracy of SARCM might strongly depend on the rotation-curvature model constants. Suggesting values for those constants are given, after some trials and errors.
文摘Calculation grid and turbulence model for numerical simulating pressure fluctuations in a high-speed train tunnel are studied through the comparison analysis of numerical simulation and moving model test.Compared the waveforms and peak-peak values of pressure fluctuations between numerical simulation and moving model test,the structured grid and the SST k-ωturbulence model are selected for numerical simulating the process of high-speed train passing through the tunnel.The largest value of pressure wave amplitudes of numerical simulation and moving model test meet each other.And the locations of the largest value of the initial compression and expansion wave amplitude of numerical simulation are in agreement with that of moving model test.The calculated pressure at the measurement point fully conforms to the propagation law of compression and expansion waves in the tunnel.
基金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.
基金supported by the National Numerical Windtunnel Project,China(No.NNW2019ZT1-A03)the National Natural Science Foundation of China(No.11721202)。
文摘The Reynolds Averaged Navier-Stokes(RANS) models are still the workhorse in current engineering applications due to its high efficiency and robustness. However, the closure coefficients of RANS turbulence models are determined by model builders according to some simple fundamental flows, and the suggested values may not be applicable to complex flows, especially supersonic jet interaction flow. In this work, the Bayesian method is employed to recalibrate the closure coefficients of Spalart-Allmaras(SA) turbulence model to improve its performance in supersonic jet interaction problem and quantify the uncertainty of wall pressure and separation length. The embedded model error approach is applied to the Bayesian uncertainty analysis. Firstly, the total Sobol index is calculated by non-intrusive polynomial chaos method to represent the sensitivity of wall pressure and separation length to model parameters. Then, the pressure data and the separation length are respectively served as calibration data to get the posterior uncertainty of model parameters and Quantities of Interests(Qo Is). The results show that the relative error of the wall pressure predicted by the SA turbulence model can be reduced from 14.99% to 2.95% through effective Bayesian parameter estimation. Besides, the calibration effects of four likelihood functions are systematically evaluated. The posterior uncertainties of wall pressure and separation length estimated by different likelihood functions are significantly discrepant, and the Maximum a Posteriori(MAP) values of parameters inferred by all functions show better performance than the nominal values. Finally, the closure coefficients are also estimated at different jet total pressures. The similar posterior distributions of model parameters are obtained in different cases, and the MAP values of parameters calibrated in one case are also applicable to other cases.
基金National Basic Research Program of China(No.2009CB724104)the Innovation Foundation of BUAA for Ph.D.Graduates(No.300521)
文摘Entropy represents the dissipation rate of energy. Through direct numerical simulation （DNS） of supersonic compression ramp flow, we find the value of entropy is monotonously decreasing along the wall-normal direction no matter in the attached or the separated region. Based on this feature, a new version of Baldwin-Lomax turbulence model （BL-entropy） is proposed in this paper. The supersonic compression ramp and cavity-ramp flows in which the original Baldwin-Lomax model fails to get convergent solutions are chosen to evaluate the performance of this model. Results from one-equation Spalart-Allmaras model （SA） and two-equation Wilcox k-x model are also included to compare with available experimental and DNS data. It is shown that BLentropy could conquer the essential deficiency of the original version by providing a more physically meaningful length scale in the complex flows. Moreover, this method is simple, computationally efficient and general, making it applicable to other models related with the supersonic boundary layer.
文摘The flow field of liquid phase (water) of agitated extraction columns is simulated with the help of computational fluid dynamics (CFD). Four kinds of Reynolds-averaged turbulence models, i.e. the standard k-ε model, the RNG (renormalization group) k-s model, the realizable k-ε model and the Reynolds stress model, are compared in detail in order to judge which is the best model in terms of the accuracy, less CPU time and memory required. The performance of the realizable k-s model is obviously improved by reducing the model constant from C2 = 1.90 to C2 = 1.61. It is concluded that the improved realizable k-e model is the optimal model.
文摘From the controlling equations of atmosphere motion, Prandtl's mixing length theory is used to derive the atmospheric turbulence models, such as Burgers equation model and Burgers-KdV equation model. And then the projective Riccati equations are applied to solve these atmospheric turbulence models, where much more patterns are obtained, including solitary wave pattern, singular pattern, and so on.
基金This work was supported by the National Natural Science Foundation of China(91852108,11872230 and 92152301).
文摘Data-driven turbulence modeling studies have reached such a stage that the basic framework is settled,but several essential issues remain that strongly affect the performance.Two problems are studied in the current research:(1)the processing of the Reynolds stress tensor and(2)the coupling method between the machine learning model and flow solver.For the Reynolds stress processing issue,we perform the theoretical derivation to extend the relevant tensor arguments of Reynolds stress.Then,the tensor representation theorem is employed to give the complete irreducible invariants and integrity basis.An adaptive regularization term is employed to enhance the representation performance.For the coupling issue,an iterative coupling framework with consistent convergence is proposed and then applied to a canonical separated flow.The results have high consistency with the direct numerical simulation true values,which proves the validity of the current approach.
基金German Research Foundation (Deutsche Forschungs-gemeinschaft-DFG) Sonderforschungsbereich Transregio 40
文摘The flows behind the base of a generic rocket, at both hypersonic and subsonic flow conditions, are numerically studied. The main concerns are addressed to the evaluation of turbulence models and the using of grid adaptation techniques. The investigation focuses on two configurations, related to hypersonic and subsonic experiments. The applicability tests of different turbu- lence models are conducted on the level of two-equation models calculating the steady state solution of the Reynolds-averaged Navier-Stokes（RANS） equations. All used models, the original Wilcox k-co, the Menter shear-stress transport （SST） and the ex- plicit algebraic Reynolds stress model（EARSM） formulation, predict an asymmetric base flow in both cases caused by the support of the models. A comparison with preliminary experimental results indicates a preference for the SST and EARSM results over the results from the older k-co model. Sensitivity studies show no significant influence of the grid topology or the location of the laminar to turbulent transition on the base flow field, but a strong influence of even small angles of attack is reported from the related experiments.
文摘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.
基金National Natural Science Foundation of China(Nos.11675177,11875256)the Anhui Province Scientific and Technological Project(No.1604a0902145).
文摘Five turbulence models of Reynolds average Navier-Stokes(RANS),including the standard k-ω model,the RNG k-e model taking into account the low Reynolds number effect,the realizable k-ω model,the SST k-ω model,and the Reynolds stress model(RSM),are employed in the numerical simulations of direct current(DC)arc plasma torches in the range of arc current from 80 A to 240 A and air gas flow rate from 10 m^3 h^-1 to 50 m^3 h^-1.The calculated voltage,electric field intensity,and the heat loss in the arc chamber are compared with the experiments.The results indicate that the arc voltage,the electric field,and the heat loss in the arc chamber calculated by using the standard k-ω model,the RNG k-ωmodel taking into account the low Reynolds number effect,and the realizable k-ω model are much larger than those in the experiments.The RSM predicts relatively close results to the experiments,but fails in the trend of heat loss varying with the gas flow rate.The calculated results of the SST k-ω model are in the best agreement with the experiments,which may be attributed to the reasonable predictions of the turbulence as well as its distribution.
基金supported by the National Natural Science Foundation of China(Grant Nos.92152301,and 91852115)the National Numerical Wind tunnel Project(Grand No.NNW2018-ZT1B01).
文摘With the rapid development of artificial intelligence techniques such as neural networks,data-driven machine learning methods are popular in improving and constructing turbulence models.For high Reynolds number turbulence in aerodynamics,our previous work built a data-driven model applicable to subsonic airfoil flows with different free stream conditions.The results calculated by the proposed model are encouraging.In this work,we aim to model the turbulence of transonic wing flows with fully connected deep neural networks,where there is less research at present.The proposed model is driven by two flow cases of the ONERA(Office National d'Etudes et de Recherches Aerospatiales)wing and coupled with the Navier-Stokes equation solver.Four subcritical and transonic benchmark cases of different wings are used to evaluate the model performance.The iteration process is stable,and final convergence is achieved.The proposed model can be used to surrogate the traditional Reynolds averaged Navier-Stokes turbulence model.Compared with the data calculated by the Spallart-Allmaras model,the results show that the proposed model can be well generalized to the test cases.The mean relative error of the drag coefficient at different sections is below 4%for each case.This work demonstrates that modeling turbulence by data-driven methods is feasible and that our modeling pattern is effective.
文摘A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion terms were calculated using the second-order AUSM(advection upstream splitting method) scheme and the second-order central difference scheme,respectively.The Runge-Kutta time marching method was employed to solve the governing equations for steady state solutions.Significant flow separation-region which indicates highly non-isotropic turbulence structure has been found in the present work due to intensity interaction under the 24° compression corner.Comparisons between the calculated results and experimental data have been carried out,including surface pressure distribution,boundary-layer static pressure profiles and mean velocity profiles.The numerical results agree well with the experimental values,which indicate Gao-Yong compressible turbulence model is suitable for the prediction of shock wave turbulent boundary layer interaction in two-dimensional compression corner flows.
基金supported by the National Natural Science Foundation of China(Nos.51976006,51790513)the Aeronautical Science Foundation of China(No.2018ZB51013)+1 种基金the National Science and Technology Major Project,China(2017-II-003-0015)the Open Fund from State Key Laboratory of Aerodynamics,China(No.SKLA2019A0101).
文摘The Rotation and Curvature(RC)correction is an important turbulence model modifi-cation approach,and the Spalart-Allmaras model with the RC correction(SA-RC)has been exten-sively studied and used.As a multiplier of the modelling equation’s production term,the rotation function f_(r1)should have a cautiously designed value range,but its limit varies in different models and flow solvers.Therefore,the need of restriction is discussed theoretically,and the common range of f_(r1)is explored in Burgers vortexes.Afterwards,the SA-RC model with different limits is tested numerically.Negative f_(r1)always appears in the SA-RC model,and the difference between simula-tion results brought by the limits is not negligible.A lower limit of 0 enhances turbulence produc-tion,and therefore the vortex structures are dissipated faster and shrink in size,while an upper limit plays an opposite role.Considering that the lower limit of 0 usually promotes the simulation accu-racy and fixes the numerical defect,whereas the upper limit worsens the predictive performance in most cases,it is recommended to limit f_(r1)non-negative while utilizing the SA-RC model.In addi-tion,the RC-corrected model has a better prediction of the attached flow near curved walls,while the SA-Helicity model largely improves the simulation accuracy of three-dimensional large-scale vortices.The model combining both corrections has the potential to become more adaptive and more accurate.
