Large eddy simulation of a hot, high speed coflowing jet

Computations are made of a short cowl coflowing jet nozzle with a bypass ratio 8 : 1. The core flow is heated, making the inlet conditions reminiscent of those for a real engine. A large eddy resolving approach is used with a 12 × 106 cell mesh. Since the code being used tends towards being dissipative the sub-grid scale (SGS) model is abandoned giving what can be termed Numerical Large Eddy Simulation (NLES). To overcome near wall modelling problems a hybrid NLES-RANS (Reynolds Averaged Navier-Stokes) related method is used. For y+ ≤ 60 a κ-l model is used. Blending between the two regions makes use of the differential Hamilton-Jabobi (HJ) equation, an extension of the eikonal equation. Results show encouraging agreement with existing measurements of other workers. The eikonal equation is also used for acoustic ray tracing to explore the effect of the mean flow on acoustic ray trajectories, thus yielding a coherent solution strategy. Copyright © 2011 by ASME.

Large-eddy simulations and measurements of a small-scale high-speed coflowing jet

Measurements and predictions are made of a short-cowl coflowing jet with a bypass ratio of 8:1. The Reynolds number is 300,000, and the inlet Mach numbers are representative of aeroengine conditions. The low Reynolds number of the measurements makes the case well suited to the assessment of large-eddy-simulation-related strategies. The nozzle concentricity is carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both laser Doppler anemometry and particle image velocimetry. The simulations are completed on 6× 106, 12× 106, and 50 × 106 cell meshes. To overcome near-wall modeling problems, a hybrid large-eddy-simulation-Reynolds-averaged-Navier-Stokesrelated method is used. The near-wall Reynolds-averaged-Navier-Stokes layer is helpful in preventing nonphysical separation from the nozzle wall.Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Comparison of les to LDA and PIV measurements of a small scale high speed coflowing jet

Measurements and predictions are made of a short cowl co-flowing jet with a bypass ratio of 8:1. The Reynolds number for computations and measurements are matched at 300,000 and the Mach numbers representative of realistic jet conditions with core and co flow velocities of 240m/s and 216m/s respectively. The low Reynolds number of the measurements makes the case well suited to the assessment of large eddy resolving computational strategies. Also, the nozzle concentricity was carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both two dimensional coincident LDA and PIV systems. The computational simulations are completed on a modest 12×106 mesh. The simulation is now being run on a 50×106 mesh using hybrid RANSNLES (Numerical Large Eddy Simulation). Close to the nozzle wall a k-l RANS model is used. For an axisymmetric jet, comparison is made between simulations which use NLES, RANSNLES and also a simple imposed velocity profile where the nozzle is not modeled. The use of a near wall RANS model is shown to be beneficial. When compared with the measurements the NLES results are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Linear and nonlinear model large-eddy simulations of a plane jet

The nonlinear Kosovic, and mixed Leray and α subgrid scale models are contrasted with linear Smagorinsky and Yoshizawa Large Eddy Simulations for a Re = 4000 plane jet simulation. Comparisons are made with Direct Numerical Simulation data and measurements. Global properties of the jet such as the spreading and centreline velocity decay rates are investigated. The mean-flow and turbulence parameters in the self-similar region are also studied. All models generally give encouraging agreement with the Direct Numerical Simulation data and reliable measurements. Solution differences for the models are relatively minor, none giving clear improvements for all data comparisons.

Computation of transient and steady turbulent isothermal flows with an impinging rectangular jet

The computation of both transient and steady turbulent incompressible isothermal flows is studied. The flow is very complex, having streamline curvature, large vortex structures and stagnation resulting from an impinging rectangular jet. For transient computations, the standard k-ε model is adopted. For steady flows, the k-ε, high and low Reynolds number k-l and mixing length models are tried. Zonal approaches combining the above turbulence models are also investigated. None of the models are found to give satisfactory agreement with velocity measurements.

Semi-implicit second order finite volume jet stream computations

Semi-implicit, second order temporal and spatial finite volume computations of the flow in a differentially heated rotating annulus are presented. For the regime considered, three cyclones and anticyclones separated by a relatively fast moving jet of fluid or "jet stream" are predicted. Two second order methods are compared with, first order spatial predictions, and experimental measurements. Velocity vector plots are used to illustrate the predicted flow structure. Computations made using second order central differences are shown to agree best with experimental measurements, and to be stable for integrations over long time periods (> 1000s). No periodic smoothing is required to prevent divergence.

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.

Shroud heat transfer measurements from a rotating cavity with an axial throughflow of air

The paper discusses measurements of heat transfer obtained from the inside surface of the peripheral shroud. The experiments were carried out on a rotating cavity, comprising two 0.985-m-dia disks, separated by an axial gap of 0.065 m and bounded at the circumference by a carbon fiber shroud. Tests were conducted with a heated shroud and either unheated or heated disks. When heated, the disks had the same temperature level and surface temperature distribution. Two different temperature distributions were tested; the surface temperature either increased, or decreased with radius. The effects of disk, shroud, and air temperature levels were also studied. Tests were carried out for the range of axial throughflow rates and speeds: 0.0025 ≤ m ≤ 0.2 kg/s and 12.5 ≤ Ω ≤ 125 rad/s, respectively. Measurements were also made of the temperature of the air inside the cavity. The shroud Nusselt numbers are found to depend on a Grashof number, which is defined using the centripetal acceleration. Providing the correct reference temperature is used, the measured Nusselt numbers also show similarity to those predicted by an established correlation for a horizontal plate in air. The heat transfer from the shroud is only weakly affected by the disk surface temperature distribution and temperature level. The heat transfer from the shroud appears to be affected by the Rossby number. A significant enhancement to the rotationally induced free convection occurs in the regions 2 ≤ Ro ≤ 4 and Ro ≥ 20. The first of these corresponds to a region where vortex breakdown has been observed. In the second region, the Rossby number may be sufficiently large for the central throughflow to affect the shroud heat transfer directly. Heating the shroud does not appear to affect the heat transfer from the disks significantly.