Numerical simulation of two‐dimensional laminar slot‐jet impingement flows confined by a parallel wall

Two‐dimensional laminar incompressible impinging slot‐jet is simulated numerically to gain insight into flow characteristics.Computations are done for vertically downward‐directed slot‐jets impinging on a plate at the bottom and confined by a parallel surface on top. The behaviour of the jet with respect to aspect ratio (AR) and Reynolds number (Re) are described in detail. The computed flow patterns for various AR (2–5) and for a range of jet‐exit Reynolds numbers (100–500) are analysed to understand the flow characteristics. The transient development of the flow is also simulated for AR = 4 and Re = 300. It is found that the reattachment length is dependent on both AR and Reynolds number for the range considered. The correlation for reattachment length is suggested. The maximum resultant velocity V_rmax and its trajectory is reported. A detailed study of horizontal velocity profile at different downstream locations is reported. It is found that the effect of Reynolds number and AR is significant to the bottom wall vorticity in the impingement and wall jet regions. Copyright © 2007 John Wiley & Sons, Ltd.     

Prediction of unobstructed flow for co‐rotating multi‐disk drive in an enclosure

The flow structure and isotherms of hard disk drives (HDD) are investigated using a finite element method (FEM). The governing equations are based on the three‐dimensional axisymmetric Navier–Stokes partial differential equations (PDEs) with Galerkin FE formulation. Co‐rotating models are selected that include the non‐ventilated configuration within an enclosure. With various operating conditions for the disk system, the following important parameters are considered: disk number (n), rotational speed (Ro), and wall temperature. The flow structure changes rapidly when the rotational Reynolds number (Re_ϕ) is increased. The flow has a greater tendency to flow radially outwards and the swirling velocity tends to be more vertically orientated, especially for high Re_ϕ values. The isotherms only have small varying regions near the rotating axis, forming outward arcs near the wall and inward arcs near the end gap of the disk. Different from the case without the enclosure, the vorticities exist along the outer disk ends. Both the swirling velocity and isotherms indicate nearly symmetrical characteristics, as expected. A higher temperature gradient occurs near the right outer disk ends, which implies the characteristic of higher heat flux. A commercial computational fluid dynamic (CFD) code, CFX‐5, was chosen to validate the results. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical investigation of a laminar pulsating flow in a rectangular duct

A pulsating laminar flow of a viscous, incompressible liquid in a rectangular duct has been studied. The motion is induced under an imposed pulsating pressure difference. The problem is solved numerically. Different flow regimes are characterized by a non‐dimensional parameter based on the frequency (ω) of the imposed pressure gradient oscillations and the width of the duct (h). This, in fact, is the Reynolds number of the problem at hand. The induced velocity has a phase lag (shift) with respect to the imposed pressure oscillations, which varies from zero at very slow oscillations, to 90° at fast oscillations. The influence of the aspect ratio of the rectangular duct and the pulsating pressure gradient frequency on the phase lag, the amplitude of the induced oscillating velocity, and the wall shear were analyzed. Copyright © 1999 John Wiley & Sons, Ltd.     

Turbulence in a three‐dimensional wall‐bounded shear flow

A new turbulent flow with distinct three‐dimensional characteristics has been designed in order to study the impact of mean‐flow skewing on the turbulent coherent vortices and Reynolds‐averaged statistics. The skewing of a unidirectional plane Couette flow was achieved by means of a spanwise pressure gradient. Direct numerical simulations of the statistically steady Couette–Poiseuille flow enabled in‐depth explorations of the turbulence field in the skewed flow. The imposition of a modest spanwise gradient turned the mean flow about 8° away from the original Couette flow direction and this turning angle remained nearly the same over the entire cross section. Nevertheless, a substantial non‐alignment between the turbulent shear stress angle and the mean velocity gradient angle was observed. The structure parameter turned out to slightly exceed that in the pure Couette flow, contrary to the observations made in some other three‐dimensional shear flows. Coherent flow structures, which are known to be associated with the Reynolds shear stress in near‐wall regions, were identified by the λ₂‐criterion. Instantaneous and ensemble‐averaged vortices resembled those found in the unidirectional Couette flow. In the skewed flow, however, the vortex structures were turned to align with the local mean‐flow direction. The conventional symmetry between Case 1 and Case 2 vortices was broken due to the mean‐flow three‐dimensionality. The turning of the coherent vortices and the accompanying symmetry‐breaking gave rise to secondary and tertiary turbulent shear stress components. By averaging the already ensemble‐averaged shear stresses associated with Case 1 and Case 2 vortices in the homogeneous directions, a direct link between the educed near‐wall structures and the Reynolds‐averaged turbulent stresses was established. These observations provide evidence in support of the hypothesis that the structural model proposed for two‐dimensional turbulent boundary layers remains valid also in flows with moderate mean three‐dimensionality. Copyright © 2009 John Wiley & Sons, Ltd.     

A second-order dynamic subgrid-scale stress model

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎＬＥＳ,ｔｈｅｌａｒｇｅｓｃａｌｅｓｉｎｔｈｅｆｌｏｗａｒｅｃｏｍｐｕｔｅｄｅｘｐｌｉｃｉｔｌｙａｎｄｔｈｅｓｕｂｇｒｉｄ＿ｓｃａｌｅｓ(ＳＧＳ),ｗｈｉｃｈａｒｅｆｉｌｔｅｒｅｄｏｕｔｂｙａｎａｖｅｒａｇｅｐｒｏｃｅｓｓｏｐｅｒａｔｉｏｎｏｎｔｈｅＮａｖｉｅｒ＿Ｓｔｏｋｅｓｅｑｕａｔｉｏｎｓ,ａｒｅｍｏｄｅｌｅｄ.ＳｉｎｃｅＬＥＳｒｅｑｕｉｒｅｓｌｅｓｓｃｏｍｐｕｔｅｒｔｉｍｅｔｈａｎｔｈｅｄｉｒｅｃｔｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｓａｎｄｕｓｅｓｓｉｍｐｌｅｒａ…     

A hybrid RANS/LES simulation of high-Reynolds-number channel flow using additional filtering at the interface

A hybrid method combining large eddy simulation (LES) with the Reynolds-averaged Navier-Stokes (RANS) equation is used to simulate a turbulent channel flow at high Reynolds number. It is known that the mean velocity profile has a mismatch between the RANS and LES regions in hybrid simulations of a channel flow. The velocity mismatch is reproduced and its dependence on the location of the RANS/LES interface and on the type of RANS model is examined in order to better understand its properties. To remove the mismatch and to obtain better velocity profiles, additional filtering is applied to the velocity components in the wall-parallel planes near the interface. The additional filtering was previously introduced to simulate a channel flow at low Reynolds number. It is shown that the filtering is effective in reducing the mismatch even at high Reynolds number. Profiles of the velocity fluctuations of runs with and without the additional filtering are examined to help understand the reason for the mismatch. Due to the additional filtering, the wall-normal velocity fluctuation increases at the bottom of the LES region. The resulting velocity field creates the grid-scale shear stress more efficiently, and an overestimate of the velocity gradient is removed. The dependence of the velocity profile on the grid point number is also investigated. It is found that the velocity gradient in the core region is underestimated in the case of a coarse grid. Attention should be paid not only to the velocity mismatch near the interface but also to the velocity profile in the core region in hybrid simulations of a channel flow at high Reynolds number. PACS47.27.Eq; 47.27.Nz; 47.60.+i     

Near‐wall turbulence modelling with enhanced dissipation

A low‐Reynolds number k–ε turbulence model is proposed that incorporates diffusion terms and modified C_ε(1,2) coefficients to amplify the level of dissipation in non‐equilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. Unlike the conventional k–ε model, it requires no wall function/distance parameter that bridges the near‐wall integration. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Copyright © 2003 John Wiley & Sons, Ltd.     

An eddy viscosity model with near‐wall modifications

An extended version of the isotropic k–ε model is proposed that accounts for the distinct effects of low‐Reynolds number (LRN) and wall proximity. It incorporates a near‐wall correction term to amplify the level of dissipation in nonequilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and the Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Comparisons indicate that the present model is a significant improvement over the standard eddy viscosity formulation. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Bulk flow parameters in turbulent wall boundary layers near separation

Experiments were conducted in a turbulent boundary layer near separation along a flat plate. The pressure gradient in flow direction was varied such that three significant boundary layer configurations could be maintained. The flow in the test section thus had simultaneously a region of favourable pressure gradient, a region of strong adverse pressure gradient with boundary layer separation and a region of reattached boundary layer. Specially designed fine probes facilitated the measurements of skin friction and velocity distribution very close to the wall. Bulk flow parameters such as skin friction coefficient C _f, Reynold's number Re_δ2 and shape factors H and G, which are significant characteristics of wall boundary layers were evaluated. The dependence of these parameters on the Reynolds number and along the test section was explored and the values were compared with other empirical and analytical formulae known in the literature.     

Experimental Investigation of the Log-Law for an Adverse Pressure Gradient Turbulent Boundary Layer Flow at Re θ = 10000

We present an experimental investigation of a turbulent boundary layer flow at a significant adverse pressure gradient at Reynolds number Re _θ ?=?10000 using large field PIV. The testcase is designed to start from a zero pressure gradient flow at Re _θ ?=?8000 with a distinct log-law region following a slowly rising adverse pressure gradient. This allows to reveal a breakdown of the log-law under the effect of the adverse pressure gradient. The region described by the log-law is progressively reduced in terms of y ^?+? and then joins into a modified log-law which gives a good fit to the data up to at least y/δ ₉₉?≈?0.2. The scaling in the overlap region is demonstrated using the mean velocity slope diagnostic function, enabled due to the high quality of the PIV data. Locally, the velocity profile is measured down to the wall using long-range microscopic PIV with particle tracking velocimetry to determine the wall shear stress directly in the adverse pressure gradient region.     

Direct Numerical Simulation and Theory of a Wall-Bounded Flow with Zero Skin Friction

We study turbulent plane Couette-Poiseuille (CP) flows in which the conditions (relative wall velocity ΔU _w ≡ 2U _w , pressure gradient dP/dx and viscosity ν) are adjusted to produce zero mean skin friction on one of the walls, denoted by APG for adverse pressure gradient. The other wall, FPG for favorable pressure gradient, provides the friction velocity u _τ , and h is the half-height of the channel. This leads to a one-parameter family of one-dimensional flows of varying Reynolds number Re ≡ U _w h/ν. We apply three codes, and cover three Reynolds numbers stepping by a factor of two each time. The agreement between codes is very good, and the Reynolds-number range is sizable. The theoretical questions revolve around Reynolds-number independence in both the core region (free of local viscous effects) and the two wall regions. The core region follows Townsend’s hypothesis of universal behavior for the velocity and shear stress, when they are normalized with u _τ and h; on the other hand universality is not observed for all the Reynolds stresses, any more than it is in Poiseuille flow or boundary layers. The FPG wall region obeys the classical law of the wall, again for velocity and shear stress. For the APG wall region, Stratford conjectured universal behavior when normalized with the pressure gradient, leading to a square-root law for the velocity. The literature, also covering other flows with zero skin friction, is ambiguous. Our results are very consistent with both of Stratford’s conjectures, suggesting that at least in this idealized flow turbulence theory is successful like it was for the classical logarithmic law of the wall. We appear to know the constants of the law within a 10% bracket. On the other hand, that again does not extend to Reynolds stresses other than the shear stress, but these stresses are passive in the momentum equation.     

Mixed convective flow of a hybrid nanofluid between two parallel inclined plates under wall-slip condition

The mixed convection flow of a hybrid nanofluid in an inclined channel with top wall-slip due to wall stripe and constant heat flux conditions is studied. Explicit analytical solutions are given to the flow velocity, temperature, as well as the pressure in non-dimensional forms. The flow regime domain, the velocity and temperature distributions, and the dependence of various physical parameters such as the hybrid nanoparticle volume fractions, the wall-slip, the Grashof number, the Reynolds numb...     

A numerical study of turbulent flow over a two‐dimensional hill

Calculations of mean velocities and Reynolds stresses are reported for the recirculating flow established in the wake of two‐dimensional polynomial‐shaped obstacles that are symmetrical about a vertical axis and mounted in the water channel downstream of a fully developed channel flow for Re=6×10⁴. The study involves calculations of mean and fluctuating flow properties in the streamwise and spanwise directions and include comparisons with experimental data [Almeida GP, Durão DFG, Heitor MV. Wake flows behind two‐dimensional model hills. Experimental Thermal and Fluid Science 1993; 7: 87–101] for flow around a single obstacle with data resulting from the interaction of consecutive obstacles, using two versions of the low‐Reynolds number differential second‐moment (DSM) closure model. The results include analysis of the turbulent stresses in local flow co‐ordinates and reveal flow structure qualitatively similar to that found in other turbulent flows with a reattachment zone. It is found that the standard isotropization of production model (IPM), based on that proposed by Gibson and Launder [Ground effects on pressure fluctuations in the atmospheric boundary layer. Journal of Fluid Mechanics 1978; 86(3): 191–511], with the incorporation of the wall reflection model of Craft and Launder [New wall‐reflection model applied to the turbulent impinging jet. AIAA Journal 1992; 32(12): 2970–2972] predicts the mean velocities quite well, but underestimates the size of the recirculation region and turbulent quantities in the shear layer. These inadequacies are circumvented by adopting a new cubic Reynolds stress closure scheme based on that more recently developed by Craft and Launder [A Reynolds stress closure designed for complex geometries. International Journal of Heat and Fluid Flow 1996; 17: 245–254] which satisfies the two component limit (TCL) of turbulence. In this model the geometry‐specific quantities, such as the wall‐normal vector or wall distance, are replaced by invariant dimensionless gradient indicators. Also, the model captures the diverse behaviour of the different components of the stress dissipation, ε_ij, near the wall and uses a novel decomposition for the fluctuating pressure terms. Copyright © 2001 John Wiley & Sons, Ltd.     

Near-Wall Stress Balance in Front of a Wall-Mounted Cylinder

The stress balance in the near-wall flow in front of a cylinder mounted on a flat plate at moderate Reynolds number is investigated by applying highly resolved Large-Eddy Simulation (LES). The flow around wall-mounted bluff bodies is subject of research due to its wide relevance for engineering applications. However, the structure of the vortex system in front of such a bluff body is complex, bears strong velocity and pressure gradients in each spatial direction and has rich dynamics. Furthermore, the vortex system is located close to the investigated flat bottom wall (Dargahi, Exp. Fluids 8(1-2):1–12, 1989; Devenport and Simpson, J. Fluid Mech. 210:23–55, 1990). Thus, classical models for the treatment of the near-wall flow based on the logarithmic law of the wall or a power law cannot be expected to suffice in such kind of flow (Pope 2011). This paper assesses which contributors to the stress balance have significant influence on the balances residual and thus have to be considered by an approach to model the investigated near-wall flow. To do so, the momentum equation in streamwise direction is integrated in wall-normal direction and applied to the results gained from the LES. The evaluation of the stress balance along four selected wall-normal profiles indicates that the significance of each single term depends on where the profile is located. Outside the viscous layer, no term except the viscous stresses can be neglected in general. The amplitude of the pressure gradient as well as horizontal gradients of mean and fluctuating velocity are multiples of the estimated wall shear stress. Wall models not including a spatial approach are therefore most likely to fail in such kind of flow.     

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.     

Numerical simulation of turbulent impinging jet on a rotating disk

The calculations of quasi‐three‐dimensional momentum equations were carried out to study the influence of wall rotation on the characteristics of an impinging jet. The pressure coefficient, the mean velocity distributions and the components of Reynolds stress are calculated. The flow is assumed to be steady, incompressible and turbulent. The finite volume scheme is used to solve the continuity equation, momentum equations and k–ε model equations. The flow characteristics were studied by varying rotation speed ω for 0?ω?167.6 rad/s, the distance from nozzle to disk (H/d) was (3, 5, 8 and 10) and the Reynolds number Re base on V_J and d was 1.45 × 10⁴. The results showed that, the radial velocity and turbulence intensity increase by increasing the rotation speed and decrease in the impingement zone as nozzle to disk spacing increases. When the centrifugal force increases, the radial normal stresses and shear stresses increase. The location of maximum radial velocity decreases as the local velocity ratio (α) increases. The pressure coefficient depends on the centrifugal force and it decreases as the distance from nozzle to plate increases. In impingement zone and radial wall jet, the spread of flow increases as the angular velocity decreases The numerical results give good agreement with the experiment data of Minagawa and Obi (Int. J. of Heat and Fluid Flow 2004; 25 :759–766). Copyright © 2006 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.     

Experimental Investigation on Turbulent Structure of Drag Reducing Channel Flow with Blowing Polymer Solution from the Wall

Turbulent drag reducing flow with blowing polymer solution from the channel wall was investigated experimentally using particle image velocimetry (PIV). Experiments were carried out with varying conditions of blowing polymer solution (e.g. weight concentration of polymer solution). Reynolds number based on the channel height and mean velocity was set to 20000 and 40000. When the polymer solution was blown from the channel wall, streamwise velocity fluctuation little increased, but wall-normal velocity fluctuation, Reynolds shear stress and correlation coefficient decreased significantly only near the blower wall. This behavior corresponds to the decrease of the ejection and sweep in the near-wall region observed by the investigation of instantaneous velocity map. On the contrary, this characteristic behavior was not observed at a position away from the blower wall (y/(H/2) > 0.4) and the scatter plot was almost the same as that of the water flow in this region. These results suggest that there are two regions in the drag reducing flow with blowing polymer solution from the wall; one is a non-Newtonian region which exists near the blower wall, and the other is a Newtonian region at a distance from the wall. The non-Newtonian region plays a key role in the drag reduction by the blowing polymer solution.     

Analysis of flow in the plate‐spiral of a reaction turbine using a streamline upwind Petrov–Galerkin method

The prediction of the flow field in a novel spiral casing has been accomplished. Hydraulic turbine manufacturers are considering the potential of using a special type of spiral casing because of the easier manufacturing process involved in its fabrication. These special spiral casings are known as plate‐spirals. Numerical simulation of complex three‐dimensional flow through such spiral casings has been accomplished using a finite element method (FEM). An explicit Eulerian velocity correction scheme has been deployed to solve the Reynolds‐average Navier–Stokes equations. The simulation has been performed to describe the flow in high Reynolds number (10⁶) regimes. For spatial discretization, a streamline upwind Petrov–Galerkin (SUPG) technique has been used. The velocity field and the pressure distribution inside the spiral casing reveal meaningful results. Copyright © 2000 John Wiley & Sons, Ltd.