APS/123-QED On the missing link between pressure drop, viscous dissipation and the turbulent energy spectrum A. Badillo∗ and O. Matary Department of Chemical Engineering, Imperial College London, UK. (Dated: April 27, 2021) We present convincing evidence of a direct connection between the pressure drop, viscous dissipa- tion, and the turbulent energy spectrum. We use this finding to explain Nikuradse's experimental results of pressure drop for turbulent flows in rough pipes, in terms of a modified Kolmogorov length scale that varies with the surface roughness. Furthermore, we use Laufer's measurements of turbu- lent energy spectra in pipe flow to calculate|to a good approximation|the turbulent component of the pressure drop directly from an averaged turbulent energy spectrum. We also show that the incompressibility assumption, leads to the conclusion that viscous dissipation in fully-developed (laminar and turbulent) pipe flow, cannot increase the temperature of the fluid through viscous heating. PACS numbers: $ For many decades, turbulence has remained stub- τ = µ ru + ruT the viscous stress tensor expressed bornly one of the most challenging open problems in in terms of the dynamic viscosity µ and the strain rate. classical physics. Achieving complete understanding of In our analysis, we consider a portion of a pipe, where turbulent flows is hindered by the their multi-scale na- the flow is adiabatic, incompressible and fully developed ture, and the complex mechanisms underlying the dis- (Fig. 1). Hence, all the quantities are statistically invari- tribution of energy to the large number of harmonics, ant along the stream-wise direction, except for the fluid which constitute the fluctuating velocity field. It is well pressure. accepted that vortex stretching is one of the main mecha- nisms responsible for transferring energy from low to high wavenumbers. This mechanism, however, only exists in viscous, three-dimensional turbulence, with no analogue in quantum [1], and two-dimensional turbulence [2]. The pictorial energy cascade proposed by Richardson [3], and then recovered theoretically first by Kolmogorov [4, 5], and then, independently, by Onsager [6], does not reveal much about the energy transfer between consecu- tive harmonics; it only represents an `equilibrium' state of energy distribution arising from a balance between the fluxes of energy going up and down the cascade. Despite the fact that viscous dissipation is irrelevant for quantum turbulence, the scaling laws dominating its energy cas- FIG. 1: Schematic of fully developed pipe flow. The pipe is cade are the same as those in viscous, three-dimensional thermally insulated (adiabatic) and the fluctuations of any turbulence [7{9]. What, then, is the true role of viscous flow quantity are statistically identical at the inlet and outlet dissipation in the formation of the energy cascade, and of the portion of the analyzed pipe. We also present how, what is its connection to macroscopic quantities such as hypothetically, the molecular radial distribution varies along the pipe (wavy lines in the upper part of the figure). The the pressure drop? We address the second question in red circle surrounded by black ones represents, schematically, the present letter. the local arrangement of atoms/molecules in the fluid. At We start our discussion from the conservation of linear $ inlet, the pressure is higher than that at the outlet, leading momentum @tp + r · (up) = −∇· σ +Σfext, with p = to a slightly higher local degree of order. This higher degree arXiv:2104.12110v1 [physics.flu-dyn] 25 Apr 2021 ρu the linear momentum, ρ the density, u the instanta- of short range order is observed as a larger amplitude of the $ second and subsequent peaks in the radial distribution. neous velocity field, σ the stress tensor, and Σfext the summation of all external forces. We further split the $ $ $ By applying a time average to the conservation of mo- stress tensor as σ= P I − τ , with P the pressure and mentum, we arrive to the following expression for the pressure drop 2L ∆P = − ρu2 (1) ∗Electronic address: [email protected] R τ yElectronic address: [email protected] p with R the pipe radius, L the pipe length, uτ = τw/ρ 2 the friction velocity, τw = (µ∂U=@r)w the wall shear indicate that the local range order (a measure of how stress, and U the stream-wise component of the time- molecules are arranged in the liquid) undergoes impor- averaged velocity field. tant changes as pressure increases. Using the MD data An equivalent expression can be obtained from the from Mahoney and Jorgensen [11], we can estimate the 2 conservation of kinetic energy ek = u =2 (per unit derivative deu=dP . Applying a linear fit to their MD $ $ data yields a value of de =dP ∼ −0:826 (see supplemen- mass) @ (ρe ) + r · uρe + σ ·u =σ : ru + Σf · u. u t k k ext tal material). Since the fluid considered in the MD study By applying a time average to the conservation law for is slightly compressible [11], we would expect a value kinetic energy and invoking incompressible and fully de- slightly different from −1. These molecular studies, show veloped pipe flow conditions, the pressure drop reads the strong degree of correlation between internal energy and the molecular arrangement of the liquid. This would 1 Z D$ E ∆P = − τ : ru dΩ; (2) imply that for strictly incompressible fluids, the energy VA Ω of the flow (pressure + kinetic) is dissipated only into the potential component of the internal energy follow- with h·i a time average, VA = R u · dA the volumetric A ing a reversible path, and not, as commonly believed, flow rate, A the pipe cross section, and V the average into the micro-kinetic component of the internal energy cross section velocity (see supplemental material for a (temperature). For a liquid flowing in a pipe, the higher detailed derivation of conservation laws and the two ex- pressure at the inlet will cause a slightly higher degree of pressions for the pressure drop). Equation (2) was also molecular order compared to that at the outlet (see Fig. presented by Pope [10]. Since the conservation law for 1). kinetic energy is not independent from that of momen- tum, Eqs. (1) and (2) are completely equivalent. Thus, To reveal the connection between Eq. (2) and the tur- we can express the wall shear stress in terms of the total bulent the energy spectrum, we follow the work of Gioia and Chakraborty[12]. They investigated the relationship viscous dissipation. 2 between the friction factor, f = 8τw/ρV and the phe- 1 Z D$ E nomenological turbulent energy spectrum. Considering τ = − τ : ru dΩ; (3) w 2πRLV the Darcy-Weisbach equation, we can express the pres- Ω sure drop as We need to emphasize that Eq. (3) is only valid for fully L ρV 2 developed (laminar or turbulent) and incompressible pipe ∆P = −f ; (4) flow. For laminar flow, it is possible to confirm analyti- DH 2 cally the validity of Eqs. (2) and (3). The first interest- ing conclusion from Eq. (2), is that viscous dissipation where DH the hydraulic pipe diameter. The negative cannot increase the temperature of a fluid at the incom- sign is to enforce a decrease in the pressure along the pressibility limit. The conservation of enthalpy per unit pipe. By replacing Eq. (2) into (4) and introducing the $ Reynolds decomposition u = U + u0, with U = hui, we mass, h, implies that ρD h + r · q = D P + τ : ru (see t t can express the friction factor in terms of two components supplemental material), with D = @ +u·r the material t t of the viscous dissipation, derivative and q the heat flux. If the pipe is adiabatic and the flow incompressible, the only possibility for the en- thalpy to change, is through viscous heating (dissipation Z Z 0 64 1 0 term). After applying a time average to the conserva- f = dΩ + 0 dΩ (5) Re 8πLV 2 ν ν tion law of enthalpy, the integration of the R.H.S. over Ω Ω the pipe volume, vanishes by virtue of Eq. (2) and, thus, R where ν is the kinematic viscosity, and hρD hi dΩ = 0. Since there is no energy exchange with D$ E Ω t ρ−1 τ : ru = + 0 is defined as the total time- the environment (adiabatic flow), this requires either a 0 0 T constant fluid temperature or that the local viscous dissi- averaged dissipation, with 0/ν = rU: rU+rU : rU 0 0 0 0T 0 pation increase and decrease the local fluid temperature and 0/ν = ru : ru + ru : ru the mean and in a way that the net change of total enthalpy is null. Be- turbulent components of the viscous dissipation respec- cause the second option would be in violation with the tively. The Reynolds number, Re = VDH /ν, is defined first law of thermodynamics, the fluid temperature must in terms of the hydraulic pipe diameter, DH . be constant at the incompressibility limit. The main difficulty faced by Gioia and Chakraborty For molecular fluids like water, the internal energy can [12], was the determination of the wall shear stress, and suffer important variations by changing the molecular ar- how to relate it to the turbulent energy spectrum and the rangement, keeping an almost constant temperature and surface roughness. They postulated that the wall shear density. A small change in the enthalpy per unit volume, stress can be expressed in terms of a momentum transfer ~ ~ h, is given by dh = deu + dP , with eu being the inter- between the bulk flow and the flow near the wall, charac- nal energy per unit volume. If the enthalpy of the fluid terized by a velocity scale us.