A detailed derivation, analysis, and verification is given for the non-orthogonal, plane-marching Parabolized Stability Equations (PSE) approach. In applying the approach to a flow distorted by a medium-amplitude crossflow vortex, we determine its linear secondary instability mechanisms. We show that converged solutions can be achieved for a broad frequency range with an existing stabilization method for the line-marching PSE approach. We verify that 1) solutions converge versus grid size in all dimensions, 2) primary disturbance solutions agree with line-marching PSE results, and 3) secondary disturbance solutions match amplitude and growth-rate evolution of reference Direct Numerical Simulation (DNS) results. We show how and why the type-II instability displays a delayed neutral point when modeled with the plane-marching approach versus the considered local stability approaches, whether the streamwise evolution of the distorted base flow is accounted for or not. This may explain why the type-II disturbance is scarcely captured by DNS in the literature.

A comprehensive set of accurate physical models and numerical simulation methods for transcritical dual fuel combustion systems is presented. The method combines multiphase real-fluid physical properties modeling, flamelet-based chemistry reduction, and large-eddy simulation (LES). It allows a fully integrated simulation of the multiphase region and the combustion zones in the turbulent flames. A central challenge in dual-fuel combustion applications is the successful ignition of low-reactivity ambient fuel by the high-reactivity pilot fuel jet. The interaction of the radical pools associated with both fuels is known to be an important characteristic. This study captures the evolution of this interaction across both the low-temperature combustion (LTC) and high-temperature combustion (HTC) regimes. Three cases are investigated, each involving the injection of an n-dodecane jet into a 6 MPa, 1000 K ambient. The cases differ only in ambient composition: (1) a single-fuel (SF) baseline case with a mixture of air and combustion products, and two dual-fuel (DF) cases with methane and a 90:10 (by volume) methane/ethane blend added to the ambient. Insights on ignition and flame structure are gained from homogeneous reactor simulations, transcritical counterflow diffusion flames, and LES. Combustion modes are defined by threshold levels of key species and temperature. The flame volume is subdivided according to the modes and for each mode the ignition delay and the temporal evolution of heat release rate (HRR) are obtained. The differences in the radical pools of SF and DF combustion are demonstrated. LTC ignition is shown to start in the multiphase region near the base of the jet. The high accuracy of the proposed methods provides a firm basis for modeling prospective dual-fuel systems using alternative fuels.

High-speed supersonic radial compressors are a critical enabling technology for meeting the requirements of future aviation-propulsion and thermal-management systems. These turbomachines must be designed to be both efficient and robust on the widest possible operating range. Flow instabilities in the form of rotating stall and surge are therefore phenomena that must be accurately predicted early in the design process. Unsteady full-Annulus computational fluid dynamics (CFD) can be used to get accurate information about the onset of instabilities, but at the expense of costly simulations. As a result, the design of new compressors continues to rely on existing correlations for the prediction of the critical mass flowrate. This approach, however, leads to suboptimal compressor designs. This article provides a review of the numerical methodologies that can be used for the accurate prediction of the critical mass flowrate in high-speed centrifugal compressors. Methods of different fidelity level and computational cost are described. Two particularly promising models, namely, those proposed by Spakovszky and Sun, are subsequently examined in more detail. Exemplary applications of these two models are finally discussed.

Leading-edge protuberances on airfoils have been shown to soften the onset of aerodynamic stall and to increase lift in the post-stall regime. The present study examines the effect of tubercles during dynamic stall. Pitching airfoils with tubercles of different amplitudes are studied by wind-tunnel experiments, where the three-dimensional time-resolved velocity field is determined using large-scale particle-tracking velocimetry. Computational fluid dynamics simulations are carried out that complement the experimental observations providing pressure distribution and aerodynamic forces. The dynamic stall is dominated by a vortex formed at the leading edge; we characterize the vorticity, circulation, and advection path of this dynamic-stall vortex (DSV). The presence of the tubercles profoundly modifies the boundary layer from the leading edge. The roll-up of the vorticity sheet is significantly delayed compared to a conventional airfoil, resulting in a weaker DSV. The vortex formation is shifted downstream, with the overall effect of a weaker and shorter lift overshoot, in turn enabling a quicker transition to deep stall. Regions of flow separation (stall cells) are visibly compartmentalized with a stable spacing of two tubercles wavelengths.

Simulations of reacting multiphase flows tend to display an inhomogeneously distributed computational intensity over the spatial and temporal domains. The time-to-solution of chemical reaction rates can span multiple orders of magnitude due to the emergence of combustible kernels and thin turbulent reaction zones. Similarly, the time to solve the equation of state (EoS) for non-ideal fluid mixtures deviates substantially between the grid cells. These effects result in a performance profile that is unbalanced and rapidly changing for transient simulations, and therefore beyond the capabilities of traditional (quasi-)static mesh partitioning methods. We analyse this loss of parallel efficiency for large-eddy simulations of the ECN Spray-A benchmark with the multi-physics solver INCA and propose to mitigate the problem by introducing two independent repartitioning stages in addition to the classic domain decomposition for fluid transport: one for the EoS and one for chemical reactions. We explore various scalable repartitioning strategies in this context and observe that rebalancing computational load yields a significant speedup that is robust for various mesh resolutions and process numbers. The dynamic multistage load-balancing thus effectively removes obstacles towards good parallel scaling of INCA and similar solvers for reacting and/or multiphase flows.

Direct numerical simulations (DNS) are conducted for reactants-to-products counterflow configurations at turbulent conditions to understand how strain affects the structure and NOx emissions of lean premixed hydrogen flames. Two nominal equivalence ratio conditions, 0.5 and 0.7, are investigated. Under unstretched conditions, the Markstein length is negative for the former and slightly positive for the latter, indicating distinct responses of heat release rate and flame consumption speed to strain in each case. For each equivalence ratio condition, three levels of applied strain rate are considered, resulting in a total of six DNS. Results indicate that overall NOx emissions decrease with increasing strain at turbulent conditions, consistent with recent results for laminar conditions presented in Porcarelli et al. (2024). However, the relative decrease of NOx with strain is faster under turbulent conditions because turbulent mixing limits the occurrence of super-adiabatic temperatures. Moreover, the decrease of NOx is strongly correlated only to the mean applied tangential strain rate, while local fluctuations of strain due to vortices exhibit more stochastic behaviour. The detailed analysis presented in this article indicates that the applied strain can be used to substantially decrease NOx emissions in premixed hydrogen flames under practical conditions. Novelty and Significance statement: This work examines for the first time in detail the coupled effects of strain and turbulence in hydrogen flames, for various conditions spanning different signs of the Markstein length and increasing applied strain levels. In particular, it clarifies the different roles of applied strain, turbulence-driven strain, and curvature on both flame structure and NOx generation. Results further show for the first time that both in-flame and post-flame NOx can be suppressed at high strain levels under turbulent conditions. This result is of paramount importance as it implies that NOx can be suppressed at combustor-relevant conditions by straining the flame.

We present novel observations from direct numerical simulations of transitional Mach 8 flow over a 15° compression ramp ablator. Heating streaks over the ramp are seen to undergo a half-wavelength shift near the location of transition from laminar to turbulent boundary layer flow. This phenomenon leads to an intriguing pattern of ablation grooves on the surface. Our analysis shows that the underlying mechanism is driven by the baroclinic torque in the strongly stratified near-wall region. We discuss the impact of this baroclinic shift for a surface undergoing ablative recession and assess its sensitivity to different thermal boundary conditions and perturbation amplitudes.

The nacelle of aircraft engines is coated with acoustic liners to reduce engine noise emissions. An undesirable side effect of acoustic liners is that they increase aerodynamic drag. For the first time, the authors study this drag penalty through pore-resolved direct numerical simulation (DNS) of a flat-plate zero pressure gradient turbulent boundary layer at friction Reynolds number Re_τ ≈ 850–2600, which is high enough to be representative of liners in operating conditions. In the configuration under scrutiny, the turbulent boundary layer experiences a step change in surface topography passing from a smooth wall to an acoustic liner array, allowing one to study the streamwise adaptation length of the boundary layer. It is found that the mean velocity profile adjusts to the new surface condition in a nearly negligible distance (less than 10 local boundary-layer thicknesses), whereas turbulent fluctuations take much longer. DNS is also performed with external acoustic noise in the form of planar monocromatic waves grazing the boundary layer with an amplitude of 150 dB. In agreement with some earlier studies, it is found that sound waves do not affect aerodynamic drag at these flow conditions.

This study introduces a new numerical framework for the accurate simulation of transcritical reacting sprays using a multiphase, real-fluid, flamelet-based model. The transcritical flamelet library is combined with large-eddy simulations (LES) and rapid vapor–liquid equilibrium calculations in the context of a modern multiphase thermodynamic approach to explore vaporization dynamics, ignition characteristics, and soot formation. Current applications focus on the combustion of polyoxymethylene dimethyl ethers (OMEs), which are carbon-neutral e-fuels, in transcritical high-pressure configurations. Validation against experimental data shows a strong match in ignition delay and penetration lengths. The analysis of three OME₃– n-dodecane fuel blends reveals differences in evaporation, ignition, and soot production. Adding OME₃ to n-dodecane reduces soot production and shortens the liquid penetration length and ignition delay time. The findings highlight the importance of further investigation into the effects of transcritical states and fuel composition on combustion performance and emissions. Novelty and significance This work introduces a modeling technique for the use of transcritical counterflow flames in flamelet modeling, expanding the capabilities of large-eddy simulations with multiphase thermodynamics (LES-MT) to accurately modeling transcritical combustion. By incorporating real-fluid effects and two-phase interactions, the transcritical flamelet library provides a high-fidelity representation of the complex behaviors in high-pressure multiphase autoignition scenarios. This calibration-free approach can significantly improve our understanding of the transcritical combustion of emerging fuels such as OME₃ or their combination with traditional fuels such as n-dodecane.

Stationary velocity-perturbation streaks have recently been identified in laminar swept-wing boundary-layer flow interacting with a surface forward-facing step. Streaky structures at the step promote early laminar-turbulent transition under certain conditions. This work utilizes direct numerical simulations to explore the mechanisms of growth of stationary streaks at the step and provides insight into their origin, nature, and spatial organization. The analysis is mainly focused on, but not restricted to, incoming perturbations in the form of stationary crossflow instability. Stationary streaky structures are found to be universal to swept forward-facing-step flow subjected to three-dimensional perturbations in the incoming boundary layer. The streaks at the step are primarily ascribed to the lift-up effect. They appear as a linear perturbation response of the highly sheared step flow to the cross-stream pattern of incoming perturbations. A mechanism of base-flow deceleration additionally contributes to feeding growth to the streaks. Linear stability analysis carried out through the harmonic Navier-Stokes method confirms that the streaks are a linear perturbation phenomenon. Effects of spanwise perturbation wavelength and effective sweep angle on the mechanisms of the streaks are also assessed.

High-speed supersonic radial compressors are a critical enabling technology for meeting the requirements of future aviation-propulsion and thermal-management systems. These turbomachines must be designed to be both efficient and robust on the widest possible operating range. Flow instabilities in the form of rotating stall and surge are therefore phenomena that must be accurately predicted early in the design process. Unsteady full-annulus computational fluid dynamics can be used to get accurate information about the onset of instabilities, but at the expense of costly simulations. As a result, the design of new compressors continues to rely on existing correlations for the prediction of the critical mass flow rate. This approach, however, leads to sub-optimal compressor designs. This article provides a review of the numerical methodologies that can be used for the accurate prediction of the critical mass flow rate in high-speed centrifugal compressors. Methods of different fidelity level and computational cost are described. Two particularly promising models, namely those proposed by Spakovszky and Sun, are subsequently examined in more detail. Exemplary applications of these two models are finally discussed.

Pore-resolved direct numerical simulations (DNS) of turbulent flows grazing over acoustic liners with aerodynamically and/or acoustically optimized orifice configurations are presented. The DNS explore a large parameter space, studying various families of orifice geometries, including the influence of orifice shape, orientation, and the number of orifices. All flow cases show an increase in drag compared to the smooth wall. However, the added drag can be reduced by as much as approximately 55% as compared to conventional acoustic liners by simply altering the shape of the orifice or its orientation, in the case of a noncircular orifice. Complementary acoustic simulations demonstrate that this reduced drag may be achieved while maintaining the same noise reduction properties over a wide range of frequencies.

We revisit the origin of low-frequency unsteadiness in turbulent recirculation bubbles (TRBs), and, in particular, the hypothesis of a dynamic feedback mechanism between unconstrained separation and reattachment locations. To this end, we conduct wall-resolved large-eddy simulations of a novel experimental configuration where a shock-induced TRB forms over a backward-facing step that is intended to intercept the hypothesized dynamic feedback. Our results demonstrate, for the first time, effective suppression of the low-frequency characteristics of the TRB without reducing its size, strongly supporting our hypothesis.

We propose several enhancements to improve the accuracy and performance of the digital filter turbulent inflow generation technique and assess their efficacy in the context of wall-resolved large-eddy simulations of a compressible turbulent boundary layer. Improvements of accuracy include a more realistic correlation function for the transversal directions, target length scales that vary with wall-distance, and a counter-intuitive approach that involves the suppression of streamwise velocity fluctuations at the inflow. For improving the computational performance, we propose to generate the inflow data in parallel in single precision and at a prescribed time interval based on the turbulence time scale, and not at every time-step of the simulation. Based on the results of 7 wall-resolved large-eddy simulations, we find that the new correlation functions and the considered performance improvements are beneficial and therefore desired. Suppressing streamwise velocity fluctuations at the inflow leads to the fastest relaxation of the pressure fluctuations; however, this approach increases the adaptation length defined in terms of compliance with the von Kármán integral equation. The adaptation length can be shortened by artificially increasing the wall-normal Reynolds stresses, thereby preserving the desired turbulence kinetic energy level. A detailed inspection of the Reynolds stress transport budgets reveals that the observed spurious spatial transients are largely driven by pressure-related terms. For instance, increased values of u^′p^′¯ are found throughout the computational domain when a physical Reynolds stress distribution is prescribed at the inflow. Therefore, efforts to enhance digital filter techniques should aim at modeling pressure fluctuations as well as their correlation with the velocity components.

We present pore-resolved direct numerical simulations (DNS) of turbulent flows grazing over acoustic liners with aerodynamically and/or acoustically optimised orifice configurations. Our DNS explore a large parameter space, studying different families of orifice geometries including the influence of orifice shape, orientation, and number. All flow cases show an increase in drag compared to the smooth wall. However, the added drag can be reduced by as much as ∼60%, as compared to conventional acoustic liners by simply changing the shape of the orifice or the orientation in the case of a non-circular orifice. Complementary acoustic simulations show that this drag reduction can be achieved while retaining the same noise reduction properties.

Wall-resolved large-eddy simulations (LES) are performed to investigate Reynolds number effects in supersonic turbulent boundary layers (TBLs) at Mach 2.0. The resulting database covers more than a decade of friction Reynolds number Re_τ, from 242 to 5554, which considerably extends the parameter range of current high-fidelity numerical studies. Reynolds number trends are identified on a variety of statistics for skin-friction, velocity and thermodynamic variables. The efficacy of recent scaling laws as well as compressibility effects are also assessed. In particular, we observe the breakdown of Morkovin's hypothesis for third-order velocity statistics, in agreement with previous observations for variable-property flows at low Mach number. Special attention is also placed on the size and topology of the turbulent structures populating the TBL, with an emphasis on the outer-layer motions at high Reynolds number. The corresponding streamwise spectra of streamwise velocity fluctuations show a clear separation between inner and outer scales, where energetic peaks are found at streamwise wavelengths of λ_x⁺≈700 and λ_x/δ₀≈6. The spanwise spacing of the outer-layer structures, in turn, is found to be insensitive to the Reynolds number and equal to ∼0.7δ₀. It is also found that the integral length-scales in spanwise direction for the temperature, streamwise and spanwise velocity fields appear to progressively collapse with increasing Reynolds number. The modulating influence that the outer-layer structures exert on the near-wall turbulence is also clearly visible in many of the metrics discussed. In addition, the present LES data is further exploited to assess the Re_τ-sensitivity of uniform momentum regions in the flow. We find that the resulting probability density function of the number of zones as well as its evolution with Re_τ agrees well with incompressible data. This suggests that uniform zones, which have been associated with outer-layer dynamics, are not strongly influenced by compressibility at the considered Mach number.

A novel passive flow-control method for shock-wave/turbulent boundary-layer interactions (STBLI) is investigated. The method relies on a structured roughness pattern constituted by streamwise-aligned ridges. Its effectiveness is assessed with wall-resolved large-eddy simulations of the interaction of a Mach 2 turbulent boundary layer flow with an oblique impinging shock with shock angle 40^∘. The structured roughness pattern, which is fully resolved by a cut-cell based immersed boundary method, covers the entire computational domain. Results show that this rough surface induces large-scale secondary streamwise vortices, which energize the boundary layer by transporting high-speed fluid closer to the wall. A parametric study is performed to investigate the effect of the spacing between the ridges. This investigation is further substantiated through spectral analysis and sparsity-promoting dynamic mode decomposition. We find that ridges with small spacing effectively mitigate the low-frequency unsteadiness of STBLI and slightly reduce total-pressure loss.

A novel mechanism is identified, through which a spanwise-invariant surface feature (a two-dimensional forward-facing step) significantly stabilizes the stationary crossflow instability of a three-dimensional boundary layer. The mechanism is termed here as reverse lift-up effect, inasmuch as it acts reversely to the classic lift-up effect; that is, kinetic energy of an already-existing shear-flow instability is transferred to the underlying laminar flow through the action of cross-stream perturbations. To characterize corresponding energy-transfer mechanisms, a theoretical framework is presented, which is applicable to generic three-dimensional flows and surface features of arbitrary shape with one invariant spatial direction. The identification of a passive geometry-induced effect responsible for dampening stationary crossflow vortices is a promising finding for laminar flow control applications.

We investigate Reynolds number effects in strong shock-wave/turbulent boundary-layer interactions (STBLI) by leveraging a new database of wall-resolved and long-integrated large-eddy simulations. The database encompasses STBLI with massive boundary-layer separation at Mach 2.0, impinging-shock angle 40^◦ and friction Reynolds numbers Re_τ 355, 1226 and 5118. Our analysis shows that the shape of the reverse-flow bubble is notably different at low and high Reynolds number, while the mean-flow separation length, separation-shock angle and incipient plateau pressure are rather insensitive to Reynolds number variations. Velocity statistics reveal a shift in the peak location of the streamwise Reynolds stress from the separation-shock foot to the core of the detached shear layer at high Reynolds number, which we attribute to increased pressure transport in the separation-shock excursion domain. Additionally, in the high Reynolds case, the separation shock originates deep within the turbulent boundary, resulting in intensified wall-pressure fluctuations and spanwise variations associated with the passage of coherent velocity structures. Temporal spectra of various signals show energetic low-frequency content in all cases, along with a distinct peak in the bubble-volume spectra at a separation-length-based Strouhal number St_Lsep ≈ 0.1. The separation shock is also found to lag behind bubble-volume variations, consistent with the acoustic propagation time from reattachment to separation and a downstream mechanism driving the shock motion. Finally, dynamic mode decomposition of three-dimensional fields suggests a Reynolds-independent statistical link among separation-shock excursions, velocity streaks and large-scale vortices at low frequencies.