Showing papers by "Sergio Pirozzoli published in 2017"

Stability and modal analysis of shock/boundary layer interactions

Abstract: The dynamics of oblique shock wave/turbulent boundary layer interactions is analyzed by mining a large-eddy simulation (LES) database for various strengths of the incoming shock. The flow dynamics is first analyzed by means of dynamic mode decomposition (DMD), which highlights the simultaneous occurrence of two types of flow modes, namely a low-frequency type associated with breathing motion of the separation bubble, accompanied by flapping motion of the reflected shock, and a high-frequency type associated with the propagation of instability waves past the interaction zone. Global linear stability analysis performed on the mean LES flow fields yields a single unstable zero-frequency mode, plus a variety of marginally stable low-frequency modes whose stability margin decreases with the strength of the interaction. The least stable linear modes are grouped into two classes, one of which bears striking resemblance to the breathing mode recovered from DMD and another class associated with revolving motion within the separation bubble. The results of the modal and linear stability analysis support the notion that low-frequency dynamics is intrinsic to the interaction zone, but some continuous forcing from the upstream boundary layer may be required to keep the system near a limit cycle. This can be modeled as a weakly damped oscillator with forcing, as in the early empirical model by Plotkin (AIAA J 13:1036–1040, 1975).

Mixed convection in turbulent channels with unstable stratification

Abstract: We study turbulent flows in pressure-driven planar channels with imposed unstable thermal stratification, using direct numerical simulations in a wide range of Reynolds and Rayleigh numbers and reaching flow conditions which are representative of fully developed turbulence. The combined effect of forced and free convection produces a peculiar pattern of quasi-streamwise rollers occupying the full channel thickness, with aspect ratio considerably higher than unity; it has been observed that they have an important redistributing effect on temperature and momentum, providing for a substantial fraction of the heat and momentum flux at bulk Richardson numbers larger than . The mean values and the variances of the flow variables do not appear to follow Prandtl’s scaling in the free-convection regime, except for the temperature and vertical velocity fluctuations, which are more directly affected by wall-attached turbulent plumes. We find that the Monin–Obukhov theory nevertheless yields a useful representation of the main flow features. In particular, the widely used Businger–Dyer flux-profile relationships are found to provide a convenient way of approximately accounting for the bulk effects of friction and buoyancy, although the individual profiles may have wide scatter from the alleged trends. Significant deviations are found in direct numerical simulations with respect to the commonly used parametrization of the momentum flux in the light-wind regime, which may have important practical impact in wall models of atmospheric dynamics. Finally, for modelling purposes, we devise a set of empirical predictive formulae for the heat flux and friction coefficients, which are within approximately standard deviation from the numerical results in a wide range of flow parameters.

Turbulence and secondary motions in square duct flow

Abstract: We study turbulent flows in pressure-driven ducts with square cross-section through direct numerical simulation in a wide enough range of Reynolds number to reach flow conditions which are representative of fully developed turbulence. Numerical simulations are carried out over extremely long integration times to get adequate convergence of the flow statistics, and specifically high-fidelity representation of the secondary motions which arise. The intensity of the latter is found to be in the order of 1-2% of the bulk velocity, and unaffected by Reynolds number variations. The smallness of the mean convection terms in the streamwise vorticity equation points to a simple characterization of the secondary flows, which in the asymptotic high-Re regime are found to be approximated with good accuracy by eigenfunctions of the Laplace operator. Despite their effect of redistributing the wall shear stress along the duct perimeter, we find that secondary motions do not have large influence on the mean velocity field, which can be characterized with good accuracy as that resulting from the concurrent effect of four independent flat walls, each controlling a quarter of the flow domain. As a consequence, we find that parametrizations based on the hydraulic diameter concept, and modifications thereof, are successful in predicting the duct friction coefficient.

Direct Numerical Simulation and Theory of a Wall-Bounded Flow with Zero Skin Friction

Abstract: We study turbulent plane Couette-Poiseuille (CP) flows in which the conditions (relative wall velocity ΔU w ≡ 2U w , pressure gradient dP/dx and viscosity ν) are adjusted to produce zero mean skin friction on one of the walls, denoted by APG for adverse pressure gradient. The other wall, FPG for favorable pressure gradient, provides the friction velocity u τ , and h is the half-height of the channel. This leads to a one-parameter family of one-dimensional flows of varying Reynolds number Re ≡ U w h/ν. We apply three codes, and cover three Reynolds numbers stepping by a factor of two each time. The agreement between codes is very good, and the Reynolds-number range is sizable. The theoretical questions revolve around Reynolds-number independence in both the core region (free of local viscous effects) and the two wall regions. The core region follows Townsend’s hypothesis of universal behavior for the velocity and shear stress, when they are normalized with u τ and h; on the other hand universality is not observed for all the Reynolds stresses, any more than it is in Poiseuille flow or boundary layers. The FPG wall region obeys the classical law of the wall, again for velocity and shear stress. For the APG wall region, Stratford conjectured universal behavior when normalized with the pressure gradient, leading to a square-root law for the velocity. The literature, also covering other flows with zero skin friction, is ambiguous. Our results are very consistent with both of Stratford’s conjectures, suggesting that at least in this idealized flow turbulence theory is successful like it was for the classical logarithmic law of the wall. We appear to know the constants of the law within a 10% bracket. On the other hand, that again does not extend to Reynolds stresses other than the shear stress, but these stresses are passive in the momentum equation.

Transitional natural convection with conjugate heat transfer over smooth and rough walls

Abstract: We study turbulent natural convection in enclosures with conjugate heat transfer. The simplest way to increase the heat transfer in this flow is through rough surfaces. In numerical simulations often the constant temperature is assigned at the walls in contact with the fluid, which is unrealistic in laboratory experiments. The DNS (Direct Numerical Simulation), to be of help to experimentalists, should consider the heat conduction in the solid walls together with the turbulent flow between the hot and the cold walls. Here the cold wall, $0.5h$ thick (where $h$ is the channel half-height) is smooth, and the hot wall has two- and three-dimensional elements of thickness $0.2h$ above a solid layer $0.3h$ thick. The independence of the results on the box size has been verified. A bi-periodic domain $4h$ wide allows to have a sufficient resolution with a limited number of grid points. It has been found that, among the different kind of surfaces at a Rayleigh number $Ra \approx 2 \cdot 10^6$, the one with staggered wedges has the highest heat transfer. A large number of simulations varying the $Ra$ from $10^3$ to $10^7$ were performed to find the different ranges of the Nusselt number ($Nu$) relationship as a function of $Ra$. Flow visualizations allow to explain the differences in the $Nu(Ra)$ relationship. Two values of the thermal conductivity were chosen, one corresponding to copper and the other ten times higher. It has been found that the Nusselt number behaves as $Nu=\alpha Ra^\gamma$, with $\alpha$ and $\gamma$ independent on the solid conductivity, and dependent on the roughness shape.

A Minimal Flow Unit for Turbulence, Combustion, and Astrophysics

Abstract: The present research is aimed at understanding the subtle (and not yet fully understood) relationship between the complex nonlinear dynamics of fluid turbulence and Kolmogorov power-law scaling in wavenumber space. Time evolution is often overlooked in DNS of turbulent flows, hence investigation of a suitably simple minimal flow unit (MFU) can help to understand the passage from a vortical-dominated stage to a turbulent stage having all the ingredients of turbulent flows. In particular, we aim at clarifying the physical phenomena associated with the formation of a finite-time singularity (FTS) in the Euler equations and of Kolmogorov’s k−5∕3 scaling in the Navier–Stokes equations. For that purpose, high-resolution simulations of the Euler and Navier–Stokes detection are carried out and analyzed by means of state-of-the-art detection techniques to isolate the contribution of tube-like and sheet-like structures. Equipping the MFU with passive scalars (relevant in turbulent combustion) further helps understanding why passive scalar spectra have a different behavior than the velocity field spectra close to the FTS, but they also attain a k−5∕3 power spectrum at subsequent times, in the presence of finite viscosity. By adding a passive vector (relevant in MHD flows), dynamical differences with respect to the vorticity field can also be established.