Exascale High Performance Computing (HPC) represents a tremendous opportunity to push the boundaries of Computational Fluid Dynamics (CFD), but despite the consolidated trend towards the use of Graphics Processing Units (GPUs), programmability is still an issue. STREAmS-2 (Bernardini et al. Comput. Phys. Commun. 285 (2023) 108644) is a compressible solver for canonical wall-bounded turbulent flows capable of harvesting the potential of NVIDIA GPUs. Here we extend the already available CUDA Fortran backend with a novel HIP backend targeting AMD GPU architectures. The main implementation strategies are discussed along with a novel Python tool that can generate the HIP and CPU code versions allowing developers to focus their attention only on the CUDA Fortran backend. Single GPU performance is analyzed focusing on NVIDIA A100 and AMD MI250x cards which are currently at the core of several HPC clusters. The gap between peak GPU performance and STREAmS-2 performance is found to be generally smaller for NVIDIA cards. Roofline analysis allows tracing this behavior to unexpectedly different computational intensities of the same kernel using the two cards. Additional single-GPU comparisons are performed to assess the impact of grid size, number of parallelized loops, thread masking and thread divergence. Parallel performance is measured on the two largest EuroHPC pre-exascale systems, LUMI (AMD GPUs) and Leonardo (NVIDIA GPUs). Strong scalability reveals more than 80% efficiency up to 16 nodes for Leonardo and up to 32 for LUMI. Weak scalability shows an impressive efficiency of over 95% up to the maximum number of nodes tested (256 for LUMI and 512 for Leonardo). This analysis shows that STREAmS-2 is the perfect candidate to fully exploit the power of current pre-exascale HPC systems in Europe, allowing users to simulate flows with over a trillion mesh points, thus reducing the gap between the Reynolds numbers achievable in high-fidelity simulations and those of real engineering applications.

@en

We develop a theoretical framework for predicting friction and heat transfer coefficients in variable-properties forced-air convection, as typical of turbine blade cooling. To do this, we borrow concepts from high-speed wall turbulence, also featuring large temperature and density variations. Using the mean momentum balance and mean thermal balance equations we develop integral transformations that account for the effect of the variable fluid properties, and apply the inverse transformations to calculate the friction and heat transfer coefficients. The proposed theory is validated using a direct numerical simulation dataset spanning both heating and cooling conditions, and the predicted friction and heat transfer coefficients match DNS data with 1–2% accuracy.

@en

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

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

We derive explicit formulas for the mean profiles of temperature (modeled as a passive scalar) in forced turbulent convection, as a function of the Reynolds and Prandtl numbers. The derivation leverages on the observed universality of the inner-layer thermal eddy diffusivity with respect to Reynolds and Prandtl number variations and across different flows, and on universality of the passive scalar defect in the core flow. Matching of the inner- and outer-layer expression yields a smooth compound mean temperature profile. We find excellent agreement of the analytical profile with data from direct numerical simulations of pipe and channel flows under various thermal forcing conditions, and over a wide range of Reynolds and Prandtl numbers.

@en

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.

@en

We present pore-resolved compressible direct numerical simulations of turbulent flows grazing over perforated plates, that closely resemble the acoustic liners found on aircraft engines. Our direct numerical simulations explore a large parameter space including the effects of porosity, thickness and viscous-scaled diameter of the perforated plates, at friction Reynolds numbers, which allows us to develop a robust theory for estimating the added drag induced by acoustic liners. We find that acoustic liners can be regarded as porous surfaces with a wall-normal permeability and that the relevant length scale characterizing their added drag is the inverse of the wall-normal Forchheimer coefficient. Unlike other types of porous surfaces featuring Darcian velocities inside the pores, the flow inside the orifices of acoustic liners is fully turbulent, with a magnitude of the wall-normal velocity fluctuations comparable to the peak in the near-wall cycle. We provide clear evidence of a fully rough regime for acoustic liners, also confirmed by the increasing relevance of pressure drag. Once the fully rough asymptote is reached, canonical acoustic liners provide an added drag comparable to that of sand-grain roughness with viscous-scaled height matching the inverse of the viscous-scaled Forchheimer permeability of the plate.

@en

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.

@en

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.

@en

This paper is associated with a video winner of a 2022 American Physical Society's Division of Fluid Dynamics (DFD) Gallery of Fluid Motion (GFM) Award for work presented at the DFD Gallery of Fluid Motion. The original video is available online at the Gallery of Fluid Motion, https://doi.org/10.1103/APS.DFD.2022.GFM.V0037

@en

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.

@en

We study turbulent flow in open channels with a free surface and rectangular cross-section, for various Reynolds numbers and duct aspect ratios. Direct numerical simulations are used to obtain accurate characterization of the secondary motions, which are found to be more intense than in closed ducts, and to scale with the bulk, rather than with the friction velocity. A notable feature of open-duct flows is the presence of a velocity dip, namely the peak velocity is achieved at some depth underneath the free surface. We find that the depth of the velocity peak increases with the Reynolds number, and correspondingly the flow becomes more symmetric with respect to the horizontal midplane. This is also confirmed from the change of the topology of the secondary motions, which exhibit a strong corner circulation at the free-surface/wall corners at low Reynolds number, which, however, weakens at higher. The structure of the mean velocity field is such that the log law applies with good approximation in the direction normal to the nearest wall, which allows us to explain why predictive friction formulae based on the hydraulic diameter concept are successful. Additional analysis shows that the secondary motions account for a large fraction of the frictional drag (up to %).

@en

The nacelle of aircraft engines is coated with acoustic liners to reduce engine noise. An undesirable effect of these liners is that they increase aerodynamic drag. We study this drag penalty by performing Direct Numerical Simulations of a turbulent boundary layer over an acoustic liner array at friction Reynolds number, Re _τ ≈ 850–2500. We use this simulation to confirm several findings that we recently brought forward using a simpler channel flow setup [1]. We show that acoustic liners lead to high wall-normal velocity fluctuations that can be directly correlated with a modulation of the classical near-wall cycle and to an increase in drag. We also confirm that the acoustic liners act as permeable surface roughness and the non-linear Forchheimer coefficient is the relevant permeability parameter for scaling the drag increase.

@en

We study the flow above non-optimal riblets, specifically large drag-increasing and two-scale trapezoidal riblets. In order to reach large Reynolds numbers and large scale separation while retaining access to flow details, we employ a combination of boundary-layer hot-wire measurements and direct numerical simulation (DNS) in minimal-span channels. Although the outer Reynolds numbers differ, we observe fair agreement between experiments and DNS at matched viscous-friction-scaled riblet spacings in the overlapping physical and spectral regions, providing confidence that both data sets are valid. We find that hot-wire velocity spectra above very large riblets with are depleted of near-wall energy at scales that are (much) greater than. Large-scale energy likely bypasses the turbulence cascade and is transferred directly to secondary flows of size, which we observe to grow in strength with increasing riblet size. Furthermore, the present very large riblets reduce the von Kármán constant of the spanwise uniform mean velocity in a logarithmic layer and, thus, reduce the accuracy of the roughness-function concept, which we link to the near-wall damping of large flow structures. Half-height riblets in the groove, which we use as a model of imperfectly repeated (spanwise-varying) riblets, impede in-groove turbulence. We show how to scale the drag optimum of imperfectly repeated riblets based on representative measurements of the true geometry by solving inexpensive Poisson equations.

@en

In order to reduce the noise emitted by aircraft engines, the nacelle is coated with acoustic liners. An undesirable effect of these surfaces is that they increase the aerodynamic drag. In the present work, we characterize this type of surface roughness by performing Direct Numerical Simulations of fully resolved acoustic liner geometries. We find evidence of a fully rough regime, whose onset is determined by the value of the viscous-scaled Forchheimer coefficient. Moreover, the intensity of the wall-normal velocity fluctuations at the wall also scales with the viscous-scaled wall-normal permeability, leading to a relation between fluctuations and added drag.

@en

This work presents a transpacific airliner designed for minimal climate impact, incorporating several novel design features. These include open rotor engines, sustainable aviation fuels, natural laminar flow airfoils, and riblets. The design’s configuration and mission have been optimised simultaneously using a combination of standard preliminary techniques, experimental data, a multi-point mission analysis, and a model of average temperature response. It is demonstrated that, on an 8000 km mission, the design offers an 89.8% reduction in average temperature response relative to an Airbus A330-200, at the expense of a 7.3% increase in direct operating cost. The sensitivity of these results is investigated by comparing the performance over a range of operating conditions. In addition, several alternative designs incorporating only some of the above-mentioned features are analysed, allowing for an assessment of their individual contribution. Finally, a life-cycle average temperature response analysis is presented to place the climate impact of operation, manufacturing and end-of-life procedures in context@en

We investigate the Reynolds analogy over riblets, namely the analogy between the fractional increase in Stanton number Ch and the fractional increase in the skin-friction coefficient Cf, relative to a smooth surface. We investigate the direct numerical simulation data of Endrikat et al. (Flow Turbul. Combust., vol. 107, 2021, pp. 1–29). The riblet groove shapes are isosceles triangles with tip angles α=30∘,60∘,90∘, a trapezoid, a rectangle and a right triangle. The viscous-scaled riblet spacing varies between s+≈10 to 60. The global Reynolds analogy is primarily influenced by Kelvin–Helmholtz rollers and secondary flows. Kelvin–Helmholtz rollers locally break the Reynolds analogy favourably, i.e. cause a locally larger fractional increase in Ch than in Cf. These rollers induce negative wall shear stress patches which have no analogue in wall heat fluxes. Secondary flows at the riblets’ crests are associated with local unfavourable breaking of the Reynolds analogy, i.e. locally larger fractional increase in Cf than in Ch. Only the triangular riblets with α=30∘ trigger strong Kelvin–Helmholtz rollers without appreciable secondary flows. This riblet shape globally preserves the Reynolds analogy from s+=21 to 33. However, the other riblet shapes have weak or non-existent Kelvin–Helmholtz rollers, yet persistent secondary flows. These riblet shapes behave similarly to rough surfaces. They unfavourably break the global Reynolds analogy, and do so to a greater extent as s+ increases.@en

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.

@en

We carry out direct numerical simulation (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 Re_τ ≈ 2000 are considered which are much higher than in previous studies, and which yield clear scale separation between inner- and outer-layer dynamics. Close similarity between the behavior of the temperature and the streamwise velocity fields is confirmed as in previous studies related to plane channels and pipes. We find good agreement between the Nusselt number of square duct and circular pipe flow when the Reynolds number based on the hydraulic diameter is used, thus corroborating the common engineering practice. Popular engineering correlations for the heat transfer reveal deviations up to 5% with respect to DNS data, which are nicely fitted by a power law.

@en

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.

@en