Simulation of the thermal method for reducing turbulent friction

The results of calculating a supersonic turbulent boundary layer on a heated surface on the basis of the algebraic two-parameter (k-ε) and four-parameter (k-ε-θ ²-ε ₆) models of turbulence are compared with experimental data. Emphasis is placed on the ability of the models to predict the behavior of the friction and heat-transfer coefficients on a heated surface. The optimal model of turbulence is chosen. The possibility of improving the efficiency of viscous drag reduction by localizing the regions of heat addition to the boundary layer is demonstrated on the basis of numerical calculations. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–68, January–February, 1998. This research was carried out with financial support from the International Scientific and Technological Center (project No. 199).     

An Attempt to Combine Large Eddy Simulation with the k-ε Model in a Channel-Flow Calculation

Large eddy simulation (LES) is combined with the Reynolds-averaged Navier–Stokes (RANS) equation in a turbulent channel-flow calculation. A one-equation subgrid-scale model is solved in a three-dimensional grid in the near-wall region whereas the standard k–ε model is solved in a one-dimensional grid in the outer region away from the wall. The two grid systems are overlapped to connect the two models smoothly. A turbulent channel flow is calculated at Reynolds numbers higher than typical LES and several statistical quantities are examined. The mean velocity profile is in good agreement with the logarithmic law. The profile of the turbulent kinetic energy in the near-wall region is smoothly connected with that of the turbulent energy for the k–ε model in the outer region. Turbulence statistics show that the solution in the near-wall region is as accurate as a usual LES. The present approach is different from wall modeling in LES that uses a RANS model near the wall. The former is not as efficient as the latter for calculating high-Reynolds-number flows. Nevertheless, the present method of combining the two models is expected to pave the way for constructing a unified turbulence model that is useful for many purposes including wall modeling. Received 11 June 1999 and accepted 15 December 2000     

Testing of Model Equations for the Mean Dissipation using Kolmogorov Flows

Direct Numerical Simulations (DNS) of Kolmogorov flows are performed at three different Reynolds numbers Re _λ between 110 and 190 by imposing a mean velocity profile in y-direction of the form U(y) = F sin(y) in a periodic box of volume (2π)³. After a few integral times the turbulent flow turns out to be statistically steady. Profiles of mean quantities are then obtained by averaging over planes at constant y. Based on these profiles two different model equations for the mean dissipation ε in the context of two-equation RANS (Reynolds Averaged Navier–Stokes) modelling of turbulence are compared to each other. The high Reynolds number version of the k-ε-model (Jones and Launder, Int J Heat Mass Transfer 15:301–314, 1972), to be called the standard model and a new model by Menter et al. (2006), to be called the Menter–Egorov model, are tested against the DNS results. Both models are solved numerically and it is found that the standard model does not provide a steady solution for the present case, while the Menter–Egorov model does. In addition a fairly good quantitative agreement of the model solution and the DNS data is found for the averaged profiles of the kinetic energy k and the dissipation ε. Furthermore, an analysis based on flow-inherent geometries, called dissipation elements (Wang and Peters, J Fluid Mech 608:113–138, 2008), is used to examine the Menter–Egorov ε model equation. An expression for the evolution of ε is derived by taking appropriate moments of the equation for the evolution of the probability density function (pdf) of the length of dissipation elements. A term-by-term comparison with the model equation allows a prediction of the constants, which with increasing Reynolds number approach the empirical values.     

A numerical study of non-isothermal turbulent coaxial jets

In this work, we propose to study non isothermal air–air coaxial jets with two different approaches: parabolic and elliptic approaches. The standard k−ε model and the RSM model were applied in this study. The numerical resolution of the equations governing this flow type was carried out for: the parabolic approach, by a “home-made” CFD code based on a finite difference method, and the elliptic approach by an industrial code (FLUENT) based on a finite volume method. In forced convection mode (Fr = ∞), the two turbulence models are valid for the prediction of the mean flow. But for turbulent sizes, k−ε model gives results closer to those achieved in experiments compared to RSM Model. Concerning the limit of validity of the parabolic and elliptic approaches, we showed that for velocities ratio r lower than 1, the results of the two approaches were satisfactory. On the other hand, for r > 1, the difference between the results became increasingly significant. In mixed convection mode (Fr ≅ 20), the results obtained by the two turbulence models for the mean axial velocity were very different even in the plume region. For the temperature and the turbulent sizes the two models give satisfactory results which agree well with the correlations suggested by the experimenters for X ≥ 20. Thus, the second order model with σ _t = 0.85 is more effective for a coaxial jet study in a mixed convection mode.     

Numerical prediction of turbulent flow in a conical diffuser usingk-ε model

Fully developed incompressible turbulent flow in a conical diffuser having a total divergence angle of 8° and an area ratio of 4∶1 has been simulated by ak-ε turbulence model with high Reynolds number and adverse pressure gradient. The research has been done for pipe entry Reynolds numbers of 1.16×10⁵ and 2.93×10⁵. The mean flow velocity and turbulence energy are predicted successfully and the advantage of Boundary Fit Coordinates approach is discussed. Furthermore, thek-ε turbulence model is applied to a flow in a conical diffuser having a total divergence angle of 30° with a perforated screen. A simplified mathematical model, where only the pressure drop is considered, has been used for describing the effect of the perforated screen. The optimum combination of the resistance coefficient and the location of the perforated screen is predicted for high diffuser efficiency or the uniform velocity distribution.     

Prediction of Flow and Heat Transfer through Stationary and Rotating Ribbed Ducts Using a Non-linear <Emphasis Type="Italic">k</Emphasis>−<Emphasis Type="Italic">ε</Emphasis> Model

The present paper deals with the prediction of three-dimensional fluid flow and heat transfer in rib-roughened ducts of square cross-section, which are either stationary, or rotate in orthogonal mode. The main objective is to assess how a recently developed variant of a cubic non-linear k−ε model (proposed by Craft et al. Flow Turbul Combust 63:59–80, 1999) can predict three-dimensional flow and heat transfer characteristics through stationary and rotating ribbed ducts. The present paper discusses turbulent air flow and heat transfer through two different configurations, namely: (I) a stationary square duct with “in-line” normal and (II) a square duct with normal ribs in a “staggered” arrangement under stationary and rotating conditions, with the axis of rotation normal to the flow direction and parallel to the ribs. In this paper the flow and thermal predictions of the linear k−ε model (EVM) are also included, as a set of baseline predictions. The mean flow predictions show that both linear and non-linear k−ε models can successfully reproduce most of the measured data for stream-wise and cross-stream velocity components. Moreover, the non-linear model is able to produce better results for the turbulent stresses. The heat transfer predictions show that both EVM and NLEVM2, the more recent variant of the non-linear k−ε, with the algebraic length-scale correction term, overestimate the measured Nusselt numbers for both geometries examined. While the EVM with the differential length-scale correction term underestimates heat transfer levels, the Nusselt number predictions with the NLEVM2 and the ‘NYP’ term are in close agreements with the measured data. Comparisons with our earlier work, Iacovides and Raisee (Int J Heat Fluid Flow, 20:320–328, 1999), show that the NLEVM2 thermal predictions are of similar quality to those of a second-moment closure.     

Development of an Analytical <Emphasis Type="Italic">β</Emphasis>-Function PDF Integration Algorithm for Simulation of Non-premixed Turbulent Combustion

Among many presumed-shape pdf approaches for modeling non-premixed turbulent combustion, the presumed β-function pdf is widely used in the literature. However, numerical integration of the β-function pdf may encounter singularity difficulties at mixture fraction values of Z = 0 or 1. To date, this issue has been addressed by few publications. The present study proposes the Piecewise Integration Method (PIM), an efficient, robust and accurate algorithm to overcome these numerical difficulties with the added benefit of improving computational efficiency. Comparison of this method to the existing numerical integration methods shows that the PIM exhibits better accuracy and greatly increases computational efficiency. The PIM treatment of the β-function pdf integration is first applied to the Burke–Schumann solution in conjunction with the k − ε turbulence model to simulate a CH₄/H₂ bluff-body turbulent flame. The proposed new method is then applied to the same flow using a more complex combustion model, the laminar flamelet model. Numerical predictions obtained by using the proposed β-function pdf integration method are compared to experimental values of the velocity field, temperature and species mass fractions to illustrate the efficiency and accuracy of the present method.     

Development of an Intermittency Equation for the Modeling of the Supersonic/Hypersonic Boundary Layer Flow Transition

An intermittency transport equation is developed in this study to model the laminar-turbulence boundary layer transition at supersonic and hypersonic conditions. The model takes into account the effects of different instability modes associated with the variations in Mach numbers. The model equation is based on the intermittency factor γ concept and couples with the well-known SST k–ω eddy-viscosity model in the solution procedures. The particular features of the present model approach are that: (1) the fluctuating kinetic energy k includes the non-turbulent, as well as turbulent fluctuations; (2) the proposed transport equation for the intermittency factor γ triggers the transition onset through a source term; (3) through the introduction of a new length scale normal to wall, the present model employs the local variables only avoiding the use of the integral parameters, like the boundary layer thickness δ, which are often cost-ineffective with the modern CFD methods; (4) in the fully turbulent region, the model retreats to SST model. This model is validated with a number of available experiments on boundary layer transition including the incompressible, supersonic and hypersonic flows past flat plates, straight/flared cones at zero incidences, etc. It is demonstrated that the present model can be successfully applied to the engineering calculations of a variety of aerodynamic flow transition with a reasonably wide range of Mach numbers.     

Computation of conjugate heat transfer in the turbulent mixed convection regime in a vertical channel with multiple heat sources

This paper presents the results of a comprehensive numerical study to analyze conjugate, turbulent mixed convection heat transfer from a vertical channel with four heat sources, uniformly flush-mounted to one of the channel walls. The results are presented to study the effect of various parameters like thermal conductivity of wall material (k _s), thermal conductivity of flush-mounted discrete heat source (k _c), Reynolds number of fluid flow (Re _s), modified Richardson number (Ri ⁺) and aspect ratio (AR) of the channel. The standard k-ε turbulence model, modified by including buoyancy effects with physical boundary conditions, i.e. without wall functions, has been used for the analysis. Semi-staggered, non-uniform grids are used to discretise the two dimensional governing equations, using finite volume method. A correlation, encompassing a wide range of parameters, is developed for the non-dimensional maximum temperature (T ^*) using the asymptotic computational fluid dynamics (ACFD) technique.     

Buoyancy modelling with incompressible SPH for laminar and turbulent flows

This work aims to model buoyant, laminar or turbulent flows, using a two‐dimensional incompressible smoothed particle hydrodynamics model with accurate wall boundary conditions. The buoyancy effects are modelled through the Boussinesq approximation coupled to a heat equation, which makes it possible to apply an incompressible algorithm to compute the pressure field from a Poisson equation. Based on our previous work [1], we extend the unified semi‐analytical wall boundary conditions to the present model. The latter is also combined to a Reynolds‐averaged Navier–Stokes approach to treat turbulent flows. The k ? ? turbulence model is used, where buoyancy is modelled through an additional term in the k ? ? equations like in mesh‐based methods. We propose a unified framework to prescribe isothermal (Dirichlet) or to impose heat flux (Neumann) wall boundary conditions in incompressible smoothed particle hydrodynamics. To illustrate this, a theoretical case is presented (laminar heated Poiseuille flow), where excellent agreement with the theoretical solution is obtained. Several benchmark cases are then proposed: a lock‐exchange flow, two laminar and one turbulent flow in differentially heated cavities, and finally a turbulent heated Poiseuille flow. Comparisons are provided with a finite volume approach using an open‐source industrial code. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical study of turbulent flows over wavy boundaries

The fully developed turbulent flows over wavy boundaries are investigated by means of thek-ε model. Predicted flow characteristics over rigid wavy walls are in good agreement with the vailable experimental data. Moreover drag reduction has been found in a 2-dimensional channel with periodical wavy walls. The energy input from turbulent wind to regular waves is also studied in the paper by the same turbulence model with carefully posed boundary conditions at wind-wave interface. Better agreement has been obtained in the predication of the growth rates of wind waves as compared with the previous theoretical and numerical results. The project supported by the National Natural Science Foundation of China.     

Numerical simulation of 3-D turbulent flows over dredged trenches

A 3-D free surface flow in open channels based on the Reynolds equations with thek-ε turbulence closure model is presented in this paper. Insted of the “rigid lid” approximation, the solution of the free surface equation is implemented in the velocity—pressure iterative procedure on the basis of the conventional SIMPLE method. This model was used to compute the flow in rectangular channels with trenches dredged across the bottom. The velocity, eddy viscosity coefficient, turbulent shear stress, turbulent kinetic energy and elevation of the free surface can be obtained. The computed results are in good agreement with previous experimental data.     

A New Low-Reynolds-Number One-Equation Model of Turbulence

In this study, we propose a new Low-Reynolds-Number (LRN)one-equation model, which is derived from an LRN two-equation(k-ε) model. The derivation of the transport equation, in principle, is based on the assumption that the turbulent structure parameter remains constant. However, the relation for the turbulent structure parameter a ₁(=|− |/k) is modified to account for near-wall turbulence. As a result, the present one-equation model contains a term which takes the near-wall limiting behavior explicitly into account. Thus, the present model provides the correct wall-limiting behavior of turbulence in the vicinity of the wall and can be applied to the analysis of heat transfer. The validity of the present model is tested in channel flows, boundary layer flows with and without pressure gradient, plane wall jet, and flow with separation and reattachment. The calculated results showed good agreement with the direct numerical simulation (DNS) and experimental data. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical investigation of corner angle and wing number effects on fluid flow characteristics of a turbulent stellar jet

In this research the fluid dynamics characteristics of a stellar turbulent jet flow is studied numerically and the results of three dimensional jet issued from a stellar nozzle are presented. A numerical method based on control volume approach with collocated grid arrangement is employed. The turbulent stresses are approximated using k–ε and k–ω models with four different inlet conditions. The velocity field is presented and the rate of decay at jet centerline is noted. Special attention is drawn on the influence of corner angle and number of wings on mixing in stellar cross section jets. Stellar jets with three; four and five wings and 15–65° corner angles are studied. Also the effect of Reynolds number (based on hydraulic diameter) as well as the inflow conditions on the evolution of the stellar jet is studied. The Numerical results show that the jet entrains more with corner angle 65° and five wings number. The jet is close to a converged state for high Reynolds numbers. Also the influence of the inflow conditions on the jet characteristics is so strong.     

Large eddy simulation of fuel injection and mixing process in a diesel engine

The large eddy simulation(LES) approach implemented in the KIVA-3V code and based on one-equation sub-grid turbulent kinetic energy model are employed for numerical computation of diesel sprays in a constant volume vessel and in a Caterpillar 3400 series diesel engine.Computational results are compared with those obtained by an RANS(RNG k-ε) model as well as with experimental data.The sensitivity of the LES results to mesh resolution is also discussed.The results show that LES generally provides flow and spray characteristics in better agreement with experimental data than RANS;and that small-scale random vortical structures of the in-cylinder turbulent spray field can be captured by LES.Furthermore,the penetrations of fuel droplets and vapors calculated by LES are larger than the RANS result,and the sub-grid turbulent kinetic energy and sub-grid turbulent viscosity provided by the LES model are evidently less than those calculated by the RANS model.Finally,it is found that the initial swirl significantly affects the spray penetration and the distribution of fuel vapor within the combustion chamber.     

Deduction of Saint-Venant"s Principle for the Asymptotic Residue of Elastic Rods

Let u(ε) be a rescaled 3-dimensional displacement field solution of the linear elastic model for a free prismatic rod Ω^ε having cross section with diameter of order ε, and let u ⁽⁰⁾ –Bernoulli–Navier displacement – and u ⁽²⁾ be the two first terms derived from the asymptotic method. We analyze the residue r(ε) = u(ε) − (u ⁽⁰⁾ + ε² u ⁽²⁾) and if the cross section is star-shaped, we prove such residue presents a Saint-Venant"s phenomenon near the ends of the rod. This revised version was published online in August 2006 with corrections to the Cover Date.     

Predictions of Axial and Transverse Injection into Supersonic Flow

CFD calculations are performed on a swept ramp injecting fuelaxially, and on a two-hole transverse fuel injection downstream of abackward-facing step, both into a supersonic turbulent flow. Theresulting complex flowfields are predicted using a cubick–ε turbulence closure. Comparisons with experimental data show very good agreement. A discussion of the main flow features is presented. The fast computational convergence demonstrates the readiness of the method for the design cycle. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Turbulence Model Sensitivity Study for CH4/H2 Bluff-Body Stabilized Flames

The objective of this work is to assess the performances of different turbulence models in predicting turbulent diffusion flames in conjunction with the flamelet model.The k– model, the Explicit Algebraic Stress Model (EASM) and the k– model withvaried anisotropy parameter C (LEA k– model)are first applied to the inert turbulent flow over a backward-facing step, demonstrating the quality of the turbulence models. Following this, theyare used to simulate the CH₄/H₂ bluff-body flame studied by the University of Sydney/Sandia.The numerical results are compared to experimental values of the mixture fraction, velocity field, temperature and constituent mass fractions.The comparisons show that the overall result depends on the turbulence model used, and indicate that theEASM and the LEA k– models perform better than the k– model and mimic most of the significant flow features.     

Spinodal Decomposition for the¶Cahn-Hilliard Equation in Higher Dimensions:¶Nonlinear Dynamics

This paper addresses the phenomenon of spinodal decomposition for the Cahn-Hilliard equation

Interaction of turbulent natural convection and surface thermal radiation in inclined square enclosures

Two-dimensional numerical studies of flow and temperature fields for turbulent natural convection and surface radiation in inclined differentially heated enclosures are performed. Investigations are carried out over a wide range of Rayleigh numbers from 10⁸ to 10¹², with the angle of inclination varying between 0° and 90°. Turbulence is modeled with a novel variant of the k–ε closure model. The predicted results are validated against experimental and numerical results reported in literature. The effect of the inclination of the enclosure on pure turbulent natural convection and the latter’s interaction with surface radiation are brought out. Profiles of turbulent kinetic energy and effective viscosity are studied to observe the net effect on the intensity of turbulence caused by the interaction of natural convection and surface radiation. The variations of local Nusselt number and average Nusselt number are presented for various inclination angles. Marked change in the convective Nusselt number is found with the orientation of enclosure. Also analyzed is the influence of change in emissivity on the flow and heat transfer. A correlation relevant to practical applications in the form of average Nusselt number, as a function of Rayleigh number, Ra, radiation convection parameter, N _RC and inclination angle of the enclosure, φ is proposed.