Numerical simulation of the interaction of a transverse jet with a supersonic flow using different turbulence models

This paper presents a numerical simulation of the flow resulting from transverse jet injection into a supersonic flow through a slot nozzle at different pressures in the injected jet and the crossflow. Calculations on grids with different resolutions use the Spalart–Allmaras turbulence model, the k–ε model, the k–ω model, and the SST model. Based on a comparison of the calculated and experimental data on the wall pressure distribution, the length of the recirculation area, and the depth of jet penetration into the supersonic flow, conclusions are made on the accuracy of the calculation results for the different turbulence models and the applicability of these models to similar problems.     

Measurements of Turbulent Jet Mixing in a Turbulent Co-Flow Including the Influence of Periodic Forcing and Heating

In this work, the turbulent mixing of a confined coaxial jet in air is investigated by means of simultaneous particle image velocimetry and planar laser induced fluorescence of the acetone seeded flow injection. The jet is injected into a turbulent duct flow at atmospheric pressure through a 90 ^° pipe bend. Measurements are conducted in a small scale windtunnel at constant mass flow rates and three modes of operation: isothermal steady jet injection at a Dean number of 20000 (R e _d =32000), pulsed isothermal injection at a Womersley number of 65 and steady injection at elevated jet temperatures of ΔT=50 K and ΔT=100 K. The experiment is aimed at providing statistically converged quantities of velocity, mass fraction, turbulent fluctuations and turbulent mass flux at several downstream locations. Stochastic error convergence over the number of samples is assessed within the outer turbulent shear layer. From 3000 samples the statistical error of time-averaged velocity and mass fraction is below 1 % while the error of Reynolds shear stress and turbulent mass flux components is in the of range 5-6 %. Profiles of axial velocity and turbulence intensity immediately downstream of the bend exit are in good agreement with hot-wire measurements from literature. During pulsed jet injection strong asymmetric growing of shear layer vortices lead to a skewed mass fraction profile in comparison with steady injection. Phase averaging of single shot PLIF-PIV measurements allows to track the asymmetric shear layer vortex evolvement and flow breakdown during a pulsation cycle with a resolution of 10^°. Steady injection with increased jet temperature supports mixing downstream from 6 nozzle diameters onward.     

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.     

Slow gas motions and the negative drag of a strongly heated spherical particle

In the slow flows of a strongly and nonuniformly heated gas, in the continuum regime (Kn → 0) thermal stresses may be present. The theory of slow nonisothermal continuum gas flows with account for thermal stresses was developed in 1969–1974. The action of the thermal stresses on the gas results in certain paradoxical effects, including the reversal of the direction of the force exerted on a spherical particle in Stokes flow. The propulsion force effect is manifested at large but finite temperature differences between the particle and the gas. This study is devoted to the thermal-stress effect on the drag of a strongly heated spherical particle traveling slowly in a gas for small Knudsen numbers (M ～ Kn → 0), small but finite Reynolds numbers (Re ≤ 1), a linear temperature dependence of the transport coefficients µ ∝ T, and large but finite temperature differences ((T _w ? T _∞ )/T _M8 ～ 1). Two different systems of equations are solved numerically: the simplified Navier-Stokes equations and the modified Navier-Stokes equations with account for the thermal stresses.     

Effect of surface permeability on the structure of a separated turbulent flow and heat transfer behind a backward-facing step

The structure and heat transfer in a turbulent separated flow in a suddenly expanding channel with injection (suction) through a porous wall are numerically simulated with the use of two-dimensional averaged Navier–Stokes equations, energy equations, and v ²–f turbulence model. It is shown that enhancement of the intensity of the transverse mass flux on the wall reduces the separation region length in the case of suction and increases the separation region length in the case of injection up to complete boundary layer displacement. The maximum heat transfer coefficient as a function of permeability is accurately described by the asymptotic theory of a turbulent boundary layer.     

Separation zones in the flow near a small-angle convex corner

On the basis of an asymptotic analysis of the Navier-Stokes system of equations for large Reynolds numbers (Re → ∞), the plane incompressible fluid flow near a surface having a convex corner with a small angle 2θ* is investigated. It is shown that for θ* = O(Re^?1/4), in addition to the known solution that describes a separated flow completely localized in a thin “viscous” sublayer of the interaction region near the corner point, another solution corresponding to a flow with a developed separation zone is possible. For θ ₀ = Re^1/4 θ* = O(1), the longitudinal dimension of this zone varies from finite values up to values of the order of Re^?3/8. The nonuniqueness of the solution is established on a certain range of variation of the parameter θ ₀. The dependence of the drag coefficient on the angle θ* is found.     

Effect of Multilateral Jet Mixing on Stability and Structure of Turbulent Partially-Premixed Flames

An experimental study was conducted to investigate the effects of multilateral jet mixing, using both three and four side-jets, on the structure and stability of turbulent partially-premixed flames. Particle Image Velocimetry and OH*-chemiluminescence were used to study the effects of geometry and operating conditions on the resulting flow-field and reaction zone structures, respectively. These effects were compared under varying ratios of side-jet to primary flow momentum, whilst keeping the bulk flow constant. It was found that the mixing regimes upstream of the nozzle exit affect the flame characteristics, i.e. an impinging regime is likely to generate a lifted flame whilst a backflow regime is likely to generate an attached flame. Unlike the 4 side-jets cases, the OH* images and v _{r m s} profiles for the 3 side-jets cases show distinct asymmetry, with intense OH* and low velocity fluctuations on the opposite sides of the fuel injection. It was also found that the flow and scalar fields become independent of the upstream conditions, for both 3 and 4 side-jets, after one diameter downstream of the nozzle exit.     

Numerical Study of Turbulent Jet Ignition in a Lean Premixed Configuration

Direct numerical simulations (DNS) of a hot combustion product jet interacting with a lean premixed hydrogen-air coflow are conducted to fundamentally investigate turbulent jet ignition (TJI) in a three-dimensional configuration. TJI is an efficient method for initiating and controlling combustion in ultra-lean combustion systems. Fully compressible gas dynamics and species equations are solved with high order finite difference methods. The hydrogen-air reaction is simulated with a reliable detailed chemical kinetics mechanism. The physical processes involved in the TJI-assisted combustion are investigated by considering the flame heat release, temperature, species concentrations, vorticity, and Baroclinc torque. The complex turbulent flame and flow structures are delineated in three main: i) hot product jet, ii) burned-mixed, and iii) flame zones. In the TJI-assisted combustion, the flow structures and the flame features such as flame speed, temperature, and species distribution are found to be quite different than those in “standard” turbulent premixed combustion due to the existence of a high energy turbulent hot product jet. The flow structures and statistics are also found to be different than those normally seen in non-isothermal non-reacting jets.     

A RANS simulation toward the effect of turbulence and cavitation on spray propagation and combustion characteristics

A multidimensional computational fluid dynamic code was developed and integrated with probability density function combustion model to give the detailed account of multiphase fluid flow. The vapor phase within injector domain is treated with Reynolds-averaged Navier–Stokes technique. A new parameter is proposed which is an index of plane-cut spray propagation and takes into account two parameters of spray penetration length and cone angle at the same time. It was found that spray propagation factor (SPI) tends to increase at lower r/d ratios, although the spray penetration tends to decrease. The results of SPI obtained by empirical correlation of Hay and Jones were compared with the simulation computation as a function of respective r/d ratio. Based on the results of this study, the spray distribution on plane area has proportional correlation with heat release amount, NO _x emission mass fraction, and soot concentration reduction. Higher cavitation is attributed to the sharp edge of nozzle entrance, yielding better liquid jet disintegration and smaller spray droplet that reduces soot mass fraction of late combustion process. In order to have better insight of cavitation phenomenon, turbulence magnitude in nozzle and combustion chamber was acquired and depicted along with spray velocity.     

Manipulation of Large-Scale Vortical Structures and Mixing in a Coaxial Jet with Miniature Jet Actuators Toward Active Combustion Control

Manipulation of large-scale vortical structures and associated mixing in a methane-air coaxial jet is carried out by using miniature jet actuators installed on the inner surface of the annular nozzle. The periodic radial miniature jet injections are achieved with a rapid-response servo-valve. The spatio-temporal primary jet structures are investigated through phase-locked 2C-PIV (2 Component Particle Image Velocimetry) and stereoscopic-PIV. It is found that intense ring-like vortices are produced perfectly in phase with the periodic miniature jet injections regardless of the valve-driven frequency f_v examined. When the Strouhal number St_v, which is defined with f_v, is larger than unity, the ring-like vortices are densely formed and thus methane/air mixing is prompted with low periodic fluctuation. The diameter of the vortices becomes small as St_v is increased, so that the transport range of the inner methane and outer air fluids can be controlled by changing St_v. In addition, the evolution of counter-rotating vortex pair is also observed in the cross-sectional plane. These streamwise vortices are directly formed as a result of the radial miniature jet injection, which leads to entrainment of the ambient fluid near the primary jet shear layer, and they also contribute to the mixing enhancement. Moreover, it is demonstrated that coaxial jet flame characteristics such as carbon monoxide (CO) emission and flame holding can be drastically improved under different equivalence ratios by the present jet control scheme.     

Biorthogonal stretching and shearing of an impermeable surface in a uniformly rotating fluid system

The flow induced by an impermeable flat surface executing orthogonal stretching and orthogonal shearing in a rotating fluid system is investigated. Both the stretching and shearing are linear in the coordinates. An exact similarity reduction of the Navier–Stokes equations gives rise to a pair of nonlinearly-coupled ordinary differential equations governed by three parameters. In this study we set one parameter and analyze the problem which leads to flow for an impermeable surface with shearing and stretching due to velocity u along the x-axis of equal strength a while the shearing and stretching due to velocity v along the y-axis of equal strength b. These solutions depend on two parameters—a Coriolis (rotation) parameter \(\sigma = \Omega /a\) and a stretching/shearing ratio \(\lambda =b/a\). A symmetry in solutions is found for \(\lambda = 1\). The exact solution for \(\sigma = 0\) and the asymptotic behavior of solutions for \(|\sigma | \rightarrow \infty\) are determined and compared with numerical results. Oscillatory solutions are found whose strength increases with increasing values of \(|\sigma |\). It is shown that these solutions tend to the well-known Ekman solution as \(|\sigma | \rightarrow \infty\).     

The stress state of a plate subjected to a piecewise homogeneous load

We consider the stress-strain state of a plate having a doubly connected domain S bounded from the outside by a circle of radius R and from the inside by an ellipse with two rectilinear cuts. The cuts lie symmetrically on the x-axis. The plate is subjected to various forces: the hole contour (the ellipse) is under the action of uniformly distributed forces of intensity q, and the cut shores are free of loads; at the points ±ib of the imaginary axis, the plate is under the action of a lumped force P.The solution of the problem is reduced to determining two analytic functions φ(z) and ψ(z) satisfying certain boundary conditions (depending on the type of the acting loads).We use the Kolosov-Muskhelishvili method to reduce the problem to a system of linear algebraic equations for the coefficients in the expansions of the functions φ(z) and ψ(z). The solution thus obtained is illustrated by numerical examples.     

A method of solving the diffusion problem for a vortex layer in a viscoplastic half-plane

A method is proposed to reduce the classical formulation of the problem to a system of two functional equations whose solution can be found numerically. A number of assertions that characterize the behavior of a rigid zone are proved. In particular, the lower estimate h ⁰(t) = 2b√t for the boundary motion is obtained; an explicit expression for b is given as a boundary stress function.     

Spontaneous rotation in the exact solution of magnetohydrodynamic equations for flow between two stationary impermeable disks

Magnetohydrodynamic (MHD) flow of a viscous electrically conducting incompressible fluid between two stationary impermeable disks is considered. A homogeneous electric current density vector normal to the surface is specified on the upper disk, and the lower disk is nonconducting. The exact von Karman solution of the complete system of MHD equations is studied in which the axial velocity and the magnetic field depend only on the axial coordinate. The problem contains two dimensionless parameters: the electric current density on the upper plate Y and the Batchelor number (magnetic Prandtl number). It is assumed that there is no external source that produces an axial magnetic field. The problem is solved for a Batchelor number of 0–2. Fluid flow is caused by the electric current. It is shown that for small values of Y, the fluid velocity vector has only axial and radial components. The velocity of motion increases with increasing Y, and at a critical value of Y, there is a bifurcation of the new steady flow regime with fluid rotation, while the flow without rotation becomes unstable. A feature of the obtained new exact solution is the absence of an axial magnetic field necessary for the occurrence of an azimuthal component of the ponderomotive force, as is the case in the MHD dynamo. A new mechanism for the bifurcation of rotation in MHD flow is found.     

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.     

Monte Carlo Assessment of the Impact of Oscillatory and Pulsating Boundary Conditions on the Flow Through Porous Media

The linear stability analysis of vertical throughflow of power law fluid for double-diffusive convection with Soret effect in a porous channel is investigated in this study. The upper and lower boundaries are assumed to be permeable, isothermal and isosolutal. The linear stability of vertical through flow is influenced by the interactions among the non-Newtonian Rayleigh number (Ra), Buoyancy ratio (N), Lewis number (Le), Péclet number (Pe), Soret parameter (Sr) and power law index (n). The results indicate that the Soret parameter has a significant influence on convective instability of power law fluid. It has also been noticed that buoyancy ratio has a dual effect on the instability of fluid flow. Further, it is noticed that the basic temperature and concentration profiles have singularities at \(Pe = 0\) and \(Le = 1\), the convective instability is looked into for the limiting case of \(Pe\rightarrow 0\) and \(Le \rightarrow 1\). For the case of pure thermal convection with no vertical throughflow, the present numerical results coincide with the solution of standard Horton–Rogers–Lapwood problem. The present results for critical Rayleigh number obtained using bvp4c and two-term Galerkin approximation are compared with those available in the literature and are tabulated.     

Numerical Analysis of the Impact of the Interior Nozzle Geometry on Low Mach Number Jet Acoustics

The flow and acoustic fields of subsonic turbulent hot jets exhausting from three divergent nozzles at a Mach number M=0.12 based on the nozzle exit velocity are conducted using a hybrid CFD-CAA method. The flow field is computed by highly resolved large-eddy simulations (LES) and the acoustic field is computed by solving the acoustic perturbation equations (APE) whose acoustic source terms are determined by the LES. The LES of the computational domain includes the interior of the nozzle geometry. Synthetic turbulence is prescribed at the inlet of the nozzle to mimic the exit conditions downstream of the last turbine stage. The LES is based on hierarchically refined Cartesian meshes, where the nozzle wall boundaries are resolved by a conservative cut-cell method. The APE solution is determined on a block structured mesh. Three nozzle geometries of increasing complexity are considered, i.e., the flow and acoustic fields of a clean geometry without any built-in components, a nozzle with a centerbody, and a nozzle with a centerbody plus struts are computed. Spectral distributions of the LES based turbulent fluctuated quantities inside the nozzle and further downstream are analyzed in detail. The noise sources in the near field are noticeably influenced by the nozzle built-in components. The centerbody nozzle increases the overall sound pressure level (OASPL) in the near field with respect to the clean nozzle and the centerbody-plus-strut nozzle reduces it compared to the centerbody nozzle due to the increased turbulent mixing. The centerbody perturbed nozzle configurations generate a remarkable spectral peak at S t=0.56 which also occurs in the APE findings in the near field region. This tone is generated by large scale vortical structures shed from the centerbody. The analysis of the individual noise sources shows that the entropy term possesses the highest acoustic contribution in the sideline direction whereas the vortex sound source dominates the downstream acoustics.     

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.     

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.     

Integral boundary layer relations in the theory of wave flows for capillary liquid films

A generalized method of deriving the model equations is considered for wave flow regimes in falling liquid films. The viscous liquid equations are used on the basis of integral boundary layer relations with weight functions. A family of systems of evolution differential equations is proposed. The integer parameter n of these systems specifies the number of a weight function. The case n = 0 corresponds to the classical IBL (Integral Boundary Layer) model. The case n ≥ 1 corresponds to its modifications called the WIBL (Weighted Integral Boundary Layer) models. The numerical results obtained in the linear and nonlinear approximations for n = 0, 1, 2 are discussed. The numerical solutions to the original hydrodynamic differential equations are compared with experimental data. This comparison leads us to the following conclusions: as a rule, the most accurate solutions are obtained for n = 0 in the case of film flows on vertical and inclined solid surfaces and the accuracy of solutions decreases with increasing n. Hence, the classical IBL model has an advantage over the WIBL models.