Direct Numerical Simulation of Subsonic Round Turbulent Jet

Direct numerical simulation(DNS) of spatially developing round turbulent jet flow with Reynolds number 4,700 was carried out. Over 20 million grid points were used in this simulation. Fully compressible three-dimensional Navier–Stokes equations were solved. High order explicit spatial difference schemes and Runge–Kutta time integration scheme were used to calculate derivatives and time marching, respectively. Non-reflecting boundary conditions and exit zone techniques were adopted. Some refined computational grids were used in order to capture the smallest turbulent structures near the centerline of the jet. Low level disturbance were imposed on the jet inflow velocity to trigger the developing of turbulence. Turbulent statistics such as mean velocity, Reynolds stresses, third order velocity moments were obtained and compared with experimental data. One-dimensional velocity autospectra was also calculated. The inertial region where the spectra decays according to the k^− 5/3 was observed. The quantitative profiles of mean velocity and all of the third order velocity moments which were difficult to measure via experimental techniques were presented here in detail. The jet flow was proven to be close to fully self-similar around 19 jet diameters downstream of jet exit. The statistic data and revealed flow feature obtained in this paper can provide valuable reference for round turbulent jet research.     

孔板空调风口送风射流的数值模拟

介绍Ｎ点风口模型用于数值模拟室内空气流动时描述孔板类送风口的入流边界条件．然后采用该风口模型对不同的孔板风口出流条件算例进行数值计算，并就轴心速度衰减、射流扩展角以及断面流速分布等射流特性与实验数据进行了对比．比较结果表明，Ｎ点风口模型用于描述数值模拟室内空气流动的孔板类风口入流边界条件，可以获得工程上足够满意的结果．     

Large-eddy simulation of a turbulent forced plume

This paper reports on an application of large-eddy simulation (LES) to a spatially-developing round turbulent buoyant jet. The numerical method used is based on a low-Mach-number version of the governing equations for compressible flow which can account for density variations. The second-order centre-difference scheme is used for spatial discretization and an Adams–Bashforth scheme for temporal discretization. Comparisons are made between LES results, experimental measurements and plume theory for the forced plume under moderate Reynolds number and good agreement has been achieved. It is found that the plume spreading and the centerline maximum mean velocity strongly depend on the forcing conditions imposed on the inflow plane. The helical mode of instability leads to a larger spreading rate as compared to an axisymmetric mode. The enhanced entrainment is directly related to the strong turbulent momentum and energy transports between the plume and surrounding fluid induced by vortex dynamics. The entrainment ratio is about 0.09 and falls into the range of experimentally determined values. Budgets of the mean momentum and energy equations are analyzed. It is found that the radial turbulent transport nearly balances the streamwise convection and the buoyancy force in the axial momentum equation. Also, the radial turbulent stress is balanced by the streamwise convection in the energy equation. The energy-spectrum for the axial velocity fluctuations shows a −5/3 power law of the Kolmogorov decay, while the power spectrum for the temperature fluctuations shows both −5/3 and −3 power laws in the inertial-convective and inertial-diffusive ranges, respectively.     

Influence of the inlet velocity profiles on the prediction of velocity distribution inside an electrostatic precipitator

The influence of the velocity profile at the inlet boundary on the simulation of air velocity distribution inside an electrostatic precipitator is presented in this study. Measurements and simulations were performed in a duct and an electrostatic precipitator (ESP). A four-hole cobra probe was used for the measurement of velocity distribution. The flow simulation was performed by using the computational fluid dynamics (CFD) code FLUENT. Numerical calculations for the air flow were carried out by solving the Reynolds-averaged Navier–Stokes equations coupled with the realizable k-ε turbulence model equations. Simulations were performed with two different velocity profiles at the inlet boundary – one with a uniform (ideal) velocity profile and the other with a non-uniform (real) velocity profile to demonstrate the effect of velocity inlet boundary condition on the flow simulation results inside an ESP. The real velocity profile was obtained from the velocity measured at different points of the inlet boundary whereas the ideal velocity profile was obtained by calculating the mean value of the measured data. Simulation with the real velocity profile at the inlet boundary was found to predict better the velocity distribution inside the ESP suggesting that an experimentally measured velocity profile could be used as velocity inlet boundary condition for an accurate numerical simulation of the ESP.     

Large Eddy Simulation of the sound field of a round turbulent jet

In this paper we will use Large Eddy Simulation (LES) to obtain the flow field of a turbulent round jet at a Reynolds number based on the jet orifice velocity of 11000. In the simulations it is assumed that the flow field is incompressible. The acoustic field of the jet is calculated with help of the Lighthill acoustic analogy. The coupling between the flow solver and the acoustic solver is discussed in detail. The Mach number used in the acoustic calculation was equal to 0.6. It is shown that the decay of the jet centerline velocity and centerline rms are in good agreement with experimental data of [12]. Furthermore, it is shown that the influence of the LES modeling on the acoustic field is very small, if the dynamic subgrid model is used.     

The structure of particle-laden,underexpanded free jets

Underexpanded, supersonic gas-particle jets were experimentally studied using the shadowgraph technique in order to examine the influence of the dispersed particles on the shape of the free jet and the structure of the imbedded shock waves. The particle mass loading at the nozzle exit was varied between zero and one, and two sizes of particles (i.e. spherical glass beads) with mean number diameters of 26 and 45 m were used. It was found that the Mach-disc moves upstream towards the orifice with increasing particle loading. The laser light sheet technique was also used to visualize the particle concentration distribution within the particle jet and the spreading rate of the particle jet. Furthermore, the particle velocity along the jet centerline was measured with a modified laser-Doppler anemometer. These measurements revealed that the particles move considerably slower than the gas flow at the nozzle exit. This is mainly the result of the particle inertia, whereby the particles are not accelerated to sonic speed in the converging part of the nozzle.In order to further explore the particle behavior in the free jet, numerical studies were performed by a combined Eulerian/Lagrangian approach for the gas and particle phases, including full coupling between the two phases. The numerical results showed that the application of different particle velocities at the nozzle exit as the inlet conditions, which were below the sonic speed of the gas phase has a significant influence on the free jet shape and the configuration of the shock waves. These results demonstrate that the assumption of equilibrium flow (i.e. zero slip between the phases) at the nozzle exit which has been applied in most of the previous numerical studies is not justified in most cases. Furthermore, the numerical calculations of the free jet shape and the particle velocity along the jet axis were compared with the measurements. Although correlations for rarefaction and compressibility effects in the drag coefficient were taken into account, the particle velocity along the center line was considerably overpredicted.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Numerical investigation of turbulent free jet flows issuing from rectangular nozzles: the influence of small aspect ratio

In this research, the fluid and thermal characteristics of a rectangular turbulent jet flow is studied numerically. The results of three-dimensional jet issued from a rectangular nozzle are presented. A numerical method employing control volume approach with collocated grid arrangement was employed. Velocity and pressure fields are coupled with SIMPLEC algorithm. The turbulent stresses are approximated using k–e{\varepsilon} model with two different inlet conditions. The velocity and temperature fields are presented and the rates of their decay at the jet centerline are noted. The velocity vectors of the main flow and the secondary flow are illustrated. Also, effect of aspect ratio on mixing in rectangular cross-section jets is considered. The aspect ratios that were considered for this work were 1:1 to 1:4. The results showed that the jet entrains more with smaller AR. Special attention has been drawn to the influence of the Reynolds number (based on hydraulic diameter) as well as the inflow conditions on the evolution of the rectangular jet. An influence on the jet evolution is found for smaller Re, but the jet is close to a converged state for higher Reynolds numbers. The inflow conditions have considerable influence on the jet characteristics.     

Large eddy simulation of incompressible free round jet with discontinuous Galerkin method

In the paper, discontinuous Galerkin method is applied to simulation of incompressible free round turbulent jet using large eddy simulation with eddy viscosity approach. The solution algorithm is based on the classical projection method, but instead of the solution of the Poisson equation, a parabolic equation is advanced in pseudo‐time, which provides the pressure field ensuring the proper pressure–velocity coupling. For time and pseudo‐time integration, explicit Runge–Kutta method is employed. The computational meshes consist of hexahedral elements with flat faces. Within a given finite element, all flow variables are expressed with modal expansions of the same order (including velocity and pressure). Discretisation of the viscous terms in the Navier–Stokes equations and Laplacian in the Poisson equation is stabilised with mixed finite element approach. The correctness of the solution algorithm is verified in a commonly used test case of laminar flow in 3D lid‐driven cavity. The results of computations of the free jet are compared with experimental and numerical reference data, the latter obtained from the high‐order pseudospectral code. The statistics of centerline flow velocity – mean velocity and its fluctuations – show satisfactory agreement with the reference data. Copyright © 2015 John Wiley & Sons, Ltd.     

Reynolds number effects on the behavior of a non-buoyant round jet

The behavior of a non-buoyant circular water jet discharged from a contraction nozzle was experimentally investigated. In this experiment, the Reynolds number of the jet, based on the mean velocity results obtained by particle image velocimetry (PIV), ranged from 177 to 5,142. From the experimental results, we found that the cross-sectional profile of the axial velocity for a laminar flow near the nozzle did not show a top-hat distribution, whereas the profiles with Reynolds number higher than 437 were almost top-hat. The length of the zone of flow establishment (ZFE) was found to decrease with increasing Reynolds number. The measured centerline velocity decayed more rapidly and, consequently, approached the theoretical equation earlier near the nozzle as the Reynolds number increased. The decay constant for the centerline velocity of the turbulent cases was relatively lower than that discovered in theory. It is assumed that this probably resulted from the use of the contraction nozzle. Verifying the similarity of the lateral velocity profiles demonstrated that the Gaussian curve was properly approximated only for the turbulent jets and not for the laminar or transitional flows. The jet half width seldom grew for the laminar or transitional flows, whereas it grew with increasing axial distance for the turbulent flows. The spreading rates for the turbulent flows gradually decreased with increasing Reynolds number. The normalized turbulence intensity along the jet centerline increased more rapidly with the axial distance as the Reynolds number increased, and tended to the constant values proposed by previous investigators. The Reynolds shear stress levels were also found to increase as the Reynolds number increased for the turbulent jets.     

Moderate Reynolds number axisymmetric jet development downstream an extended conical diffuser: Influence of extension length

Low-velocity (bulk velocity of 4.4 m/s) and moderate Reynolds (7350) axisymmetrical jet development is studied by hot-film single sensor anemometry. The jet issues from a conical convergent-divergent diffuser with uniform extension (diameter 25 mm). Decreasing the length-to-diameter ratio of the extension tube from 20 down to 0.4 is shown to alter severely the mean velocity profile at the tube outlet from Blasius to top-hat whereas turbulence intensities increases from 6 up to 50%. Next, the influence of the initial velocity profile at the tube outlet on axisymmetrical jet development is assessed. The velocity development exhibits a self-similar far field characteristic for axisymmetrical jet development. Although, the jet centerline decay constant increases and the jet spreading rate decreases as length-to-diameter ratios are increased from 0.4 up to 7.2 for which the initial centerline velocity decreases. Therefore, scaling of the centerline decay constant and inverse scaling of the spreading rate with initial centerline velocity U₀ or initial velocity Reynolds number Re₀ reported for moderate Reynolds numbers and low initial turbulence levels [18,22] does not hold as the turbulence level exceeds a threshold value in the range 12–27%. In addition, the influence of initial conditions on near and far field turbulence properties is shown. A transition in near field behaviour is observed for length-to-diameters around 3.6. Flow and geometrical configurations under study are relevant to e.g. upper airway flow.     

Nanoscale Poiseuille flow and effects of modified Lennard–Jones potential function

Numerical simulation of Poiseuille flow of liquid Argon in a nanochannel using the non-equilibrium molecular dynamics simulation (NEMD) is performed. The nanochannel is a three-dimensional rectangular prism geometry where the concerned numbers of Argon atoms are 2,700, 2,550 and 2,400 at 102, 108 and 120 K. Poiseuille flow is simulated by embedding the fluid particles in a uniform force field. An external driving force, ranging from 1 to 11 PN (Pico Newton), is applied along the flow direction to inlet fluid particles during the simulation. To obtain a more uniform temperature distribution across the channel, local thermostating near the wall are used. Also, the effect of other mixing rules (Lorenthz–Berthelot and Waldman–Kugler rules) on the interface structure are examined by comparing the density profiles near the liquid/solid interfaces for wall temperatures 108 and 133 K for an external force of 7 PN. Using Kong and Waldman–Kugler rules, the molecules near the solid walls were more randomly distributed compared to Lorenthz–Berthelot rule. These mean that the attraction between solid–fluid atoms was weakened by using Kong rule and Waldman–Kugler rule rather than the Lorenthz–Berthelot rule. Also, results show that the mean axial velocity has symmetrical distribution near the channel centerline and an increase in external driving force can increase maximum and average velocity values of fluid. Furthermore, the slip length and slip velocity are functions of the driving forces and they show an arising trend with an increase in inlet driving force and no slip boundary condition is satisfied at very low external force (<1 PN).     

入口速度比对射流拟序结构的影响

为了深入了解湍流流动机理以及湍流拟序结构发现过程的影响因素，本文采用大涡模拟方法对不同入口射流伴流速度比的平面湍射流流动进行了数值模拟。采用分步投影法求解动量方程，亚格子项采用标准Smagorinsky亚格子模式模拟，压力泊松方程采用修正的循环消去法快速求解，空间方程采用二阶精度的差分格式，在时间方向上采用二阶精度的显式差分格式。模拟结果给出了平面射流中湍流拟序结构的瞬态发展演变过程，分析了入口速度比对射流拟序结构发展演化过程及宏观流场形态的影响。为进一步研究射流拟序结构及其在湍流流动中的作用提供了基础。     

自激振荡脉冲射流喷嘴装置系统频率特性理论研究

根据相似系统原理和流体网络理论建立了自激振荡脉冲射流喷嘴装置的等效网络模型,利用系统传递函数推导了系统频率特性方程并进行了数值计算。结果表明:喷嘴装置的固有频率主要由喷嘴形状、结构参数、入口流速、射流中压力扰动波波速决定;自激振荡腔腔径、自激振荡腔腔长、上喷嘴直径、下喷嘴直径都对系统频率特性影响很大。提出了相应的自激振荡脉冲射流喷嘴设计准则,即喷嘴装置在最佳阻尼比下产生谐波共振。     

On the effect of inflow conditions in simulation of a turbulent round jet

This paper investigates the impact of the inflow conditions on simulations of a round jet discharging from a wall into a large space. The fluid dynamic characteristics of a round jet are studied numerically. A numerical method based on the control volume approach with collocated grid arrangement is employed. The k-e{k-\varepsilon} model is utilized to approximate turbulent stresses by considering six different inlet conditions. The velocity field is presented, and the rate of decay at the jet centerline is determined. The results showed that inflow conditions had a strong influence on the jet characteristics. This paper also investigates both sharp-edged and contoured nozzles. The effects of velocity, turbulence intensity, turbulence kinetic energy, and turbulence dissipation rate on flow field characteristics are examined. Results showed that the present simulations in both types of nozzles are in good agreement with experiments when considering the appropriate inflow conditions.     

The advantages of using high‐order finite differences schemes in laminar–turbulent transition studies

This paper presents various finite difference schemes and compare their ability to simulate instability waves in a given flow field. The governing equations for two‐dimensional, incompressible flows were solved in vorticity–velocity formulation. Four different space discretization schemes were tested, namely, a second‐order central differences, a fourth‐order central differences, a fourth‐order compact scheme and a sixth‐order compact scheme. A classic fourth‐order Runge–Kutta scheme was used in time. The influence of grid refinement in the streamwise and wall normal directions were evaluated. The results were compared with linear stability theory for the evolution of small‐amplitude Tollmien–Schlichting waves in a plane Poiseuille flow. Both the amplification rate and the wavenumber were considered as verification parameters, showing the degree of dissipation and dispersion introduced by the different numerical schemes. The results confirmed that high‐order schemes are necessary for studying hydrodynamic instability problems by direct numerical simulation. Copyright © 2005 John Wiley & Sons, Ltd.     

One-dimensional turbulence modeling for cylindrical and spherical flows: model formulation and application

The one-dimensional turbulence (ODT) model resolves a full range of time and length scales and is computationally efficient. ODT has been applied to a wide range of complex multi-scale flows, such as turbulent combustion. Previous ODT comparisons to experimental data have focused mainly on planar flows. Applications to cylindrical flows, such as round jets, have been based on rough analogies, e.g., by exploiting the fortuitous consistency of the similarity scalings of temporally developing planar jets and spatially developing round jets. To obtain a more systematic treatment, a new formulation of the ODT model in cylindrical and spherical coordinates is presented here. The model is written in terms of a geometric factor so that planar, cylindrical, and spherical configurations are represented in the same way. Temporal and spatial versions of the model are presented. A Lagrangian finite-volume implementation is used with a dynamically adaptive mesh. The adaptive mesh facilitates the implementation of cylindrical and spherical versions of the triplet map, which is used to model turbulent advection (eddy events) in the one-dimensional flow coordinate. In cylindrical and spherical coordinates, geometric stretching of the three triplet map images occurs due to the radial dependence of volume, with the stretching being strongest near the centerline. Two triplet map variants, TMA and TMB, are presented. In TMA, the three map images have the same volume, but different radial segment lengths. In TMB, the three map images have the same radial segment lengths, but different segment volumes. Cylindrical results are presented for temporal pipe flow, a spatial nonreacting jet, and a spatial nonreacting jet flame. These results compare very well to direct numerical simulation for the pipe flow, and to experimental data for the jets. The nonreacting jet treatment overpredicts velocity fluctuations near the centerline, due to the geometric stretching of the triplet maps and its effect on the eddy event rate distribution. TMB performs better than TMA. A hybrid planar-TMB (PTMB) approach is also presented, which further improves the results. TMA, TMB, and PTMB are nearly identical in the pipe flow where the key dynamics occur near the wall away from the centerline. The jet flame illustrates effects of variable density and viscosity, including dilatational effects.     

Numerical study of unsteady flow and heat transfer of circular tangential direction jets flowing over the inner cylinder surface in the annular chamber

In order to understand the dynamics of vortices on heat transfer, the unsteady flow field of tangential direction jets flowing in the annular chamber is numerically investigated by scale-adaptive simulation (SAS). The jet Reynolds number is 332,000 based on the jet’s diameter and inflow velocity for a specific geometric model. The analogy theory is used to obtain the convective heat transfer coefficient distribution on the hub surface. Spectral analysis via fast Fourier transform (FFT) is used to analyze frequency information that flows inside the chamber. The proper orthogonal decomposition (POD) method is performed on the velocity field in the chamber and the convective heat transfer coefficient on the hub surface using a snapshot method. The fast Fourier transform helps find the dominant frequency of the unsteady flow in the chamber. The time sequence of velocity fields on the radial plane shows the presence of cyclic flapping of the jet. The proper orthogonal decomposition analysis indicates that the unsteady periodic flow phenomenon in the chamber and unsteady heat transfer on the hub surface are mainly related to the dynamics of the counter-rotating vortices caused by the jet.     

A study of the flow characteristics in micro-abrasive jets

An experimental study of particle velocities in micro-abrasive jets by using the particle image velocimetry (PIV) technique is presented. It has been found that the particle jet flow has a nearly linear expansion downstream. The particle velocities increase with air pressure, and the increasing rate increases with nozzle diameter within the range considered. The instantaneous velocity profile of the particle flow field in terms of the particle velocity distribution along the axial and radial directions of the jets is discussed. For the axial profile in the jet centerline downstream, there exists an extended acceleration stage, a transition stage, and a deceleration stage. For the radial velocity profiles, a relatively flat shape is observed at a jet cross-section near the nozzle exit. Mathematical models for the particle velocities in the air jet are then developed. It is shown that the results from the models agree well with experimental data in both the variation trend and magnitude.     

Three-dimensional numerical study of jet-in-crossflow characteristics at low Reynolds number

A numerical study of a square jet in a cross flow is carried out at a Reynolds number of 100. The flow field and heat transfer characteristic downstream of the jet have been explored by solving three-dimensional unsteady Navier–Stokes equations and energy equation using higher order spatial and temporal discretization. The projection of vortical structure on a plane is seen to give the component of vortex normal to the plane. Four combinations of velocity profile namely (1) uniform crossflow and uniform jet, (2) laminar boundary layer crossflow and uniform jet, (3) uniform crossflow and parabolic jet profile, and (4) laminar boundary layer crossflow and parabolic jet are compared at same phase to see their effect on the flow field and heat transfer characteristic. All the four cases are seen to exhibit unsteadiness but the jet with parabolic profile is seen to show stronger unsteadiness. The instantaneous vortical structures of all the cases at the same phase show that the structures are more complex for the jet with parabolic velocity profile. The temperature field is seen to be correlated with the vortical structures. Comparison of the time averaged flow field reveals that the jet penetration is the highest for the jet having parabolic profile and boundary layer crossflow. The adiabatic effectiveness is observed to be more for the jet with uniform velocity profile and uniform crossflow and was least for the jet with parabolic velocity profile and boundary layer crossflow.     

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.