Channel Flow Research Articles

A local synergy angle of velocity, temperature and pressure fields for enhancement of comprehensive heat transfer performance

The field synergy principles for convective heat transfer and flow resistance are integrated in this study. A local three-field synergy angle of velocity, temperature and pressure fields is proposed to evaluate the local comprehensive performance of heat transfer and pumping power consumption. The angle is based on the polar angle of a plot with dot products of velocity and enthalpy/total pressure gradients, respectively, as coordinates. The polar angle is then shifted by the observation that the enhancement of heat transfer rate at the cost of pumping power decreases from the second quadrant to the forth quadrant. The application of the synergy angle is demonstrated by examples of flow and heat transfer in channels with different fin structures. Inspired by the distribution of the synergy angle, Airfoil-Rect fin and Rhom-Rect fin structures are designed, which can improve the heat transfer rate with constraints of pumping power consumption. The present study provides a possible approach for heat transfer enhancement.

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.

On hydromagnetic two-phase gas-liquid flow in ciliary channel: An application of a metachronal rhythm

On hydromagnetic two-phase gas-liquid flow in ciliary channel: An application of a metachronal rhythm

Optimization of heat transfer in a channel with stretching walls using a magnetized tetra-hybrid nanofluid

Nanoparticles have numerous applications and are used frequently in different cooling, heating, treatment of cancer cells and manufacturing processes. The current investigation covers the utilization of tetra hybrid nanofluid (aluminium oxide, iron dioxide, titanium dioxide and copper) for Crossflow model over a vertical disk by considering the shape effects (bricks, cylindrical and platelet) of nanoparticles, electro-magneto-hydrodynamic effect and quadratic thermal radiation. The present study is devoted to present the mathematical formulation of the tetra-hybrid nanofluid flow in a porous channel with stretching/shrinking walls. The present inspection model is derived from given partial differential equations (PDEs) and then transformed into a system of ordinary differential equations (ODEs) by incorporating similarity variables. The transformed ODEs are solved using the bvp4c methodology, which yields numerical results. From the obtained results it is observed that the maximum amount of [Formula: see text] happens when there is slightly increase in [Formula: see text] with stretched wall that is [Formula: see text]. A high radiation value indicates a considerable contribution from radiative heat transfer, whereas stretching and shrinking increase and reduce the size of the boundary. The combined impact caused an increase in the Nusselt number. The graphical and tabular representations are used to explain the physical behaviour of various model parameters. Previous outcomes are also contrasted with the current outcomes.

Impact of dual-fracture location on heat extraction from Enhanced geothermal system in low-permeability reservoirs

As environmental and economic issues exacerbated by traditional fossil fuels intensify, Dry Hot Rock geothermal energy has garnered significant attention due to its immense potential. Current research predominantly focuses on the heat extraction effectiveness in complex fractures within Enhanced Geothermal Systems, with relatively less emphasis on fluid behavior between large fractures. This study employs finite element software to establish a multi-horizontal well model in a low permeability reservoir (5 × 10-17 m2), analyzing the impact of different spatial configurations of two fractures on fluid flow, heat transfer, and heat extraction efficiency. The results indicate that fractures are the primary channels for fluid flow in low permeability reservoirs. Optimal heat extraction occurs when fractures are parallel, spaced 100 m apart, and perpendicular to the horizontal well, with significant threshold effects of fracture spacing and angle on temperature distribution. The highest heat extraction is achieved at a 37° intersection angle of fractures at the production well. Over a 60-year production period, reservoir extraction degrees for parallel or intersecting fractures range from 8 % to 15 %, while that for fracture communication between injection and production well horizontal segments is only 1.91 % to 3.01 %. Under the optimal injection scenario, cumulative 60-year heat production for parallel, intersecting, and communicating fractures is 2.13 × 1018J, 2.20 × 1018J, and 2.22 × 1018J, respectively, with the best heat extraction efficiency when fractures communicate. Additionally, under constant flow rate and inlet temperature, outlet temperature and system thermal power show near-linear declines, while reservoir extraction degree rises linearly. This study provides crucial theoretical support and practical guidance for the efficient extraction of geothermal energy.

Reversible to irreversible transitions for ac driven skyrmions on periodic substrates

Abstract Using atomistic simulations, we investigate the dynamical behavior of magnetic skyrmions in dimer and trimer molecular crystal arrangements, as well as bipartite lattices at 3/2 and 5/2 fillings, under ac driving over a square array of anisotropy defects. For low ac amplitudes, at all fillings reversible motion appears in which the skyrmions return to their original positions at the end of each ac drive cycle and the diffusion is zero. We also identify two distinct irreversible regimes. The first is a translating regime in which the skyrmions form channels of flow in opposing directions and translate by one substrate lattice constant per ac drive cycle. The translating state appears in the dimer and trimer arrangements, and produces pronounced peaks in the diffusivity in the direction perpendicular to the external drive. For larger ac amplitudes, we find chaotic irreversible motion in which the skyrmions can randomly exchange places with each other over time, producing long-time diffusive behavior both parallel and perpendicular to the ac driving direction.

Dual-functional adsorptive membranes for PFAS removal: Mechanism, CFD simulation, and selective enrichment

The remediation of emerging water contaminants, particularly per- and polyfluoroalkyl substances (PFAS), presents challenges due to their refractory nature and the presence of competing substances. Dual-functional adsorptive membranes, with hydrophobic backbone and quaternary ammonium moieties, were thereby designed to selectively intercept organic competitors while enrich PFAS. A 96.8 % removal of perfluorooctanoic acid (PFOA) was achieved and this effective removal (>90 %) maintained across five reuse cycles with a total treatment capacity of 650 L m−2. Rather than the rejection mechanism of nanofiltration process, these adsorptive membranes utilize synergistic electrostatic attraction and hydrophobic interactions, leading to a greater enrichment factor of 18.5 (PFOA over humic acid) and a permeability of 34.6 L m-2h−1 bar−1 (1.9- and 4.5-fold higher than reported NF 270 membranes, respectively). Furthermore, computational fluid dynamics (CFD) modeling revealed that the sponge-like matrix effectively prevents channeling flow and enhance access to adsorption sites. Sensitivity analysis and the high Damkohler number indicated that the adsorption process is mass transfer-controlled, with the key parameters ranked in order of significance: residence time > fluid viscosity > intrinsic adsorption rate. With consistent removal performance with co-existing competitors, efficient regeneration, and reusability, the dual-functional adsorptive membranes offer promising practical efficacy for PFAS remediation.

An open source electrochemical channel flow cell setup for kinetics studies. Application to investigations on oxygen electrocatalysis.

An open source electrochemical channel flow cell setup for kinetics studies. Application to investigations on oxygen electrocatalysis.

Modeling of non-equilibrium suspended sediment transport process in open channel flows with vegetation

A mathematical model based on advection-diffusion theory is established to study the non-equilibrium sediment transport process in vegetated channels. The effects of vegetation on velocity distribution and sediment diffusion coefficients were considered, respectively. Validation against experimental data from flume studies confirms the model's ability to accurately predict the longitudinal sediment deposition rate and the vertical distribution of suspended sediment concentration (SSC). A comparative analysis of three sediment diffusion coefficient formulations indicates that the linear-exponential formula provides a more precise estimate of εsz, and the linear-exponential formula performs well in predicting the turbulent diffusion coefficients of both rigid and flexible vegetation when gently swaying. Moreover, the distance required for SSC to regain equilibrium is influenced by the submergence level of the vegetation canopy. At lower submergence levels, the canopy shear vortices significantly affect the vertical exchange of sediment, and the sediment diffusion coefficients exhibit pronounced stratification near the vegetation canopy. An increase in vegetation density at these lower submergence levels intensifies the shear vortices, thereby extending the distance needed for SSC to reach equilibrium. At higher submergence levels, the impact of canopy shear vortices is lessened, which reduces sediment diffusion coefficient stratification characteristics, and the flow is similar to rough boundary layer flow. An increase in vegetation density increases flow resistance, which shortens the distance required for SSC to attain equilibrium. However, further efforts are required to explore turbulent characteristics with highly flexible vegetation motion and the grain size distribution of non-uniform sediments in vegetated flows.

Numerical investigation of turbulent flow over a randomly packed sediment bed using a variable porosity continuum model

A large eddy simulation (LES) is performed for a turbulent open channel flow over a porous sediment bed at permeability Reynolds number of ReK∼2.56 (Reτ = 270) representative of aquatic systems. A continuum approach based on the upscaled, volume-averaged Navier−Stokes (VaNS) equations is used by defining smoothly varying porosity across the sediment water interface (SWI) and modeling the drag force in the porous bed using a modified Ergun equation with Forchheimer corrections for inertial terms. The results from the continuum approach are compared with a pore-resolved direct numerical simulation (PR-DNS) in which turbulent flow over a randomly packed sediment bed of monodispersed particles is investigated [Karra et al.,J. Fluid Mech. 971, A23 (2023)] A spatially varying porosity profile generated from the pore-resolved DNS is used in the continuum approach. Mean flow, Reynolds stress statistics, and net momentum exchange between the freestream and the porous bed are compared between the two studies, showing reasonably good agreement. Small deviations within the transitional region between the sediment bed and the freestream as compared to the PR-DNS results are attributed to the local protrusions of particles in a randomly packed bed that are absent in the continuum approach but are present in the PR-DNS. A better representation of the effective permeability in the top transition layer that accounts for roughness effect of exposed particles is necessary. The continuum approach significantly reduces the computational cost, thereby making it suitable to study hyporheic exchange of mass and momentum in large scale aquatic domains with combined influence of bedform and bed roughness.

Sediment transport on rippled beds

We conduct an Euler-Lagrange, direct numerical simulation of a turbulent channel flow at a shear Reynolds number of Reτ=180 over an erodible particle bed. The particle bed consists of approximately 1.3 × 106 monodisperse particles, resulting in a bed thickness of around 12–13 particles. The particle density and size are chosen to achieve a ratio of 4 for the Shields stress to the critical Shields stress necessary for incipient motion such that particle transport occurs primarily as bedload. The simulation is run long enough for ripples to form. We track the temporal evolution of the particle flux and excess Shields stress for the entire bed as well as for the four regions of a ripple, namely, the crest, trough, lee side, and stoss side. We find that the particle flux and excess Shields stress closely match the Wong and Parker correlation when the particle bed is featureless at early time but diverge from the correlation when ripples form. This deviation primarily arises from particle transport in the trough and lee side regions. Conversely, particle transport in the crest and stoss side regions remains largely consistent with the Wong and Parker correlation. A root mean square-based correction for the bed is proposed to be used in conjunction with the Wong and Parker correlation. Additionally, ripples attain a self-similar profile in the shape and near-bed shear stress when they are sufficiently distant from their upstream neighbor. Any departure from self-similarity occurs when the upstream neighbor gets within close proximity.

Wall temperature effects in high-enthalpy supersonic turbulent channel flows considering air dissociation

Wall temperature effects in high-enthalpy supersonic turbulent channel flows considering air dissociation

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.

Vertical concentration distribution of fine settling particles in a pulsatile laminar open channel flow

Vertical concentration distribution of fine settling particles in a pulsatile laminar open channel flow

Wake dynamics of side-by-side hydrokinetic turbines in open channel flows

Lateral placement of hydrokinetic turbines is an interesting topic, as the blockage effect can increase the flow speed and increase the power coefficient (CP) for neighboring turbines. This study investigates wake dynamics in hydrokinetic turbine arrays with single- (1T), double- (2T), and triple-turbine (3T) configurations under various tip speed ratios (λ = 3.5, 5.8, and 7.1) using large eddy simulation coupled with the actuator line (AL) model. Results indicate that CP increases as lateral spacing decreases, which highlights the advantages of tighter lateral placement. The CP of the 3T-S turbine (the side turbine in the 3T configuration) is larger than those of the other configurations, following the trend CP,3T−S>CP,3T−M>CP,2T>CP,1T, which reflects a growing blockage effect with more turbines. Wake dynamics are analyzed using time-averaged and instantaneous methods. In 3T scenarios, blockage enhances turbulence kinetic energy, facilitating faster wake recovery, aided by turbine interference. Mean kinetic energy budget analysis shows that 3T-S wakes recover fastest due to increased turbulent convection. For instantaneous analysis, pre-multiplied power spectral density reveals vertical meandering begins at approximately 3D (D is the rotor diameter) and horizontal meandering starts near 4D, with a dominant frequency of St=0.28. Integral length scales show an initial increase followed by a downstream decrease, with minima marking the onset of wake meandering. Dynamic mode decomposition analysis reveals that high-frequency disturbance amplitudes increase with the number of turbines. At the optimal λ, wake effects dominate over inflow effects.

Interaction between hydrodynamic and morphodynamic processes at the confluence of pipe and channel

Confluences between a pipeline and a channel are commonly found in urban drainage networks. The mechanisms of sediment transport and the coherent structure within the confluence zone remain ambiguous. The hydro-morpho-sedimentary processes at pipe and channel confluences, which are influenced by discharge ratios, are investigated using a laboratory-scale confluence. The pipe current compresses and diverts downstream upon entering the main channel, creating a scour hole that carries sediment downstream. Coarser sediment deposits first, forming the slope of the downstream boundary of the scour hole (D50 = 60.31 μm) and the inner bar (D50 = 62.67 μm). Finer sediments are then deposited in the outer bar (D50 = 53.12 μm), resulting in a bed morphology consisting of the scour hole, outer and inner bars, and corridor. The penetration of the pipe current causes redistribution of the main channel flow, creating two shear layers. The height of the bar is directly related to the intensity of turbulence. Within the range of discharge ratios from 0.06 to 0.09, the shear layers and bars are displaced toward the outer bank by a distance of approximately one pipe diameter (0.05 m). Turbulence in this area is not isotropic, with power-law relations for streamwise and lateral velocity having an exponent of approximately −5/3. Conversely, the separation zone deviates from this pattern. Overall, this study highlights the significant influence of the bed morphology in modifying the flow structure compared to channel confluences and jets into crossflow, which are important in drainage engineering designs and fluvial studies.

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.

Analysing Turbulence Patterns in Nature‐Like Fishways: An Experimental Approach

ABSTRACTNature‐like fishways are an important measure for restoring fish passage, mitigating the impacts of barriers, and providing valuable habitat benefits in a more natural and effective way than traditional technical fishways. Their ecological advantages make them a preferred solution in many river restoration projects. The turbulence at the fishway's pool is a key challenge in the design of fishways. This study evaluated the hydraulic performance of the nature‐like fishway section of the Zongyang Fishway Project in China. This utilized a trapezoidal cross‐section with a 1:242 bottom slope, 1:2 side slopes, and an operating water depth of 1–3 m. The bottom and sides were paved with 0.2–0.5 m diameter pebbles to create suitable habitat for bottom‐dwelling fish. The target design velocity was 0.7–0.9 m s−1. A physical model was constructed to assess the hydraulic characteristics of the nature‐like fishway, including flow velocities within the slots and turbulence levels in the pools. The results showed that the slot velocities were within the target range, the head loss was parallel to the bottom elevation, and the flow field near the bottom mimicked technical fishway conditions, whereas the surface flow resembled open channel flow. Turbulence intensity remained below 60% of the design velocity. These findings provide valuable insights into the hydraulic performance and suitability of this nature‐like fishway design for facilitating fish passage of the target freshwater species along with directly aiding to SDGs 6 and 14.

An experimental investigation of boundary layer over permeable interfaces in Hele-Shaw micromodels

This study experimentally investigates boundary layer development over permeable interfaces using Hele-Shaw micromodels and high-resolution micro-particle image velocimetry (micro-PIV). Velocity vectors, captured at a 5 μm scale, reveal the flow behavior at the interface between free-flow and porous media with ordered structures and porosities ranging from 50% to 85%. The results show that the boundary layer streamline alignment decreases with increasing porosity, while lower permeability fosters more uniform and parallel flow near the interface. Flow channeling occurs along paths of the least resistance, with more flow directed through the Hele-Shaw free-flow region as the solid fraction of the porous material increases. The Reynolds number (0.14–0.94), based on the Hele-Shaw hydraulic diameter, has a minimal effect on the normalized velocity distribution. Furthermore, an analytical solution for the external boundary layer thickness exhibited good agreement with experimental data, confirming a thickness of 2–4 times the square root of the free-flow Hele-Shaw permeability. Additionally, a Q-criterion analysis identified, for the first time, distinct zones within the external boundary layer, capturing the balance between rotational and deformation components as a function of permeability. These findings offer insight into flow dynamics in porous media systems, with implications for both natural and industrial applications, and contribute to the improved modeling of fluid dynamics and momentum transport in coupled free-flow and porous media environments.

Thermo-hydraulic performance evaluation of lattice structures with triply periodic minimal surfaces for latent heat storage devices

Triply periodic minimal surface (TPMS) structures exhibit significant potential in latent heat storage (LHS) applications. In this study, numerical investigations were conducted to study the thermal storage performance of four TPMS lattice types: Primitive, IWP, FRD, and Gyroid lattices, used as metallic skeletons in LHS devices. The volumes of both the TPMS skeletons and the phase change material (PCM) were kept constant, with PCM embedded either inside or outside the lattices. The melting process of PCM, the heat storage energy density, and the power density of the LHS units for different configurations were examined. Additionally, a performance evaluation coefficient was introduced to assess the heat transfer and flow characteristics of heat transfer fluid (HTF) flowing on both sides of the lattices. The results indicated that when using Primitive and FRD lattices, embedding the PCM within the internal space of the lattices resulted in a 20 % and 9.8 % higher energy storage compared to embedding it outside the lattices. The internal space of IWP and Gyroid lattices was more suitable for HTF flow channels, resulting in performance evaluation coefficients (PECs) that improved by 27.7 % and 86.3 %, respectively, when HTF flowed through the internal channels of these lattices. Moreover, the Gyroid lattice exhibited the most balanced overall heat transfer performance due to its minimal flow resistance, whereas the FRD lattice had the lowest flow performance due to its complex structure. Therefore, the Primitive and Gyroid lattices could be considered as new heat transfer structures for LHS devices, while the FRD lattice was more suitable for replacing traditional fin structures in heat sink devices.