Numerical simulations and analysis of the turbulent flow field in a practical gas turbine engine combustor

The turbulent flow field in a practical gas turbine combustor is very complex because of the interactions between various flows resulting from components like multiple types of swirlers, dilution holes, and liner effusion cooling holes. Numerical simulations of flows in such complex combustor configurations are challenging. The challenges result from (a) the complexities of the interfaces between multiple three-dimensional shear layers, (b) the need for proper treatment of a large number of tiny effusion holes with multiple angles, and (c) the requirements for fast turnaround times in support of engineering design optimization. Both the Reynolds averaged Navier–Stokes simulation (RANS) and the large eddy simulation (LES) for the practical combustor geometry are considered. An autonomous meshing using the cut-cell Cartesian method and adaptive mesh refinement (AMR) is demonstrated for the first time to simulate the flow in a practical combustor geometry. The numerical studies include a set of computations of flows under a prescribed pressure drop across the passage of interest and another set of computations with all passages open with a specified total flow rate at the plenum inlet and the pressure at the exit. For both sets, the results of the RANS and the LES flow computations agree with each other and with the corresponding measurements. The results from the high-resolution LES simulations are utilized to gain fundamental insights into the complex turbulent flow field by examining the profiles of the velocity, the vorticity, and the turbulent kinetic energy. The dynamics of the turbulent structures are well captured in the results of the LES simulations.

Secondary Flow in Smooth and Rough Turbulent Circular Pipes: Turbulence Kinetic Energy Budgets

Direct Numerical Simulations have been performed for turbulent flow in circular pipes with smooth and corrugated walls. The numerical method, based on second-order finite discretization together with the immersed boundary technique, was validated and applied to various types of flows. The analysis is focused on the turbulence kinetic energy and its budget. Large differences have been found in the near-wall region at low Reynolds number. The change in the near-wall turbulent structures is responsible for increase of drag and turbulence kinetic energy. To investigatselinae the effects of wall corrugations, the velocity fields have been decomposed so as to isolate coherent and incoherent motions. For corrugated walls, we find that coherent motions are strongest for walls covered with square bars aligned with the flow direction. In particular, the coherent contribution is substantial when the bars are spaced apart by a distance larger than their height. Detailed analysis of the turbulence kinetic energy budget shows for this set-up a very different behavior than for the other types of corrugations.

Vortex Flow on the Surface Generated by the Onset of a Buoyancy-Induced Non-Boussinesq Convection in the Bulk of a Normal Liquid Helium

The onset of the Rayleigh–Benard convection (RBC) in a heated from above normal He-I layer in a cylindrical vessel in the temperature range Tλ < T ≤ Tm (RBC in non-Oberbeck–Boussinesq approximation) is attended by the emergence of a number of vortices on the free liquid surface. Here, Tλ = 2.1768 K is the temperature of the superfluid He-II–normal He-I phase transition, and the liquid density passes through a well-pronounced maximum at Tm ≈ Tλ + 6 mK. The inner vessel diameter was D = 12.4 cm, and the helium layer thickness was h ≈ 2.5 cm. The mutual interaction of the vortices between each other and their interaction with turbulent structures appeared in the layer volume during the RBC development gave rise to the formation of a vortex dipole (two large-scale vortices) on the surface. Characteristic sizes of the vortices were limited by the vessel diameter. The formation of large-scale vortices with characteristic sizes twice larger than the layer thickness can be attributed to the arising an inverse vortex cascade on the two-dimensional layer surface. Moreover, when the layer temperature exceeds Tm, convective flows in the volume decay. In the absence of the energy pumping from the bulk, the total energy of the vortex system on the surface decreases with time according to a power law.

Effect of Local Turbulent Structures on the Hot Jet in Transitional Flow

Turbulent parts localized in flow direction may arise in a pipe with transitional regime of the stable laminar Poiseuille flow. A key condition for occurrence of such structures is a pipe with rather long length relative to its diameter. Our paper presents numerical modelling of the hot air jet flowing from the long pipe into the cold open volume at Re=2426. The modelling was performed in OpenFOAM software on the basis of the large eddy simulation (LES) method. The WALE (Wall-adapting local eddy-viscosity) model was used for closure of Navier-Stokes equations on subgrid scales. We demonstrated that local turbulent structures have a weak effect on the hot jet at flowing into the cold open volume.

Destructiveness of pyroclastic surges controlled by turbulent fluctuations

Pyroclastic surges are lethal hazards from volcanoes that exhibit enormous destructiveness through dynamic pressures of 100–102 kPa inside flows capable of obliterating reinforced buildings. However, to date, there are no measurements inside these currents to quantify the dynamics of this important hazard process. Here we show, through large-scale experiments and the first field measurement of pressure inside pyroclastic surges, that dynamic pressure energy is mostly carried by large-scale coherent turbulent structures and gravity waves. These perpetuate as low-frequency high-pressure pulses downcurrent, form maxima in the flow energy spectra and drive a turbulent energy cascade. The pressure maxima exceed mean values, which are traditionally estimated for hazard assessments, manifold. The frequency of the most energetic coherent turbulent structures is bounded by a critical Strouhal number of ~0.3, allowing quantitative predictions. This explains the destructiveness of real-world flows through the development of c. 1–20 successive high-pressure pulses per minute. This discovery, which is also applicable to powder snow avalanches, necessitates a re-evaluation of hazard models that aim to forecast and mitigate volcanic hazard impacts globally.

Formation of the radial electric field profile in the WEST tokamak

Sheared flows are known to reduce turbulent transport by decreasing the correlation length and/or intensity of turbulent structures. The transport barrier that takes place at the edge during improved regimes such as H mode, corresponds to the establishment of a large shear of the radial electric field. In this context, the radial shape of the radial electric field or more exactly of the perpendicular $E\times B$ velocity appears as a key element in accessing improved confinement regimes. In this paper, we present the radial profile of the perpendicular velocity measured using Doppler back-scattering system at the edge of the plasma, dominated by the $E\times B$ velocity, during the first campaigns of the WEST tokamak. It is found that the radial velocity profile is clearly more sheared in LSN than in USN configuration for ohmic and low current plasmas ($B=3.7T$ and $q_{95}=4.7$), consistently with the expectation comparing respectively “favourable” versus “unfavourable” configuration. Interestingly, this tendency is sensitive to the plasma current and to the amount of additional heating power leading to plasma conditions in which the $E\times B$ velocity exhibits a deeper well in USN configuration. For example, while the velocity profile exhibits a clear and deep well just inside the separatrix concomitant with the formation of a density pedestal during L-H transitions observed in LSN configuration, deeper $E_r$ wells are observed in USN configuration during similar transitions with less pronounced density pedestal.

Turbulent Characteristics and Air Entrainment Patterns in Breaking Surge Waves

Breaking surge waves are highly turbulent three-dimensional (3D) flows, which occur when the water flow encounters a sudden change in depth or velocity. The 3D turbulent structures across a breaking surge are induced by the velocity gradient across the surge and phase discontinuity at the front. This paper examined the turbulent structures in breaking surge waves with Froude numbers of 1.71 and 2.13 by investigating the air entrainment and perturbation patterns across the surge front. A combination of the Volume Of Fluid (VOF) method and Large Eddy Simulation (LES) was utilized to capture air entrainment and turbulent structures simultaneously. The 3D nature of the vortical structures was simulated by implementing a spanwise periodic boundary. The water surface perturbation and air concentration profiles were extracted, and the averaged air concentration profiles obtained from the numerical simulations were consistent with laboratory observations reported in the literature. The linkage between turbulent kinetic energy distribution and air entrainment was also explored in this paper. Finally, using quadrant analysis and the Q-criterion, this paper examined the role of the spanwise perturbations in the development of turbulent structures in the surge front.

Dynamic LES of the magnetohydrodynamic flow in a square duct with the varied wall conductance parameters

The current study is focused on the magnetohydrodynamics and demonstrates how electrical conductivity of the wall can affect the turbulent flow in the square duct. Different variations of the boundary walls have been considered including arbitrary conductive walls. The Large Eddy Simulations method with the dynamic Smagorinsky sub-grid scale model have been used for the turbulent structures resolving. Results show the significant impact of the wall conductance parameters for both Hartmann and side walls.