Friction Reynolds Number Research Articles

Direct numerical simulation of turbulent heat transfer over surfaces with hemisphere protrusions

Direct numerical simulations (DNSs) of turbulent heat transfer over walls with regularly distributed hemispheric protrusions were conducted to explore how the arrangement and number density of hemispheres affect the turbulent heat transfer. For the rough walls, we considered eight rough surfaces, in which the number density, i.e., the number of hemispheres per unit area, and the distances between two neighboring hemispheres in the streamwise and spanwise directions were systematically varied. The friction Reynolds number was fixed at 660, and we considered an incompressible airflow with a Prandtl number of 0.71, neglecting the buoyancy effects. The results showed that the velocity roughness function strongly depends on the hemisphere arrangement and number density. The spanwise-aligned hemisphere array yields the larger velocity roughness function than its streamwise-aligned hemisphere counterpart, whereas the temperature roughness function depends only on number density. The Reynolds analogy factor decreases with increasing inner-scaled equivalent sand grain roughness, and the decreasing trend is weakly affected by the number density. We analyzed the momentum and heat transfer budgets and found that the pressure drag, which dominates the momentum transfer near the wall, is strongly affected by the hemisphere arrangement and number density. In contrast, the roughness-induced wall heat transfer term, which is the dominant contributor to near-wall heat transfer, depends only on the number density. This difference is the driving mechanism that causes dissimilar trends in momentum and heat transfer over walls with hemispheric protrusions.

Direct numerical simulation of sodium in vertical channel flow: From forced convection to natural convection at friction Reynolds number 180

The buoyancy-aided sodium flow in a vertical channel is investigated using direct numerical simulation (DNS) to study turbulent flow and heat transfer at six different Richardson numbers (Ri = 0, Ri = 0.025, Ri = 0.25, Ri = 2.5, Ri = 7.5, and Ri = 15) with a fixed friction Reynolds number (Reτ = 180). The results reveal that the velocity profile shows an “M” shape under buoyancy effect and reverses at the center under strong buoyancy. Additionally, the temperature profile exhibits a thicker boundary layer compared to the velocity profile. Global coefficients, such as the skin friction coefficient and the Nusselt number, are analyzed using Fukagata, Iwamoto, and Kasai (FIK) decomposition to elucidate their respective contributions. Furthermore, anisotropy analysis indicates that buoyancy makes the turbulence more isotropic, and the buoyancy also makes the turbulent Prandtl number (Prt) unpredictable; however, a comparison among the molecular heat flux, the definition of turbulent heat flux, and the calculation of the standard gradient diffusion hypothesis (SGDH) model suggests that the turbulent heat flux can be neglected without significant influence in this study. Finally, the turbulent structures in the viscous layer are presented, and the quadrant analysis is performed to quantitatively analyze the influence of buoyancy on the turbulent structure.

The isolated effect of yield stress in viscoplastic turbulent flow

Turbulent flows of viscoplastic fluids are present in several industrial and natural applications. The effects of yield stress on this problem have always been studied as a part of a larger physical context, because real viscoplastic materials have many properties that cannot be easily isolated. Direct numerical simulations have recently emerged as a viable tool for investigating non-Newtonian fluid flow in turbulent regimes. In the present work, we solve the turbulent flow of an ideal Bingham fluid, focusing on the isolated effect of yield stress. A numerical scheme for viscoplastic flows was implemented based on the lattice Boltzmann method. An outstanding characteristic of this scheme is the possibility of representing infinite viscosity by setting the relaxation frequency to zero, enabling the representation of the Bingham constitutive equation without artifacts, and producing a more accurate representation of the yield surfaces. In the turbulent channel flow simulations, the friction Reynolds number was fixed at 180, while the Bingham number varied from 0 (Newtonian) to 0.15. It is shown that unyielded portions of material travel along with the flow near the centerline. These unyielded spots do not disappear quickly, but rather have a significant lifetime. Another interesting outcome is that the yield stress increases the turbulence anisotropy, by lowering the spanwise and normal velocity fluctuations, while the streamwise component becomes higher. Reynolds stresses and budgets of turbulent kinetic energy have been analyzed regarding the increased bulk velocities that were found by increasing the yield stress.

On the interplay between fluid flow characteristics and small particle deposition in turbulent wall bounded flows

This study investigates the transport and deposition of small particles (500 nm≤dp≤10 μm) in a fully developed turbulent channel flow, focusing on two fluid friction Reynolds numbers: Reτ=180 and Reτ=1000. Using the point particle–direct numerical simulation method under the assumption of one-way coupling, we study how fluid flow (carrier phase) characteristics influence particle deposition. Our findings suggest that changes in flow conditions can significantly alter the deposition behavior of particles with the same size and properties. Furthermore, we show for the first time that gravity has minimal impact on deposition dynamics only at high Reynolds numbers. This research enhances our understanding of small particle deposition and transport in turbulent flows at high Reynolds numbers, which is crucial for various industrial applications.

High-fidelity simulations of microramp-controlled shock wave/boundary layer interaction

Microvortex generators (MVGs) are a promising solution to control shock wave/turbulent boundary layer interactions (SBLIs), especially in supersonic inlets. In this study, we examine the effects of a microramp vortex generator on an SBLI generated by an oblique shock wave and a turbulent boundary layer using direct numerical simulations (DNSs). Two cases, with and without the presence of a microramp, are compared in terms of their mean and unsteady flow features at free-stream Mach number equal to 2 and friction Reynolds number at the inviscid shock impingement equal to 600. The long integration period allows us to assess how microramps affect the typical low-frequency unsteadiness observed in SBLIs, and the data generated may serve as a reference for simulations of lower fidelity or reduced order models. The analysis shows that the three-dimensional microramp wake alters the interaction region dramatically, inducing a significant spanwise modulation and topology change of the separation. For example, tornado-like structures redistribute the flow in both the spanwise and wall-normal directions inside the recirculation region. The increase in momentum close to the wall by the ramp vortices effectively delays the onset of the separation and, thus, the separation length, but at the same time leads to a significant increase in the intensity of the wall-pressure fluctuations. We then characterise the mutual interaction between the arch-like vortices around the ramp wake and the SBLI. The specific spanwise vorticity shows that these vortices follow the edge of the separation and their intensity, apart from mean compressibility effects, is not affected by the shocks. The shocks, instead, are deformed in shape by the periodic impingement of the vortices, although the spectral analysis did not reveal any significant trace of their shedding frequency in the separation region. These Kelvin–Helmholtz vortices, however, may be relevant in the closure of the separation bubble. Fourier analysis also shows a constant increase, in both value and magnitude, in the low-frequency peak all along the span, suggesting that the motion of the separation shock remains coherent while being disturbed by the arch-like vortices and oscillating at a higher frequency in absolute terms.

Physics informed data-driven near-wall modelling for lattice Boltzmann simulation of high Reynolds number turbulent flows

Data-driven approaches offer novel opportunities for improving the performance of turbulent flow simulations, which are critical to wide-ranging applications from wind farms and aerodynamic designs to weather and climate forecasting. However, current methods for these simulations often require large amounts of data and computational resources. While data-driven methods have been extensively applied to the continuum Navier-Stokes equations, limited work has been done to integrate these methods with the highly scalable lattice Boltzmann method. Here, we present a physics-informed neural network framework for improving lattice Boltzmann-based simulations of near-wall turbulent flow. Using a small amount of data and integrating physical constraints, our model accurately predicts flow behaviour at a wide range of friction Reynolds numbers up to 1.0 × 106. In contradistinction with other models that use direct numerical simulation datasets, this approach reduces data requirements by three orders of magnitude and allows for sparse grid configurations. Our work broadens the scope of lattice Boltzmann applications, enabling efficient large-scale simulations of turbulent flow in diverse contexts.

Investigation of the impact of transient pressure gradients on turbulent channel flow dynamics

AbstractThe effect of transient pressure gradients, or transient pumping in turbulent channel flow configuration, is investigated. Employing a cost‐efficient reduced‐order stochastic method known as one‐dimensional turbulence (ODT) modeling, simulations explore different signal shapes for the modulation of the prescribed pressure gradient forcing, including sinusoidal modulation, step‐like modulation, and piecewise sinusoidal beating. Various cycle periods and active pumping times are investigated. The simulations are conducted at a frictional Reynolds number of and a molecular Prandtl number of . The study adopts a passive scalar formulation to investigate heat transfer properties. The study quantifies the effects of transient pressure gradients on heat transfer rate, drag, and pumping power. Preliminary ODT predictions suggest that all transient cases exhibit lower heat transfer rates and a higher pumping power requirement than the constant pressure gradient case, with the step‐like modulation yields superior skin‐friction drag and heat transfer rate reductions relative to other signals.

Streamwise energy-transfer mechanisms in zero- and adverse-pressure-gradient turbulent boundary layers

The present study investigates streamwise ( $\overline {u^2}$ ) energy-transfer mechanisms in the inner and outer regions of turbulent boundary layers (TBLs). Particular focus is placed on the $\overline {u^2}$ production, its inter-component and wall-normal transport as well as dissipation, all of which become statistically significant in the outer region with increasing friction Reynolds number ( $Re_{\tau }$ ). These properties are analysed using published data sets of zero, weak and moderately strong adverse-pressure-gradient (APG) TBLs across a decade of $Re_{\tau }$ , revealing similarity in energy-transfer pathways for all these TBLs. It is found that both the inner and outer peaks of $\overline {u^2}$ are always associated with local maxima in the $\overline {u^2}$ production and its inter-component transport, and the regions below/above each of these peaks are always dominated by wall-ward/away-from-wall transport of $\overline {u^2}$ , thereby classifying the $\overline {u^2}$ profiles into four distinct regimes. This classification reveals existence of phenomenologically similar energy-transfer mechanisms in the ‘inner’ and ‘outer’ regions of moderately strong APG TBLs, which meet at an intermediate location coinciding with the minimum in $\overline {u^2}$ profiles. Conditional averaging suggests existence of similar phenomena even in low $Re_{\tau }$ canonical and/or weak APG TBLs, albeit with the outer-region mechanisms weaker than those in the inner region. This explains the absence of their $\overline {u^2}$ outer peak and the dominance of $\overline {u^2}$ wall-normal transport away from the wall, which potentially originates from the inner region. Given that the wall-ward/away-from-wall transport of $\overline {u^2}$ is governed by the $Q_4$ (sweeps)/ $Q_2$ (ejections) quadrants of the Reynolds shear stress, it is argued that the emergence of the $\overline {u^2}$ outer peak corresponds with the statistical dominance of $Q_4$ events in the outer region. Besides unravelling the dynamical significance of $Q_2$ and $Q_4$ events in the outer region of TBLs, the present analysis also proposes new phenomenological arguments for testing on canonical wall-turbulence data at very high $Re_{\tau }$ .

Multiscale analysis of Reynolds stresses and its dissipation rates for premixed flame–wall interaction

Direct numerical simulation (DNS) data of flame–wall interaction (FWI) has been utilized to analyze the multiscale nature of turbulent Reynolds stresses and dissipation rate tensor anisotropies within turbulent reacting flow boundary layers across a broad range of scales. The DNS data of head-on quenching of premixed flames propagating through turbulent boundary layers, representative of friction Reynolds numbers Reτ of 110 and 180, has been explicitly filtered using both two- and three-dimensional Gaussian filter kernels for the purpose of multiscale analysis. The low-pass filter results demonstrate the transition from a 2-component limit to a 1-component limit near the wall with increasing filter width, accompanied by a decrease in isotropy, suggesting a significant alteration in dominant flow patterns and a diminishing tendency toward isotropy. The high-pass filter results indicate a progressive increase in anisotropy with the progress of FWI at the channel center, emphasizing the anisotropy of the large scales with the progress of FWI. Furthermore, behaviors of the second and third invariants of the Reynolds stress tensor remain qualitatively similar to that of the dissipation rate tensor at all stages of FWI, suggesting a link between viscous dissipation and Reynolds stress distributions; notably, there is a stronger isotropic tendency in the dissipation rate tensor when the flame is away from the wall, intensifying with an increase in Reynolds numbers. However, as FWI progresses, the shift in the trend toward the 1-component limit indicates an increase in anisotropy within the turbulent reacting flow for the region near the center of the channel.

Predicting Heat Transfer Enhancement with Twisted Tape Inserts Using Fuzzy Logic Techniques in Heat Exchangers

Fuzzy logic, introduced by Lotfi Zadeh in 1965, is a powerful method for modelling complex experiments. This study utilizes fuzzy logic to simulate and predict heat transfer in a double-pipe heat exchanger equipped with wavy inserts. The inserts, in the form of twisted tapes, have varying twist ratios (TR=9, 7, 6). The study investigates a range of Reynolds numbers (Re) from 6000 to 18000, with friction factors ranging from 0.03620 to 0.08231, and Nusselt numbers (Nu) between 66.13 and 253.28. The results for different twist ratios are compared to the ideal case. The experimental results indicate that the highest heat transfer occurs with a twist ratio of 6, leading to a significant increase of 162% in the Nusselt number and a 36.21% rise in the friction factor compared to the ideal scenario. In the fuzzy logic framework, the input variables are the twist ratio (Tr), temperature, and Reynolds number (Re), while the output variables are the friction factor (f) and Nusselt number (Nu). The study demonstrates that the Mamdani fuzzy inference system is an exceptionally effective tool for predicting experimental outcomes, given its low error rate. Upon analysing the data, it is observed that the graphs plotting the Nusselt number versus Reynolds number and friction factor versus Reynolds number, derived from both experimental data and the fuzzy logic model, exhibit nearly identical trends with a margin of error of just 3%. This high level of accuracy underscores the reliability of the fuzzy logic model in replicating the experimental results.

Flame dynamic insights into pollutant production characteristics of NH[formula omitted]/H[formula omitted]/air combustion during turbulent flame–wall interaction

A turbulent premixed NH3/H2/air flame with Damköhler number Da = 0.004 has been simulated using direct numerical simulations (DNS) with detailed chemistry. The flame is propagating in a fully-developed turbulent channel flow with friction Reynolds number Reτ = 313. The global and local production speeds of NO, N2O and NH3 have been analyzed. Special attention is paid on the flame dynamic insights into the pollutant production characteristics during the flame–wall, flame–turbulence and flame–flame interaction process. Finally, important findings have been concluded: (1) The unburned NH3 is reduced in the turbulent case due to the smaller quenching distance; (2) The N2O consumption speed would decrease just before flame quenching on the wall in the turbulent case, due to the thickened flame thickness, which is contradictory to the laminar case; (3) Effect of flame curvature on pollutant production speeds should be combined with the effect of flame surface density; (4) Low Damköhler number flame tends to produce more N2O and less NO, due to the more frequent flame–flame interaction events. These findings would help to better understand the pollutant production characteristics of turbulent premixed NH3/H2/air flames in lab-scale burners with strong flame–wall interactions.Novelty and Significance StatementThe novelty of this research is summarized into 3 points: (1) It is, to our knowledge the first time, pollutant production characteristics are investigated in detail for the premixed NH3/H2/air flame propagating in a fully developed turbulent boundary layer; (2) The effects of flame–wall, flame–turbulence and flame–flame interactions, the underlying reasons and controlling mechanisms are investigated in detail, which has not yet been (comprehensively) done in previous studies; (3) Novel and important findings have been concluded for the wall effects, turbulence effects (including flame curvature and surface density function effects) and flame–flame interaction effects. Some of these findings have never been reported before, but important for the comprehensive understanding of the pollutant production characteristics. It is significant because understanding the turbulence and wall effects on NO, N2O and NH3 production characteristics in the premixed NH3/H2/air flames is crucial for the future clean combustion energy system design.

Large-scale control in turbulent flows over surface riblets

The drag reduction efficacy of a large-scale flow control over a rough surface is studied via direct numerical simulations of turbulent channels (at friction Reynolds numbers Reτ=180) by combining together wall riblets and streamwise counter-rotating swirls. In particular, the height of triangular riblets is h+≈10 (+indicating wall units), while the number of riblets (NRib in the range 1–56) along the periodic spanwise direction is varied to find the optimum. The swirls are generated by the spanwise opposed wall-jet forcing (SOJF) in the Navier–Stokes equation, whose controlling parameters follow the optimal ones as for the smooth wall. In total, 12 cases of combined SOJF and riblets are performed to investigate the coupling effects between the two methods. We find a range of NRib=7–14 (with the spanwise width z+≈140−280) yields the largest drag reduction (up to 20%) for Reτ=180, much higher than riblets control only (about 3%). Compared to SOJF control only, riblets suppress the secondary swirls of SOJF hence decreasing drag, while the lateral and down washing motions of SOJF impinging on riblets would increase drag—the opposite two effects thus giving rise to an optimal. Through examinations on coherent structures, we elucidate that the attenuation of both large-scale coherent motions and small-scale random fluctuations leads to the net drag reduction. We conclude that large-scale control is a robust approach in the cases of rough surfaces, and the parameters can be selected for maximum drag reduction in each particular situation.

An improved viscoelastic k- ε -v2-f turbulence model for intermediate and high drag reduction regimes

In this work, an improved anisotropic k-ε-v2-f model based on the finite extensible nonlinear elastic model with the Peterlin approximation for viscoelastic channel flows is proposed. This model is tested using direct numerical simulation (DNS) data for friction Reynolds numbers (Reτ) in the range of 120–1000, friction Wiesenberg numbers (Wiτ) in the range of 25–116, viscosity ratios (β) in the range of 0.6–0.9, and maximum polymer extensibility values (L2) in the range of 900–14 400. The flow characteristics of viscoelastic fluids with various parameters obtained from the new model agree well with existing DNS results. By adding closures for the flow, shear, and transverse components, the incomplete prediction of nonlinear terms from interactions between the fluctuating components of the conformation and velocity gradient tensors is improved. Compared with DNS results, these closures can fully obtain each component and significantly improve the accuracy of the flow direction component in the intermediate and high drag reduction regimes. Furthermore, the model in this paper retains the advantages of an anisotropic model, does not require a damping function, is simple to construct, and is easily extended to a variety of bounded flows.

Harnessing an elastic flow instability to improve the kinetic performance of chromatographic columns

Despite decades of research and development, the optimal efficiency of slurry-packed HPLC columns is still hindered by inherent long-range flow heterogeneity from the wall to the central bulk region of these columns. Here, we show an example of how this issue can be addressed through the straightforward addition of a semidilute amount (500 ppm) of a large, flexible, synthetic polymer (18 MDa partially hydrolyzed polyacrylamide, HPAM) to the mobile phase (1% NaCl aqueous solution, hereafter referred to as “brine”) during operation of a 4.6 mm × 300 mm column packed with 10μm BEHTM 125 Å particles. Addition of the polymer imparts elasticity to the mobile phase, causing the flow in the interparticle pore space to become unstable above a threshold flow rate. We verify the development of this elastic flow instability using pressure drop measurements of the friction factor versus Reynolds number. In prior work, we showed that this flow instability is characterized by large spatiotemporal fluctuations in the pore-scale flow velocities that may promote analyte dispersion across the column. Axial dispersion measurements of the quasi non-retained tracer thiourea confirm this possibility: they reveal that operating above the onset of the instability improves column efficiency by greater than 100%. These experiments thereby suggest that elastic flow instabilities can be harnessed to mitigate the negative impact of trans-column flow heterogeneities on the efficiency of slurry-packed HPLC columns. While this approach has its own inherent limitations and constraints, our results lay the groundwork for future targeted development of polymers that can impart elasticity when dissolved in commonly used liquid chromatography mobile phases, and can thereby generate elastic flow instabilities to help improve the resolution of HPLC columns.

Effects of steepness on turbulent heat transfer over sinusoidal rough surfaces

We conducted a direct numerical simulation (DNS) study to investigate the impact of surface undulation steepness on rough wall turbulent heat transfer. The flow geometry was turbulent open-channel flow over three-dimensional sinusoidal rough surfaces. To examine the effects of steepness, we systematically varied the streamwise and spanwise wavelengths of the sinusoidal roughness while keeping the roughness height constant. The friction Reynolds number ranged from 180 to 600, and we considered a passive scalar with the fluid Prandtl number was 0.7, assuming air flow conditions. In the fully rough regime, the velocity roughness function is expressed as a function of the inner-scaled equivalent sand grain roughness ks+ independent of steepness, whereas the steeper surfaces with shorter wavelengths result in larger temperature roughness functions at the same ks+ value. Analysis of the physical mechanisms that increases the roughness function shows that the pressure drag primarily contributes to the increase in the velocity roughness function, while the temperature roughness function is mainly augmented by the roughness-induced wall heat transfer term, correlating with the steepness of the surface undulations. It is also suggested that the effective slope, which quantifies the steepness of rough surfaces, could improve the predictive accuracy of existing correlations for the temperature roughness function.

Direct numerical simulation study on wall-modeling of turbulent water channel flows with temperature-dependent viscosity

We discuss wall-modeling of turbulent heat transfer of water channel flows with temperature-dependent viscosity via direct numerical simulations. We considered a top-cooled wall (293[K]) and a bottom-heated wall (353[K]) and varied the friction Reynolds numbers from 300 to 1000. The fluid viscosity varied depending on the local fluid temperature, whereas the other physical properties were assumed to be constant. The results show that semi-local scaling based on the local viscosity and wall friction velocity reasonably accounts for the effects of variable viscosity on turbulent flows, except in the vicinity of the wall, where wall cooling intensifies the turbulent vortical motion, leading to increased semi-locally scaled eddy diffusivities compared with those near the heated wall. In the vicinity of the cooled wall, turbulent transport is enhanced by increased viscous transport, which transfers more turbulent kinetic energy toward the cooled wall. The effectiveness of semi-local scaling for wall-modeling was validated by performing a wall-modeled large-eddy simulation at Reτ=1000, where we incorporated the semi-local viscous length scale into the classical mixing-length model. The modified mixing-length model reasonably reproduced the effects of variable viscosity on turbulent flows.

Large-eddy simulation of shock train in convergent-divergent nozzles with isothermal walls

ABSTRACT High-order accurate Large-eddy simulations (LES) have been carried out for compressible, turbulent flow through two bell-shaped convergent-divergent nozzles- one with a short divergent portion and the other with a longer one. An explicit filtering approach based on approximate deconvolution method is used for the LES. The nozzles have isothermal walls and circular cross-sections. The incoming flow is a turbulent, subsonic fully developed pipe flow at centreline Mach number of around 0.4 and friction Reynolds number 216. An adaptive filter is used to capture the shocks appearing in the divergent portion of the nozzles. The complicated behaviour of the mean flow, Reynolds stresses, production and pressure-strain terms of the axial Reynolds stress budget in the shock train region is shown. The broadband nature of the shock train oscillations having multiple dominant frequencies in these nozzles is observed in the present LES.

Uniform momentum zones in turbulent channel flow

The characteristics of uniform momentum zones (UMZs) and UMZ interfaces in turbulent channel flow (TCF) at friction Reynolds number ranging from 600 to 2000 are investigated. The focus is on the Reynolds number effect and the difference among the turbulent channel flow, turbulent pipe flow (TPF), and turbulent boundary layer (TBL). The UMZ is identified using the approach of probability density histograms, a discrete approximation of the probability density function. We find that the average number of UMZs increases with the Reynolds number, similar to the findings in the TBL. The first group of UMZ in the TCF, namely the quiescent core, has a thicker thickness at a higher Reynolds number. Therefore, the UMZ closer to the channel centre significantly influences high Reynolds number wall turbulence by occupying a large portion of the channel. With increasing Reynolds number, the intensity of spanwise vorticity at the UMZ interface increases, while the thickness of the local maximum decreases. The dominance of prograde vorticity at the interface also increases with the Reynolds number. Moreover, the contortion of the UMZ interface is significantly enhanced at high Reynolds numbers than at low Reynolds numbers in both the streamwise and spanwise directions. Furthermore, the number of UMZ for a given snapshot is related to the overall velocity fluctuations in the subregion, increasing with large-scale Q2 events and decreasing with large-scale Q4 events. Comparing the three flow configurations, UMZs show similarities in their trends, along with two notable differences. First, the TPF has more UMZs compared to the TCF and TBL. Second, the innermost group in the TPF is considerably thicker than that of the higher ranks, a pattern absent in the TCF and TBL. The present findings advance the understanding of UMZs in TCF and demonstrate the similarities and differences between the cases of TCF and TBL.

Direct numerical simulation of open-channel flow over a heterogeneous particle bed at low relative submergence

This study investigates turbulent open-channel flows over beds of irregularly arranged particles, using direct numerical simulations at a friction Reynolds number of Reτ=300. Two distinct cases are examined: a polydisperse bed (P800) composed of multiple layers of randomly distributed spheres of varying sizes, and a monodisperse bed (M1015) formed by a random distribution of uniform sized spheres, with a bottommost single layer of varied-sized particles to introduce realistic randomness. Our investigation unveils a rich network of low- and high-speed streaks within the flow field, exhibiting distinctive behaviors in different bed configurations. The P800 case presents a poorly organized flow pattern induced by the varied particle sizes and arrangements, while the M1015 case shows a more regular flow pattern, marked by larger streaks. We also observe that total wall shear stress is substantially influenced by surface roughness-induced drag, extending beyond the effects documented in existing studies of open-channel flows. The present study reveals intricate secondary flow patterns over irregular particle beds. Large-scale circulations are discerned around particle crests in the P800 case and localized circulations with increased turbulence in the M1015 case. Furthermore, analysis of Reynolds stress tensor components indicates that roughness disrupts coherent turbulent eddies, consequently mitigating peak stress. We quantify correlations between drag force and local fluid velocity fluctuations. Notably, a larger deviation in drag is observed in the P800 case compared to M1015, accentuating the influence of particle size and distribution on fluid–particle interactions.

Turbulence modulation by charged inertial particles in channel flow

Large amounts of small inertial particles embedded in a turbulent flow are known to modify the turbulent statistics and structures, a phenomenon referred to as turbulence modulation. While particle electrification is ubiquitous in particle-laden turbulence and significantly alters particle behaviour, the effects of inter-particle electrostatic forces on turbulence modulation and the underlying physical mechanisms remain unclear. To fill this gap, we perform a series of point-particle direct numerical simulations of turbulent channel flows at a friction Reynolds number of approximately 540, laden with uncharged and charged bidisperse particles. The results demonstrate that, compared to flows laden with uncharged particles, the presence of inter-particle electrostatic forces leads to substantial changes in both turbulent intensities and structures. In particular, the inner-scaled mean streamwise fluid velocity is found to shift towards lower values, indicating a noticeable increase in fluid friction velocity. Turbulent intensities appear to be further suppressed through facilitating the particles to extract momentum from the fluid phase and increasing extra turbulent kinetic dissipation by particles. Importantly, the overall drag is enhanced by indirectly strengthening the contribution of particle stress, even though the contribution of the total fluid stress is decreased. On the other hand, the magnitude of the large-scale motions is weakened by simultaneously reducing turbulent production and increasing particle feedback around the scales of the large-scale motions. Meanwhile, the average streaky fluid structures in the streamwise–spanwise planes and inclined fluid structures in the streamwise–wall-normal planes become expanded and flattened, respectively.