Production Of Turbulent Kinetic Energy Research Articles

The Vertical Structure of Turbulence Kinetic Energy Near the Arctic Sea‐Ice Surface

AbstractAtmospheric turbulence over the Arctic sea‐ice surface has been understudied due to the lack of observational data. In this study, we focus on the turbulence kinetic energy (TKE) over sea ice and distinguish its two different vertical structures, the “Surface” type and the “Elevated” type, using observations during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition (MOSAiC). The “Surface” type has the maximum TKE near the surface (at 2 m), while the “Elevated” type has the maximum TKE at a higher level (6 m). The TKE budget analysis indicates that the “Elevated” type is caused by the increased shear production of TKE at 6 m. In addition, spectral analysis reveals that the contribution to TKE by horizontal large eddies is enhanced in the “Elevated” type. Finally, how the vertical structure of TKE affects the parameterization of turbulent momentum flux is discussed.

Study on power losses and pressure fluctuations of diffuser mixed flow pump as turbine on different power generation speeds based on energy power models

Variable speed regulation represents an effective means to ensure the efficient operation of the diffuser mixed flow pump as turbine (DMFPAT). This paper employed computational fluid dynamics-based energy power models and high-frequency water pressure fluctuation tests to investigate the power losses and pressure fluctuations of DMFPAT under different power generation speeds. Analysis of power losses reveals that internal power losses in DMFPAT predominantly originate from the turbulent dissipation power (PEPD′) in the entropy production and the turbulent kinetic energy production power (PEBL3) in the energy balance equation, indicating that fully developed turbulence significantly contributes to substantial power losses in DMFPAT. A strong correlation exists between power generation speed and power losses, and decreasing rotational speed results in gradual reduction of power losses in DMFPAT. Additionally, a comprehensive analysis of the pressure fluctuation signals of DMFPAT under variable speed mode was conducted, revealing that the main frequency at each monitoring point is the rotational frequency and its harmonics. As speed decreases, the pressure fluctuation amplitude diminishes. The research presented demonstrates that adjusting the speed not only markedly reduces power losses but also eliminates adverse pressure fluctuations, thereby providing a theoretical basis for the efficient, stable and safe operation of DMFPAT.

Optimization of Low-Pressure Turbine Blade by Means of Fine Inspection of Loss Production Mechanisms

Abstract In this study, Proper Orthogonal Decomposition (POD) was applied to Large Eddy Simulations (LES) of two high-loaded low-pressure turbine cascades under unsteady inflow conditions to investigate entropy production within different parts of the blade passage. The terms for Turbulent Kinetic Energy (TKE) production, diffusion, and dissipation from the stagnation pressure transport equation were integrated across the computational domain. The POD-based method enables the decomposition of contributions from various coherent flow dynamics to TKE production and dissipation, such as the migration, bowing, tilting, and reorientation of incoming wake filaments, as well as the breakdown of streaky structures in the blade boundary layer and the formation of Von Karman vortices in the blade wake. This approach helps designers identify the dominant POD modes (i.e., turbulent flow structures) responsible for loss generation, understand their dynamics and locate where they primarily act. Using this optimization strategy, a new low-pressure turbine profile was designed. LES calculations on the optimized geometry showed that it is possible to design a higher-loaded profile with a lower global loss coefficient. The new loading distribution, characterized by stronger acceleration in the early part of the blade passage, results in reduced upstream wake migration losses and lower TKE production in the trailing edge wake zone due to early suction side boundary layer transition.

Parameterization of Vertical Turbulent Transport in the Inner Core of Tropical Cyclones and Its Impact on Storm Intensification. Part I: Sensitivity to Turbulent Mixing Length

Abstract In the inner core of a tropical cyclone, turbulence not only exists in the boundary layer (BL) but can also be generated above the BL by eyewall and rainband clouds. Thus, the treatment of vertical turbulent mixing must go beyond the conventional scope of the BL. The turbulence schemes formulated based on the turbulent kinetic energy (TKE) are attractive as they are applicable to both deep and shallow convection regimes in the tropical cyclone (TC) inner core provided that the TKE production and dissipation can be appropriately determined. However, TKE schemes are not self-closed. They must be closed by an empirically prescribed vertical profile of mixing length. This motivates this study to investigate the sensitivity of the simulated TC intensification to the sloping curvature and asymptotic length scale of mixing length, the two parameters that determine the vertical distribution of a prescribed mixing length. To tackle the problem, both idealized and real-case TC simulations are performed. The results show that the simulated TC intensification is sensitive to the sloping curvature of mixing length but only exhibits marginal sensitivity to the asymptotic length scale. The underlying reasons for such sensitivities are explored analytically based on the Mellor and Yamada level-2 turbulence model and the analyses of azimuthal-mean tangential wind budget. The results highlight the uncertainty and importance of mixing length in the numerical prediction of TCs and suggest that future research should focus on searching for physical constraints on mixing length, particularly in the low- to midtroposphere, using observations and large-eddy simulations. Significance Statement The parametric representation of subgrid-scale turbulent mixing is one of the major sources of uncertainty in numerical predictions of tropical cyclones (TCs). This study investigates how the numerical prediction of TC intensification is affected by the turbulent mixing length, a length scale that is required to close a turbulence scheme formulated based on the turbulent kinetic energy (TKE). The research highlights the uncertainty and importance of mixing length in numerical prediction of TCs and suggests that future research should focus on searching for physical constraints on the mixing length, particularly in the low- to midtroposphere, using observations and large-eddy simulations.

Flow Field Characteristics at the Spanwise Edge of a Vegetative Canopy Model

This study investigates the complex flow field across a spanwise vegetative model canopy edge focusing on turbulent transport processes. Utilizing stereoscopic particle image velocimetry, the three velocity components were measured in wall-parallel planes at various elevations within canopy across the spanwise canopy edge. Conventional ensemble averaged results were contrasted with those obtained by conditionally averaged flow properties across instantaneous internal interfaces in the flow to understand their contribution to the ensemble average. The conditional average captured the strong gradients in mean velocities, Reynolds stresses, vorticity, swirling strength, and turbulent kinetic energy production across the dynamically changing instantaneous interface. In contrast, the conventional ensemble average smeared out the strong gradients. Small magnitudes of advective terms in the turbulent kinetic energy transport equation suggested weak secondary transverse flows in the present model canopy. The turbulent flow structure across the spanwise canopy edge was further investigated using Quadrant-Hole analysis for both averaging approaches. Conventional ensemble averaged results indicated a shift from sweep to ejection dominance when moving from canopy into the open patch, while the conditional average showed only sweep dominated transport. In contrast to a homogeneous canopy layout, below canopy height at the canopy edge, sweeps and ejections lose their dominance in vertical turbulent transport. The present results show that the dynamics of internal interfaces govern the ensemble averaged results and a possible implementation into existing models is proposed. The present results are expected to increase understanding of spanwise turbulent transport and aid in developing strategies to mitigate desertification.

Particle image velocimetry in the impeller of a centrifugal pump: Relationship between turbulent flow and energy loss

Turbulent flows play a dominant role in the operation of centrifugal pumps, which find widespread use in industrial settings and various aspects of human life. The dissipation rate of turbulent kinetic energy emerges as a key parameter within these devices, with its local values exerting a significant influence on centrifugal pump performance. Recent advances in particle image velocimetry (PIV) techniques have expanded the ability to analyze complex turbulent flows across a broad spectrum of scales. In this context, this paper aims to deepen our understanding of the turbulent flow field and its correlation with energy loss in centrifugal pump impellers. To achieve this, experiments were conducted using PIV on a transparent pump operating under different conditions. Statistics of the turbulent flow were then obtained from phase-ensemble averages of velocities, vorticity, turbulence production, and local dissipation of turbulent kinetic energy. To overcome the limited spatial resolution constraint of PIV, the large-eddy PIV (LES-PIV) method was employed to estimate the local dissipation rate. In this method, it is assumed that the motion of larger scales is measured by the PIV technique, while the smaller scales (unresolved scales) are modeled by a sub-grid scale model, calculated from the strain rate tensors obtained from the measured fields. Energy losses in the impeller were studied using two methodologies: (i) a conventional method based on power measurements, and (ii) an alternative approach based on the budget of turbulent kinetic energy. Our results reveal that turbulent loss caused by turbulence production is the main source of energy loss in the pump impeller, and it is particularly pronounced in low-flow operating conditions characterized by large-scale structures. On the other hand, in situations where flow rates exceed the best efficiency point (BEP) condition, the predominant flow structures are marked by small-scale features, mainly attributed to local dissipation of turbulence. Our findings clarify the characteristics of energy losses in centrifugal pump impellers and their relationship with the turbulent flow field, and, in addition, providing a methodology for calculating the local turbulent dissipation rate and its limitations when derived from PIV measurements.

Numerical Simulations of Scalar Transport on Rough Surfaces

Numerical simulations provide unfettered access to details of the flow where experimental measurements are difficult to obtain. This paper summarises the progress achieved in the study of passive scalars in flows over rough surfaces thanks to recent numerical simulations. Townsend’s similarity applies to various scalar statistics, implying the differences due to roughness are limited to the roughness sublayer (RSL). The scalar field exhibits a diffusive sublayer that increasingly conforms to the roughness surface as ks+ or Pr increase. The scalar wall flux is enhanced on the windward slopes of the roughness, where the analogy between momentum and scalar holds well; the momentum and scalar fields, however, have very different behaviours downwind of the roughness elements, due to recirculation, which reduces the scalar wall flux. Roughness causes breakdown of the Reynolds analogy: any increase in St is accompanied by a larger increase in cf. A flattening trend for the scalar roughness function, ΔΘ+, is observed as ks+ increases, suggesting the possibility of a scalar fully rough regime, different from the velocity one. The form-induced (FI) production of scalar fluctuations becomes dominant inside the RSL and is significantly different from the FI production of turbulent kinetic energy, resulting in notable differences between the scalar and velocity fluctuations. Several key questions remain open, in particular regarding the existence of a fully rough scalar regime and its characteristics. With the increase in Re and Pr, various quantities such as scalar roughness function, the dispersive fluxes, FI wall flux, etc., appear to trend towards saturation. However, the limited range of Re and Pr achieved by numerical simulations only allows us to speculate regarding such asymptotic behaviour. Beyond extending the range of Re and Pr, systematic coverage of different roughness types and topologies is needed, as the scalar appears to remain sensitive to the geometrical details.

A numerical study of rainfall effects on wind turbine wakes

Rainfall has significant impacts on the operations of wind turbines due to erosion of turbine blades and alterations in aerodynamic performance of wind turbines. However, little is understood in the influences of rainfall on wind turbine wakes. In this study, large-eddy simulation (LES) coupled with actuator-disk model with rotation (ADM-R) is used to investigate impacts of rainfall on wind turbine wakes. A double Euler method is employed in ADM-R to simulate the rainfall injection. The results of the model indicate that rainfall reduces wind speed in the sweep area and increases wind speed in the outer region of the top tip level by 2.1 %. Furthermore, rainfall reduces turbulence intensity in the near wake and outer wake region of the top tip (up to 2.0 %), with the influence extending up to 10 d (diameters of the wind turbine) in the downstream. These changes have a positive correlation with the rainfall intensity and an inverse correlation with the wind speed. The mean kinetic energy (MKE) and turbulent kinetic energy (TKE) budget analysis reveals that (1) the turbulent radial transport variation of MKE is the primary reason for the rainfall-induced wind speed change; (2) the change in shear production of TKE is responsible for the rainfall-induced turbulent intensity change; (3) the rainfall-induced Reynolds stress −u′w′‾ is the dominant factor of this dynamics.

An innovative modification to the Menter shear-stress transport turbulence model employing the symbolic regression approach

Drawing from the non-equilibrium link between the production Pk and dissipation ε of turbulent kinetic energy (TKE), we advocate for the introduction of a limiter to modulate the TKE production term within the Menter shear-stress transport (SST) model. The original SST model is made more sensitive to the adverse pressure gradient (APG) by Bradshaw's assumption. Bradshaw's assumption introduces the equilibrium condition Pk/ε = 1 in most regions of the turbulent boundary layer. In the APG flows with Pk≫ε, the equilibrium condition suppresses the magnitude of TKE (k) within the boundary layer, resulting in an early separation problem. To address this issue, we employ the symbolic regression (SR) to scrutinize the physical correlation between Pk/ε and local turbulence parameters, obtaining an approximate function FSR that encapsulates the relationship between Pk/ε, Sk/ε, and y+ in the APG flow. Following its incorporation into the original SST model in the form of a limiter, the FSR evolves into the SST-Symbolic Regression Evolution model. The SST-SRE is then applied to four cases with APGs. The modification leads to an increase in the skin-friction coefficient Cf in the APGs region and causes a downstream shift in the separation location, improving the consistency with high-accuracy data and experimental results. It is demonstrated that this correction can improve the early separation problem in the Menter SST turbulence model.

Turbulent kinetic energy budget of sediment-laden open-channel flows: bedload-induced wall-roughness similarity

New experiments in highly turbulent, steady, subcritical and uniform water open-channel flows have been carried out to measure the mean turbulent kinetic energy (TKE) budget of sediment-laden boundary layer flows with two sizes (dp = 3 mm and 1 mm) of Plexiglas particles $({\rm relative\ density}\ = 1.192)$ . The experiments covered energetic sediment transport conditions (Shields number of $0.35 < \theta < 1.2$ ) ranging from non-capacity to full-capacity flows in bedload-to-suspension-dominated transport modes (suspension number of $0.5 < {w_s}/{u_\ast } < 1.3$ where ${w_s}$ is the settling velocity and $u_*$ is the friction velocity) and for weakly to highly inertial, finite size turbulence-particle conditions (Stokes number of 0.1 < St < 3.5 and dp/η > 10 where $\eta$ is the Kolmogorov length scale). It was shown that the effects of sediments on the TKE budget are very pronounced in all large particle experiments for which a bedload layer of several grain diameter thickness is developed above the channel bed. When compared with the corresponding reference clear-water flows, the TKE shear-production rate for the 3 mm particle flows is strongly reduced in the wall region corresponding to the bedload layer. This turbulence damping is seen to increase with sediment load until full capacity for flows with constant Shields value, as well as with Shields number value. Inside this damped TKE shear-production zone, a distinct peak of maximal turbulence production appears to coincide with the upper edge of the bedload layer delimited by a sharp gradient in mean sediment concentration. This vertically upshifted peak of TKE production is accompanied by an enhanced net downward oriented TKE flux when compared with the reference clear-water flows. The downward diffused TKE is found to act in the bedload layer as a local energy source in reasonable balance with the sediment transport term. The mechanism behind this downward TKE transport was further analysed on the basis of coherent flow structure dynamics controlled by ejection- and sweep-type events. The agreement between the height of downward directed mean TKE flux and the height below which sweep-type events dominate the Reynolds shear-stress contribution over ejections, revealed the leading role played by sweeps in mean TKE transport. This agreement holds for all reference clear-water flows supporting the well-known wall-roughness-induced dominance of the sweep contribution in turbulent, rough clear-water boundary layer flows. Furthermore, for all 3 mm particle flows, the two referred to transition levels were significantly and similarly upshifted to the upper edge of the bedload layer. Only for these sediment-laden flows, the bedload layer thickness is seen to exceed the wall-roughness sublayer of the reference clear-water flows. This supports a strong analogy between wall-roughness effects in clear-water flows and bedload layer effects in sediment-laden flows, on the mean TKE budget induced by a similarly modified coherent flow structure dynamics. The bedload layer-controlled wall roughness is finally confirmed by the good prediction of the wall-roughness parameter ks of the logarithmic velocity distribution. An empirical formulation fitting the presented measurements is presented, valid over the range of Shields number values covered herein.

Investigation on the turbulent structures in combined wave-current boundary layers

The problem of turbulent vortices and vortex dynamics under the influence of ocean waves is a subject of continuing research in the field of wave-current interaction. Thus, this study investigates the evolution of turbulence structures and vortex interactions in combined wave-current boundary layers, where water waves propagate with currents. Experiments were performed in two glass-walled laboratory channels with water depths of 160 and 400 mm. A two-dimensional two-component (2D2C) Particle Image Velocimetry (PIV) system was used to measure the velocities. Periodic variations in turbulence are observed for combined wave-current flows. Pairs of counter-rotating vortices are identified in the outer region when the combined flow is reversed, leading to destructive vortex interactions. This process leads to a reduction in the turbulent kinetic energy (TKE) and TKE production and dissipation during the passage of the wave trough. A conceptual model of hairpin vortices in combined wave-current flows is proposed in this study. The results of this study provide novel insights into the hydrodynamics of combined wave-current flows.

A modification of production term for the shear stress transport turbulence model based on a new wall-law for adverse pressure gradient

This paper presents an attempt to improve the predictive accuracy of the Menter shear stress transport (SST) turbulence model for adverse pressure gradient (APG) flows. The foundation of this work lies in the newly modified law of the wall and turbulent kinetic energy (TKE) transport equation within the APG turbulent boundary layer. The modification of the law of the wall is achieved through a newly proposed variable λ, which is derived theoretically and calibrated using direct numerical simulation data. Based on the new law of the wall, we analyze the TKE transport equation of the SST model at APG and demonstrate that the analytical solution satisfying the law of the wall cannot achieve equilibrium in the current SST model. Further theoretical analysis reveals that this imbalance is caused by the TKE production term Pk, and the decomposition of the wall friction coefficient Cf indicates that Pk directly affects the prediction of flow separation. Based on the analysis and the new law of the wall, an additional term for Pk is suggested to correct Pk and improve the model's predictive accuracy for APG flows. Numerical validation results show that the correction leads to an increase in Pk within the APG turbulent boundary layer, resulting in more accurate predictions of Cf and the mean velocity profile in the inner layer of the turbulent boundary layer.

Spacing Ratio Effects on the Evolution of the Flow Structure of Two Tandem Circular Cylinders in Proximity to a Wall

The flow around two tandem circular cylinders in proximity to a wall is investigated using particle image velocimetry (PIV) for Re = 2 × 103. The spacing ratios L/D are 1, 2, and 5, and the gap ratios G/D are 0.3, 0.6, and 1. The proper orthogonal decomposition (POD) method and λci vortex identification method are used to investigate the evolution of flow structure, and the influences of L/D and G/D on flow physics are shown. At L/D = 2 and G/D = 0.3, a “pairing” process occurs between the wall shear layer and the upstream cylinder’s lower shear layer, resulting in a small separation bubble behind the upstream cylinder. At L/D = 1, the Strouhal number (St) increases with decreasing G/D. At three gap ratios, the St gradually decreases as L/D increases. At G/D = 0.3, there is nearly a 49.98% decrease from St = 0.3295 at L/D = 1 to St = 0.1648 at L/D = 5, which is larger than the reductions in cases of G/D = 0.6 and G/D = 1. The effects of L/D on the evolution of flow structure at G/D = 0.6 are revealed in detail. At L/D = 1, the vortex shedding resembles that of the single cylinder. As L/D increases to 2, a squarish flow structure is formed between two cylinders, and a small secondary vortex is formed due to induction of the lower shear layer of the upstream cylinder. At L/D = 5, there is a vortex merging process between the upper wake vortices of the upstream and downstream cylinders, and the lower wake vortex of the upstream cylinder directly impinges the downstream cylinder. In addition, the shear layers and wake vortices of the upstream cylinder interact with the wake of the downstream cylinder as L/D increases, resulting in reductions in velocity fluctuations, and the production and turbulent diffusion of turbulent kinetic energy are decreased behind the downstream cylinder.

Reynolds number effect on the parameters of turbulent flows over open channels

Abstract Experimental investigations were conducted to analyze the effect of Reynolds numbers on turbulent flow properties in a nonuniform sand bed channel. Steady flow simulations were performed over the nonuniform sand bed channel, considering five Reynolds numbers within the range of 36500–53886. This article endeavors to delineate the influence of Reynolds number on turbulent flow properties through meticulous laboratory studies. Observations revealed that higher Reynolds numbers corresponded to increased longitudinal velocity. As the Reynolds number increases by 10 to 47%, various turbulent flow properties exhibit distinct trends. Specifically, the longitudinal velocity, longitudinal turbulent intensity, vertical turbulent intensity, turbulent kinetic energy, Reynolds shear stress, and Taylor scale show increases ranging from 5 to 30%, 15 to 25%, 15 to 20%, 25 to 60%, 20 to 40%, and 35 to 45%, respectively. Taylor scale analysis indicated higher magnitudes associated with higher Reynolds numbers. In-depth examinations of turbulent anisotropy, third-order moments of velocity fluctuations, kurtosis, turbulent kinetic energy production, and dissipation provided additional insights into flow behavior across different Reynolds numbers. This study contributes to a more comprehensive understanding of flow dynamics in nonuniform sand bed channels under varying Reynolds number conditions, bridging the gap between laboratory studies and real-world scenarios.

Experimental assessment of square-wave spatial spanwise forcing of a turbulent boundary layer

We present an experimental realisation of spatial spanwise forcing in a turbulent boundary layer flow, aimed at reducing the frictional drag. The forcing is achieved by a series of spanwise running belts, running in alternating spanwise direction, thereby generating a steady spatial square-wave forcing. Stereoscopic particle image velocimetry in the streamwise–wall-normal plane is used to investigate the impact of actuation on the flow in terms of turbulence statistics, drag performance characteristics, and spanwise velocity profiles, for a non-dimensional wavelength of λx+=397\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\lambda _x^+ = 397$$\\end{document}. In line with reported numerical studies, we confirm that a significant flow control effect can be realised with this type of forcing. The scalar fields of the higher-order turbulence statistics show a strong attenuation of stresses and production of turbulence kinetic energy over the first belt already, followed by a more gradual decrease to a steady-state energy response over the second belt. The streamwise velocity in the near-wall region is reduced, indicative of a drag-reduced flow state. The profiles of the higher-order turbulence statistics are attenuated up to a wall-normal height of y+≈100\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$y^+ \\approx 100$$\\end{document}, with a maximum streamwise stress reduction of 45% and a reduction of integral turbulence kinetic energy production of 39%, for a non-dimensional actuation amplitude of A+=12.7\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$A^+ = 12.7$$\\end{document}. An extension of the classical laminar Stokes layer theory is introduced, based on the linear superposition of Fourier modes, to describe the non-sinusoidal boundary condition that corresponds to the current case. The experimentally obtained spanwise velocity profiles show good agreement with this extended theoretical model. The drag reduction was estimated from a linear fit in the viscous sublayer in the range 2≤y+≤5\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2 \\le y^+\\le 5$$\\end{document}. The results are found to be in good qualitative agreement with the numerical implementations of Viotti et al. (Phys Fluids 21, 2009), matching the drag reduction trend with A+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$A^+$$\\end{document}, and reaching a maximum of 20%.Graphical abstract

Laminarization of Low Reynolds Number Turbulent Flow in Heated Rotating Pipe

Numerical investigation of laminarization of low Reynolds number turbulent flow in a heated rotating vertical pipe is carried out using a rotation-sensitized k-ω Shear Stress Transport (SST) model with modifications to include low Reynolds number flow effects. The results show that rotation accelerates laminarization in heated pipes and leads to as much as 72% deterioration in heat transfer for a typical inlet Reynolds number of 5000 and at a rotation number of 2. Rotation causes a reduction of the velocity gradients in the near-wall region. The lower gradients result in a reduction in the rate of production of turbulence kinetic energy, which causes the flow to laminarize. A laminarization map that relates the nondimensional heat flux as a function of inlet Reynolds number and rotation number is presented in this study. Correlation to predict the nondimensional wall heat flux required to laminarize the flow as a function of inlet Reynolds number and rotation number is proposed. The proposed laminarization map and the correlation which predicts the nondimensional heat flux required for laminarization to within an accuracy of ±8.1% accuracy for flow of air in a heated rotating vertical pipe would be of value to operators and designers of heat exchangers using air as a working medium to maintain a heat flux or a rotation number below the threshold value to prevent flow laminarization.

Tracking the Nearfield Evolution of an Initially Shallow, Neutrally Buoyant Plane Jet Over a Sloping Bottom Boundary

AbstractUnderstanding coastal plane jets which occur when a body of water discharges into an ocean or a lake through a channel or outlet is important, since they play a significant role in sediment, nutrient, and pollutant exchange. This study investigates the nearfield of initially shallow, neutrally buoyant plane jets, bounded by a free surface and a sloping bottom (Sloping Bottom Jet; SBJ) that issue into a laterally unconfined quiescent ambient both experimentally and numerically, and compares them with a plane jet flowing over a horizontal bottom (Horizontal Bottom Jet; HBJ). Results revealed that, different from the HBJ, the width and centerline velocity of SBJ decrease near the mouth. The SBJ width gradually increases after that as the transverse longitudinal velocity profile progressively transforms from a “top‐hat” into a Gaussian distribution. Once the Gaussian distribution is established, both jets diverge and centerline velocity decreases. Shear layers are generated on the sides of both jets with Kelvin Helmholtz‐type Coherent Structures (KHCS) developing inside. KHCS produce periodic velocity fluctuations with a Strouhal number of ∼0.079 and contribute significantly to momentum exchange and turbulent kinetic energy production. Since the thickness of the SBJ increases longitudinally, the vertical extent of KHCS also increases. When the two shear layers meet and merge at the centerline, they cause a flapping motion of the jet. This location is closer to the jet mouth for SBJs than for the HBJ. These findings demonstrate that a sloping bottom modifies the flow field from quasi‐2D for the HBJ to strongly 3D for SBJs.

Assessment of Laws of the Wall During Flame–Wall Interaction of Premixed Flames Within Turbulent Boundary Layers

The validity of the usual laws of the wall for Favre mean values of the streamwise velocity component and temperature for non-reacting flows has been assessed for turbulent premixed flame-wall interaction using Direct Numerical Simulation (DNS) data. Two different DNS databases corresponding to friction velocity-based Reynolds number of 110 and 180 representing unsteady head-on quenching of statistically planar flames within turbulent boundary layers have been considered. The usual log-law based expressions for the Favre mean values of the streamwise velocity and temperature for the inertial layer have been found to be inadequate at capturing the corresponding variations obtained from DNS data. The underlying assumptions of constant shear stress and the equilibrium of production and dissipation of turbulent kinetic energy underpinning the derivation of the usual log-law for the mean streamwise velocity have been found to be rendered invalid within the usual inertial layer during flame-wall interaction for both cases considered here. The heat flux does not remain constant within the usual inertial layer, and the turbulent flux of temperature exhibits counter-gradient transport within the so-called inertial layer for the cases considered in this work. These render the assumptions behind the derivation of the usual log-law for temperature to be invalid for application to turbulent flame-wall interaction. It has been found that previously proposed empirical modifications to the existing laws of the wall, which account for density and kinematic viscosity variations with temperature, do not significantly improve the agreement with the corresponding DNS data in the inertial layer and the inaccurate approximations for the kinematic viscosity compensated wall normal distance and the density compensated streamwise velocity component contribute to this disagreement. The DNS data has been utilised here to propose new expressions for the kinematic viscosity compensated wall normal distance and the density compensated streamwise velocity component, which upon using in the empirically modified law of wall expressions have been demonstrated to provide reasonable agreement with DNS data.

Direct numerical simulations of supersonic flat-plate turbulent boundary layers with uniform blowing

The effect of uniform blowing on a spatially developing flat-plate turbulent boundary layer at Mach 2.25 is investigated using direct numerical simulations. Two values of the wall blowing ratio are considered, corresponding to low and high blowing rates. Uniform blowing is found to significantly reduce the near-wall turbulence anisotropy, although the turbulent kinetic energy still exhibits near-wall asymptotic behavior and the Reynolds analogy is relatively insensitive to changes in the blowing ratio. The pre-multiplied spectra of turbulent kinetic energy production demonstrate that increasing the blowing ratio significantly energizes the large-scale structures in the outer region, while suppressing the inner small-scale structures. An increase in the blowing ratio also has a strong influence on the behavior of the fluctuating wall pressure, amplifying the fluctuation intensity and reducing the dominant frequencies in the power spectrum. Two-point space–time correlations indicate that the characteristic length scale of the pressure fluctuations increases with increasing blowing ratio, whereas the convection velocity exhibits the opposite trend. Analysis of the reduced mean wall heat flux reveals that it is dominated by the relative balance between the work of the Reynolds stress and the turbulent transport of heat, but is insensitive to uniform blowing. Importantly, bidimensional empirical mode decomposition of the turbulent structures highlights the increasingly dominant contributions related to the significantly energized outer large-scale structures in the blowing region.

Direct numerical simulations of high-enthalpy supersonic turbulent channel flows including finite-rate reactions

Direct numerical simulations of temporally evolving high-enthalpy supersonic turbulent channel flows are conducted at a Mach number of 3.0 and Reynolds number of 4880 under isothermal wall conditions. Air is assumed to behave as a five-species mixture, and chemical non-equilibrium and equilibrium assumptions are adopted to investigate the influence of finite-rate reactions on the turbulent statistics and large-scale structures. The two wall temperatures of 1733.2 and 3500 K are such that the mixture components undergo strong dissociation and recombination reactions along the channel. Investigation shows that the turbulent intensity is weakened and the mean and fluctuating temperatures are smaller when finite-rate reactions are considered. The mean dissociation degree is a quadratic function of the normal position, and its curvature under the chemical non-equilibrium assumption is greater than that under the chemical equilibrium assumption. The fluctuating mass fractions of the generated species seem to decrease slightly in the near-wall region, and their distributions are obviously different from those of the fluctuating velocity and fluctuating temperature. Finite-rate reactions increase the proportion of turbulent kinetic energy production in the skin friction decomposition, especially when the wall temperature is 3500 K. The large-scale structures visualized by the cross correlation between temperature and species mass fraction become stronger in the normal direction. The turbulent Schmidt number and several velocity–temperature correlations, including the recovery enthalpy and strong Reynolds analogy, are insensitive to the chemical reaction rate and wall temperature.