Numerical simulation of turbulent duct flows with constant power input

The numerical simulation of a flow through a duct requires an externally specified forcing that makes the fluid flow against viscous friction. To this aim, it is customary to enforce a constant value for either the flow rate (CFR) or the pressure gradient (CPG). When comparing a laminar duct flow before and after a geometrical modification that induces a change of the viscous drag, both approaches (CFR and CPG) lead to a change of the power input across the comparison. Similarly, when carrying out the (DNS and LES) numerical simulation of unsteady turbulent flows, the power input is not constant over time. Carrying out a simulation at constant power input (CPI) is thus a further physically sound option, that becomes particularly appealing in the context of flow control, where a comparison between control-on and control-off conditions has to be made. We describe how to carry out a CPI simulation, and start with defining a new power-related Reynolds number, whose velocity scale is the bulk flow that can be attained with a given pumping power in the laminar regime. Under the CPI condition, we derive a relation that is equivalent to the Fukagata--Iwamoto--Kasagi relation valid for CFR (and to its extension valid for CPG), that presents the additional advantage of natively including the required control power. The implementation of the CPI approach is then exemplified in the standard case of a plane turbulent channel flow, and then further applied to a flow control case, where the spanwise-oscillating wall is used for skin friction drag reduction. For this low-Reynolds number flow, using 90% of the available power for the pumping system and the remaining 10% for the control system is found to be the optimum share that yields the largest increase of the flow rate above the reference case, where 100% of the power goes to the pump.Comment: Accepted for publication in J. Fluid Mec

Does the choice of the forcing term affect flow statistics in DNS of turbulent channel flow?

We seek possible statistical consequences of the way a forcing term is added to the Navier--Stokes equations in the Direct Numerical Simulation (DNS) of incompressible channel flow. Simulations driven by constant flow rate, constant pressure gradient and constant power input are used to build large databases, and in particular to store the complete temporal trace of the wall-shear stress for later analysis. As these approaches correspond to different dynamical systems, it can in principle be envisaged that these differences are reflect by certain statistics of the turbulent flow field. The instantaneous realizations of the flow in the various simulations are obviously different, but, as expected, the usual one-point, one-time statistics do not show any appreciable difference. However, the PDF for the fluctuations of the streamwise component of wall friction reveals that the simulation with constant flow rate presents lower probabilities for extreme events of large positive friction. The low probability value of such events explains their negligible contribution to the commonly computed statistics; however, the very existence of a difference in the PDF demonstrates that the forcing term is not entirely uninfluential. Other statistics for wall-based quantities (the two components of friction and pressure) are examined; in particular spatio-temporal autocorrelations show small differences at large temporal separations, where unfortunately the residual statistical uncertainty is still of the same order of the observed difference. Hence we suggest that the specific choice of the forcing term does not produce important statistical consequences, unless one is interested in the strongest events of high wall friction, that are underestimated by a simulation run at constant flow rate

A modified Parametric Forcing Approach for modelling of roughness

Surface roughness in turbulent channel flow is effectively modelled using a modified version of the Parametric Forcing Approach introduced by Busse and Sandham (2012). In this modified approach, the model functions are determined based on the surface geometry and two model constants, whose value can be fine tuned. In addition to a quadratic forcing term, accounting for the effect of form drag due to roughness, a linear forcing term, analogous to the Darcy term in the context of porous media, is employed in order to represent the viscous drag. Comparison of the results with full-geometry resolved Direct Numerical Simulation (DNS) data for the case of dense roughness (frontal solidity ≅0.4) shows a satisfactory prediction of mean velocity profile, and hence the friction factor, by the model. The model is found to be able to reproduce the trends of friction factor with morphological properties of surface such as skewness of the surface height probability density function and coefficient of variation of the peak heights

Money versus Time: Evaluation of Flow Control in Terms of Energy Consumption and Convenience

Flow control with the goal of reducing the skin friction drag on the fluid-solid interface is an active fundamental research area, motivated by its potential for significant energy savings and reduced emissions in the transport sector. Customarily, the performance of drag reduction techniques in internal flows is evaluated under two alternative flow conditions, i.e. at constant mass flow rate or constant pressure gradient. Successful control leads to reduction of drag and pumping power within the former approach, whereas the latter leads to an increase of the mass flow rate and pumping power. In practical applications, however, money and time define the flow control challenge: a compromise between the energy expenditure (money) and the corresponding convenience (flow rate) achieved with that amount of energy has to be reached so as to accomplish a goal which in general depends on the specific application. Based on this idea, we derive two dimensionless parameters which quantify the total energy consumption and the required time (convenience) for transporting a given volume of fluid through a given duct. Performances of existing drag reduction strategies as well as the influence of wall roughness are re-evaluated within the present framework; how to achieve the (application-dependent) optimum balance between energy consumption and convenience is addressed. It is also shown that these considerations can be extended to external flows

DNS of Turbulent Impinging JETS on rough surfaces using a parametric forcing Approach

This work presents direct numerical simulations (DNS) of a circular turbulent jet impinging on rough plates. A parametric forcing approach (PFA) accounts for surface roughness effects by applying a forcing term into the Navier-Stokes equations within a thin layer in the near-wall region. The application of the PFA in the context of spatially developing flows is the essential aspect of the investigation. The method is known to produce accurate predictions of the velocity field in fully developed turbulent flows while avoiding the demanding grid resolution required by an immersed boundary method (IBM) approach. The comparison between PFA results and IBM-resolved roughness DNS allows addressing the advantages and limitations of the PFA in the context of spatially developing flows

Heat Transfer of an Impinging JET : Sensitivity Towards Inflow Conditions

The wall-heat transfer in the impingement region of a turbulent incompressible jet impinging onto a flat plate is assessed using direct numerical simulations. For impinging jets, the mean wall-heat transfer distribution along the plate features a global maximum close to the jet axis and, for specific configurations, a secondary peak further downstream of the impingement region. While the appearance of the latter maximum has been the focus of many investigations, the occurrence of the former and its links to the inflow conditions of the jet have only been alluded to by existing studies. The present research considers two different inflow conditions: a fully developed turbulent pipe flow and a slightly convergent nozzle. The study shows that turbulent fluctuations, characterising the impingement region of the fully-developed turbulent inflow case, are significant in determining events of positive wall-heat transfer rates very close to the jet axis. On the contrary, similar events are not observed in the stagnation region of the convergent nozzle case, and the associated mean wall-heat transfer displays a markedly different distribution compared to the fully- developed turbulent inflow case