Numerical calculation of Reynolds stresses in a square duct with secondary flow

A new model for the Reynolds stress equations is presented. This model is used to obtain a theoretical solution for the problem of fully developed turbulent flow in a square duct. Nine governing equations for the axial velocity, lateral vorticity, lateral stream function and six components of the Reynolds stresses are simultaneously solved, by a finite-difference technique. To ensure numerical stability of the solution a special linearised implicit representation of the source terms is proposed, and simultaneous solution of the equations at each.mesh point is obtained. Near the wall a special procedure is used, by which the Reynolds stress equations are assumed to be in local equilibrium, and the velocity profile is assumed to be logarithmic. However, due to the secondary motion the logarithmic velocity profile is inclined to the axial direction. The results bear reasonable agreement with experimental data. Computer time requirements are moderate.     

Experiments in Turbulent Pipe Flow with Polymer Additives at Maximum Drag Reduction

In this paper we report on (two-component) LDV experiments in a fully developed turbulent pipe flow with a drag-reducing polymer (partially hydrolyzed polyacrylamide) dissolved in water. The Reynolds number based on the mean velocity, the pipe diameter and the local viscosity at the wall is approximately 10000. We have used polymer solutions with three different concentrations which have been chosen such that maximum drag reduction occurs. The amount of drag reduction found is 60–70%. Our experimental results are compared with results obtained with water and with a very dilute solution which exhibits only a small amount of drag reduction. We have focused on the observation of turbulence statistics (mean velocities and turbulence intensities) and on the various contributions to the total shear stress. The latter consists of a turbulent, a solvent (viscous) and a polymeric part. The polymers are found to contribute significantly to the total stress. With respect to the mean velocity profile we find a thickening of the buffer layer and an increase in the slope of the logarithmic profile. With respect to the turbulence statistics we find for the streamwise velocity fluctuations an increase of the root mean square at low polymer concentration but a return to values comparable to those for water at higher concentrations. The root mean square of the normal velocity fluctuations shows a strong decrease. Also the Reynolds (turbulent) shear stress and the correlation coefficient between the stream wise and the normal components are drastically reduced over the entire pipe diameter. In all cases the Reynolds stress stays definitely non-zero at maximum drag reduction. The consequence of the drop of the Reynolds stress is a large polymer stress, which can be 60% of the total stress. The kinetic-energy balance of the mean flow shows a large transfer of energy directly to the polymers instead of the route by turbulence. The kinetic energy of the turbulence suggests a possibly negative polymeric dissipation of turbulent energy. This revised version was published online in July 2006 with corrections to the Cover Date.     

Singularities of the boundary layer equations and the structure of the flow in the vicinity of the convergence plane on conical bodies

The singularities of the boundary layer equations and the laminar viscous gas flow structure in the vicinity of the convergence plane on sharp conical bodies at incidence are analyzed. In the outer part of the boundary layer the singularities are obtained in explicit form. It is shown that in the vicinity of a singularity a boundary domain, in which the flow is governed by the shortened Navier-Stokes equations, is formed; their regular solutions are obtained. The viscous-inviscid interaction effect predominates in a region whose extent is of the order of the square root of the boundary layer thickness, in which the flow is described by a two-layer model, namely, the Euler equations in the slender-body approximation for the outer region and the three-dimensional boundary layer equations; the pressure is determined from the interaction conditions. On the basis of an analysis of the solutions for the outer part of the boundary layer it is shown that interaction leads to attenuation of the singularities and the dependence of the nature of the flow on the longitudinal coordinate, but does not make it possible to eliminate the singularities completely.     

Asymptotic theory of the Boltzmann system,for a steady flow of a slightly rarefied gas with a finite Mach number: General theory

A steady rarefied gas flow with Mach number of the order of unity around a body or bodies is considered. The general behaviour of the gas for small Knudsen numbers is studied by asymptotic analysis of the boundary-value problem of the Boltzmann equation for a general domain. The effect of gas rarefaction (or Knudsen number) is expressed as a power series of the square root of the Knudsen number of the system. A series of fluid-dynamic type equations and their associated boundary conditions that determine the component functions of the expansion of the density, flow velocity, and temperature of the gas is obtained by the analysis. The equations up to the order of the square root of the Knudsen number do not contain non-Navier–Stokes stress and heat flow, which differs from the claim by Darrozes (in Rarefied Gas Dynamics, Academic Press, New York, 1969). The contributions up to this order, except in the Knudsen layer, are included in the system of the Navier–Stokes equations and the slip boundary conditions consisting of tangential velocity slip due to the shear of flow and temperature jump due to the temperature gradient normal to the boundary.     

Damage analysis of thermopiezoelectric properties: Part I - crack tip singularities

This is Part I of the work on a two-dimensional analysis of thermal and electric fields of a thermopiezoelectric solid damaged by cracks. It deals with finding the singular crack tip behavior for the temperature, heat flow, displacements, electric potential, stresses and electric displacements. By application of Fourier transformations and the extended Stroh formalism, the problem is reduced to a pair of dual integral equations for the temperature field with the aid of an auxiliary function. The electroelastic field is governed by another pair of dual integral equations. The inverse square root singularity is found for the heat flow field while the logarithmic singularity prevailed for the electroelastic field regardless of whether the crack lies in a homogeneous piezoelectric solid or at an interface of two dissimilar piezoelectric materials. Results are given for the energy release rate and a finite length crack oriented at an arbitrarily angle with reference to the external disturbances. Part II of this paper considers the modelling of a piezoelectric material containing microcracks. A representative cracked area element is used to obtain the effective conductivity and electroelastic modulus. Numerical results are given for a peizoelectric Bati O₃ ceramic with cracks.     

Exact solution of the Navier-Stokes equations in a fluid layer between the moving parallel plates

This paper studies exact solutions of the Navier-Stokes equations for a layer between parallel plates the distance between which increases proportionally to the square root of time. A countable set of exact solutions and their derived countable set of continuous families of exact solutions are obtained. It is shown that certain intervals of the Reynolds parameter have two solutions and some of them one solution.     

Natural convection flow in a square cavity revisited: Laminar and turbulent models with wall functions

Numerical simulations have been undertaken for the benchmark problem of natural convection flow in a square cavity. The control volume method is used to solve the conservation equations for laminar and turbulent flows for a series of Rayleigh numbers (Ra) reaching values up to 10¹⁰. The k-? model has been used for turbulence modelling with and without logarithmic wall functions. Uniform and non-uniform (stretched) grids have been employed with increasing density to guarantee accurate solutions, especially near the walls for high Ra-values. ADI and SIP solvers are implemented to accelerate convergence. Excellent agreement is obtained with previous numerical solutions, while some discrepancies with others for high Ra-values may be due to a possibly different implementation of the wall functions. Comparisons with experimental data for heat transfer (Nusselt number) clearly demonstrates the limitations of the standard k-? model with logarithmic wall functions, which gives significant overpredictions.     

Rotation correction of wall function and velocity structure under the couple effects of buoyancy force and rotating effects at the entry section of the smooth straight channel

This paper studies the effect of rotation on the turbulent boundary-layer flow in a rotating duct with a square cross section by using hot-wire. The experiments were conducted with the Reynolds numbers, based on the duct's hydraulic diameter (D = 80 mm) equaling 19,000. The rotation numbers (Ro) studied ranged from 0 to 0.362. Hot-wire measurements of the flow field were made at four cross sections of the rotating duct. The effects of rotation on velocity profile, semi-logarithmic mean velocity profile, and wall shear stress are discussed in this paper. Results obtained show the velocity deficit about the leading surface of the rotating duct, created by the secondary flows induced by the Coriolis force, to not increase monotonically with the increase in the Rotation number. Results obtained also show the effects of rotation to penetrate into the logarithm region, and the flow near the leading surface tends to laminarize. In this study, a correction factor is developed for logarithmic law to account for the effects of rotation, which can be used in CFD studies of rotating ducts that use wall functions.     

Computation of three‐dimensional flow around square and circular piers

The aim of the present study is to investigate, by numerical simulation, the three‐dimensional turbulent flow field around square and circular piers. The numerical model employs a finite volume method based on MacCormack's explicit predictor–corrector scheme to solve weakly compressible hydrodynamic equations for turbulent flow. Computed results are compared with Dargahi's experimental measurements to assess the validity of the proposed model. Very good agreements are obtained. The results of flow simulation indicate that near the upstream face of the pier there exists a downflow, which joins the separated flow to form the horseshoe vortex stretched around the pier. This horseshoe vortex interacts with the wake vortex to create the upflow behind the pier. These phenomena appear to be very important to the mechanism of scouring around the pier. In general, the flow patterns for the square and circular piers are similar. However, the strengths of the downflow and horseshoe vortex are greater in the case of the square pier. The position of the horseshoe vortex around the circular pier is closer to the front face than that around the square pier. In the meantime, the domain of the wake flow in the case of the square pier is greater than that in the case of the circular one. Copyright © 2000 John Wiley & Sons, Ltd.     

Investigation of trajectories of inviscid fluid particles in two-dimensional rotating boxes

The particle trajectories of inviscid fluid flow within two-dimensional rotating (elliptic, triangular, and square) boxes are numerically investigated. The source panel method is employed to represent the instantaneous potential interior flow field, and the Runge–Kutta method is used to track the fluid particles. The analytic solutions for the fluid trajectories for the elliptic box are used to verify the numerical accuracy of the method. The numerical error can be reduced to the level of the round-off error if the panels are properly configured and an appropriate number of panels is used. The stagnation of the particles at the corners of the triangular box is successfully predicted with this method. The corner of the square box is found to be a singularity. A logarithmic complex potential is proposed to account for the singularity, using which the stagnation of the particles at the corner in the square box is also captured. The natural frequency of the particles in the rotating elliptic box is constant throughout the flow domain, and the fluid trajectories are epitrochoidal curves. In the triangular box and the square box, the natural frequency strongly depends on the particle position, and the particle trajectories are similar to epitrochoidal curves. In general, the trajectory patterns depend only on the box rotating frequency and the natural frequency of the fluid particle motion.      

Parameters analysis of fluid diffusion

The influences of fluid density, diffusivity, viscosity, width of the flow channel, travel distance, and flow velocity on fluid diffusion are analyzed theoretically and numerically. Concentration boundary layer is taken as the quantitative index of fluid diffusion in this work. The results show that diffusion is a function of travel distance, diffusivity, fluid density, and flow velocity. Diffusion is independent of the width of the channel. Viscous effect determines the velocity gradient and does not affect diffusion directly. The usually used Péclet number uL/D cannot govern the full condition of fluid diffusion. For two-fluids co-flowing in a two-dimensional straight channel with relative low viscous effect, diffusion is proportional to the square root of travel distance and diffusivity, and is inversely proportional to the square root of flow velocity.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Aerodynamic force evaluation for ice shedding phenomenon using vortex in cell scheme,penalisation and level set approaches

In this work, we propose a formulation to evaluate aerodynamic forces for flow solutions based on Cartesian grids, penalisation and level set functions. The formulation enables the evaluation of forces on closed bodies moving at different velocities. The use of Cartesian grids bypasses the meshing issues, and penalisation is an efficient alternative to explicitly impose boundary conditions so that the body fitted meshes can be avoided. Penalisation enables ice shedding simulations that take into account ice piece effects on the flow. Level set functions describe the geometry in a non-parametric way so that geometrical and topological changes resulting from physics, and particularly shed ice pieces, are straightforward to follow. The results obtained with the present force formulation are validated against other numerical formulations for circular and square cylinder in laminar flow. The capabilities of the proposed formulation are demonstrated on ice trajectory calculations for highly separated flow behind a bluff body, representative of inflight aircraft ice shedding.     

Large eddy simulation of hot and cold fluids mixing in a T-junction for predicting thermal fluctuations

Temperature fluctuations in a mixing T-junction have been simulated on the FLUENT platform using the large eddy simulation （LES） turbulent flow model and a sub-grid scale Smagorinsky-Lilly model. The normalized mean and root mean square temperatures for describing time-averaged temperature and temperature fluctuation intensity, and the velocity are obtained. The power spectrum densities of temperature fluctuations, which are key parameters for thermal fatigue analysis and lifetime evaluation, are analyzed. Simulation results are consistent with experimental data published in the literature, showing that the LES is reliable. Several mixing processes under different conditions are simulated in order to analyze the effects of varying Reynolds number and Richardson number on the mixing course and thermal fluctuations.     

A one-dimensional moving grid solution for the coupled non-linear equations governing multiphase flow in porous media. 2: Example simulations and sensitivity analysis

This paper presents numerical examples for the moving grid finite element algorithm derived in Part Ito solve the non-linear coupled set of PDEs governing immiscible multiphase flow in porous media in one dimension. Examples include single- and double-front simulations for two- and three-phase flow regimes and incorporating a mass sink. The modelling approach is shown to achieve significant savings in computation time and memory allocation when compared with fixed grid solutions of equivalent accuracy. This work includes sensitivity analyses for the parameters which are incorporated in the grid adaptation method, including the curvature weights, artificial viscosity and artificial repulsive force. It is found that the curvature weights are exponential functions of the negative ratio of the square root of the domain length to the number of discrete nodes. These weighting parameters are also shown to depend upon the shape of the front. On the basis of the examined simulations, it is recommended that artificial viscosity be neglected in the solution of the coupled non-linear set of PDEs governing multiphase flow in porous media. Similarly, use of a repulsive force is found to be unnecessary in simulations involving the migration of two liquid phases. For multiphase flows incorporating a gas phase it is recommended to use a non-zero value for the repulslive force to avoid development of an ill-conditioned nodal distribution matrix. An equation to evaluate the repulsive force under these circumstances is suggested.     

Comparison of Symmetric and Asymmetric Schemes with Arithmetic and Harmonic Averaging for Fracture Flow on Cartesian Grids

Performance of four finite-difference schemes for fluid flow in rough-walled fractures on regular Cartesian grids is evaluated numerically. The four schemes are an asymmetric scheme with arithmetic averaging, an asymmetric scheme with harmonic averaging, a symmetric scheme with arithmetic averaging, and a symmetric scheme with harmonic averaging. The schemes are compared with respect to their simulated hydraulic aperture and the mass balance error. 1320 flow simulations with different grid sizes, mean fracture aperture and root mean square (RMS)/mean aperture ratio are completed. The asymmetric scheme with arithmetic averaging arises naturally, without any extra assumptions about the correct transmissivity averaging procedure, when one uses second-order finite differences to approximate the generalized Laplace operator expanded as a derivative of a product. Hydraulic apertures obtained with harmonic averaging are found to usually be smaller than those obtained with arithmetic averaging, especially when the ratio of aperture RMS to the mean aperture is larger. The traditionally used asymmetric schemes are found to be superior to symmetric schemes in terms of mass balance accuracy.

     

An explicit transient algorithm for predicting incompressible viscous flows in Arbitrary Geometry

A numerical method for predicting viscous flows in complex geometries has been presented. Integral mass and momentum conservation equations are deploved and these are discretized into algebraic form through numerical quadrature. The physical domain is divided into a number of non-orthogonal control volumes which are isoparametrically mapped on to standard rectangular cells. Numerical integration for unsteady mementum equations is performed over such non-orthogonal cells. The explicitly advanced velocity components obtained from unsteady momentum equations may not necessarily satisfy the mass conservation condition in each cell. Compliance of the mass conservation equation and the consequent evolution of correct pressure distribution are accomplished through an iterative correction of pressure and velocity till divergence-free condition is obtained in each cell. The algorithm is applied on a few test problems, namely, lid-driven square and oblique cavities, developing flow in a rectangular channel and flow over square and circular cylinders placed in rectangular channels. The results exhibit good accuracy and justify the applicability of the algorithm. This Explicit Transient Algorithm for Flows in Arbitrary Geometry is given a generic name EXTRAFLAG.     

Mixed convection boundary layer flow over a horizontal plate with thermal radiation

The effects of thermal radiation and thermal buoyancy on the steady, laminar boundary layer flow over a horizontal plate is investigated. The plate temperature is assumed to be inversely proportional to the square root of the distance from the leading edge. The set of similarity equations is solved numerically, and the solutions are given for some values of the radiation and buoyancy parameters for Prandtl number unity. It is found that dual solutions exist for negative values of the buoyancy parameter, up to certain critical values. Beyond these values, the solution does no longer exist. Moreover, it is found that there is no local heat transfer at the surface except in the singular point at the leading edge. The radiation parameter is found to increase the local Stanton number.