Numerical prediction of turbulent flow structure generated by a synthetic cross‐jet into a turbulent boundary layer

A detailed numerical study using large‐eddy simulation (LES) and unsteady Reynolds‐averaged Navier–Stokes (URANS) was undertaken to investigate physical processes that are engendered in the injection of a circular synthetic (zero‐net mass flux) jet in a zero pressure gradient turbulent boundary layer. A complementary study was carried out and was verified by comparisons with the available experimental data that were obtained at corresponding conditions with the aim of achieving an improved understanding of fluid dynamics of the studied processes. The computations were conducted by OpenFOAM C++, and the physical realism of the incoming turbulent boundary layer was secured by employing random field generation algorithm. The cavity was computed with a sinusoidal transpiration boundary condition on its floor. The results from URANS computation and LES were compared and described qualitatively and quantitatively. There is a particular interest for acquiring the turbulent structures from the present numerical data. The numerical methods can capture vortical structures including a hairpin (primary) vortex and secondary structures. However, the present computations confirmed that URANS and LES are capable of predicting current flow field with a more detailed structure presented by LES data as expected. Copyright © 2011 John Wiley & Sons, Ltd.     

On the possible formation of a vortex flow in a supersonic inlet with three-dimensional compression

The structure of the flow in an inlet with three-dimensional compression and an adjoining channel of rectangular cross-section is experimentally investigated for the external flow parameters M_∞ = 4 and Re ≈ 52 · 10⁶ 1/m. The phenomenon of vortex generation in the air-intake and its breakdown at the channel entry, where it encounters an elevated pressure gradient, is studied     

Three‐dimensional vortex simulation of unsteady flow in hydraulic turbines

On the basis of the Helmholtz decomposition, a grid‐free numerical scheme is provided for the solution of unsteady flow in hydraulic turbines. The Lagrangian vortex method is utilized to evaluate the convection and stretch of the vorticity, and the BEM is used to solve the Neumann problem to define the potential flow. The no‐slip boundary condition is satisfied by generating vortex sticks at the solid surface. A semi‐analytical regularization technique is applied to evaluate the singular boundary surface integrals of the potential velocity and its gradients accurately. The fast multipole method was extended to evaluate the velocity and velocity gradients induced by the discretized vortex blobs in the Lagrangian vortex method. The successful simulation for the unsteady flow through a hydraulic turbine's runner has manifested the effectiveness of the proposed method. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical investigation of flow over a two-dimensional square cylinder with a synthetic jet generated by a bi-frequency signal

Applied Mathematics and Mechanics - The flow around a square cylinder with a synthetic jet positioned at the rear surface is numerically investigated with the unsteady Reynolds-averaged...     

Numerical simulation of magnetohydrodynamic flow in a toroidal duct of square cross-section

We present numerical simulation results of the quasi-static magnetohydrodynamic (MHD) flow in a toroidal duct of square cross-section with insulating Hartmann walls and conducting side walls. Both laminar and turbulent flows are considered. In the case of steady flows, we present a comprehensive analysis of the secondary flow. It consists of two counter-rotating vortex cells, with additional side wall vortices emerging at sufficiently high Hartmann number. Our results agree well with existing asymptotic analysis. In the turbulent regime, we make a comparison between hydrodynamic and MHD flows. We find that the curvature induces an asymmetry between the inner and outer side of the duct, with higher turbulence intensities occurring at the outer side wall. The magnetic field is seen to stabilize the flow so that only the outer side layer remains unstable. These features are illustrated both by a study of statistically averaged quantities and by a visualization of (instantaneous) coherent vortices.     

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.     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

A three‐dimensional vortex particle‐in‐cell method for vortex motions in the vicinity of a wall

A new vortex particle‐in‐cell method for the simulation of three‐dimensional unsteady incompressible viscous flow is presented. The projection of the vortex strengths onto the mesh is based on volume interpolation. The convection of vorticity is treated as a Lagrangian move operation but one where the velocity of each particle is interpolated from an Eulerian mesh solution of velocity–Poisson equations. The change in vorticity due to diffusion is also computed on the Eulerian mesh and projected back to the particles. Where diffusive fluxes cause vorticity to enter a cell not already containing any particles new particles are created. The surface vorticity and the cancellation of tangential velocity at the plate are related by the Neumann conditions. The basic framework for implementation of the procedure is also introduced where the solution update comprises a sequence of two fractional steps. The method is applied to a problem where an unsteady boundary layer develops under the impact of a vortex ring and comparison is made with the experimental and numerical literature. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of quasi-streamwise hairpin-like vortex generation in turbulent boundary layer

A mechanism for generation of near wall quasi-streamwise hairpin-like vortex （QHV） and secondary quasi-streamwise vortices （SQV） is presented. The conceptual model of resonant triad in the theory of hydrodynamic instability and direct numerical simulation of a turbulent boundary layer were applied to reveal the formation of QHV and SQV. The generation procedures and the characteristics of the vortex structures are obtained, which share some similarities with previous numerical simulations. The research using resonant triad conceptual model and numerical simulation provides a possibility for investigating and controling the vortex structures, which play a dominant role in the evolution of coherent structures in the near-wall region.     

Suppression of the unsteady vortex wakes of a circular cylinder pair by a doublet‐like counter‐rotation

The present investigation examines the suppression of unsteady, two‐dimensional wake instabilities of a pair of identical circular cylinders, placed side‐by‐side normal to freestream at a low Reynolds number of 150. It is found that when the cylinders are counter‐rotated, unsteady vortex wakes can be completely suppressed. At fast enough rotational speeds, a virtual elliptic body is produced by a closed streamline, strongly resembling a doublet potential flow. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical analysis of the developing fluid flow in a circular duct rotating steadily about a parallel axis

A numerical analysis of the flow pattern in the inlet region of a circular pipe rotating steadily about an axis parallel to its own is presented. Both finite cell and finite element methods are used to analyse the problem and they give qualitatively similar results which show that a swirling fluid motion is induced in the pipe inlet region. The analyses show that the direction of swirl is opposite to that of the pipe rotation when viewed along the flow axis and that its magnitude depends on the speed of pipe rotation and throughflow Reynolds number. Neither numerical analysis predicts the marked upturn in friction factor (or pressure drop) which has been observed experimentally. However, a dependence on the pipe inlet boundary conditions is demonstrated.     

微型射流涡流器对湍流边界层控制的数值研究

运用数值方法,模拟出展向分布的同向倾斜微型射流列与平板湍流边界层相互作用形成流向涡列的流场结构,验证了利用其来对湍流边界层进行控制的可能性.随射流间距减小,流向涡列控制作用流向渗透能力增强,但作用区域减小;随射流速度提高,流向涡列控制作用增强,但过大的射流速度反而会导致流向涡列在局部区域内控制作用的下降;随射流俯仰角减小、倾斜角增大,流向涡列初始控制作用增强,但过小的俯仰角、过大的倾斜角会导致流向涡列流向控制区域明显缩小.要保证流向涡列具有较强的湍流边界层控制作用,必须通过合理配置射流列各主要参数,在保证各流向涡具有一定强度的同时,还要确保各流向涡在形成时部分嵌入边界层内部.     

Study of a plane‐free jet exhausting from a channel by vortex‐in‐cell method

A two‐dimensional simulation of a plane jet exhausting from a channel has been performed using the vortex in cell algorithm in the Reynolds number range of 100–900. The vorticity is generated on the wall of the entrance channel whose length has been fixed in order to obtain a fully developed velocity profile at the entry of the jet. The transient behaviour of the velocity field starting from rest has been observed until reaching a quasi steady regime. The mean value of the velocity field is compared with the results of a finite volume computation on the same mesh. The velocity fluctuations obtained using this method are analysed. Their effect on the mean flow is estimated to be smaller than the viscous effect. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation of an isolated vortex and shear flow interaction

A numerical solution for the Navier-Stokes equations in the unbounded region is considered for the interaction of an isolated vortex and shear flow. A Chebyshey collocation method in space and finite-difference method for temporal discretization are used. The results of the numerical experiments for the interaction are discussed. It is shown that shear flow can both increase and decrease the vortex dissipation rate.     

Numerical simulation of jet-forced flow in a circular reservoir using discrete and random vortex methods

This paper describes a Biot–Savart discrete vortex model for simulating the flow patterns which occur when a single high-velocity inflow jet is used to stir the fluid within a circular container. The first stage of the model consists of conformally mapping the circular perimeter of the container onto a rectangle by means of a Schwarz–Christoffel transformation. A potential flow solution is then obtained for the flow inside the rectangle and this is transformed to give the potential flow inside the circle. In the second stage of the simulation, discrete vortices are added at the inlet of the physical system in order to model the inflow shear layers. Velocity components resulting from the discrete vortices and their images in the walls of the cylinder are superimposed on the uniform potential flow solution. The positions of the vortices are updated using a Lagrangian tracking procedure. Viscous effects are incorporated through the use of random walks. From the results it is shown that the discrete vortex method does predict qualitatively the important features of jet-forced reservoir flow.     

The solution of viscous incompressible jet flows using non‐staggered boundary fitted co‐ordinate methods

A new approach for the solution of the steady incompressible Navier–Stokes equations in a domain bounded in part by a free surface is presented. The procedure is based on the finite difference technique, with the non‐staggered grid fractional step method used to solve the flow equations written in terms of primitive variables. The physical domain is transformed to a rectangle by means of a numerical mapping technique. In order to design an effective free solution scheme, we distinguish between flows dominated by surface tension and those dominated by inertia and viscosity. When the surface tension effect is insignificant we used the kinematic condition to update the surface; whereas, in the opposite case, we used the normal stress condition to obtain the free surface boundary. Results obtained with the improved boundary conditions for a plane Newtonian jet are found to compare well with the available two‐dimensional numerical solutions for Reynolds numbers, up to Re=100, and Capillary numbers in the range of 0≤Ca<1000. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of trajectory and deformation of bubble in tip vortex

According to the behaviors of a bubble in the ship wake flow,the numerical simulation is divided into two stages,quasi-spherical motion and non-spherical motion,based on whether the bubble is captured by the vortex or not.The one-way coupled particle tracking method(PTM) and the boundary element method(BEM) are adopted to simulate these two stages,respectively.Meanwhile,the initial condition of the second stage is taken as the output of the first one,and the entire simulation is connected and completed.Based on the numerical results and the published experimental data,the cavitation inception is studied,and the wake bubble is tracked.Besides,the split of the bubble captured by the vortex and the following sub-bubbles are simulated,including motion,deformation,and collapse.The results provide some insights into the control on wake bubbles and optimization of the wake flow.     

Interaction of a heavy fluid flow in a channel with a transverse jet barrier

An experimental and numerical analysis of the interaction between a plane horizontal water flow in a rectangular channel (free water current) and a plane thin water jet (water jet curtain) is presented; the jet flows out vertically from either a slot nozzle in the bottom of the channel or the crest of a rigid spillway at a velocity appreciably (several times) greater than the water velocity in the channel. Numerical calculations were carried out using the STAR-CD software package preliminarily tested against the experimental data obtained. The dependence of the water level in the channel at a certain distance ahead of the jet barrier on the main jet parameters and the water flow rate in the horizontal channel is studied. It is found that in the region of the interface between the flows both steady and unsteady (self-oscillatory) flow patterns can be realized. Steady stream/jet interaction patterns of the “ejection” and “ejection-spillway” types are distinguished and a criterion separating these regimes is obtained. The notion of a rigid spillway equivalent to a jet curtain is introduced and an approximate dependence of its height on the relevant parameters of the problem is derived. The possibility of effectively controlling the water level ahead of a rigid spillway with a sharp edge by means of a plane water jet flowing from its crest is investigated. The boundary of transition to self-oscillation interaction patterns in the region of the flow interface is determined. The structure of these flows and a possible mechanism of their generation are described. Within the framework of the inviscid incompressible fluid model in the approximate formulation for a “thin” jet, an analytical dependence of the greatest possible depth of a reservoir filled with a heavy fluid at rest and screened by a vertical jet barrier on the jet parameters is obtained.