On turbulent friction in straight ducts with complex cross section: the wall law and the hydraulic diameter

We develop predictive formulas for friction resistance in ducts with complex cross-sectional shape based on the use of the log law and neglect of wall shear stress nonuniformities. The traditional hydraulic diameter naturally emerges from the analysis as the controlling length scale for common duct shapes as triangles and regular polygons. The analysis also suggests that a new effective diameter should be used in more general cases, yielding corrections of a few percent to friction estimates based on the traditional hydraulic diameter. Fair but consistent predictive improvement is shown for duct geometries of practical relevance, including rectangular and annular ducts, and circular rod bundles

Direct numerical simulation of supersonic pipe flow at moderate Reynolds number

We study compressible turbulent flow in a circular pipe, at computationally high Reynolds number. Classical related issues are addressed and discussed in light of the DNS data, including validity of compressibility transformations, velocity/temperature relations, passive scalar statistics, and size of turbulent eddies.Regarding velocity statistics, we find that Huang's transformation yields excellent universality of the scaled Reynolds stresses distributions, whereas the transformation proposed by Trettel and Larsson (2016) yields better representation of the effects of strong variation of density and viscosity occurring in the buffer layer on the mean velocity distribution. A clear logarithmic layer is recovered in terms of transformed velocity and wall distance coordinates at the higher Reynolds number under scrutiny (\Rey_{\tau} \approx 1000), whereas the core part of the flow is found to be characterized by a universal parabolic velocity profile. Based on formal similarity between the streamwise velocity and the passive scalar transport equations, we further propose an extension of the above compressibility transformations to also achieve universality of passive scalar statistics. Analysis of the velocity/temperature relationship provides evidence for quadratic dependence which is very well approximated by the thermal analogy proposed by Zhang et Al.(2014). The azimuthal velocity and scalar spectra show an organization very similar to canonical incompressible flow, with a bump-shaped distribution across the flow scales, whose peak increases with the wall distance. We find that the size growth effect is well accounted for through an effective length scale accounting for the local friction velocity and for the local mean shear

Vortex shedding in a two-dimensional diffuser: theory and simulation of separation control by periodic mass injection

We develop a reduced-order model for large-scale unsteadiness (vortex shedding) in a two-dimensional diffuser and use the model to show how periodic mass injection near the separation point reduces stagnation pressure loss. The model estimates the characteristic frequency of vortex shedding and stagnation pressure loss by accounting for the accumulated circulation due to the vorticity flux into the separated region. The stagnation pressure loss consists of two parts: a steady part associated with the time-averaged static pressure distribution on the wall, and an unsteady part caused by vortex shedding. To validate the model, we perform numerical simulations of compressible unsteady laminar diffuser flows in two dimensions. The model and simulation show good agreement as we vary the Mach number and the area ratio of the diffuser. To investigate the effects of periodic mass injection near the separation point, we also perform simulations over a range of the injection frequencies. Periodic mass injection causes vortices to be pinched off with a smaller size as observed in experiments. Consequently, their convective velocity is increased, absorption of circulation from the wall is enhanced, and the reattached point is shifted upstream. Thus, in accordance with the model, the stagnation pressure loss, particularly the unsteady part, is substantially reduced even though the separation point is nearly unchanged. This study helps explain experimental results of separation control using unsteady mass injection in diffusers and on airfoils

Passive scalars in turbulent channel flow at high Reynolds number

We study passive scalars in turbulent plane channels at computationally high Reynolds number, thus allowing us to observe previously unnoticed effects. The mean scalar profiles are found to obey a generalized logarithmic law which includes a linear correction term in the whole lower half-channel, and they follow a universal parabolic defect profile in the core region. This is consistent with recent findings regarding the mean velocity profiles in channel flow. The scalar variances also exhibit a near universal parabolic distribution in the core flow and hints of a sizeable log layer, unlike the velocity variances. The energy spectra highlight the formation of large scalar-bearing eddies with size proportional to the channel height which are caused by a local production excess over dissipation, and which are clearly visible in the flow visualizations. Close correspondence of the momentum and scalar eddies is observed, with the main difference being that the latter tend to form sharper gradients, which translates into higher scalar dissipation. Another notable Reynolds number effect is the decreased correlation of the passive scalar field with the vertical velocity field, which is traced to the reduced effectiveness of ejection event

High-Reynolds-number effects on turbulent scalings in compressible channel flow

The effect of the Reynolds number in a supersonic isothermal channel flow is studied using a direct numerical simulation (DNS). The bulk Mach number based on the wall temperature is 1.5, and the bulk Reynolds number is increased up to Reτ ≈ 1000. The use of van Driest velocity transformation in the presence of heated walls has been questioned due to the poor accuracy at low Reynolds number. For this reason alternative transformations of the velocity profile and turbulence statistics have been proposed, as, for instance, semi-local scalings. We show that the van Driest transformation recovers its accuracy as the Reynolds number is increased. The Reynolds stresses collapse on the incompressible ones, when properly scaled with density, and very good agreement with the incompressible stresses is found in the outer layer

On algebraic TVD-VOF methods for tracking material interfaces

We revisit simple algebraic VOF methods for advection of material interfaces based of the well established TVD paradigm. We show that greatly improved representation of contact discontinuities is obtained through use of a novel CFL-dependent limiter whereby the classical TVD bounds are exceeded. Perfectly crisp numerical interfaces are obtained with very limited numerical atomization (flotsam and jetsam) as compared to previous SLIC schemes. Comparison of the algorithm with accurate geometrical VOF shows larger error at given mesh resolution, but comparable efficiency when the reduced computational cost is accounted for

On the role of secondary motions in turbulent square duct flow

We use a direct numerical simulations (DNS) database for turbulent flow in a square duct up to bulk Reynolds number \Rey_b=40000, to quantitatively analyze the role of secondary motions on the mean flow structure. For that purpose we derive a generalized form of the identity of Fukagata, Iwamoto and Kasagi (FIK), which allows to quantify the effect of cross-stream convection on the mean streamwise velocity, wall shear stress and bulk friction coefficient. Secondary motions are found to contribute for about $6\%$ of total friction, and to act as a self-regulating mechanism of turbulence whereby wall shear stress nonuniformities induced by corners are equalized, and universality of the wall-normal velocity profiles is established. We also carry out numerical experiments whereby the secondary motions are artificially suppressed, in which case their equalizing role is partially taken by the turbulent stresses

An aerothermodynamic design optimization framework for hypersonic vehicles

In the aviation field great interest is growing in passengers transportation at hypersonic speed. This requires, however, careful study of the enabling technologies necessary for the optimal design of hypersonic vehicles. In this framework, the present work reports on a highly integrated design environment that has been developed in order to provide an optimization loop for vehicle aerothermodynamic design. It includes modules for geometrical parametrization, automated data transfer between tools, automated execution of computational analysis codes, and design optimization methods. This optimization environment is exploited for the aerodynamic design of an unmanned hypersonic cruiser flying at M∞=8 and 30 km altitude. The original contribution of this work is mainly found in the capability of the developed optimization environment of working simultaneously on shape and topology of the aircraft. The results reported and discussed highlight interesting design capabilities, and promise extension to more challenging and realistic integrated aerothermodynamic design problems

Turbulence and secondary motions in square duct flow

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