Source author record

Matteo Bernardini

Matteo Bernardini appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn physics.comp-ph

Catalog footprint

What is connected

6works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

On the Reynolds analogy for high-speed rough-wall flows: implications for wall modelling

We study the validity of the generalized Reynolds analogy (GRA) in compressible turbulent boundary layers over prism-shaped roughness by mining direct numerical simulation data of Mach 2 and Mach 4 compressible turbulent boundary layers with adiabatic and cooled surfaces. Although the direct influence of roughness strongly disrupts the near-wall coupling between momentum and energy, we show that this breakdown is confined to the roughness sublayer. Above this layer, the enthalpy and velocity fields recover a smooth-wall-like similarity, and the GRA becomes asymptotically valid by naturally accounting for roughness-enhanced wall shear stress and heat flux. Building on these results, we propose a GRA-based wall model for predicting heat transfer over rough surfaces, which is coupled with a drag-predictive physics-based method developed for prism-shaped roughness by means of compressibility transformations.

preprint2020arXiv

Passive flow-separation control in a dual-bell rocket nozzle

A delayed detached eddy simulation of a sub-scale cold-gas dual-bell nozzle flow at high Reynolds number and in sea-level mode is carried out at nozzle pressure ratio NPR=45.7. In this regime the over-expanded flow exhibits a symmetric and controlled flow separation at the inflection point, that is the junction between the two bells, leading to the generation of a low content of aerodynamic side loads with respect to conventional bell nozzles. The nozzle wall-pressure signature is analysed in the frequency domain and compared with the experimental data available in the literature for the same geometry and flow conditions. The Fourier spectra in time and space (azimuthal wavenumber) show the presence of a persistent tone associated to the symmetric shock movement. Asymmetric modes are only slightly excited by the shock and the turbulent structures. The low mean value of the side-loads magnitude is in good agreement with the experiments and confirms that the inflection point dampens the aero-acoustic interaction between the separation-shock and the detached shear layer.

preprint2020arXiv

STREAmS: a high-fidelity accelerated solver for direct numerical simulation of compressible turbulent flow

We present STREAmS, an in-house high-fidelity solver for large-scale, massively parallel direct numerical simulations (DNS) of compressible turbulent flows on graphical processing units (GPUs). STREAmS is written in the Fortran 90 language and it is tailored to carry out DNS of canonical compressible wall-bounded flows, namely turbulent plane channel, zero-pressure gradient turbulent boundary layer and supersonic oblique shock-wave/boundary layer interactions. 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. The use of cuf automatic kernels allowed an easy and efficient porting on the GPU architecture minimizing the changes to the original CPU code, which is also maintained. We discuss a memory allocation strategy based on duplicated arrays for host and device which carefully minimizes the memory usage making the solver suitable for large scale computations on the latest GPU cards. Comparison between different CPUs and GPUs architectures strongly favor the latter, and executing the solver on a single NVIDIA Tesla P100 corresponds to using approximately 330 Intel Knights Landing CPU cores. STREAmS shows very good strong scalability and essentially ideal weak scalability up to 2048 GPUs, paving the way to simulations in the genuine high-Reynolds number regime, possibly at friction Reynolds number $Re_τ > 10^4$. The solver is released open source under GPLv3 license and is available at https://github.com/matteobernardini/STREAmS.

preprint2016arXiv

Detached eddy simulation of shock unsteadiness in an over-expanded planar nozzle

This work investigates the self-excited shock wave oscillations in a three-dimensional planar over-expanded nozzle turbulent flow by means of Detached Eddy Simulations. Time resolved wall pressure measurements are used as primary diagnostics. The statistical analysis reveals that the shock unsteadiness has common features in terms of the root mean square of the pressure fluctuations with other classical shock wave/boundary layer interactions, like compression ramps and incident shocks on a flat plate. The Fourier transform and the continuous wavelet transform are used to conduct the spectral analysis. The results of the former indicate that the pressure in the shock region is characterized by a broad low-frequency content, without any resonant tone. The wavelet analysis, which is well suited to study non stationary process, reveals that the pressure signal is characterized by an amplitude and a frequency modulation in time.

preprint2016arXiv

Heat transfer and wall temperature effects in shock wave turbulent boundary layer interactions

Direct numerical simulations are carried out to investigate the effect of the wall temperature on the behavior of oblique shock-wave/turbulent boundary layer interactions at freestream Mach number $2.28$ and shock angle of the wedge generator $φ= 8^{\circ}$. Five values of the wall-to-recovery-temperature ratio ($T_w/T_r$) are considered, corresponding to cold, adiabatic and hot wall thermal conditions. We show that the main effect of cooling is to decrease the characteristic scales of the interaction in terms of upstream influence and extent of the separation bubble. The opposite behavior is observed in the case of heating, that produces a marked dilatation of the interaction region. The distribution of the Stanton number shows that a strong amplification of the heat transfer occurs across the interaction, and the maximum values of thermal and dynamic loads are found in the case of cold wall. The analysis reveals that the fluctuating heat flux exhibits a strong intermittent behavior, characterized by scattered spots with extremely high values compared to the mean. Furthermore, the analogy between momentum and heat transfer, typical of compressible, wall-bounded, equilibrium turbulent flows does not apply for most part of the interaction domain. The pre-multiplied spectra of the wall heat flux do not show any evidence of the influence of the low-frequency shock motion, and the primary mechanism for the generation of peak heating is found to be linked with the turbulence amplification in the interaction region.

preprint2015arXiv

DNS of laminar-turbulent boundary layer transition induced by solid obstacles

Results of numerical simulations obtained by a staggered finite difference scheme together with an efficient immersed boundary method are presented to understand the effects of the shape of three-dimensional obstacles on the transition of a boundary layer from a laminar to a turbulent regime. Fully resolved Direct Numerical Simulations (DNS), highlight that the closer to the obstacle the symmetry is disrupted the smaller is the transitional Reynolds number. It has been also found that the transition can not be related to the critical roughness Reynolds number used in the past. The simulations highlight the differences between wake and inflectional instabilities, proving that two-dimensional tripping devices are more efficient in promoting the transition. Simulations at high Reynolds number demonstrate that the reproduction of a real experiment with a solid obstacle at the inlet is an efficient tool to generate numerical data bases for understanding the physics of boundary layers. The quality of the numerical method to fully resolve the small scales, that is one ingredient for a DNS was shown by a comparison of the exponential range of the velocity spectra, in Kolmogorov units, with those for isotropic turbulence. The good comparison reinforces the idea of local isotropy at the smallest scales.