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 1024³ 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. Program summary: Program Title: CUBic Equation of state Navier-Stokes (CUBENS) CPC Library link to program files: https://doi.org/10.17632/6jfy758gyv.1 Developer's repository link: https://github.com/pcboldini/CUBENS Licensing provisions: MIT Programming language: Fortran 90, OpenACC, MPI, Python, MATLAB Nature of problem: This code solves the three-dimensional Navier-Stokes equations for non-ideal gas flows in a Cartesian domain, applicable to boundary layers and channels. Solution method: This code uses high-order central finite-differences with split-convective form, preserving kinetic energy and entropy (KEEP) and pressure-equilibrium-preserving (PEP) property, for spatial discretization. The time advancement is performed with a third-order Total Variation Diminishing low-storage Runge-Kutta scheme. Flow non-ideality is accounted for by cubic equations of state and complex transport-properties models. Alongside MPI parallelization, the solver is GPU-accelerated using OpenACC for computation offloading and CPU-GPU data transfer, along with GPU-aware MPI for GPU-GPU communication.

@en

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)0022-112010.1017/jfm.2019.348]. 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.

@en

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.@en

Ren et al. (2019) recently studied the stability of the boundary layer flow over a flat plate for supercritical CO2. While only one unstable mode usually exists for boundary layer flows, the authors found an additional unstable mode, whose origin has so far not been identified. In the present work, we carry out a stability analysis in the general case of a fluid following the Van der Waals equation of state and flowing over a heated flat plate in the limit of zero Eckert number. In this framework, the second unstable mode is also recovered, ruling out an acoustic origin. From the Rayleigh equation derived in the presence of density gradients, a generalised inflection point (GIP) criterion of instability exists, similar to that of fully compressible flows. Inviscid stability calculations confirm the existence of an unstable mode in the presence of a GIP, which is linked to the additional second mode found at finite Reynolds numbers. A theoretical analysis is then carried out by approximating the momentum equation for a base flow exhibiting strong gradients of dynamic viscosity. It is shown that the origin of the GIP, and hence the additional unstable mode, is associated with a minimum of kinematic viscosity at the Widom line. The universality of this result beyond supercritical fluids is eventually discussed.

@en