Multi-scale Equations for Incompressible Turbulent Flows

The short-range property of interactions between scales in incompressible turbulent flow was examined. Some formulae for the short-range eddy stress were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress were also derived. Multi-scale equations for the incompressible turbulent flows were proposed. Key words turbulence - incompressible flow - interactions between scales - multi-scale equations MSC 2000 76F70

Optimal convergence rates for three-dimensional turbulent flow equations

In this paper, the convergence turbulent flow equations are considered. By rates of solutions to the three-dimensional combining the LP-Lq estimate for the linearized equations and an elaborate energy method, the convergence rates are obtained in various norms for the solution to the equilibrium state in the whole space when the initial perturbation of the equilibrium state is small in the H3-framework. More precisely, the optimal convergence rates of the solutions and their first-order derivatives in the L2-norm are obtained when the LP-norm of the perturbation is bounded for some p ε [1, 6）.

Multi-scale Equations for Compressible Turbulent Flows

The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfer etc. were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress and heat transfer etc. were also given. Multi-scale equations for the compressible turbulent flows were presented. The multi-scale equations are approximately closed and do not contain any empirical constants. The compressibility effects on turbulence are determined by the Farve averaged variables and the nonlinear relationships between the Farve- and physical-averaged variables.

Solution of Self-similar Equations of the k-ε Model in the Shear Turbulent Mixing Problem and Its Numerical Simulation

The paper presents the k-ε model equations of turbulence with a single set of constants chosen by the authors, which is appropriate to simulate a wide range of turbulent flows. The model validation has been performed for a number of flows and its main results are given in the paper. The turbulent mixing of flow with shear in the tangential velocity component is discussed in details. An analytical solution to the system of ordinary differential equations of the k-ε model of turbulent mixing has been found for the self-similar regime of flow. The model coefficients were chosen using simulation results for some simplest turbulent flows. The solution can be used for the verification of codes. The numerical simulation of the problem has been performed by the 2D code EGAK using this model. A good agreement of the numerical simulation results with the self-similar solution, 3D DNS results and known experimental data has been achieved. This allows stating that the k-ε model constants chosen by the authors are acceptable for the considered flow.

Solution of general dynamic equation for nanoparticles in turbulent flow considering fluctuating coagulation

A new averaged general dynamic equation （GDE） for nanoparticles in the turbulent flow is derived by considering the combined effect of convection, Brownian diffusion, turbulent diffusion, turbulent coagulation, and fluctuating coagulation. The equation is solved with the Taylor-series expansion moment method in a turbulent pipe flow. The experiments are performed. The numerical results of particle size distribu- tion correlate well with the experimental data. The results show that, for a turbulent nanoparticulate flow, a fluctuating coagulation term should be included in the averaged particle GDE. The larger the Schmidt number is and the lower the Reynolds number is, the smaller the value of ratio of particle diameter at the outlet to that at the inlet is. At the outlet, the particle number concentration increases from the near-wall region to the near-center region. The larger the Schmidt number is and the higher the Reynolds num- ber is, the larger the difference in particle number concentration between the near-wall region and near-center region is. Particle polydispersity increases from the near-center region to the near-wall region. The particles with a smaller Schmidt number and the flow with a higher Reynolds number show a higher polydispersity. The degree of particle polydispersity is higher considering fluctuating coagulation than that without considering fluctuating coagulation.

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.

Heat Transfer and Flow Characteristics Predictions with a Refined k-ε-f_u Turbulent Model in Impinging Jet

Local heat transfer and flow characteristics in a round turbulent impinging jet for Re≈23 000 is predicted numerically with the RANS approach and a k-ε-fu turbulence model. The heat transfer predictions and turbulence parameters are verified against the axis-symmetric free jet impingement measurements and compared with previous other turbulence models, and results show the k-ε-fu model has a good performance in predictions of the local wall heat transfer coefficient, and in agreement with measurements in mean velocity profiles at different radial positions as well. The numerical model is further used to examine the effect of the fully confined impingement jet on the local Nusselt number. Local Nusselt profiles in x and y-centerlines for the target plate over three separation distances are predicted. Compared with the experimental data, the numerical results are accurate in the central domain around the stagnation region and present a consistent structure distribution.

Rate Equation of Small Particle-Air Bubble Attachment in Turbulent Flow of Flotation Cells

A rate equation of small particle-air bubble attachment in the turbulent now of flotation cells has beenderived. The equation, integrating both the collision probability and adhesion probability together, represents theprobability of attachment between particle and bubble in the turbulent flow. 'Capture efficiency' f(a) is introducedinto the rate equation to reflect the influence of energy hairier on the attachment rate. Three typical situations of particle-bubble interaction in flotation process have been discussed. For a completely hydrophobic particle-bubble system,f(a) = 1. This means that all collision leads to attachment. Whereas for hydrophilic particle-bubble systems, .f(a) =0. Thus no adhesion of particle on bubble occurs at all. In real notation circumstances, however, there always existsa certain energy barrier between the particle and the bubble. Therefore, f(a) = 0～1. In such cases, not all collisionsresult in particle-bubble attachment.

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.

Structural ensemble dynamics based closure model for wall-bounded turbulent flow

Wall-bounded turbulent flow involves the development of multi-scale turbulent eddies, as well as a sharply varying boundary layer. Its theoretical descriptions are yet phenomenological. We present here a new framework called structural ensemble dynamics （SED）, which aims at using systematically all relevant statistical properties of turbulent structures for a quantitative description of ensemble means. A new set of closure equations based on the SED approach for a turbulent channel flow is presented. SED order functions are defined, and numerically determined from data of direct numerical simulations （DNS）. Computational results show that the new closure model reproduces accurately the solution of the original Navier-Stokes simulation, including the mean velocity profile, the kinetic energy of the streamwise velocity component, and every term in the energy budget equation. It is suggested that the SED-based studies of turbulent structure builds a bridge between the studies of physical mechanisms of turbulence and the development of accurate model equations for engineering predictions.

Assessment of k–ε models using tetrahedral grids to describe the turbulent flow field of a PBT impeller and validation through the PIV technique

In turbulence modeling, the RNG and Realizable models have important improvements in the turbulent production and dissipation terms in comparison to the Standard. The selection of the appropriate turbulence model has an impact on the convergence and solution in STRs, and they are used in mixing, multiphase modeling or as starting solution of transient models as DES and LES. Although there are several studies with the pitched blade turbine(PBT) impeller, most of them used the Standard model as representative of all k–ε models, using structured hexahedral grids composed of low number of cells, and in some cases under axial symmetry assumptions.Accordingly, in this work the assessment of the Standard, RNG and Realizable models to describe the turbulent flow field of this impeller, using the Multiple Reference Frame(MRF) and Sliding Mesh(SM) approaches with tetrahedral domains in dense grids, is presented. This kind of cell elements is especially suitable to reproduce complex geometries. Flow velocities and turbulent parameters were verified experimentally by PIV and torque measurements. The three models were capable of predicting fairly the pumping number, the power number based on torque, and velocities. Although the RNG improved the predictions of the turbulent kinetic energy and dissipation rate, the Realizable model presented better performance for both approaches. All models failed in the prediction of the total dissipation rate, and a dependence of its value on the number of cells for the MRF was found.

Three-dimensional turbulent model of heat transfer and fluid flow in GTAW process

A two-equation K-ε turbulent fluid flow model is built to model the heat transfer and fluid flow in gas tungsten arc welding （GTAW） process of stainless steel S US310 and S US316. This model combines the buoyancy force, lorentz force and marangni force as the driving forces of thefluidflow in the weld pool. The material properties are functions of temperature in this model. The simulated results show that the molten metal flowing outward is mainly caused by the marangoni convection, which makes the weld pool become wider and shallower. The comparison of the weld pool shape of SUS310 and SUS316 shows that the slight differences of the value of thermal conductivity mainly attributes to the difference of the weld pool shape and the distinction of heat transport in laminar and turbulent model makes large diversity in the simulated results.

Numerical Studies of Turbulent Flows in Channels with Abrupt Cross-sectional Changes

The turbulent flows through the channels with abrupt cross-sectional changes are common and importantphysical process in nature.For a better prediction of the mean flow and turbulent characteristics for this problem,atwo-dimensional depth-averaged numerical model is developed.The model is robust and accurate in reproducing therecirculation flow behind a groyne and turbulent flows in channels with abrupt cross-sectional changes,when com-pared to the available experimental data of mean velocities and turbulence kinetic energy.Our results reveal that theabrupt cross-sectional change of a channel can affect the flow pattern significantly and introduces the complex turbu-lence characteristics.In particular,when the channel has an abrupt expansion,the mean flow pattern is mainly in lon-gitudinal direction with rather small transverse component.Meanwhile,a recirculating region forms behind the expan-sion position and the turbulence has very strong intensity within this region.For the flow in the channel with an ab-rupt contraction,the longitudinal component of the flow is decreased by the obstruction on one side and accelerated onthe other side,whereas the transverse velocity is small.The turbulence is extraordinarily strong in the regions adja-cent to the contraction wall in the narrow channel.In both cases of abrupt cross-sectional changes,the TKE is genera-ted dominantly by the shear of the longitudinal velocities.

Simulation of Non-Isothermal Turbulent Flows Through Circular Rings of Steel

This article is intended to examine the fluid flow patterns and heat transfer in a rectangular channel embedded with three semi-circular cylinders comprised of steel at the boundaries.Such an organization is used to generate the heat exchangers with tube and shell because of the production of more turbulence due to zigzag path which is in favor of rapid heat transformation.Because of little maintenance,the heat exchanger of such type is extensively used.Here,we generate simulation of flow and heat transfer using nonisothermal flow interface in the Comsol multiphysics 5.4 which executes the Reynolds averaged Navier stokes equation(RANS)model of the turbulent flow together with heat equation.Simulation is tested with Prandtl number(Pr=0.7)with inlet velocity magnitude in the range from 1 to 2 m/sec which generates the Reynolds number in the range of 2.2×10^(5) to 4.4×10^(5) with turbulence kinetic energy and the dissipation rate in ranges(3.75×10^(−3) to 1.5×10^(−2))and(3.73×10^(−3)−3×10^(−2))respectively.Two correlations available in the literature are used in order to check validity.The results are displayed through streamlines,surface plots,contour plots,isothermal lines,and graphs.It is concluded that by retaining such an arrangement a quick distribution of the temperature over the domain can be seen and also the velocity magnitude is increasing from 333.15%to a maximum of 514%.The temperature at the middle shows the consistency in value but declines immediately at the end.This process becomes faster with the decrease in inlet velocity magnitude.

Secondary Instability of Large Scale Structure in Free Turbulent Shear Layer

The secondary instability theory is used to study the behavior of spatially growingdisturbance in free turbulent shear layer.The numerical results indicate that secondaryinstability of subharmonic mode shows a strong choice of spanwise wavenumber andthe maximum growth two dimensional case.In contrast to thatsecondary instabilities of the fundamental mode occur in a wide scope of spanwisewavenumber.We have found so called translative atβ=0 and bifurcationphenomenon for an amplitude of the KH wave larger than 0.06.Dey words instability,large scale structure,bifurcation

Investigation of turbulence and skin friction modification in particle-laden channel flow using lattice Boltzmann method

Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.

NUMERICAL SIMULATION OF WING-BODY JUNCTION TURBULENCE FLOW

Wing-body junction turbulence flow is simulated by using RANS equation and boundary fitted coordinate technique. Three order differential scheme is used in the computation of convection term and two layers turbulence model are employed in the calculation.

2－D SIMULATION OF OPEN CHANNEL TURBULENT FLOW WITH HIGH REYNOLDS NUMBER BY EMPLOYING N－S EQUATIONS AND IEM

Since 1986, we have been studying the topic of simulating open channel turbulent flow, especially the back flow in extended section, by directly employing N-S equations. Now a breakthrough may be obtained. The simulated turbulent flow pattern is realistic and lively, and has a good agreement with physical model test results.

A Pressure-Correction Scheme for Rotational Navier-Stokes Equations and Its Application to Rotating Turbulent Flows

The rotational incremental pressure-correction(RIPC)scheme,described in Timmermans et al.[Int.J.Numer.Methods.Fluids.,22(1996)]and Shen et al.[Math.Comput.,73(2003)]for non-rotational Navier-Stokes equations,is extended to rotating incompressible flows.The method is implemented in the context of a pseudo Fourier-spectral code and applied to several rotating laminar and turbulent flows.The performance of the scheme and the computational results are compared to the socalled diagonalization method(DM)developed by Morinishi et al.[Int.J.Heat.Fluid.Flow.,22(2001)].The RIPC predictions are in excellent agreement with the DM predictions,while being simpler to implement and computationally more efficient.The RIPC scheme is not in anyway limited to implementation in a pseudo-spectral code or periodic boundary conditions,and can be used in complex geometries and with other suitable boundary conditions.

基于新的描述湍流耗散方程的k-ζ两方程湍流模型的数值算法研究

对Navier-Stokes方程进行雷诺平均后出现的各关联项建模,通过新的描述湍流脉动耗散的变量构造耗散方程,建立k-ζ两方程湍流模型,研究了k-ζ两方程湍流模型的数值求解方法。通过求解有限体积法离散的RANS流动控制方程,数值模拟了平板,翼型,机翼等不同湍流流场,并与理论解、实验值及SSTk-ω模型进行比较,全面考察了k-ζ两方程湍流模型在湍流流场计算中的准确性及适用性。数值计算表明,通过建立新的耗散方程研究湍流的方法是可行的,目前的k-ζ两方程湍流模型具有良好的数值稳定性,并且计算结果要优于或者至少与传统的两方程模型精度相当。