Turbulence in a skewed three‐dimensional wall‐bounded shear flow: effect of mean vorticity on structure modification

Mean‐flow three‐dimensionalities affect both the turbulence level and the coherent flow structures in wall‐bounded shear flows. A tailor‐made flow configuration was designed to enable a thorough investigation of moderately and severely skewed channel flows. A unidirectional shear‐driven plane Couette flow was skewed by means of an imposed spanwise pressure gradient. Three different cases with 8°, 34°and 52°skewing were simulated numerically and the results compared with data from a purely two‐dimensional plane Couette flow. The resulting three‐dimensional flow field became statistically stationary and homogeneous in the streamwise and spanwise directions while the mean velocity vector V and the mean vorticity vector Ω remained parallel with the walls. Mean flow profiles were presented together with all components of the Reynolds stress tensor. The mean shear rate in the core region gradually increased with increasing skewing whereas the velocity fluctuations were enhanced in the spanwise direction and reduced in the streamwise direction. The Reynolds shear stress is known to be closely related to the coherent flow structures in the near‐wall region. The instantaneous and ensemble‐averaged flow structures were turned by the skewed mean flow. We demonstrated for the medium‐skewed case that the coherent structures should be examined in a coordinate system aligned with V to enable a sound interpretation of 3D effects. The conventional symmetry between Case 1 and Case 2 vortices was broken and Case 1 vortices turned out to be stronger than Case 2. This observation is in conflict with the common understanding on the basis of the spanwise (secondary) mean shear rate. A refined model was proposed to interpret the structure modifications in three‐dimensional wall‐flows. What matters is the orientation of the mean vorticity vector Ω relative to the vortex vorticity vector ω _v, that is, the sign of Ω · ω _v. In the present situation, Ω · ω _v > 0 for the Case 1 vortices causing a strengthening relative to the Case 2 vortices. Copyright © 2011 John Wiley & Sons, Ltd.     

Turbulence in plane Couette flow

Experimental results for fully developed turbulent plane Couette flow are reported and compared to earlier experimental and numerical results. In addition some turbulent statistics not previously reported are shown.     

Similar solution for three‐dimensional flow in an Oldroyd‐B fluid over a stretching surface

The present investigation deals with the three‐dimensional flow of an Oldroyd‐B fluid over a stretching surface. The governing equations for the three‐dimensional flow are developed. Similarity transformations are invoked for the conversion of nonlinear partial differential equations into the coupled system of ordinary differential equations. Computations for the series solution are presented through implementation of homotopy analysis method. The salient features of the involved parameters have been pointed out. Copyright © 2011 John Wiley & Sons, Ltd.     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.     

Performance of very‐high‐order upwind schemes for DNS of compressible wall‐turbulence

The purpose of the present paper is to evaluate very‐high‐order upwind schemes for the direct numerical simulation (DNS ) of compressible wall‐turbulence. We study upwind‐biased (UW ) and weighted essentially nonoscillatory (WENO ) schemes of increasingly higher order‐of‐accuracy (J. Comp. Phys. 2000; 160 :405–452), extended up to WENO 17 (AIAA Paper 2009‐1612, 2009). Analysis of the advection–diffusion equation, both as Δx→0 (consistency), and for fixed finite cell‐Reynolds‐number Re_Δx (grid‐resolution), indicates that the very‐high‐order upwind schemes have satisfactory resolution in terms of points‐per‐wavelength (PPW ). Computational results for compressible channel flow (Re∈[180, 230]; M?_CL ∈[0.35, 1.5]) are examined to assess the influence of the spatial order of accuracy and the computational grid‐resolution on predicted turbulence statistics, by comparison with existing compressible and incompressible DNS databases. Despite the use of baseline O(Δt²) time‐integration and O(Δx²) discretization of the viscous terms, comparative studies of various orders‐of‐accuracy for the convective terms demonstrate that very‐high‐order upwind schemes can reproduce all the DNS details obtained by pseudospectral schemes, on computational grids of only slightly higher density. Copyright © 2009 John Wiley & Sons, Ltd.     

Low Reynolds number k–ε model for near‐wall flow

A wall‐distance free k–ε turbulence model is developed that accounts for the near‐wall and low Reynolds number effects emanating from the physical requirements. The model coefficients/functions depend non‐linearly on both the strain rate and vorticity invariants. Included diffusion terms and modified C_ε(1,2) coefficients 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. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

Unsteady three‐dimensional boundary layer flow due to a stretching surface in a micropolar fluid

In this paper, the unsteady three‐dimensional boundary layer flow due to a stretching surface in a viscous and incompressible micropolar fluid is considered. The partial differential equations governing the unsteady laminar boundary layer flow are solved numerically using an implicit finite‐difference scheme. The numerical solutions are obtained which are uniformly valid for all dimensionless time from initial unsteady‐state flow to final steady‐state flow in the whole spatial region. The equations for the initial unsteady‐state flow are also solved analytically. It is found that there is a smooth transition from the small‐time solution to the large‐time solution. The features of the flow for different values of the governing parameters are analyzed and discussed. The solutions of interest for the skin friction coefficient with various values of the stretching parameter c and material parameter K are presented. Copyright © 2011 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

DNS of backward‐facing step flow with fully turbulent inflow

Direct numerical simulation (DNS) has been performed to study the channel flow over a backward‐facing step at a Reynolds number Re_b=5600 based on the step height h and the inflow bulk velocity U_b. A dynamic method has been used in order to generate realistic turbulent inflow conditions. The results upstream of the step compared well with the fully developed channel flow. Downstream of the step our results show excellent agreement with experimental data. Copyright © 2009 John Wiley & Sons, Ltd.     

Start‐up flow in a three‐dimensional lid‐driven cavity by means of a massively parallel direction splitting algorithm

The purpose of this paper is to validate a new highly parallelizable direction splitting algorithm. The parallelization capabilities of this algorithm are illustrated by providing a highly accurate solution for the start‐up flow in a three‐dimensional impulsively started lid‐driven cavity of aspect ratio 1 × 1 × 2 at Reynolds numbers 1000 and 5000. The computations are done in parallel (up to 1024 processors) on adapted grids of up to 2 billion nodes in three space dimensions. Velocity profiles are given at dimensionless times t = 4, 8, and 12; at least four digits are expected to be correct at Re = 1000. Copyright © 2011 John Wiley & Sons, Ltd.     

Fringe element reconstruction for front tracking for three‐dimensional incompressible flow analysis

Fringe element reconstruction technique for tracking the free surface in three‐dimensional incompressible flow analysis was developed. The flow field was calculated by the mixed formulation based on a four‐node tetrahedral element with a bubble function at the centroid (P1+/P1). Since an Eulerian approach was employed in this study, the flow front interface was advected by the flow through a fixed mesh. For accurate modelling of interfacial movement, a fringe element reconstruction method developed can provide not only an accurate treatment of material discontinuity but also surface tension across the interface. The effect of surface tension was modelled by imposing tensile stress directly on the constructed surface elements at the flow front interface. To verify the numerical approach developed, the developed algorithm was applied to two examples whose solutions are available in references. Good agreement was obtained between the simulation results and these solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

Direct Numerical Simulation of Turbulence in a Stably Stratified Fluid and Wave-Shear Interaction

Turbulence decay in a strongly stratified medium is simulated by a direct pseudo-spectral code solving the three-dimensional equations of motion under the Boussinesq approximation. The results are compared to non-stratified simulations results. We focus on the production of mean shear energy observed in the stratified case. We then simulate the decay of stratified turbulence when affected by an initial horizontal mean flow and show that this mean flow is the major component remaining at large t. Next, we give some analytical elements on wave-shear interaction by using a simple refraction calculation with WKB hypothesis. This calculation is illustrated by simulating the interaction between one monochromatic internal wave and a vertical shear profile. We conclude that the existence of singularities in the mean shear production term in the presence of internal gravity waves may be one of the possible mechanisms involved within stratified turbulent shear flows.     

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.     

Finite element solution of three‐dimensional turbulent flows applied to mold‐filling problems

This paper presents a finite element solution algorithm for three‐dimensional isothermal turbulent flows for mold‐filling applications. The problems of interest present unusual challenges for both the physical modelling and the solution algorithm. High‐Reynolds number transient turbulent flows with free surfaces have to be computed on complex three‐dimensional geometries. In this work, a segregated algorithm is used to solve the Navier–Stokes, turbulence and front‐tracking equations. The streamline–upwind/Petrov–Galerkin method is used to obtain stable solutions to convection‐dominated problems. Turbulence is modelled using either a one‐equation turbulence model or the κ–ε two‐equation model with wall functions. Turbulence equations are solved for the natural logarithm of the turbulence variables. The change of dependent variables allows for a robust solution algorithm and good predictions even on coarse meshes. This is very important in the case of large three‐dimensional applications for which highly refined meshes result in untreatable large numbers of elements. The position of the flow front in the mold cavity is computed using a level set approach. Finally, equations are integrated in time using an implicit Euler scheme. The methodology presents the robustness and cost effectiveness needed to tackle complex industrial applications. Copyright © 2000 John Wiley & Sons, Ltd.     

Object‐oriented parallel algorithms for computing three‐dimensional isopycnal flow

In this paper, we derive an object‐oriented parallel algorithm for three‐dimensional isopycnal flow simulations. The matrix formulation is central to the algorithm. It enables us to apply an efficient preconditioned conjugate gradient linear solver for the global system of equations, and leads naturally to an object‐oriented data structure design and parallel implementation. We discuss as well, in less detail, a similar algorithm based on the reduced system, suitable also for parallel computation. Favorable performances are observed on test problems. Copyright © 2002 John Wiley & Sons, Ltd.     

A hybrid vortex method for the simulation of three‐dimensional flows

This paper presents an integral vorticity method for solving three‐dimensional Navier–Stokes equations. A finite volume scheme is implemented to solve the vorticity transport equation, which is discretized on a structured hexahedral mesh. A vortex sheet algorithm is used to enforce the no‐slip boundary condition through a vorticity flux at the boundary. The Biot–Savart integral is evaluated to compute the velocity field, in conjunction with a fast algorithm based on multipole expansion. This method is applied to the simulation of uniform flow past a sphere. Copyright © 2007 John Wiley & Sons, Ltd.     

Direct numerical simulation of decaying two‐dimensional turbulence in a no‐slip square box using smoothed particle hydrodynamics

This paper explores the application of SPH to a DNS of decaying turbulence in a two‐dimensional no‐slip wall‐bounded domain. In this bounded domain, the inverse energy cascade, and a net torque exerted by the boundary, results in a spontaneous spin‐up of the fluid, leading to a typical end state of a large monopole vortex that fills the domain. The SPH simulations were compared against published results using a high‐accuracy pseudo‐spectral code. Ensemble averages of the kinetic energy, enstrophy and average vortex wavenumber compared well against the pseudo‐spectral results, as did the evolution of the total angular momentum of the fluid. However, although the pseudo‐spectral results emphasised the importance of the no‐slip boundaries as generators of long‐lived coherent vortices in the flow, no such generation was seen in the SPH results. Vorticity filaments produced at the boundary were always dissipated by the flow shortly after separating from the boundary layer. The kinetic energy spectrum of the SPH results was calculated using an SPH Fourier transform that operates directly on the disordered particles. The ensemble kinetic energy spectrum showed the expected k^?3 scaling over most of the inertial range. However, the spectrum flattened at smaller length scales (initially less than 7.5 particle spacings and growing in size over time), indicating an excess of small‐scale kinetic energy.Copyright © 2011 John Wiley & Sons, Ltd.     

Vortex particle‐in‐cell method for three‐dimensional viscous unbounded flow computations

A new vortex particle‐in‐cell (PIC) method is developed for the computation of three‐dimensional unsteady, incompressible viscous flow in an unbounded domain. The method combines the advantages of the Lagrangian particle methods for convection and the use of an Eulerian grid to compute the diffusion and vortex stretching. The velocity boundary conditions used in the method are of Dirichlet‐type, and can be calculated using the vorticity field on the grid by the Biot–Savart equation. The present results for the propagation speed of the single vortex ring are in good agreement with the Saffman's model. The applications of the method to the head‐on and head‐off collisions of the two vortex rings show good agreement with the experimental and numerical literature. Copyright © 2000 John Wiley & Sons, Ltd.