Pietro Carlo Boldini

The objective of this work is to investigate the unexplored laminar-to-turbulent transition of a heated flat-plate boundary layer with a fluid at supercritical pressure. Two temperature ranges are considered: a subcritical case, where the fluid remains entirely in the liquid-like regime, and a transcritical case, where the pseudo-critical (Widom) line is crossed and pseudo-boiling occurs. Fully compressible direct numerical simulations are used to study (i) the linear and nonlinear instabilities, (ii) the breakdown to turbulence, and (iii) the fully developed turbulent boundary layer. In the transcritical regime, two-dimensional forcing generates not only a train of billow-like structures around the Widom line, resembling Kelvin-Helmholtz instability, but also near-wall travelling regions of flow reversal. These spanwise-oriented billows dominate the early nonlinear stage. When high subharmonic three-dimensional forcing is applied, staggered Λ-vortices emerge more abruptly than in the subcritical case. However, unlike the classic H-type breakdown under zero pressure gradient observed in ideal-gas and subcritical regimes, the H-type breakdown is triggered by strong shear layers caused by flow reversals – similar to that observed in adverse-pressure-gradient boundary layers. Without oblique wave forcing, transition is only slightly delayed and follows a naturally selected fundamental breakdown (K-type) scenario. Hence, in the transcritical regime, it is possible to trigger nonlinearities and achieve transition to turbulence relatively early using only a single two-dimensional wave that strongly amplifies background noise. In the fully turbulent region, we demonstrate that variable-property scaling accurately predicts turbulent skin-friction and heat-transfer coefficients.

We present a massively parallel GPU-accelerated solver for direct numerical simulations of transitional and turbulent flat-plate boundary layers and channel flows involving fluids in non-ideal thermodynamic states. While several high-fidelity solvers are currently available as open source, all of them are restricted to the ideal-gas region. In contrast, the CUBic Equation of state Navier-Stokes solver (CUBENS) can accurately model and simulate the non-ideal thermodynamics of single-phase compressible fluids in the vicinity of the vapor-liquid saturation line or the thermodynamic critical point. By employing high-order finite-difference schemes and convective terms in split, kinetic-energy-, and entropy-preserving form, the solver is numerically stable, and robust with minimal numerical dissipation, enabling it to capture the steep variations of non-ideal thermodynamic properties. For cost-effective high-fidelity simulations, in addition to MPI parallelization, CUBENS is GPU-accelerated using OpenACC directives for computation offloading, and asynchronous GPU-aware MPI for efficient GPU-GPU communication. Moreover, CUBENS is compatible with both NVIDIA and AMD GPU architectures, achieving significant performance results while ensuring energy-efficient simulations. For instance, using 64 NVIDIA A100 GPUs compared to 8192 CPUs at the same computational cost results in a speedup of approximately 130×. In multi-node and multi-GPU configurations ranging from 2 to 128 compute nodes (8 to 512 GPUs), a strong scaling efficiency of around 52% and a weak scaling efficiency of 0.88 with 10243 points per GPU, corresponding to approximately 5 billion degrees of freedom, are achieved. The CUBENS solver is validated against selected cases from the literature, covering transitional to turbulent ideal and non-ideal flows up to the transonic regime. In particular, we demonstrate the solver’s suitability and applicability for direct numerical simulations of transitional boundary layers with fluids at supercritical pressure and with buoyancy effects. The development of this high-fidelity solver offers the potential for future fundamental research in non-ideal compressible fluid dynamics.

This paper addresses the stability of plane Couette flow in the presence of strong density and viscosity stratifications. It demonstrates the existence of a generalised inflection point that satisfies the generalised Fjørtoft criterion of instability when a minimum of kinematic viscosity is present in the base flow. The characteristic scales associated with this minimum are identified as the primary controlling parameters of the associated instability, regardless of the type of stratification. To support this finding, analytical stability models are derived in the long-wave approximation using piecewise linear base flows. Numerical stability calculations are carried out to validate these models and to provide further information on the production of disturbance vorticity. All instabilities are interpreted as arising from the interaction between two vorticity waves. Depending on the type of stratification, these two waves are produced by different physical mechanisms. When both strong density and viscosity stratifications are present, we show that they result from the concurrent action of shear and inertial baroclinic effects. The stability models developed for simple fluid models ultimately shed light on a recently observed unstable mode in supercritical fluids (Ren et al., J. Fluid Mech., vol. 871, 2019, pp. 831–864), providing a quantitative prediction of the stability diagram and identifying the dominant mechanisms at play. Furthermore, our study suggests that the minimum of kinematic viscosity reached at the Widom line in these fluids is the leading cause of their instability. The existence of similar instabilities in different fluids and flows (e.g. miscible fluids) is finally discussed.

In the region close to the thermodynamic critical point and in the proximity of the pseudoboiling (Widom) line, strong property variations substantially alter the growth of modal instabilities, as revealed in Ren et al. [J. Fluid Mech. 871, 831 (2019)]. Here, we study nonmodal disturbances in the spatial framework using an eigenvector decomposition of the linearized Navier-Stokes equations under the assumption of locally parallel flow. To account for nonideality, a new energy norm is derived. Several heat transfer scenarios at supercritical pressure are investigated, which are of practical relevance in technical applications. The boundary layers with the fluid at supercritical pressure are heated or cooled by prescribing the wall and free-stream temperatures so that the temperature profile is either entirely subcritical (liquidlike), supercritical (gaslike), or transcritical (across the Widom line). The free-stream Mach number is set to 10^−3. In the nontranscritical regimes, the resulting streamwise-independent streaks originate from the lift-up effect. Wall cooling enhances the energy amplification for both subcritical and supercritical regimes. When the temperature profile is increased beyond the Widom line, a strong suboptimal growth is observed over very short streamwise distances due to the Orr mechanism. Due to the additional presence of transcritical Mode II, the optimal energy growth at large distances is found to arise from an interplay between lift-up and Orr mechanism. As a result, optimal disturbances are streamwise-modulated streaks with strong thermal components and with a propagation angle inversely proportional to the local Reynolds number. The nonmodal growth is put in perspective with modal growth by means of an N-factor comparison. In the nontranscritical regimes, modal stability dominates regardless of a wall-temperature variation. In contrast, in the transcritical regime, nonmodal N factors are found to resemble the imposition of an adverse pressure gradient in the ideal-gas regime. When cooling beyond the Widom line, optimal growth is greatly enhanced, yet strong inviscid instability prevails. When heating beyond the Widom line, optimal growth could be sufficiently large to favor transition, particularly with a high free-stream turbulence level.