We present STREAmS-2.0, an updated version of the flow solver STREAmS, first introduced in Bernardini et al. (2021) [1]. STREAmS-2.0 has an object-oriented design which separates the physics equations from the specific back-end, making the code more suitable for future expansions, such as porting to novel computing architectures or implementation of additional flow physics. Similarly to the previous version, STREAmS-2.0 supports NVIDIA-GPU and CPU back-ends. Additionally, this version features improvements of the input/output data management, new energy and entropy preserving schemes for the discretization of the convective fluxes, recycling/rescaling inflow boundary condition, and a model for thermally perfect gases with variable specific heats. New version program summary: Program Title: STREAmS CPC Library link to program files: https://doi.org/10.17632/hdcgjpzr3y.2 Developer's repository link: https://github.com/STREAmS-CFD/STREAmS-2 Licensing provisions: GPLv3 Programming language: Fortran, CUDA Journal reference of previous version: M. Bernardini, D. Modesti, F. Salvadore, and S. Pirozzoli. STREAmS: a high-fidelity accelerated solver for direct numerical simulation of compressible turbulent flows. Comput. Phys. Commun. 263 (2021) 107906. Does the new version supersede the previous version?: Yes. Reasons for the new version: New code structure and release of new features. Summary of revisions: • The original solver [1] has been rewritten following an object-oriented design implemented through Fortran derived types that include variables and type bound procedures. The new software architecture has been designed to increase modularity and extensibility of the code, allowing users to add new back-ends and physics equations while maintaining the same code structure. This allows users to reuse portions of the code that are independent of the physics equations, the back-end, or both. The layer of computing procedures maintains a lean structure that can be highly optimized with respect to the implemented back-end. • Input handling is now based on the classic.ini format improving both user readability and input data management. • A family of new kinetic energy and entropy preserving schemes (KEEP) are now available and can be selected for stable, non-dissipative and accurate spatial discretization of the convective terms of the Navier–Stokes equations in smooth flow regions [2]. Concerning the shock-capturing flux, the improved low-dissipative WENO-Z scheme proposed by [3] is now available. • New inflow boundary conditions based on the recycling/rescale approach [4] have been implemented for the simulation of spatially evolving compressible turbulent boundary layers. Moreover, a new inflow condition based on the solution of the compressible Blasius equation is available to take into account the case of laminar boundary layers. • The constitutive relations have been generalized to take into account thermally perfect gases with variable specific heats, approximated with polynomial functions of the temperature that can be specified by the user [5]. • A new stretching function has been implemented to improve the distribution of grid nodes for the computation of wall-bounded turbulent flows. The formulation blends uniform near-wall spacing with uniform resolution in terms of Kolmogorov units in the outer wall layer, guaranteeing accuracy with higher computational efficiency [6]. Nature of problem: The code solves the compressible Navier–Stokes equations in Cartesian coordinates for a thermally perfect gas. The solver is designed for direct numerical simulation (DNS) of compressible supersonic turbulent boundary layers and various canonical configurations are supported, including turbulent channel flow, laminar and turbulent boundary layer and shock-wave/boundary layer interaction. Solution method: The equations are discretized using high-order finite difference approximations with hybrid low-dissipative/shock-capturing capabilities and the time advancement is performed using a Runge–Kutta scheme. References: [1] M. Bernardini, D. Modesti, F. Salvadore, S. Pirozzoli, STREAmS: A high-fidelity accelerated solver for direct numerical simulation of compressible turbulent flows, Comput. Phys. Commun. 263 (2021) 107906. [2] Y. Tamaki, Y. Kuya, S. Kawai, Comprehensive analysis of entropy conservation property of non-dissipative schemes for compressible flows: KEEP scheme redefined, J. Comput. Phys. 468 (2022) 111494. [3] R. Borges, M. Carmona, B. Costa, W. Don, An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, J. Comput. Phys. 227 (6) (2008) 3191–3211, https://doi.org/10.1016/j.jcp.2007.11.038 [4] S. Pirozzoli, M. Bernardini, F. Grasso, Direct numerical simulation of transonic shock/boundary layer interaction under conditions of incipient separation, J. Fluid Mech. 657 (2010) 361–393. [5] B. J. McBride, M. J. Zehe, S. Gordon, NASA Glenn coefficients for calculating thermodynamic properties of individual species, NASA/TP 211556, NASA, 2002. [6] S. Pirozzoli, P. Orlandi, Natural grid stretching for DNS of wall-bounded flows, J. Comput. Phys. 439 (2021) 110408.

Microvortex generators are passive control devices smaller than the boundary layer thickness that energise the boundary layer to prevent flow separation with limited induced drag. In this work, we use direct numerical simulations (DNS) to investigate the effect of the Reynolds number in a supersonic turbulent boundary layer over a microramp vortex generator. Three friction Reynolds numbers are considered, up to, for fixed free stream Mach number and fixed relative height of the ramp with respect to the boundary layer thickness. The high-fidelity data set sheds light on the instantaneous and highly three-dimensional organisation of both the wake and the shock waves induced by the microramp. The full access to the flow field provided by DNS allows us to develop a qualitative model of the near wake, explaining the internal convolution of the Kelvin-Helmholtz vortices around the low-momentum region behind the ramp. The overall analysis shows that numerical results agree excellently with recent experimental measurements in similar operating conditions and confirms that microramps effectively induce a significantly fuller boundary layer even far downstream of the ramp. Moreover, results highlight significant Reynolds number effects, which in general do not scale with the ramp height. Increasing Reynolds number leads to enhanced coherence of the typical vortical structures in the field, faster and stronger development of the momentum deficit region, increased upwash between the primary vortices from the sides of the ramp - and thus increased lift-up of the wake - and faster transfer of momentum towards the wall.

We carry out direct numerical simulations (DNS) of turbulent flow and heat transfer in pressure-driven plane channels, by considering cases with heating on both walls, as well as asymmetric heating limited to one of the channel walls. Friction Reynolds numbers up to are considered, and Prandtl numbers from to, the temperature field being regarded as a passive scalar. Whereas cases with symmetric heating show close similarity between the temperature and the streamwise velocity fields, with turbulent structures confined to either half of the channel, in the presence of one-sided heating the temperature field exhibits larger regions with coherent fluctuations extending beyond the channel centreline. Validity of the logarithmic law for the mean temperature is confirmed, as well as universality of the associated von Kármán constant, which we estimate to be. Deviations from the logarithmic behaviour are much clearer in cases with one-sided heating, which feature a wide outer region with parabolic mean temperature profile. The DNS data are exploited to construct a predictive formula for the heat transfer coefficient as a function of both Reynolds and Prandtl number. We find that the reduction of the thermal efficiency in the one-sided case is approximately at unit Prandtl number; however, it can become much more significant at low Prandtl number.

We perform direct numerical simulations of turbulent flow at friction Reynolds number Re_τ≈ 500 - 2000 grazing over perforates plates with moderate viscous-scaled orifice diameter d⁺≈ 40 - 160 and analyse the relation between permeability and added drag. Unlike previous studies of turbulent flows over permeable surfaces, we find that the flow inside the orifices is dominated by inertial effects, and that the relevant permeability is the Forchheimer and not the Darcy one. We find evidence of a fully rough regime where the relevant length scale is the inverse of the Forchheimer coefficient, which can be regarded as the resistance experienced by the wall-normal flow. Moreover, we show that, for low porosities, the Forchheimer coefficient can be estimated with good accuracy using a simple analytical relation.

We carry out direct numerical simulations (DNS) of flow in a turbulent square duct by focusing on heat transfer effects, considering the case of unit Prandtl number. Reynolds numbers up to are considered that are much higher than in previous studies, and that yield clear scale separation between inner- and outer-layer dynamics. Close similarity between the behaviour of the temperature and the streamwise velocity fields is confirmed as in previous studies related to plane channels and pipes. Just like the mean velocity, the mean temperature is found to exhibit logarithmic layers as a function of the nearest wall, but with a different slope. The most important practical implication is the validity of the traditional hydraulic diameter as the correct reference length for reporting heat transfer data, as we show rigorously here. Temperature and velocity fluctuations also have similar behaviour, but apparently logarithmic growth of their inner-scaled peak variances is not observed here, unlike in canonical wall-bounded flows. Analysis of the split contributions to the heat transfer coefficient shows that mean cross-stream convection associated with secondary motions is responsible for about of the total. Finally, we use the DNS database to highlight shortcomings of traditional linear closures for the turbulent heat flux, and show that substantial modelling improvement in principle may be obtained by retaining at least the three terms in the vector polynomial integrity basis expansion.

We perform direct numerical simulation of supersonic turbulent channel flow over cubical roughness elements, spanning bulk Mach numbers -, both in the transitional and fully rough regime. We propose a novel definition of roughness Reynolds number which is able to account for the viscosity variations at the roughness crest and should be used to compare rough-wall flows across different Mach numbers. As in the incompressible flow regime, the mean velocity profile shows a downward shift with respect to the baseline smooth wall cases, however, the magnitude of this velocity deficit is largely affected by the Mach number. Compressibility transformations are able to account for this effect, and data show a very good agreement with the incompressible fully rough asymptote, when the relevant roughness Reynolds number is used. Velocity statistics present outer layer similarity with the equivalent smooth wall cases, however, this does not hold for the thermal field, which is substantially affected by the roughness, even in the channel core. We show that this is a direct consequence of the quadratic temperature-velocity relation which is also valid for rough walls. Analysis of the heat transfer shows that the relative drag increase is always larger than the relative heat transfer enhancement, however, increasing the Mach number brings data closer to the Reynolds analogy line due to the rising relevance of the aerodynamic heating.

We present STREAmS, an in-house high-fidelity solver for direct numerical simulations (DNS) of canonical compressible wall-bounded flows, namely turbulent plane channel, zero-pressure gradient turbulent boundary layer and supersonic oblique shock-wave/boundary layer interaction. The solver incorporates state-of-the-art numerical algorithms, specifically designed to cope with the challenging problems associated with the solution of high-speed turbulent flows and can be used across a wide range of Mach numbers, extending from the low subsonic up to the hypersonic regime. From the computational viewpoint, STREAmS is oriented to modern HPC platforms thanks to MPI parallelization and the ability to run on multi-GPU architectures. This paper discusses the main implementation strategies, with particular reference to the CUDA paradigm, the management of a single code for traditional and multi-GPU architectures, and the optimization process to take advantage of the latest generation of NVIDIA GPUs. Performance measurements show that single-GPU optimization more than halves the computing time as compared to the baseline version. At the same time, the asynchronous patterns implemented in STREAmS for MPI communications guarantee very good parallel performance especially in the weak scaling spirit, with efficiency exceeding 97% on 1024 GPUs. For overall evaluation of STREAmS with respect to other compressible solvers, comparison with a recent GPU-enabled community solver is presented. It turns out that, although STREAmS is much more limited in terms of flow configurations that can be addressed, the advantage in terms of accuracy, computing time and memory occupation is substantial, which makes it an ideal candidate for large-scale simulations of high-Reynolds number, compressible wall-bounded turbulent flows. The solver is released open source under GPLv3 license. Program summary: Program Title: STREAmS CPC Library link to program files: https://doi.org/10.17632/hdcgjpzr3y.1 Developer's repository link: https://github.com/matteobernardini/STREAmS Code Ocean capsule: https://codeocean.com/capsule/8931507/tree/v2 Licensing provisions: GPLv3 Programming language: Fortran 90, CUDA Fortran, MPI Nature of problem: Solving the three-dimensional compressible Navier–Stokes equations for low and high Mach regimes in a Cartesian domain configured for channel, boundary layer or shock-boundary layer interaction flows. Solution method: The convective terms are discretized using a hybrid energy-conservative shock-capturing scheme in locally conservative form. Shock-capturing capabilities rely on the use of Lax–Friedrichs flux vector splitting and weighted essentially non-oscillatory (WENO) reconstruction. The system is advanced in time using a three-stage, third-order RK scheme. Two-dimensional pencil distributed MPI parallelization is implemented alongside different patterns of GPU (CUDA Fortran) accelerated routines.

We carry out a priori tests of linear and nonlinear eddy viscosity models using direct numerical simulation (DNS) data of square duct flow up to friction Reynolds number Re _τ= 1055. We focus on the ability of eddy viscosity models to reproduce the anisotropic Reynolds stress tensor components a_ij responsible for turbulent secondary flows, namely the normal stress a₂₂ and the secondary shear stress a₂₃. A priori tests on constitutive relations for a_ij are performed using the tensor polynomial expansion of Pope (J Fluid Mech 72:331–340, 1975), whereby one tensor base corresponds to the linear eddy viscosity hypothesis and five bases return exact representation of a_ij. We show that the bases subset has an important effect on the accuracy of the stresses and the best results are obtained when using tensor bases which contain both the strain rate and the rotation rate. Models performance is quantified using the mean correlation coefficient with respect to DNS data C~ _ij, which shows that the linear eddy viscosity hypothesis always returns very accurate values of the primary shear stress a₁₂ (C~ ₁₂> 0.99), whereas two bases are sufficient to achieve good accuracy of the normal stress and secondary shear stress (C~ ₂₂= 0.911 , C~ ₂₃= 0.743). Unfortunately, RANS models rely on additional assumptions and a priori analysis carried out on popular models, including k–ε and v²–f, reveals that none of them achieves ideal accuracy. The only model based on Pope’s expansion which approaches ideal performance is the quadratic correction of Spalart (Int J Heat Fluid Flow 21:252–263, 2000), which has similar accuracy to models using four or more tensor bases. Nevertheless, the best results are obtained when using the linear correction to the v²–f model developed by Pecnik and Iaccarino (AIAA Paper 2008-3852, 2008), although this is not built on the canonical tensor polynomial as the other models.