Vortical structures behind a transverse jet in a supersonic flow at high jet to crossflow pressure ratios

A three-dimensional supersonic turbulent flow with symmetric normal injection of circular jets from the channel walls is numerically simulated. The initial Favre-averaged Navier–Stokes equations closed by the k–ω turbulence model are solved by an algorithm based on an ENO scheme. The mechanism of the formation of vortical structures due to the interaction of the jet with the free stream is studied for jet to crossflow total pressure ratios ranging from 3 to 50. It is known from experiments reported in the literature that, for n ? 10, mixing of the jet with the high-velocity flow leads to the formation of a pair of vortices and of an additional separation zone near the wall behind the jet. It is demonstrated that the present numerical results are consistent with such findings and that the pressure distribution on the wall ahead of the jet in the plane of symmetry is also in reasonable agreement with available experimental data.     

Application of the γ-Reθ Transition Model to Simulations of the Flow Past a Circular Cylinder

Based on the finite volume method, the flow past a two-dimensional circular cylinder at a critical Reynolds number (Re = 8.5 × 10⁵) was simulated using the Navier-Stokes equations and the γ-Re_θ transition model coupled with the SST k ? ω turbulence model (hereinafter abbreviated as γ-Re_θ model). Considering the effect of free-stream turbulence intensity decay, the SST k ? ω turbulence model was modified according to the ambient source term method proposed by Spalart and Rumsey, and then the modified SST k ? ω turbulence model is coupled with the γ-Re_θ transition model (hereinafter abbreviated as γ-Re_θ-SR model). The flow past a circular cylinder at different inlet turbulence intensities were simulated by the γ-Re_θ-SR model. At last, the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers were each simulated by the γ-Re_θ-SR model, and the three flow states were analyzed. It was found that compared with the SST k ? ω turbulence model, the γ-Re_θ model could simulate the transition of laminar to turbulent, resulting in better consistency with experimental result. Compared with the γ-Re_θ model, for relatively high inlet turbulence intensities, the γ-Re_θ-SR model could better simulate the flow past a circular cylinder; however the improvement almost diminished for relatively low inlet turbulence intensities The γ-Re_θ-SR model could well simulate the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers.     

Experimental and Theoretical Sensitivity Analysis of Turbulent Flow Past a Square Cylinder

The Siemens SGT-800 3^rd generation DLE burner fitted to an atmospheric combustion rig has been numerically investigated. Pure methane and methane enriched by 80 vol% hydrogen flames have been considered. A URANS (Unsteady Reynolds Averaged Navier-Stokes) approach was used in this study along with the k ? ω SST and the k ? ω SST-SAS models for the turbulence transport. The chemistry is coupled to the turbulent flow simulations by the use of a laminar flamelet library combined with a presumed PDF. The effect of the mesh density in the mixing and the flame region and the effect of the turbulence model and reaction rate model constant are first investigated for the methane/air flame case. The results from the k ? ω SST-SAS along with flamelet libraries are shown to be in excellent agreement with experimental data, whereas the k ? ω SST model is too dissipative and cannot capture the unsteady motion of the flame. The k ? ω SST-SAS model is used for simulation of the 80 vol% hydrogen enriched flame case without further adjusting the model constants. The global features of the hydrogen enrichment are very well captured in the simulations using the SST-SAS model. With the hydrogen enrichment the time averaged flame front location moves upstream towards the burner exit nozzle. The results are consistent with the experimental observations. The model captures the three dominant low frequency unsteady motion observed in the experiments, indicating that the URANS/LES hybrid model indeed is capable of capturing complex, time dependent, features such as an interaction between a PVC and the flame front.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

Numerical analysis of the flow pattern and vortex breakdown over a pitching delta wing at supersonic speeds

A supersonic compressible flow over a 60° swept delta wing with a sharp leading edge undergoing pitching oscillations is computationally studied. Numerical simulations are performed by the finite volume method with the use of the k?ω turbulence model for various Mach numbers and angles of attack. Variations of flow patterns in a crossflow plane, hysteresis loops associated with the vortex core location, and vortex breakdown positions during a pitching cycle are investigated. Trends for various Mach numbers, mean angles of attack, pitching amplitudes, and pitching frequencies are illustrated.     

Combining the SSG/LRR-<Emphasis Type="Italic">ω</Emphasis> differential reynolds stress model with the detached eddy and laminar-turbulent transition models

The procedure of incorporating the detached eddy method and a model of laminar-turbulent transition into the SSG/LRR-ω turbulence model is presented. The approach proposed can be regarded as the generalization of the existing models intended to perform calculations with the SST turbulence model to the case of their use with the SSG/LRR-ω model. The advantage of the approach developed over the RANS turbulence models based on the Boussinesq hypothesis is demonstrated with respect to the problems of flow past an airfoil and cold jet outflow.     

Flow visualization of supersonic laminar flow over a backward-facing step via NPLS

An experimental study on a supersonic laminar flow over a backward-facing step of 5 mm height was undertaken in a low-noise indraft wind tunnel. To investigate the fine structures of Ma = 3.0 and 3.8 laminar flow over a backward-facing step, nanotracer planar laser scattering was adopted for flow visualization. Flow structures, including supersonic laminar boundary layer, separation, reattachment, redeveloping turbulent boundary layer, expansion wave fan and reattachment shock, were revealed in the transient flow fields. In the Ma = 3.0 BFS (backward-facing step) flow, by measuring four typical regions, it could be found that the emergence of weak shock waves was related to the K–H (Kelvin–Helmholtz) vortex which appeared in the free shear layer and that the convergence of these waves into a reattachment shock was distinct. Based on large numbers of measurements, the structure of time-averaging flow field could be gained. Reattachment occurred at the location downstream from the step, about 7–7.5 h distance. After reattachment, the recovery boundary layer developed into turbulence quickly and its thickness increased at an angle of 4.6°. At the location of X = 14h, the redeveloping boundary layer was about ten times thicker than its original thickness, but it still had not changed into fully developed turbulence. However, in the Ma = 3.8 flow, the emergence of weak shock waves could be seen seldom, due to the decrease of expansion. The reattachment point was thought to be near X = 15h according to the averaging result. The reattachment shock was not legible, which meant the expansion and compression effects were not intensive.     

Turbulent flows in channels with injection. results of large eddy simulation and the two-equation turbulence model

Turbulent flows in channels with intense distributed injection are modeled using the large eddy method and the two-equation k-? turbulence model. The calculations are carried out for different velocities of injection from the channel walls. For a channel with one-sided injection the results of large eddy simulation are in good agreement with the measured data, whereas the calculations in accordance with the k-? model give a less convex cross-sectional velocity profile and an appreciable error in determining the surface friction coefficient on the impermeable wall and also have certain other shortcomings. In the case of two-sided injection, the results of the calculations by the large eddy method and the k-? model are in good agreement with one another and the data of physical experiments.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

Theory of unsteady separation in supersonic flow around a convex corner

A theory of unsteady separation in inviscid supersonic flow around a convex corner is developed. Within the framework of the hypothesis suggested the mechanism of separationless-to-separated flow transition is explained and the forces leading to flow separation are determined as functions of the angle θ and the oncoming flow velocity. The values of the angle θ _k at which the flow is separated from the corner vertex and the stall angle θ _s determining the separated flow direction obtained previously in experiments and by numerical simulation are confirmed.     

Reynolds-Averaged,Scale-Adaptive and Large-Eddy Simulations of Premixed Bluff-Body Combustion Using the Eddy Dissipation Concept

A lean premixed propane/air bluff-body stabilized flame (Volvo test rig) is calculated using the Scale-Adaptive Simulation turbulence model (SAS) and Large-Eddy simulations (LES) as well as the conventional Reynolds-averaged approach (RAS). RAS and SAS are closed by the standard k-?? and the k-ω Shear Stress Transport (SST) turbulence models, respectively. The conventional Smagorinsky and the k-equation sub-grid scales models are used for the LES closure. Effects of the sub-grid scalar flux modeling using the classical gradient hypothesis and Clark’s tensor diffusivity closures both for the inert and reactive LES flows are discussed. The Eddy Dissipation Concept (EDC) is used for the turbulence-chemistry interaction. It assumes that molecular mixing and the subsequent combustion occur in the ’fine structures’ (smaller dissipative eddies, which are close to the Kolmogorov scales). Assuming the full turbulence energy cascade, the characteristic length and velocity scales of the ’fine structures’ are evaluated using different turbulence models (RAS, SAS and LES). The finite-rate chemical kinetics is taken into account by treating the ’fine structures’ as constant pressure and adiabatic homogeneous reactors, calculated as a system of ordinary-differential equations (ODEs) described by a Perfectly Stirred Reactor (PSR) concept. Several further enhancements to model the PSRs are proposed, including a new Livermore Solver (LSODA) for integrating stiff ODEs and a new correction to calculate the PSR time scales. All models have been implemented as a stand-alone application \(\text {edcPisoFoam}\) based on the OpenFOAM technology. Additionally, several RAS calculations were performed using the Turbulence Flame Speed Closure model in Ansys Fluent to assess effects of the heat losses by modeling the conjugate heat transfer between the bluff-body and the reactive flow. Effects of the turbulence Schmidt number on RAS results are discussed as well. Numerical results are compared with available experimental data. Reasonable consistency between experimental data and numerical results provided by RAS, SAS and LES is observed. In general, there is satisfactory agreement between present LES-EDC simulations, numerical results by other authors and measurements without any major modification to the EDC closure constants, which gives a quite reasonable indication on the adequacy and accuracy of the method and its further application for turbulent premixed combustion simulations.     

Numerical Simulation of the Negative Magnus Effect of a Two-Dimensional Spinning Circular Cylinder

Based on the finite volume method, the flow past a spinning circular cylinder at a low subcritical Reynolds number (Re =1 × 10 ⁵), high subcritical Reynolds number (Re =1.3 ×10 ⁵), and critical Reynolds number (Re =1.4 ×10 ⁵) were each simulated using the Navier-Stokes equations and the γ-Re _?? transition model coupled with the SST k?ω turbulence model. The system was solved using an implicit algorithm. The freestream turbulence intensity decay was effectively controlled by the source term method proposed by Spalart and Rumsey. The variations in the Magnus force as a function of the spin ratio, α were obtained for the three Reynolds numbers, and the flow mechanism was analyzed. The results indicate that the asymmetric transitions induced by spin affect the asymmetric separations at the top and bottom surfaces of the circular cylinder, which further affects the pressure distributions at the top and bottom surfaces of the circular cylinder and ultimately result in a negative Magnus force, whose direction is opposite to that of the classical Magnus force. This study is the first to use a numerical simulation method to predict a negative Magnus force acting on a spinning circular cylinder. At the low subcritical Reynolds number, the Magnus force remained positive for all spin ratios. At the high subcritical Reynolds number, the sign of the Magnus force changed twice over the range of the spin ratio. At the critical Reynolds number, the sign of the Magnus force changed only once over the range of the spin ratio. For relatively low spin ratios, the Magnus force significantly differed by Reynolds number; however, this variation diminished as the spin ratio increased.     

Helical Structures in the Near Field of a Turbulent Pipe Jet

We perform a finely resolved Large-eddy simulation to study coherent vortical structures populating the initial (near-nozzle) zone of a pipe jet at the Reynolds number of 5300. In contrast to ‘top-hat’ jets featured by Kelvin-Helmholtz rings with the non-dimensional frequency S t≈0.3?0.6, no high-frequency dominant mode is observed in the near field of a jet issuing from a fully-developed pipe flow. Instead, in shear layers we observe a relatively wide peak in the power spectrum within the low-frequency range (S t≈0.14) corresponding to the propagating helical waves entering with the pipe flow. This is confirmed by the Fourier transform with respect to the azimuthal angle and the Proper Orthogonal Decomposition complemented with the linear stability analysis revealing that this low-frequency motion is not connected to the Kelvin-Helmholtz instability. We demonstrate that the azimuthal wavenumbers m=1?5 contain the most of the turbulent kinetic energy and that a common form of an eigenmode is a helical vortex rotating around the axis of symmetry. Small and large timescales are identified corresponding to “fast” and “slow” rotating modes. While the “fast” modes correspond to background turbulence and stochastically switch from co- to counter-rotation, the “slow” modes are due to coherent helical structures which are long-lived and have low angular velocities, in agreement with the previously described spectral peak at low S t.     

Large-eddy and unsteady RANS simulations of a shock-accelerated heavy gas cylinder

Two-dimensional numerical simulations of the Richtmyer–Meshkov unstable “shock-jet” problem are conducted using both large-eddy simulation (LES) and unsteady Reynolds-averaged Navier–Stokes (URANS) approaches in an arbitrary Lagrangian–Eulerian hydrodynamics code. Turbulence statistics are extracted from LES by running an ensemble of simulations with multimode perturbations to the initial conditions. Detailed grid convergence studies are conducted, and LES results are found to agree well with both experiment and high-order simulations conducted by Shankar et al. (Phys Fluids 23, 024102, 2011). URANS results using a k–L approach are found to be highly sensitive to initialization of the turbulence lengthscale L and to the time at which L becomes resolved on the computational mesh. It is observed that a gradient diffusion closure for turbulent species flux is a poor approximation at early times, and a new closure based on the mass-flux velocity is proposed for low-Reynolds-number mixing.     

Assessment of Five SGS Models for Low-Rem MHD Turbulence

Assessment of three regularization-based and two eddy-viscosity-based subgrid-scale (SGS) turbulence models for large eddy simulations (LES) are carried out in the context of magnetohydrodynamic (MHD) decaying homogeneous turbulence (DHT) with a Taylor scale Reynolds number (Re_λ) of 120 and a MHD transition-to-turbulence Taylor-Green vortex (TGV) problems with a Reynolds number of 3000, through direct comparisons to direct numerical simulations (DNS). Simulations are conducted using the low-magnetic Reynolds number approximation (Re_m<<1). LES predictions using the regularization-based Leray- α,LANS- α, and Clark- α SGS models, along with the eddy viscosity-based non-dynamic Smagorinsky and the dynamic Smagorinsky models are compared to in-house DNS for DHT and previous results for TGV. With regard to the regularization models, this work represents their first application to MHD turbulence. Analyses of turbulent kinetic energy decay rates, energy spectra, and vorticity fields made between the varying magnetic field cases demonstrated that the regularization models performed poorly compared to the eddy-viscosity models for all MHD cases, but the comparisons improved with increase in magnitude of magnetic field, due to a decrease in the population of SGS eddies within the flow field.     

Effect of Inflow Turbulence on an Airfoil Flow with Laminar Separation Bubble: An LES Study

The present paper is concerned with numerical investigations on the effect of inflow turbulence on the flow around a SD7003 airfoil. At a Reynolds number Re_c =?60,000, an angle of attack α =?4^° and a low or zero turbulence intensity of the oncoming flow, the flow past the airfoil is known to be dominated by early separation, subsequent transition and reattachment leading to a laminar separation bubble with a distinctive pressure plateau. The objective of the study is to investigate the effect of inflow turbulence on the flow behavior. For this purpose, a numerical methodology relying on a wall-resolved large-eddy simulation, a synthetic turbulence inflow generator and a specific source term concept for introducing the turbulence fluctuations within the computational domain is used. The numerical technique applied allows the variation of the free-stream turbulence intensity (TI) in a wide range. In order to analyze the influence of TI on the arising instantaneous and time-averaged flow field past the airfoil, the present study evaluates the range 0% ≤ TI ≤?11.2%, which covers typical values found in atmospheric boundary layers. In accordance with experimental studies it is shown that the laminar separation bubble first shrinks and finally completely vanishes for increasing inflow turbulence. Consequently, the aerodynamic performance in terms of the lift-to-drag ratio increases. Furthermore, the effect of the time and length scales of the isotropic inflow turbulence on the development of the flow field around the airfoil is analyzed and a perceptible influence is found. Within the range of inflow scales studied decreasing scales augment the receptivity of the boundary layer promoting an earlier transition.     

Analytical and numerical solutions of the problem of gravity-induced mixing of a light layer using the <Emphasis Type="Italic">k</Emphasis>–<Emphasis Type="Italic">ε</Emphasis> turbulence model

An analytical solution for the self-similar stage in the problem of gravity-induced turbulent mixing in a light (heavy) layer is obtained on the basis of the k–ε model equations. The solution obtained is compared with the results of a numerical investigation of the problem using both three-dimensional direct numerical simulation and the k–ε model. The calculations were performed using the two- and three-dimensional versions of the EGAK method. The results of all the calculations and the available experimental data are in reasonable agreement.     

Effect of Free-Stream Turbulence on the Reduction of Surface Friction in a Turbulent Boundary Layer Behind LEBU-Devices

The effect of increased free-stream turbulence on the reduction of the surface friction coefficient c _f in a turbulent boundary layer behind large-eddy break-up (LEBU) devices is investigated using a gravimetric method. The turbulence level was ε ≈ 1.9–4.9 % and the turbulence scale L _e ≈ 40–110 mm. The boundary layer Reynolds number Re** was varied from 2300 to 7500, with the boundary layer thickness being varied on the range δ = 33–44 mm. It is shown that an increase in the turbulence level ε has almost no impact on the relative reduction of friction behind LEBU-devices, whereas, under similar conditions of elevated free-stream turbulence, for another method, namely, the use of surface riblets, the friction reduction may be more strongly expressed.     

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.