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.

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.

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.

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.

A USM-Θ two-phase turbulence model for simulating dense gas-particle flows

A second-order moment two-phase turbulence model for simulating dense gas-particle flows （USM-Θ model）, combining the unified second-order moment twophase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM-Θ model is also better than the k-ε-kp-Θ model and the k-ε-kp-εp-Θ model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.

Modeling of Fluid Turbulence Modification Using Two-time-scale Dissipation Models and Accounting for the Particle Wake Effect

Presently developed two-phase turbulence models under-predict the gas turbulent fluctuation, because their turbulence modification models cannot fully reflect the effect of particles. In this paper, a two-time-scale dis- sipation model of turbulence modification, developed for the two-phase velocity correlation and for the dissipation rate of gas turbulent kinetic energy, is proposed and used to simulate sudden-expansion and swirling gas-particle flows. The proposed two-time scale model gives better results than the single-time scale model. Besides, a gas tur- bulence augmentation model accounting for the finite-size particle wake effect in the gas Reynolds stress equation is proposed. The proposed turbulence modification models are used to simulate two-phase pipe flows. It can prop- erly predict both turbulence reduction and turbulence enhancement for a certain size of particles observed in ex- periments.

Applying resolved-scale linearly forced isotropic turbulence in rational subgrid-scale modeling

In previous attempts of rational subgrid-scale (SGS) modeling by employing the Kolmogorov equation of filtered (KEF) quantities, it was necessary to assume that the resolved-scale second-order structure function is stationary. Forced isotropic turbulence is often used as a framework for establishing and validating such SGS models based on stationary restrictions, for it generates statistical stationary samples. However, traditional forcing method at low wavenumbers cannot provide an analytic form of forcing term for a complete KEF in physical space, which has been illustrated to be essential in the modeling of such SGS models. Thus, an alternative forcing method giving an analytic forcing term in physical space is needed for rational SGS modeling. Giving an analytic linear driving term in physical space, linearly forced isotropic turbulence should be considered an ideal theoretical framework for rational SGS modeling. In this paper, we demonstrate the feasibility of establishing a rational SGS model with stationary restriction based on linearly forced isotropic turbulence. The performance of this rational SGS model is validated. We, therefore, propose the use of linearly forced isotropic turbulence as a complement to free-decaying isotropic turbulence and low-wavenumber forced isotropic turbulence for SGS model validations.

Mechanism of local scour around submarine pipelines based on numerical simulation of turbulence model

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.

Turbulence and cavitation models for time-dependent turbulent cavitating flows

Cavitation typically occurs when the fluid pressure is lower than the vapor pressure at a local thermodynamic state, and the flow is frequently unsteady and turbulent. To assess the state-of-the-art of computational capabilities for unsteady cavitating flows, different cavitation and turbulence model combinations are conducted. The selected cavitation models include several widely-used models including one based on phenomenological argument and the other utilizing interface dynamics. The k-e turbulence model with additional implementation of the filter function and density correction function are considered to reduce the eddy viscosity according to the computed turbulence length scale and local fluid density respectively. We have also blended these alternative cavitation and lustrate that the eddy viscosity turbulence treatments, to ilnear the closure region can significantly influence the capture of detached cavity. From the experimental validations regarding the force analysis, frequency, and the cavity visualization, no single model combination performs best in all aspects. Furthermore, the implications of parameters contained in different cavitation models are investigated. The phase change process is more pronounced around the detached cavity, which is better illustrated by the interfacial dynamics model. Our study provides insight to aid further modeling development.

Modeling of Fluid Turbulence Modification Using Two-time-scale Dissipation Models and Accounting for the Particle Wake Effect

Presently developed two-phase turbulence models under-predict the gas turbulent fluctuation, because their turbulence modification models cannot fully reflect the effect of particles. In this paper, a two-time-scale dissipation model of turbulence modification, developed for the two-phase velocity correlation and for the dissipation rate of gas turbulent kinetic energy, is proposed and used to simulate sudden-expansion and swirling gas-particle flows. The proposed two-time scale model gives better results than the single-time scale model. Besides, a gas turbulence augmentation model accounting for the f＇mite-size particle wake effect in the gas Reynolds stress equation is proposed. The proposed turbulence modification models are used to simulate two-phase pipe flows. It can properly predict both turbulence reduction and turbulence enhancement for a certain size of particles observed in experiments.

Influence of moderate-to-strong anisotropic non-Kolmogorov turbulence on intensity fluctuations of a Gaussian–Schell model beam in marine atmosphere

The scintillation index（SI） of a Gaussian–Schell model（GSM） beam in a moderate-to-strong anisotropic nonKolmogorov turbulent atmosphere is developed based on the extended Rytov theory. The on-axis SI in a marine atmosphere is higher than that in a terrestrial atmosphere, but the off-axis SI exhibits the opposite trend. The on-axis SI first increases and then begins to decrease and saturate as the turbulence strength increases. Turbulence inner and outer scales have different effects on the on-axis SI in different turbulent fluctuation regions. The anisotropy characteristic of atmospheric turbulence leads to the decline in the on-axis SI, and the rise in the off-axis SI. The on-axis SI can be lowered by increasing the anisotropy of turbulence, wavelength, and source partial coherence before entering the saturation attenuation region. The developed model may be useful for evaluating ship-to-ship/shore free-space optical communication system performance.

Calculation grid and turbulence model for numerical simulating pressure fluctuations in high-speed train tunnel

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.