Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

BOUNDARY TREATMENT AND AN EFFICIENT PRESSURE ALGORITHM FOR INTERNAL TURBULENT FLOWS USING THE PDF METHOD

The generalized Langevin model, which is used to model the motion of stochastic particles in the velocity–composition joint probability density function (PDF) method for reacting turbulent flows, has been extended to incorporate solid wall effects. Anisotropy of Reynolds stresses in the near-wall region has been addressed. Numerical experiments have been performed to demonstrate that the forces in the near-wall region of a turbulent flow cause the stochastic particles approachi ng a solid wall to reverse their direction of motion normal to the wall and thereby, leave the near-wall layer. This new boundary treatment has subsequently been implemented in a full-scale problem to prove its validity. The test problem considered here is that of an isothermal, non-reacting turbulent flow in a two-dimensional channel with plug inflow and a fixed back-pressure. An efficient pressure correction method, developed in the spirit of the PISO algorithm, has been implemented. The pressure correction strategy is easy to implement and is completely consistent with the time- marching scheme used for the solution of the Lagrangian momentum equations. The results show remarkable agreement with both k–ϵ and algebraic Reynolds stress model calculations for the primary velocity. The secondary flow velocity and the turbulent moments are in better agreement with the algebraic Reynolds stress model predictions than the k– ϵ predictions. © 1997 by John Wiley & Sons, Ltd.     

Calculation of turbulent fluid flow and heat transfer in ducts by a full Reynolds stress model

A computational method has been developed to predict the turbulent Reynolds stresses and turbulent heat fluxes in ducts by different turbulence models. The turbulent Reynolds stresses and other turbulent flow quantities are predicted with a full Reynolds stress model (RSM). The turbulent heat fluxes are modelled by a SED concept, the GGDH and the WET methods. Two wall functions are used, one for the velocity field and one for the temperature field. All the models are implemented for an arbitrary three‐dimensional channel. Fully developed condition is achieved by imposing cyclic boundary conditions in the main flow direction. The numerical approach is based on the finite volume technique with a non‐staggered grid arrangement. The pressure–velocity coupling is handled by using the SIMPLEC‐algorithm. The convective terms are treated by the van Leer scheme while the diffusive terms are handled by the central‐difference scheme. The hybrid scheme is used for solving the ε equation. The secondary flow generation using the RSM model is compared with a non‐linear k–ε model (non‐linear eddy viscosity model). The overall comparison between the models is presented in terms of the friction factor and Nusselt number. Copyright © 2003 John Wiley & Sons, Ltd.     

LES of the adverse-pressure gradient turbulent boundary layer

We describe large-eddy simulations (LES) of the flat-plate turbulent boundary layer in the presence of an adverse pressure gradient. The stretched-vortex subgrid-scale model is used in the domain of the flow coupled to a wall model that explicitly accounts for the presence of a finite pressure gradient. The LES are designed to match recent experiments conducted at the University of Melbourne wind tunnel where a plate section with zero pressure gradient is followed by section with constant adverse pressure gradient. First, LES are described at Reynolds numbers based on the local free-stream velocity and the local momentum thickness in the range 6560–13,900 chosen to match the experimental conditions. This is followed by a discussion of further LES at Reynolds numbers at approximately 10 times and 100 times these values, which are well out of range of present day direct numerical simulation and wall-resolved LES. For the lower Reynolds number runs, mean velocity profiles, one-point turbulent statistics of the velocity fluctuations, skin friction and the Clauser and acceleration parameters along the streamwise, adverse pressure-gradient domain are compared to the experimental measurements. For the full range of LES, the relationship of the skin-friction coefficient, in the form of the ratio of the local free-stream velocity to the local friction velocity, to both Reynolds number and the Clauser parameter is explored. At large Reynolds numbers, a region of collapse is found that is well described by a simple log-like empirical relationship over two orders of magnitude. This is expected to be useful for constant adverse-pressure gradient flows. It is concluded that the present adverse pressure gradient boundary layers are far from an equilibrium state.     

Near-wall damping in model predictions of separated flows

Different near-wall scalings are reviewed by the use of data from direct numerical simulations (DNS) of attached and separated adverse pressure gradient turbulent boundary layers. The turbulent boundary layer equation is analysed in order to extend the validity of existing wall damping functions to turbulent boundary layers under severe adverse pressure gradients. A proposed near-wall scaling is based on local quantities and the wall distance, which makes it applicable for general computational fluid dynamics (CFD) methods. It was found to have a similar behaviour as the pressure-gradient corrected analytical y^* scaling and avoids the inconsistencies present in the y⁺ scaling. The performance of the model is illustrated by model computations using explicit algebraic Reynolds stress models with near-wall damping based on different scalings.     

Quadratic programming algorithm for wall slip and free boundary pressure condition

Wall slip is often observed in a highly sheared fluid film in a solid gap. This makes a difficulty in mathematical analysis for the hydrodynamic effect because fluid velocity at the liquid–solid interfaces is not known a priori. If the gap has a convergent–divergent wedge, a free boundary pressure condition, i.e. Reynolds pressure boundary condition, is usually used in the outlet zone in numerical solution. This paper, based on finite element method and parametric quadratic programming technique, gives a numerical solution technique for a coupled boundary non‐linearity of wall slip and free boundary pressure condition. It is found that the numerical error decreases with the number of elements in a negative power law having an index larger than 2. Our method does not need an iterative process and can simultaneously gives rise to fluid film pressure distribution, wall slip velocity and surface shear stress. Wall slip always decreases the hydrodynamic pressure. Large wall slip even causes a null hydrodynamic pressure in a pure sliding solid gap. Copyright © 2005 John Wiley & Sons, Ltd.     

A bypass wake induced laminar/turbulent transition

The process of laminar to turbulent transition induced by a von Karman vortex street wake, was studied for the case of a flat plate boundary layer. The boundary layer developed under zero pressure gradient conditions. The vortex street was generated by a cylinder positioned in the free stream. An X-type hot-wire probe located in the boundary layer, measured the streamwise and normal to the wall velocity components. The measurements covered two areas; the region of transition onset and development and the region where the wake and the boundary layer merged producing a turbulent flow. The evolution of Reynolds stresses and rms-values of velocity fluctuations along the transition region are presented and discussed. From the profiles of the Reynolds stress and the mean velocity profile, a ‘negative' energy production region along the transition region, was identified. A quadrant splitting analysis was applied to the instantaneous Reynolds stress signals. The contributions of the elementary coherent structures to the total Reynolds stress were evaluated, for several x-positions of the near wall region. Distinct regions in the streamwise and normal to the wall directions were identified during the transition.     

Laminar flow through channels with porous walls and with an applied transverse magnetic field

Summary The problem of two-dimensional steady laminar flow of a viscous incompressible and electrically conducting fluid through a channel with two equally porous walls in the presence of a transverse magnetic field has been extended to include all values of Hartmann number and small suction velocity at the walls. Expressions for the velocity components, the pressure and the wall friction in terms of the Hartmann number and the suction Reynolds number are given. It is found that the pressure drop in the major flow direction and the wall friction decrease with the increase in suction and increase with the increase in the strength of the magnetic field.     

倾斜吹吸控制下湍流边界层减阻的直接数值模拟

雷诺切应力是壁湍流高摩擦阻力的重要来源, 有理论认为可以通过壁面生成负雷诺应力(数值上为正)的方式来削弱湍流流场中雷诺应力的分布, 以此获得流动减阻. 而通过对雷诺平均运动方程的法向二次积分, 可以发现壁面生成正雷诺应力(数值上为负)对壁面摩擦阻力系数才有负贡献. 文中在湍流边界层流动的控制区域下边界设置一系列倾斜狭缝, 利用该装置通过周期性吹吸的方法产生壁面生成正(负)雷诺应力, 并采用直接数值模拟方法考察和验证上文提到的减阻理论. 文中采用的湍流边界层流动模型, 其流动雷诺数(基于外流速度及动量损失厚度)从300 发展到860. 文中通过多组数值模拟算例, 考察了射流强度和频率对壁面摩擦阻力系数的影响, 并对比了壁面生成正或负雷诺应力对流动的影响. 研究表明, 壁面生成正雷诺应力控制的减阻率能达到3.26, 而壁面生成负雷诺应力控制的减阻效果较壁面生成正雷诺应力控制的要差; 壁面生成的正雷诺应力对壁面摩擦阻力有负贡献, 而壁面生成的负雷诺应力对壁面摩擦阻力有正贡献; 通过考察控制的收支比, 发现控制方案不能获得能量净收益.      

A fully coupled Navier–Stokes solver for calculation of turbulent incompressible free surface flow past a ship hull

This paper deals with the calculation of free surface flow of viscous incompressible fluid around the hull of a boat moving with rectilinear motion. An original method used to avoid a large part of the theoretical problems connected with free surface boundary conditions in three‐dimensional Navier–Stokes–Reynolds equations is proposed here. The linearised system of convective equations for velocities, pressure and free surface elevation unknowns is discretised by finite differences and two methods to solve the fully coupled resulting matrix are presented here. The non‐linear convergence of fully coupled algorithm is compared with the velocity–pressure weakly coupled algorithm SIMPLER. Turbulence is taken into account through Reynolds decomposition and k–ε or k–ω model to close the equations. These two models are implemented without wall function and numerical calculations are performed up to the viscous sub‐layer. Numerical results and comparisons with experiments are presented on the Series 60 CB=0.60 ship model for a Reynolds number Rn=4.5×10⁶ and a Froude number Fn=0.316. Copyright © 1999 John Wiley & Sons, Ltd.     

Pulsatile magneto-biofluid flow and mass transfer in a non-Darcian porous medium channel

A numerical study of pulsatile flow and mass transfer of an electrically conducting Newtonian biofluid via a channel containing porous medium is considered. The conservation equations are transformed and solved under boundary conditions prescribed at both walls of the channel, using a finite element method with two-noded line elements. The influence of magnetic field on the flow is studied using the dimensionless hydromagnetic number, Nm, which defines the ratio of magnetic (Lorentz) retarding force to the viscous hydrodynamic force. A Darcian linear impedance for low Reynolds numbers is incorporated in the transformed momentum equation and a second order drag force term for inertial (Forchheimer) effects. Velocity and concentration profiles across the channel width are plotted for various values of the Reynolds number (Re), Darcy parameter (λ), Forchheimer parameter (Nf), hydro-magnetic number (Nm), Schmidt number (Sc) and also with dimensionless time (T). Profiles of velocity varying in space and time are also provided. The conduit considered is rigid with a pulsatile pressure applied via an appropriate pressure gradient term. Increasing the hydromagnetic number (Nm) from 1 to 15 considerably depresses biofluid velocity (U) indicating that a magnetic field can be used as a flow control mechanism in, for example, medical applications. A rise in Nf from 1 to 20 strongly retards the flow development and decreases the velocity, U, across the width of the channel. The effects of other parameters on the flowfield are also discussed at length. The flow model also has applications in the analysis of electrically conducting haemotological fluids flowing through filtration media, diffusion of drug species in pharmaceutical hydromechanics, and also in general fluid dynamics of pulsatile systems.     

Effects of slip velocity and surface acoustic wave on the laminar flow stability

The stability of slip flows when a surface acoustic wave (SAW) propagating along the walls of a microchannel in the laminar flow regime is investigated. The governing equation which was derived by considering the weakly nonlinear coupling between the deformable wall and streaming slip flow is linearized and then the eigenvalue problem is solved by a numerical code together with the associated interface and slip velocity boundary conditions. The value of the critical Reynolds number was found to be near 1,441 for a Knudsen number being 0.001 (associated with a physical parameter K ₀ characterizing the SAW effect) which is much smaller than the static-wall case for conventional pressure-driven flows.     

Computational and experimental investigation of flow over a transient pitching hydrofoil

The present study is developed within the framework of marine structure design operating in transient regimes. It deals with an experimental and numerical investigation of the time–space distribution of the wall-pressure field on a NACA66 hydrofoil undergoing a transient up-and-down pitching motion from 0° to 15° at four pitching velocities and a Reynolds number Re = 0.75 × 10⁶. The experimental investigation is performed using an array of wall-pressure transducers located on the suction side and by means of time–frequency analysis and Empirical Modal Decomposition method. The numerical study is conducted for the same flow conditions. It is based on a 2D RANS code including mesh reconstruction and an ALE formulation in order to take into account the foil rotation and the tunnel walls. Due to the moderate Reynolds number, a laminar to turbulent transition model was also activated. For the operating flow conditions of the study, experimental and numerical flow analysis revealed that the flow experiences complex boundary layer events as leading-edge laminar separation bubble, laminar to turbulent transition, trailing-edge separation and flow detachment at stall. Although the flow is relatively complex, the calculated wall pressure shows a quite good agreement with the experiment provided that the mesh resolution and the temporal discretization are carefully selected depending on the pitching velocity. It is particularly shown that the general trend of the wall pressure (low frequency) is rather well predicted for the four pitching velocities with for instance a net inflection of the wall pressure when transition occurs. The inflection zone is reduced as the pitching velocity increases and tends to disappear for the highest pitching velocity. Conversely, high frequency wall-pressure fluctuations observed experimentally are not captured by the RANS model. Based on the good agreement with experiment, the model is then used to investigate the effects of the pitching velocity on boundary layer events and on hydrodynamic loadings. It is shown that increasing the pitching velocity tends to delay the laminar-to-turbulence transition and even to suppress it for the highest pitching velocity during the pitch-up motion. It induces also an increase of the stall angle (compared to quasi-static one) and an increase of the hysteresis effect during pitch-down motion resulting to a significant increase of the hydrodynamic loading.     

On the validity of the perturbation approach for the flow inside weakly modulated channels

The equations governing the flow of a viscous fluid in a two‐dimensional channel with weakly modulated walls have been solved using a perturbation approach, coupled to a variable‐step finite‐difference scheme. The solution is assumed to be a superposition of a mean and perturbed field. The perturbation results were compared to similar results from a classical finite‐volume approach to quantify the error. The influence of the wall geometry and flow Reynolds number have extensively been investigated. It was found that an explicit relation exists between the critical Reynolds number, at which the wall flow separates, and the dimensionless amplitude and wavelength of the wall modulation. Comparison of the flow shows that the perturbation method requires much less computational effort without sacrificing accuracy. The differences in predicted flow is kept well around the order of the square of the dimensionless amplitude, the order to which the regular perturbation expansion of the flow variables is carried out. Copyright © 2002 John Wiley & Sons, Ltd.     

A spectral–finite difference solution of the Navier–Stokes equations in three dimensions

A new computational code for the numerical integration of the three-dimensional Navier–Stokes equations in their non-dimensional velocity–pressure formulation is presented. The system of non-linear partial differential equations governing the time-dependent flow of a viscous incompressible fluid in a channel is managed by means of a mixed spectral–finite difference method, in which different numerical techniques are applied: Fourier decomposition is used along the homogeneous directions, second-order Crank–Nicolson algorithms are employed for the spatial derivatives in the direction orthogonal to the solid walls and a fourth-order Runge–Kutta procedure is implemented for both the calculation of the convective term and the time advancement. The pressure problem, cast in the Helmholtz form, is solved with the use of a cyclic reduction procedure. No-slip boundary conditions are used at the walls of the channel and cyclic conditions are imposed at the other boundaries of the computing domain. Results are provided for different values of the Reynolds number at several time steps of integration and are compared with results obtained by other authors. © 1998 John Wiley & Sons, Ltd.     

On the performance of sudden‐expansion pipe without/with cross‐flow injection: experimental and numerical studies

The paper explores the possibilities that different turbulence closures offer, for in‐depth analysis of a complex flow. The case under investigation is steady, turbulent flow in a pipe with sudden expansion without/with normal‐to‐wall injection through jets. This is a typical geometry where generation of turbulence energy takes place, due to sudden change in boundary conditions. This study is aimed at investigating the capability of a developed computational program by the present authors with three different turbulence models to calculate the mean flow variables. Three two‐equation models are implemented, namely the standard linear k ? ε model, the low Reynolds number k ? ε model and the cubic nonlinear eddy viscosity (NLEV) k ? ε model. The performance of the chosen turbulence models is investigated with regard to the available data in the literature including velocity profiles, turbulent kinetic energy and reattachment position in a pipe expansion. In order to further assess the reliability of the turbulence models, an experimental program was conducted by the present authors also at the fluid mechanics laboratory of Menoufiya University. Preliminary measurements, including the surface pressure along the two walls of the expansion pipe and the pressure drop without and with the presence of different arrangements of wall jets produced by symmetrical or asymmetrical fluid cross‐flow injection, are introduced. The results of the present studies demonstrate the superiority of the cubic NLEV k ? ε model in predicting the flow characteristics over the entire domain. The simple low Reynolds number k ? ε model also gives good prediction, especially when the reattachment point is concerned. The evaluation of the reattachment point and the pressure‐loss coefficient is numerically addressed in the paper using the cubic NLEV k ? ε model. The results show that the injection location can control the performance of the pipe‐expansion system. It is concluded that the introduction of flow injection can increase the energy loss in the pipe expansion. The near‐field turbulence structure is also considered in the present study and it is noticed that the turbulence level is strongly affected by the cross‐flow injection and the jet location. Copyright © 2011 John Wiley & Sons, Ltd.     

NUMERICAL PREDICTION OF SEMI-CONFINED JET IMPINGEMENT AND COMPARISON WITH EXPERIMENTAL DATA

The standard k–ε eddy viscosity model of turbulence in conjunction with the logarithmic law of the wall has been applied to the prediction of a fully developed turbulent axisymmetric jet impinging within a semi-confined space. A single geometry with a Reynolds number of 20,000 and a nozzle-to-plate spacing of two diameters has been considered with inlet boundary conditions based on measured profiles of velocity and turbulence. Velocity, turbulence and heat transfer data have been obtained using laser–Doppler anemometry and liquid crystal thermography respectively. In the developing wall jet, numerical results of heat transfer compare to within 20% of experiment where isotropy prevails and the trends in turbulent kinetic energy are predicted. However, stagnation point heat transfer is overpredicted by about 300%, which is attributed directly to the turbulence model and inapplicability of the wall function.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

Numerical simulation of turbulent free surface flow with two‐equation k–ε eddy‐viscosity models

This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k–ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k–ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k–ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k–ε model are somewhat unsatisfactory. Copyright © 2004 John Wiley & Sons, Ltd.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.