On Coupling of RANS and LES for Integrated Computations of Jet Engines

Author(s):  
Gorazd Medic ◽  
Donghyun You ◽  
Georgi Kalitzin

Large scale integrated computations of jet engines can be performed by using the unsteady RANS framework to compute the flow in turbomachinery components while using the LES framework to compute the flow in the combustor. This requires a proper coupling of the flow variables at the interfaces between the RANS and LES solvers. In this paper, a novel approach to turbulence coupling is proposed. It is based on the observation that in full operating conditions the mean flow at the interfaces is highly non-uniform and local turbulence production dominates convection effects in regions of large velocity gradients. This observation has lead to the concept of using auxilliary ducts to compute turbulence based on the mean velocity at the interface. In the case of the RANS/LES interface, turbulent fluctuations are reconstructed from an LES computation in an auxiliary three-dimensional duct using a recycling technique. For the LES/RANS interface, the turbulence variables for the RANS model are computed from an auxilliary solution of the RANS turbulence model in a quasi-2D duct. We have demonstrated the feasibility of this approach for the integrated flow simulation of a 20° sector of an entire jet engine.

2021 ◽  
Vol 929 ◽  
Author(s):  
N. Agastya Balantrapu ◽  
Christopher Hickling ◽  
W. Nathan Alexander ◽  
William Devenport

Experiments were performed over a body of revolution at a length-based Reynolds number of 1.9 million. While the lateral curvature parameters are moderate ( $\delta /r_s < 2, r_s^+>500$ , where $\delta$ is the boundary layer thickness and r s is the radius of curvature), the pressure gradient is increasingly adverse ( $\beta _{C} \in [5 \text {--} 18]$ where $\beta_{C}$ is Clauser’s pressure gradient parameter), representative of vehicle-relevant conditions. The mean flow in the outer regions of this fully attached boundary layer displays some properties of a free-shear layer, with the mean-velocity and turbulence intensity profiles attaining self-similarity with the ‘embedded shear layer’ scaling (Schatzman & Thomas, J. Fluid Mech., vol. 815, 2017, pp. 592–642). Spectral analysis of the streamwise turbulence revealed that, as the mean flow decelerates, the large-scale motions energize across the boundary layer, growing proportionally with the boundary layer thickness. When scaled with the shear layer parameters, the distribution of the energy in the low-frequency region is approximately self-similar, emphasizing the role of the embedded shear layer in the large-scale motions. The correlation structure of the boundary layer is discussed at length to supply information towards the development of turbulence and aeroacoustic models. One major finding is that the estimation of integral turbulence length scales from single-point measurements, via Taylor's hypothesis, requires significant corrections to the convection velocity in the inner 50 % of the boundary layer. The apparent convection velocity (estimated from the ratio of integral length scale to the time scale), is approximately 40 % greater than the local mean velocity, suggesting the turbulence is convected much faster than previously thought. Closer to the wall even higher corrections are required.


1999 ◽  
Vol 390 ◽  
pp. 325-348 ◽  
Author(s):  
S. NAZARENKO ◽  
N. K.-R. KEVLAHAN ◽  
B. DUBRULLE

A WKB method is used to extend RDT (rapid distortion theory) to initially inhomogeneous turbulence and unsteady mean flows. The WKB equations describe turbulence wavepackets which are transported by the mean velocity and have wavenumbers which evolve due to the mean strain. The turbulence also modifies the mean flow and generates large-scale vorticity via the averaged Reynolds stress tensor. The theory is applied to Taylor's four-roller flow in order to explain the experimentally observed reduction in the mean strain. The strain reduction occurs due to the formation of a large-scale vortex quadrupole structure from the turbulent spot confined by the four rollers. Both turbulence inhomogeneity and three-dimensionality are shown to be important for this effect. If the initially isotropic turbulence is either homogeneous in space or two-dimensional, it has no effect on the large-scale strain. Furthermore, the turbulent kinetic energy is conserved in the two-dimensional case, which has important consequences for the theory of two-dimensional turbulence. The analytical and numerical results presented here are in good qualitative agreement with experiment.


Author(s):  
Takuma Katayama ◽  
Shinsuke Mochizuki

The present experiment focuses on the vorticity diffusion in a stronger wall jet managed by a three-dimensional flat plate wing in the outer layer. Measurement of the fluctuating velocities and vorticity correlation has been carried out with 4-wire vorticity probe. The turbulent vorticity diffusion due to the large scale eddies in the outer layer is quantitatively examined by using the 4-wire vorticity probe. Quantitative relationship between vortex structure and Reynolds shear stress is revealed by means of directly measured experimental evidence which explains vorticity diffusion process and influence of the manipulating wing. It is expected that the three-dimensional outer layer manipulator contributes to keep convex profile of the mean velocity, namely, suppression of the turbulent diffusion and entrainment.


1976 ◽  
Vol 77 (3) ◽  
pp. 473-497 ◽  
Author(s):  
L. J. S. Bradbury

This paper describes an investigation into the response of both the pulsed-wire anemometer and the hot-wire anemometer in a highly turbulent flow. The first part of the paper is concerned with a theoretical study of some aspects of the response of these instruments in a highly turbulent flow. It is shown that, under normal operating conditions, the pulsed-wire anemometer should give mean velocity and longitudinal turbulent intensity estimates to an accuracy of better than 10% without any restriction on turbulence level. However, to attain this accuracy in measurements of turbulent intensities normal to the mean flow direction, there is a lower limit on the turbulent intensity of about 50%. An analysis is then carried out of the behaviour of the hot-wire anemometer in a highly turbulent flow. It is found that the large errors that are known to develop are very sensitive to the precise structure of the turbulence, so that even qualitative use of hot-wire data in such flows is not feasible. Some brief comments on the possibility of improving the accuracy of the hot-wire anemometer are then given.The second half of the paper describes some comparative measurements in the highly turbulent flow immediately downstream of a normal flat plate. It is shown that, although it is not possible to interpret the hot-wire results on their own, it is possible to calculate the hot-wire response with a surprising degree of accuracy using the results from the pulsed-wire anemometer. This provides a rather indirect but none the less welcome check on the accuracy of the pulsed-wire results, which, in this very highly turbulent flow, have a certain interest in their own right.


2012 ◽  
Vol 699 ◽  
pp. 320-351 ◽  
Author(s):  
Johan Malm ◽  
Philipp Schlatter ◽  
Dan S. Henningson

AbstractDominant frequencies and coherent structures are investigated in a turbulent, three-dimensional and separated diffuser flow at $\mathit{Re}= 10\hspace{0.167em} 000$ (based on bulk velocity and inflow-duct height), where mean flow characteristics were first studied experimentally by Cherry, Elkins and Eaton (Intl J. Heat Fluid Flow, vol. 29, 2008, pp. 803–811) and later numerically by Ohlsson et al. (J. Fluid Mech., vol. 650, 2010, pp. 307–318). Coherent structures are educed by proper orthogonal decomposition (POD) of the flow, which together with time probes located in the flow domain are used to extract frequency information. The present study shows that the flow contains multiple phenomena, well separated in frequency space. Dominant large-scale frequencies in a narrow band $\mathit{St}\equiv fh/ {u}_{b} \in [0. 0092, 0. 014] $ (where $h$ is the inflow-duct height and ${u}_{b} $ is the bulk velocity), yielding time periods ${T}^{\ensuremath{\ast} } = T{u}_{b} / h\in [70, 110] $, are deduced from the time signal probes in the upper separated part of the diffuser. The associated structures identified by the POD are large streaks arising from a sinusoidal oscillating motion in the diffuser. Their individual contributions to the total kinetic energy, dominated by the mean flow, are, however, small. The reason for the oscillating movement in this low-frequency range is concluded to be the confinement of the flow in this particular geometric set-up in combination with the high Reynolds number and the large separated zone on the top diffuser wall. Based on this analysis, it is shown that the bulk of the streamwise root mean square (r.m.s.) value arises due to large-scale motion, which in turn can explain the appearance of two or more peaks in the streamwise r.m.s. value. The weak secondary flow present in the inflow duct is shown to survive into the diffuser, where it experiences an imbalance with respect to the upper expanding corners, thereby giving rise to the asymmetry of the mean separated region in the diffuser.


1969 ◽  
Vol 36 (4) ◽  
pp. 785-803 ◽  
Author(s):  
Francis P. Bretherton

A train of internal gravity waves in a stratified liquid exerts a stress on the liquid and induces changes in the mean motion of second order in the wave amplitude. In those circumstances in which the concept of a slowly varying quasi-sinusoidal wave train is consistent, the mean velocity is almost horizontal and is determined to a first approximation irrespective of the vertical forces exerted by the waves. The sum of the mean flow kinetic energy and the wave energy is then conserved. The circulation around a horizontal circuit moving with the mean velocity is increased in the presence of waves according to a simple formula. The flow pattern is obtained around two- and three-dimensional wave packets propagating into a liquid at rest and the results are generalized for any basic state of motion in which the internal Froude number is small. Momentum can be associated with a wave packet equal to the horizontal wave-number times the wave energy divided by the intrinsic frequency.


2015 ◽  
Vol 774 ◽  
pp. 143-169 ◽  
Author(s):  
J.-H. Xie ◽  
J. Vanneste

Wind forcing of the ocean generates a spectrum of inertia–gravity waves that is sharply peaked near the local inertial (or Coriolis) frequency. The corresponding near-inertial waves (NIWs) are highly energetic and play a significant role in the slow, large-scale dynamics of the ocean. To analyse this role, we develop a new model of the non-dissipative interactions between NIWs and balanced motion. The model is derived using the generalised-Lagrangian-mean (GLM) framework (specifically, the ‘glm’ variant of Soward & Roberts, J. Fluid Mech., vol. 661, 2010, pp. 45–72), taking advantage of the time-scale separation between the two types of motion to average over the short NIW period. We combine Salmon’s (J. Fluid Mech., vol. 719, 2013, pp. 165–182) variational formulation of GLM with Whitham averaging to obtain a system of equations governing the joint evolution of NIWs and mean flow. Assuming that the mean flow is geostrophically balanced reduces this system to a simple model coupling Young & Ben Jelloul’s (J. Mar. Res., vol. 55, 1997, pp. 735–766) equation for NIWs with a modified quasi-geostrophic (QG) equation. In this coupled model, the mean flow affects the NIWs through advection and refraction; conversely, the NIWs affect the mean flow by modifying the potential-vorticity (PV) inversion – the relation between advected PV and advecting mean velocity – through a quadratic wave term, consistent with the GLM results of Bühler & McIntyre (J. Fluid Mech., vol. 354, 1998, pp. 301–343). The coupled model is Hamiltonian and its conservation laws, for wave action and energy in particular, prove illuminating: on their basis, we identify a new interaction mechanism whereby NIWs forced at large scales extract energy from the balanced flow as their horizontal scale is reduced by differential advection and refraction so that their potential energy increases. A rough estimate suggests that this mechanism could provide a significant sink of energy for mesoscale motion and play a part in the global energetics of the ocean. Idealised two-dimensional models are derived and simulated numerically to gain insight into NIW–mean-flow interaction processes. A simulation of a one-dimensional barotropic jet demonstrates how NIWs forced by wind slow down the jet as they propagate into the ocean interior. A simulation assuming plane travelling NIWs in the vertical shows how a vortex dipole is deflected by NIWs, illustrating the irreversible nature of the interactions. In both simulations energy is transferred from the mean flow to the NIWs.


1999 ◽  
Vol 394 ◽  
pp. 303-337 ◽  
Author(s):  
A. VERNET ◽  
G. A. KOPP ◽  
J. A. FERRÉ ◽  
FRANCESC GIRALT

Simultaneous velocity and temperature measurements were made with rakes of sensors that sliced a slightly heated turbulent wake in the spanwise direction, at different lateral positions 150 diameters downstream of the cylinder. A pattern recognition analysis of hotter-to-colder transitions was performed on temperature data measured at the mean velocity half-width. The velocity data from the different ‘slices’ was then conditionally averaged based on the identified temperature events. This procedure yielded the topology of the average three-dimensional large-scale structure which was visualized with iso-surfaces of negative values of the second eigenvector of [S2+Ω2]. The results indicate that the average structure of the velocity fluctuations (using a triple decomposition of the velocity field) is found to be a shear-aligned ring-shaped vortex. This vortex ring has strong outward lateral velocities in its symmetry plane which are like Grant's mixing jets. The mixing jet region extends outside the ring-like vortex and is bounded by two foci separated in the spanwise direction and an upstream saddle point. The two foci correspond to what has been previously identified in the literature as the double rollers.The ring vortex extracts energy from the mean flow by stretching in the mixing jet region just upstream of the ring boundary. The production of the small-scale (incoherent) turbulence by the coherent field and one-component energy dissipation rate occur just downstream of the saddle point within the mixing jet region. Incoherent turbulence energy is extracted from the mean flow just outside the mixing jet region, but within the core of the structure. These processes are highly three-dimensional with a spanwise extent equal to the mean velocity half-width.When a double decomposition is used, the coherent structure is found to be a tube-shaped vortex with a spanwise extent of about 2.5l0. The double roller motions are integral to this vortex in spite of its shape. Spatial averages of the coherent velocity field indicate that the mixing jet region causes a deficit of mean streamwise momentum, while the region outside the foci of the double rollers has a relatively small excess of streamwise momentum.


1982 ◽  
Vol 119 ◽  
pp. 121-153 ◽  
Author(s):  
Udo R. Müller

An experimental study of a steady, incompressible, three-dimensional turbulent boundary layer approaching separation is reported. The flow field external to the boundary layer was deflected laterally by turning vanes so that streamwise flow deceleration occurred simultaneous with cross-flow acceleration. At 21 stations profiles of the mean-velocity components and of the six Reynolds stresses were measured with single- and X-hot-wire probes, which were rotatable around their longitudinal axes. The calibration of the hot wires with respect to magnitude and direction of the velocity vector as well as the method of evaluating the Reynolds stresses from the measured data are described in a separate paper (Müller 1982, hereinafter referred to as II). At each measuring station the wall shear stress was inferred from a Preston-tube measurement as well as from a Clauser chart. With the measured profiles of the mean velocities and of the Reynolds stresses several assumptions used for turbulence modelling were checked for their validity in this flow. For example, eddy viscosities for both tangential directions and the corresponding mixing lengths as well as the ratio of resultant turbulent shear stress to turbulent kinetic energy were derived from the data.


Water ◽  
2019 ◽  
Vol 11 (2) ◽  
pp. 237
Author(s):  
Ming Chen ◽  
Haijin Huang ◽  
Xingxing Zhang ◽  
Senpeng Lv ◽  
Rengmin Li

Three-dimensional (3D) confined wall jets have various engineering applications related to efficient energy dissipation. This paper presents experimental measurements of mean flow development for a 3D rectangular wall jet confined by a vertical baffle with a fixed distance (400 mm) from its surface to the nozzle. Experiments were performed at three different Reynolds numbers of 8333, 10,000 and 11,666 based on jet exit velocity and square root of jet exit area (named as B), with water depth of 100 mm. Detailed measurements of current jet were taken using a particle image velocimetry technique. The results indicate that the confined jet seems to behave like an undisturbed jet until 16B downstream. Beyond this position, however, the mean flow development starts to be gradually affected by the baffle confinement. The baffle increases the decay and spreading of the mean flow from 16B to 23B. The decay rate of 1.11 as well as vertical and lateral growth rates of 0.04 and 0.19, respectively, were obtained for the present study, and also fell well within the range of values which correspond to the results in the radial decay region for the unconfined case. In addition, the measurements of the velocity profiles, spreading rates and velocity decay were also found to be independent of Reynolds number. Therefore, the flow field in this region appears to have fully developed at least 4B earlier than the unconfined case. Further downstream (after 23B), the confinement becomes more pronounced. The vertical spreading of current jet shows a distinct increase, while the lateral growth was found to be decreased significantly. It can be also observed that the maximum mean velocity decreases sharply close to the baffle.


Sign in / Sign up

Export Citation Format

Share Document