COMPUTATION OF TURBULENCE FLOW ON HYBRID GRIDS USING k-ω TURBULENCE MODEL AND OSHER SCHEME

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.

A sub-grid scale model for Burgers turbulence based on the artificial neural network method

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.

Artificial neural network-based one-equation model for simulation of laminar-turbulent transitional flow

A mapping function between the Reynolds-averaged Navier-Stokes mean flow variables and transition intermittency factor is constructed by fully connected artificial neural network(ANN),which replaces the governing equation of the intermittency factor in transition-predictive Spalart-Allmaras(SA)-γmodel.By taking SA-γmodel as the benchmark,the present ANN model is trained at two airfoils with various angles of attack,Mach numbers and Reynolds numbers,and tested with unseen airfoils in different flow states.The a posteriori tests manifest that the mean pressure coefficient,skin friction coefficient,size of laminar separation bubble,mean streamwise velocity,Reynolds shear stress and lift/drag/moment coefficient from the present two-way coupling ANN model almost coincide with those from the benchmark SA-γmodel.Furthermore,the ANN model proves to exhibit a higher calculation efficiency and better convergence quality than traditional SA-γmodel.

Stochastic Analysis and Modeling of Velocity Observations in Turbulent Flows

Highly turbulent water flows,often encountered near human constructions like bridge piers,spillways,and weirs,display intricate dynamics characterized by the formation of eddies and vortices.These formations,varying in sizes and lifespans,significantly influence the distribution of fluid velocities within the flow.Subsequently,the rapid velocity fluctuations in highly turbulent flows lead to elevated shear and normal stress levels.For this reason,to meticulously study these dynamics,more often than not,physical modeling is employed for studying the impact of turbulent flows on the stability and longevity of nearby structures.Despite the effectiveness of physical modeling,various monitoring challenges arise,including flow disruption,the necessity for concurrent gauging at multiple locations,and the duration of measurements.Addressing these challenges,image velocimetry emerges as an ideal method in fluid mechanics,particularly for studying turbulent flows.To account for measurement duration,a probabilistic approach utilizing a probability density function(PDF)is suggested to mitigate uncertainty in estimated average and maximum values.However,it becomes evident that deriving the PDF is not straightforward for all turbulence-induced stresses.In response,this study proposes a novel approach by combining image velocimetry with a stochastic model to provide a generic yet accurate description of flow dynamics in such applications.This integration enables an approach based on the probability of failure,facilitating a more comprehensive analysis of turbulent flows.Such an approach is essential for estimating both short-and long-term stresses on hydraulic constructions under assessment.

Artificial neural network-based subgrid-scale models for LES of compressible turbulent channel flow

Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.

Numerical simulation of gas–liquid flow in the bubble column using Wray–Agarwal turbulence model coupled with population balance model

In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..

Simulation Analysis of New Energy Vehicle Engine Cooling System Based on K-E Turbulent Flow Mathematical Model

New energy vehicles have better clean and environmental protection characteristics than traditional fuel vehicles.The new energy engine cooling technology is critical in the design of new energy vehicles.This paper used oneand three-way joint simulation methods to simulate the refrigeration system of new energy vehicles.Firstly,a k-εturbulent flow model for the cooling pump flow field is established based on the principle of computational fluid dynamics.Then,the CFD commercial fluid analysis software FLUENT is used to simulate the flow field of the cooling pump under different inlet flow conditions.This paper proposes an optimization scheme for new energy vehicle engines’“boiling”phenomenon under high temperatures and long-time climbing conditions.The simulation results show that changing the radiator’s structure and adjusting the thermostat’s parameters can solve the problem of a“boiling pot.”The optimized new energy vehicle engine can maintain a better operating temperature range.The algorithm model can reference each cryogenic system component hardware selection and control strategy in the new energy vehicle’s engine.

Numerical prediction of inner turbulent flow in conical diffuser by using a new five-point scheme and DLR k-ε turbulence model

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.

Numerical Evaluation of Two k-εTurbulence Model for Predicting Flow and Solidification in Continuous Casting Slab

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.

REVISED k-ε TURBULENCE MODEL IN ELECTROMAGNETIC CONTINUOUS CASTING OF MELT

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 K-εTWO-EQUATION TURBULENCE MODEL FOR THE SOLID-LIQUID TWO-PHASE FLOWS

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.

An assessment of k-εturbulence models for gas distribution analysis

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.

Numerical analysis of turbulent mixed convection air flow in inclined plane channel with k-εtype turbulence 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.

Transition Modelling with the k-Ω Turbulence Model and an Intermittency Transport Equation

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.

Modification of k-ω Turbulence Model for Predicting Airfoil Aerodynamic Performance

Predicting wind turbine S825 airfoil＇s aerodynamic performance is crucial to improving its energy efficiency and reducing its environmental impact. In this paper, a numerical simulation on the wind turbine S825 airfoil is con- ducted with k-to turbulence model at different attack angles. By comparing with experimental data, a new method of modifying k-to model is proposed. A modifying function is proposed to limit the production term in ω equation based on fluid rotation and deformation. This method improves turbulent viscosity and decreases separating re- gion when the airfoil works at large separating conditions. The predictive accuracy could be improved by using the modified k-to turbulence model.

COMPARISON OF K-ε AND K-ω TURBULENCE MODELS FOR SIMULATION OF SHALLOW RECIRCULATING FLOWS IN AN OPEN CHANNEL

The paper describes the comparison of commonly used k-ε and k-ω turbulence models for simulating the recirculating flows behind a sudden expansion in a straight open channel. The depth-average numerical model and the quasi three-dimensional multilayer model are presented respectively. The governing equations are split into three parts in the finite difference approach: advection, dispersion and propagation. Comparison of the results computed with the k-ε and k-ω models shows that the two turbulence models have strong robustness in simulating the reattachment length of the recirculating flows. The computed results are found in good agreement with the measured data.

An Improved k-Equation Turbulence Model

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.

Effect of RANS Turbulence Model on Aerodynamic Behavior of Trains in Crosswind

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.

New perspective in statistical modeling of wall-bounded turbulence

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.

Effect of Time Step Size and Turbulence Model on the Open Water Hydrodynamic Performance Prediction of Contra-Rotating Propellers

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.