FINITE ELEMENT SOLUTIONS OF LAMINAR AND TURBULENT MIXED CONVECTION IN A DRIVEN CAVITY

This paper presents laminar and turbulent mixed convection solutions of a driven cavity flow using the finite element method. For the laminar flow, distributions of velocity and temperature with and without the effect of buoyancy force are presented and compared. For the turbulent flow, governing partial differential equations of the thermal turbulence two-equati on model and kinetic turbulence two-equation model are used. Corresponding results such as kinetic eddy diffusivity, kinetic eddy energy, thermal eddy energy and their dissipations are presented.     

Numerical modeling of a two-dimensional vertical turbulent two-phase jet

The results of numerically modeling two-dimensional two-phase flow of the “gas-solid particles” type in a vertical turbulent jet are presented for three cases of its configuration, namely, descending, ascending, and without account of gravity. Both flow phases are modeled on the basis of the Navier-Stokes equations averaged within the framework of the Reynolds approximation and closed by an extended k-? turbulence model. The averaged two-phase flow parameters (particle and gas velocities, particle concentration, turbulent kinetic energy, and its dissipation) are described using the model of mutually-penetrating continua. The model developed allows for both the direct effect of turbulence on the motion of disperse-phase particles and the inverse effect of the particles on turbulence leading to either an increase or a decrease in the turbulent kinetic energy of the gas. The model takes account for gravity, viscous drag, and the Saffman lift. The system of equations is solved using a difference method. The calculated results are in good agreement with the corresponding experimental data which confirms the effect of solid particles on the mean and turbulent characteristics of gas jets.     

A stability condition for turbulence model： From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale （EMMS） model for gas-solid systems, is formulated to close the dynamic constraint equa- tions of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the ＂EMMS-based turbulence model＂, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room, The numerical results show that the EMMS-hased turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

A stability condition for turbulence model: From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale (EMMS) model for gas–solid systems, is formulated to close the dynamic constraint equations of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the “EMMS-based turbulence model”, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room. The numerical results show that the EMMS-based turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

MODELLING OF BYPASS TRANSITION WITH CONDITIONED NAVIER–STOKES EQUATIONS COUPLED TO AN INTERMITTENCY TRANSPORT EQUATION

A differential method is proposed to simulate bypass transition. The intermittency in the transition zone is taken into account by conditioned averages. These are averages taken during the fraction of time the flow is turbulent or laminar respectively. Starting from the Navier–Stokes equations, conditioned continuity, momentum and energy equations are derived for the laminar and turbulent parts of the intermittent flow. The turbulence is described by a classical k−ϵ model. The supplementary parameter, the intermittency factor, is determined by a transport equation applicable for zero, favourable and adverse pressure gradients. Results for these pressure gradients are given.     

Turbulent boundary-layer analysis using finite elements

The use of finite element methods for turbulent boundary-layer flow is relatively recent and of limited extent.¹ In the present study, we extend the group variable approach of Fletcher and Fleer^2,3 to treat turbulent boundary layer flows with heat transfer using a two-equation turbulence model. The main concepts in the formulations include a Dorodnitsyn-type transformation which uses a velocity component as the transverse variable, a ‘variational’ formulation for the transformed equations using special test functions and development of a two-equation turbulence model in terms of the turbulent kinetic energy and turbulence dissipation rate as additional field variables. Several numerical test cases have been examined comparing the results with finite difference calculations and comparing the two-equation turbulence model with an algebraic turbulence model.     

Prediction of the particle-laden jet with a two-equation turbulence model

A two-equation turbulence model for two-phase flows has recently been proposed by Elghobashi & Abou-Arab (1983). They derived the exact equations of the kinetic energy of turbulence and the rate of dissipation of that energy, and modeled the turbulent correlations, resulting from time-averaging, up to third order. In order to validate the proposed model, a turbulent axisymmetric gaseous jet laden with spherical uniform-size solid particles is studied here. The predictions of the mean flow properties of the two-phases and the turbulence kinetic energy and shear stress of the carrier phase show good agreement with the experimental data.     

Computations for a jet impinging obliquely on a flat surface

A SIMPLE-C algorithm and Jones-Launder k-ε two-equation turbulence model are used to simulate a two-dimensional jet impinging obliquely on a flat surface. Both the continuity and momentum equations for the unsteady state are cast into suitable finite difference equations. The pressure, velocity, turbulent kinetic energy and turbulent energy dissipation rate distributions are solved and show good agreement with various experimental data. The calculations show that the flow field structure of the jet impinging obliquely on a flat surface is strongly affected by the oblique impingement angle. The maximum pressure zone of the obliquely impinging jet flow field moves towards the left as the oblique impingement angle is decreased.     

Structure of turbulent pipe flow

The problem of turbulent flow in a straight circular pipe is solved. We consider a system consisting of the equation of motion, the equation for the turbulence energy, the expression relating the turbulence coefficient with the turbulence scale, and the integral formula for determining the turbulence scale. A numerical solution is presented for this closed system of equations for turbulent flow. The results of calculations are compared with experimental data.     

Computation of Compressible Flows using a Two-equation Turbulence Model on Unstructured Grids

This paper presents a numerical method for solving compressible turbulent flows using a k - l turbulence model on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multi-stage Runge-Kutta time stepping scheme, while the turbulence equations are advanced using a multi-stage point-implicit scheme. The positivity of turbulence variables is achieved using a simple change of dependent variables. The developed method is used to compute a variety of turbulent flow problems. The results obtained are in good agreement with theoretical and experimental data, indicating that the present method provides a viable and robust algorithm for computing turbulent flows on unstructured meshes.     

Consequences of the second law for a turbulent fluid flow

Second Law statements in thermomechanics applicable to turbulent fluid flow, in which the internal energy in a macroscopic field theory includes contributions both from molecular vibrations and from turbulent fluctuations, are discussed. In the absence of turbulence, these statements naturally reduce to the known and accepted Second Law statements for a nonturbulent medium. The usual version of the Second Law statements — which deny the existence of perpetual motion and place restrictions on the constitutive equations —is extended here in the presence of turbulence; and an additional statement is introduced associated with the tendency of turbulent fluctuations to decay in the absence of external work or the addition of thermal heat. The mathematical representations of various Second Law statements are then used to derive several restrictions on the response variables of the macroscopic turbulence theory. Examples of such variables include the rates of production and dissipation of turbulent fluctuations, the rate of thermal entropy production, internal energy (involving constitutive coefficients which may be taken to be the thermal and turbulent specific heats), turbulent viscosity coefficients and other response functions which control the degree of flow anisotropy in the medium. These Second Law restrictions are then applied to a recent theory of macroscopic turbulent flow by the present authors in which fairly general constitutive equations are presented for the dependent variables of the theory. It is found that not only is the range of values of several constitutive coefficients limited by these Second Law restrictions, but the presence of a number of terms in the constitutive equations is entirely denied.     

Numerical predictions of flows over backward-facing steps

Predictions are reported for two-dimensional, steady, incompressible flows over rearward-facing steps for both laminar and turbulent conditions. The standard k-? turbulence model was used for the turbulent flow. Attention was focused on obtaining accurate solutions to the differential equations. It is concluded that some of the serious discrepancies that have occurred between prediction and observation, and attributed in earlier studies to the inadequacy of the turbulence model, may have been due to the inaccuracy of the solution.     

Renormalization Group Model of Large-Scale Turbulence in Porous Media

Renormalization group methods are used to develop a macroscopic (large-scale) turbulence model for incompressible flow in porous media. The model accounts for the large-distance and large-time behavior of velocity correlations generated by the momentum equation for a randomly stirred, incompressible flow. Utilizing the renormalization procedure, the transport equations for the large-scale modes and expressions for effective transport coefficients are obtained. Expressions for renormalized turbulent viscosity, which accounts for the ultraviolet subrange of the turbulent kinetic energy spectrum, are also obtained.     

Direct Numerical Simulation of Self-Similar Turbulent Boundary Layers in Adverse Pressure Gradients

Direct numerical simulations of the Navier–Stokes equations have been carried out with the objective of studying turbulent boundary layers in adverse pressure gradients. The boundary layer flows concerned are of the equilibrium type which makes the analysis simpler and the results can be compared with earlier experiments and simulations. This type of turbulent boundary layers also permits an analysis of the equation of motion to predict separation. The linear analysis based on the assumption of asymptotically high Reynolds number gives results that are not applicable to finite Reynolds number flows. A different non-linear approach is presented to obtain a useful relation between the freestream variation and other mean flow parameters. Comparison of turbulent statistics from the zero pressure gradient case and two adverse pressure gradient cases shows the development of an outer peak in the turbulent energy in agreement with experiment. The turbulent flows have also been investigated using a differential Reynolds stress model. Profiles for velocity and turbulence quantities obtained from the direct numerical simulations were used as initial data. The initial transients in the model predictions vanished rapidly. The model predictions are compared with the direct simulations and low Reynolds number effects are investigated.     

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.     

Similarity solutions of vertical plane wall plume based on finite analytic method

The turbulent flow of vertical plane wall plume with concentration variation was studied with the finite analytical method. The k-epsilon model with the effect of buoyancy on turbulent kinetic energy and its dissipation rate was adopted. There were similarity solutions in the uniform environment for the system of equations including the equation of continuity, the equation of momentum along the flow direction and concentration, and equations of k, epsilon. The finite analytic method was applied to obtain the similarity solution. The calculated data of velocity, relative density difference, the kinetic energy of turbulence and its dissipation rate distribution for vertical plane plumes are in good agreement with the experimental data at the turbulent Schmidt number equal to 1.0. The variations of their maximum value along the direction of main flow were also given. It shows that the present model is good, i.e., the effect of buoyancy on turbulent kinetic energy and its dissipation rate should be taken into account, and the finite analytic method is effective.     

A k–$${\varepsilon}$$ turbulence closure model of an isothermal dry granular dense matter

The turbulent flow characteristics of an isothermal dry granular dense matter with incompressible grains are investigated by the proposed first-order k–\({\varepsilon}\) turbulence closure model. Reynolds-filter process is applied to obtain the balance equations of the mean fields with two kinematic equations describing the time evolutions of the turbulent kinetic energy and dissipation. The first and second laws of thermodynamics are used to derive the equilibrium closure relations satisfying turbulence realizability conditions, with the dynamic responses postulated by a quasi-linear theory. The established closure model is applied to analyses of a gravity-driven stationary flow down an inclined moving plane. While the mean velocity decreases monotonically from its value on the moving plane toward the free surface, the mean porosity increases exponentially; the turbulent kinetic energy and dissipation evolve, respectively, from their minimum and maximum values on the plane toward their maximum and minimum values on the free surface. The evaluated mean velocity and porosity correspond to the experimental outcomes, while the turbulent dissipation distribution demonstrates a similarity to that of Newtonian fluids in turbulent shear flows. When compared to the zero-order model, the turbulent eddy evolution tends to enhance the transfer of the turbulent kinetic energy and plane shearing across the flow layer, resulting in more intensive turbulent fluctuation in the upper part of the flow. Solid boundary as energy source and sink of the turbulent kinetic energy becomes more apparent in the established first-order model.     

Applications of a higher‐order bounded numerical scheme to turbulent flows

This paper applies the higher‐order bounded numerical scheme Weighted Average Coefficients Ensuring Boundedness (WACEB) to simulate two‐ and three‐dimensional turbulent flows. In the scheme, a weighted average formulation is used for interpolating the variables at cell faces and the weighted average coefficients are determined from a normalized variable formulation and total variation diminishing (TVD) constraints to ensure the boundedness of the solution. The scheme is applied to two turbulent flow problems: (1) two‐dimensional turbulent flow around a blunt plate; and (2) three‐dimensional turbulent flow inside a mildly curved U‐bend. In the present study, turbulence is evaluated by using a low‐Reynolds number version of the k–ω model. For the flow simulation, the QUICK scheme is applied to the momentum equations while either the WACEB scheme (Method 1) or the UPWIND scheme (Method 2) is used for the turbulence equations. The present study shows that the WACEB scheme has at least second‐order accuracy while ensuring boundedness of the solutions. The present numerical study for a pure convection problem shows that the ‘TVD’ slope ranges from 2 to 4. For the turbulent recirculating flow, two different mixed procedures (Method 1 and Method 2) produce a substantial difference for the mean velocities as well as for the turbulence kinetic energy. Method 1 predicts better results than Method 2 does, comparing the analytical solution and the experimental data. For the turbulent flow inside the mildly curved U‐bend, although the predictions of velocity distributions with two procedures are very close, a noticeable difference of turbulence kinetic energy is exhibited. It is noticed that the discrepancy exists between numerical results and the experimental data. The reason is the limit of the two‐equation turbulence model to such complex turbulent flows with extra strain‐rates. Copyright © 2001 John Wiley & Sons, Ltd.     

Buoyancy effects on turbulent swirling flows

A three-parameter model of turbulence applicable to free boundary layers has been developed and applied for the prediction of axisymmetric turbulent swirling flows in uniform and stagnant surroundings under the action of buoyancy forces. The turbulent momentum and heat fluxes appearing in the time-averaged equations for the mean motion have been determined from algebraic expressions, derived by neglecting the convection and diffusion terms in the differential transport equations for these quantities, which relate the turbulent fluxes to the kinetic energy of turbulence, k, the dissipation length scale of turbulence, L, and the temperature covariance, T ′². Differential transport equations have been used to determine these latter quantities. The governing equations have been solved using fully implicit finite difference schemes. The turbulence model is capable of reproducing the gross features of pure jet flows, buoyant flows and swirling flows for weak and moderate swirl. The behaviour of a turbulent buoyant swirling jet has been found to depend solely on exit swirl and Froude numbers. The predicted results indicate that the incorporation of buoyancy can cause significant changes in the behaviour of a swirling jet, particularly when the buoyancy strength is high. The jet exhibits similarity behaviour in the initial region for weak swirl and weak buoyancy strengths only, and the asymptotic case of a swirling jet under the action of buoyancy forces is a pure plume in the far field. The predicted results have been found to be in satisfactory agreement with the available experimental data and in good qualitative agreement with other predicted results.