Surface Flow Research Articles

The Impact of Upstream Static Deformation on Flow Past a Cylinder/Flare

Reynolds-averaged Navier–Stokes simulations are performed for supersonic turbulent flow over a cylinder/flare with upstream surface distortion representative of structural deformation induced via fluid–structural and fluid–thermal–structural behavior. Broad parametric analysis is carried out through the generation of Kriging-response surfaces from a database of general simulations. A posteriori simulations are then carried out at parametric combinations that correspond to extrema in the Kriging response surfaces to gain deeper insights into the interaction between the surface distortion and flow responses. Upstream distortions tend to decrease, rather than increase, the peak pressure and heat flux loads on the flare compared to an undeformed cylinder. Furthermore, decreases in these quantities reach up to O(10%) compared to up to O(1%) for increases. Integrated quantities over the flare are relatively insensitive to upstream distortion. The corner separation length is the most sensitive quantity to upstream distortion, with protrusions tending to increase the separation length and recessions reducing the separation length. Modifications in the separation length of up to 40% are observed. Reductions in peak loads tend to correspond to increases in the corner separation length. The movement of the surface distortion relative to the corner indicates a negligible impact beyond 1.5 distortion lengths from the corner, and the largest impact on the corner separation length occurs when distortion is directly adjacent. These results are an important step toward understanding and quantifying the impact of surface deformations on downstream components.

Effects of episodic stream dewatering on brook trout spatial population structure

Abstract Stream dewatering is expected to become more prevalent due to climate change, and we explored the potential consequences for brook trout (Salvelinus fontinalis) within a temperate forest ecosystem in eastern North America. We estimated fish density within stream pools (n = 386) from electrofishing surveys over 10 years (2012–2021) to compare a stream that exhibits episodic dewatering (Paine Run) against a stream of similar size that remains flow‐connected (Staunton River) within Shenandoah National Park, Virginia (U.S.A.). Annual surveys encompassed fluvial distances ranging from 2.6 to 4.4 km in each stream. Mean annual fish density (fish/pool m2) was not different between streams for juvenile or adult age classes, but spatial variation in density was greater in Paine Run for both age classes of fish. Paine Run also included a greater proportion of unoccupied pools than Staunton River and exhibited stronger spatial autocorrelation in fish density among nearby pools, suggesting dispersal limitation due to surface flow fragmentation. Fish density in pools increased during years with low summer precipitation, and this effect was observed in both streams but was stronger in Paine Run than Staunton River, further indicating the importance of fish movement into pools in response to low‐flow thresholds. Our results indicate the importance of pools as ecological refuges during low‐flow conditions and that episodic dewatering may affect extirpation risks for brook trout by sequestering more fish into fewer areas. Our findings also highlight the importance of hydrological variation within stream networks because downstream river gages could not predict the observed spatial heterogeneity in fish density or pool occupancy.

Shaking a container full of perfect liquid; a tractable case, a torus shell, exhibits a virtual wall

Manipulation (‘shaking’) of a rigid container filled with incompressible liquid starting from stationary generally results in some displacement, or mixing, of the liquid within it. If the liquid also has zero viscosity, a ‘perfect’ or ‘Euler’ liquid, Kelvin’s theorems dramatically simplify the flow analysis. Response is instantaneous; stop the container and all liquid motion stops. In fact an arbitrary manipulation can be considered as alternating infinitesimal translations and rotations of the container. Relative to the container, the liquid is stationary during every translation. Infinitesimal rotations (an infinitesimal vector along the rotation axis) resolve into three orthogonal components in the container frame. Each generates its own infinitesimal liquid displacement vector field. Their combined consequences are usually obscure. Rather than a volume flow, a surface flow in 3D is considerably easier, the liquid slipping freely in a shell, sandwiched between two nested closed surfaces with constant infinitesimal gap. The closedness avoids extra boundaries. The two dimensionality admits a scalar streamfunction determined by the container angular velocity vector. Manipulation of the container angular velocity at will, usually leads to an infinitely rich variety of area preserving re-configurations of the liquid. In particular, any chosen point of the liquid can be moved, in the shell container frame, to any other chosen point. However, for a torus (shell) with a small enough hole (diameter < 0.195 torus diameter), there exists a virtual wall, a hypothetical axial cylinder intersecting the torus. No matter how the torus is manipulated, liquid inside the cylinder stays inside; outside stays outside. The analysis, solving for the streamfunction, is based on the relative vorticity, and the conformal mapping of a torus to a (periodic) rectangle, which lead to a fairly simple convolution integral formula for the flow.

Parametric analysis of the performance of the SPH solutions of unsteady free-surface flow

ABSTRACT In this article, the performances of the SPH method to solve Shallow Water Equations SWEs with three investigation parameters were studied, such as the type of kernel functions, namely: cubic spline, Gaussian and quintic spline kernels, the number of particles used and the stabilization terms injected, specifically: Lax Friedrichs flux, artificial viscosity and two shocks Riemann solver. Three benchmarking tests make the subject of unsteady free surface flow in this study. It is 1D typical dam-break on wet and dry bottom; 2D partial dam-break on a dry floodplain and 2D partial dam-break on channel with 90° bend. The analysis of the different errors between the numerical and analytical solutions and/or the experimental data shows that the SPH method gives reliable values with the selected optimal parameters which are the cubic kernel function and the artificial viscosity term. The increase in the number of particles increases the precision but also the calculation time.

Characterization of the Endwall Flow in a Low-Pressure Turbine Cascade Perturbed by Periodically Incoming Wakes, Part 2: Unsteady Blade Surface Measurements Using Pressure-Sensitive Paint

Unsteady pressure-sensitive paint (i-PSP) measurements were performed at a sampling rate of 30 kHz to investigate the near-endwall blade suction surface flow inside a low-pressure turbine cascade operating at engine-relevant high-speed and low-Re conditions. The investigation focuses on the interaction of periodically incoming bar wakes at 500 Hz with the secondary flow and the blade suction surface. The results build on extensive PIV measurements presented in the first part of this two-part publication, which captured the ’negative-jet-effect’ of the wakes throughout the blade passage. The surface pressure distributions are combined with CFD to analyze the flow topology, such as the passage vortex separation line. By analyzing data from phase-locked PIV and PSP measurements, a wake-induced moving pressure gradient negative in space and positive in time is found, which is intensified in the secondary flow region by 33% with respect to midspan. Furthermore, two methods of frequency-filtering based on FFT and SPOD are compared and utilized to associate a pressure fluctuation peak around 678 Hz with separation bubble oscillation.

Tropical cyclone low-level wind speed, shear, and veer: sensitivity to the boundary layer parametrization in the Weather Research and Forecasting model

Abstract. Mesoscale modeling can be used to analyze key parameters for wind turbine load assessment in a large variety of tropical cyclones. However, the modeled wind structure of tropical cyclones is known to be sensitive to the boundary layer scheme. We analyze modeled wind speed, shear, and wind veer across a wind turbine rotor plane in the eyewall and outer cyclone. We further assess the sensitivity of wind speed, shear, and veer to the boundary layer parametrization. Three model realizations of Typhoon Megi are analyzed over the open ocean using three frequently used boundary layer schemes in the Weather Research and Forecasting (WRF) model. All three typhoon simulations reasonably reproduce the cyclone track and structure. The boundary layer parametrization causes up to 15 % differences in median wind speed at hub height between the simulations. The simulated wind speed variability also depends on the boundary layer scheme. The modeled median wind shear is smaller than or equal to 0.11 used in the current IEC (International Electrotechnical Commission) standard regardless of the boundary layer scheme for the eyewall and outer cyclone region. However, up to 43.6 % of the simulated wind profiles in the eyewall region exceed 0.11. While the surface inflow angle is sensitive to the boundary layer scheme, wind veer in the lowest 400 m of the atmospheric boundary layer is less affected by the boundary layer scheme. Simulated median wind veer reaches values up to 1.7×10-2° m−1 (1.2×10-2° m−1) in the eyewall region (outer cyclone region) and is relatively small compared to moderate-wind-speed regimes. On average, simulated wind speed shear and wind veer are highest in the eyewall region. Yet strong spatial organization of wind shear and veer along the rainbands may increase wind turbine loads due to rapid changes in the wind profile at the turbine location.

Predicting the hydraulic response of critical transport infrastructures during extreme flood events

Understanding the effects of extreme floods on critical infrastructures such as bridges is paramount for ensuring safety and resilient design in the face of climate change and extreme events. This study develops robust computational and predictive modeling tools for assessing the impacts of extreme floods on the hydraulic response and structural resilience of bridges. A computational fluid dynamic (CFD) model utilizing RANS equations and k-ω Shear Stress Transport (SST) for simulating supercritical flows is adopted to compute hydrodynamic pressures and the water levels on bridge piers of cylindrical and rectangular shapes during a flood event. The CFD model is validated based on the case study data obtained from the Haj Omran Bridge, built on the Khorramabad River in Iran. The numerical simulations consider hydrological conditions and exclude geotechnical parameters and abutment damages. The numerical results are evaluated based on well-established design guidelines. Machine learning techniques, including Extreme Gradient Boosting (XGBoost), Random Forest (RF), and Support Vector Regression (SVR), optimized with Grid Search Cross-Validation (GSCV), are adopted to enhance the accuracy of hydrodynamic pressure forecasting at bridge piers. The XGBoost model exhibits superior performance (R2 = 0.908, RMSE = 0.0279, and E = 3.41%) compared to the RF and SVR models. All the estimated pressure data by XGBoost falls within ±6 percent error lines, highlighting the model's robustness for out-of-range hydrodynamic pressure prediction. Additionally, an optimized Long Short-Term Memory (LSTM) model is adopted to effectively predict free surface flow profiles (i.e. flood depth) over the bridge (R2 = 0.937 and RMSE = 0.083), demonstrating its potential for practical applications of flood depth predictions over bridge infrastructures. The proposed methodological framework outlined in this study can facilitate modeling the impacts of extreme floods on bridges, enabling robust climate resilience assessment of critical infrastructures.

Energetic variational approaches for inviscid multiphase flow systems with surface flow and tension

We consider the governing equations for the motion of the inviscid fluids in two moving domains and an evolving surface from an energetic point of view. We employ our energetic variational approaches to derive inviscid multiphase flow systems with surface flow and tension. More precisely, we calculate the variation of the flow maps to the action integral for our model to derive both surface flow and tension. We also study the conservation and energy laws of our multiphase flow systems. The key idea of deriving the pressure of the compressible fluid on the surface is to make use of the feature of the barotropic fluid, and the key idea of deriving the pressure of the incompressible fluid on the surface is to apply a generalized Helmholtz-Weyl decomposition on a closed surface. In Appendix, we introduce one of the candidates for the viscous terms of viscous multiphase flow with a tangential compressible surface flow.

Predicting Ili River streamflow change and identifying the major drivers with a novel hybrid model

Study regionIli River, the main river in the Lake Balkhash basin. Study focusThis study introduces a novel hybrid model by coupling the Soil and Water Assessment Tool (SWAT) hydrological model with machine learning algorithms. It aims to simulate and forecast the streamflow of the Ili River, clearly delineating the roles of meteorological factors and anthropogenic factors in the process of streamflow change. New hydrological insights for the regionDuring the period 1960–2020, the contribution of climate change to streamflow variation was 104.59%-113.07%, the contribution of land use ranged from −10.75% to −4.59%, and the impact of reservoir construction was −2.27%. The predictive outcomes of the hybrid model indicate that, under the SSP2–4.5 and SSP5–8.5 scenarios, the streamflow of the Ili River is projected to increase by 12.8% and 14.3% respectively in the future period (2021–2100), in comparison to the historical period (1960–2020). The warming and humidification from 2021 to 2100 will lead to changes in the streamflow components of the Ili River, with a decrease in the proportion of groundwater flow and an increase in the proportion of surface flow. The results of this study provide a reference for the rational utilization and scientific management of water resources in the Ili River Basin.

The Navier-Stokes system with temperature and salinity for free surface flows. Numerical scheme and validation

In this paper, we propose and study a numerical scheme for the Navier-Stokes-Fourier system derived and studied by the authors in [12]. This system models hydrostatic free surface flows with density variations depending on temperature or salinity. We show that the finite volume/finite element scheme presented – based on a layer averaged formulation of the model – is well-balanced with regards to the steady state of the lake at rest and preserves the nonnegativity of the water height. A maximum principle on the density is also proved as well as a discrete entropy inequality (when the thermal and viscous effects are neglected). Some numerical validations are finally shown with comparisons to 3D analytical solutions and experiments.

Experimental heat transfer analysis of a new turbulator in a concentric heat exchanger tube: A full factorial design approach

Abstract People have been looking for alternative energy sources because current sources may be depleted. Furthermore, it is critical to utilize available energy sources as effectively as possible. Turbulators are among the topics to consider when it comes to energy efficiency. In practice, turbulators are often used in exchangers to enhance heat transfer. Putting newly constructed turbulators in determined locations of the constant surface temperature heat exchanger, velocity, temperatures, pressure, and flow measurements have been performed in the inlet and outlet. In the case of using different numbers of turbulators, their placement in different locations, their arrangement at different distances, and their use at different Reynolds, the changes in pressure drop, Nusselt number, and friction factor were studied with the full factor experimental design. In addition, efficiency, exergy loss rate, and NTU efficiency parameters have been calculated, and the relationship between them has been presented graphically. When the test data were evaluated, it was seen that as the number of turbulators increased, the thermal efficiency, friction factor and pressure loss also increased. Using turbulators in different numbers, at different positions, at different distances, and with different Reynolds numbers, the effect ratio on the pressure loss, Nusselt number, and friction factor parameters was detected. In the analyses made, the most efficient parameters on heat transfer were determined, respectively, as Reynolds number (64.55%), position of turbulators (19.73%), distance between turbulators (6.09%), and number of turbulators (5.94%).

Capillarity Constructed Open Siphon for Sustainable Drainage.

Siphon is an effective method to transfer liquid from a higher to a lower level, which has many applications in hygienic design, clinical apparatus, and hydraulic engineering. Traditional operation requires energy to overcome gravity and establish flow in a closed system. Achieving sustainable high flux siphon drainage without energy input remains a challenge due to viscous dissipation. Here, an unexpected open siphon behavior on the South American pitcher plant Heliamphora minor consisting of trichomes covered pitcher and a wedge-shaped sheath is examined. Exploiting the concept of Digital Twin, a new biomimetic research method by transforming the biological sample to a virtual 3D model is proposed and unveiled that maintained connection of wicking on sub-millimeter long trichomes due to asymmetric pressure distribution and ascending in wedge sheath under unbalanced pressure forms continuous surface flow. Exploring this mechanism, a biomimetic siphon device achieving continuous high flux exposed to ambient air is constructed. Besides, particles floating on the meniscus in the outside wedge move under a curvature gradient as water ascends, which implies a biological nutrient capture method and new dust collection manner in the drainage system. Applying the underlying principle enhances the siphon efficiency of floor drains and has the potential for other liquid transfer device design improvements.

Developing plant functional groups to identify changes in functional composition and diversity in a dryland river experiencing artificially sustained flows

Land use and climate changes are driving significant shifts in the magnitude and persistence of dryland stream surface flows. The impact of these shifts on ecological functioning is largely unknown, particularly where streams have become wetter rather than drier. This study investigated relationships between hydrologic regime (including surface water persistence, differences in groundwater depth and altered flooding dynamics) with plant traits and riverine vegetation functional composition. Our study system was a previously ephemeral creek in semi-arid northwest Australia that has received groundwater discharge from nearby mining operations for >15 years; surface flows are now persistent for ~27 km downstream of the discharge point. We aimed to (i) identify plant functional groups (FGs) associated with the creek and adjacent floodplain; and (ii) assess their distribution across hydrological gradients to predict shifts in ecological functioning in response to changing flow regimes. Seven FGs were identified using hierarchical clustering of 40 woody perennial plant species based on morphometric, phenological and physiologic traits. We then investigated how FG abundance (projective foliar cover), functional composition, and functional and taxonomic richness varied along a 14 km gradient from persistent to ephemeral flows, varying groundwater depths, and distances from the stream channel. Dominant FGs were (i) drought avoidant mesic trees that are fluvial stress tolerant, or (ii) drought tolerant xeric tall shrubs that are fluvial stress intolerant. The drought avoidant mesic tree FG was associated with shallow groundwater but exhibited lower cover in riparian areas closer to the discharge (persistent surface flows). However, there were more FGs and higher species richness closer to the discharge point, particularly on the floodplain. Our findings demonstrate that quantifying FG distribution and diversity is a significant step in both assessing the impacts of mine water discharge on riverine ecosystems and for planning for post-mining restoration.

Solute flow and particle transport in aquatic ecosystems: A review on the effect of emergent and rigid vegetation

In-channel vegetation is ubiquitous in aquatic environments and plays a critical role in the fate and transport of solutes and particles in aquatic ecosystems. Recent studies have advanced our understanding of the role of vegetation in solute flow and particle transport in aquatic ecosystems. This review summarizes these papers and discusses the impacts of emergent and rigid vegetation on the surface flow, the advection and dispersion of solutes, suspended load transport, bedload transport, and hyporheic exchange. The two competing effects of emergent vegetation on the above transport processes are discussed. On the one hand, emergent vegetation reduces mean flow velocity at the same surface slope, which reduces mass transport. On the other hand, at the same mean flow velocity, vegetation generates turbulence, which enhances mass transport. Mechanistic understanding of these two competing effects and predictive equations derived from laboratory experiments are discussed. Predictive equations for the mean flow velocity and turbulent kinetic energy inside an emergent vegetation canopy are derived based on force and energy balance. The impacts of emergent vegetation on the advection-dispersion process, the suspended load and bedload transport, and the hyporheic exchange are summarized. The impacts of other vegetation-related factors, such as vegetation morphology, submergence, and flexibility, are briefly discussed. The role of vegetation in transporting other particles, such as micro- and macro-plastics, is also briefly discussed. Finally, suggestions for future research directions are proposed to advance the understanding of the dynamic interplays among natural vegetation, flow dynamics, and sedimentary processes.

Adaptation of Wind Drag Coefficient Parameterization: Improvement of Hydrodynamic Modeling by a Wave‐Dependent Cd in Large Shallow Lakes

AbstractWind is a critical driving force in hydrodynamic and water quality modeling of large shallow lakes, and is characterized by the wind drag coefficient Cd, representing the momentum transfer at the air‐water interface. Contemporary empirical formulae for Cd estimation were derived over oceans and some of which are solely wind velocity U10 dependent. These formulae were previously found to be inadequate in inland lake models often resulting in the water velocity underestimation. To address this problem, a physical scale experiment was designed, in which Cd was measured using a wind profile and eddy covariance methodology. A new wind‐induced wave‐dependent Cd parameterization was also established and validated in two lake studies. The driving force was modified by the wave‐dependent Cd formula in a hydrodynamic model of the shallow Upper Klamath Lake (UKL), OR, USA. The experimental Cd was negatively correlated to the wind velocity up until the critical U10 = 1.6 m s−1 which was 1.0~3.1 times previous empirical extrapolations at light winds. The variation partitioning results showed that wave parameters contributed to more than 30% of Cd variation combined with wind parameters. The modified wind stress field was spatially heterogeneous and the modeled water velocity was closer to the observations at two sites. Significant main circulation and outer bank circulation were modeled accompanied by higher surface vorticity, compared to the original UKL model. Overall, the wave‐dependent Cd formula provided an improvement of the surface flow field in the UKL model and will improve the management of the lake ecosystems.

SURFACE FLOW ANALYSIS AS AN EFFORTS FLOOD MITIGATION

Problems that often arise in the city of Cirebon are flooding caused by rain, land changes not supported by adequate infrastructure, narrowing of drainage channels and sedimentation of channels. This is the location where flooding often occurs in the Pemuda Street area and its surroundings. Therefore, it is necessary to carry out research on surface flow analysis as an effort to handle floods with the aim of inundation management strategies to reduce excessive rainwater runoff in drainage areas and channels. The research methods used were problem identification, literature study, data collection, analysis and design planning. The results can show that inundation is handled in two ways, first by changing the dimensions of the channel at the initial height H to 0.45m, it still experiences inundation, this is due to the difference in elevation which causes inundation in the channel. There was a change in the flood inundation area to 23,127 ha from the original 58,958 ha, meaning that the flood free area increased from 9,388 ha to 35,831 ha. The second way is by making 378 infiltration wells spread across 34 channels that experience flooding. With these infiltration wells, the flooding can be reduced to a minimum from an area of ​​23,127 ha to 0.040 ha. The success rate for flood management efforts reached 99.908%, with failure being 0.092%. Of the area of ​​68,346 ha, the flood-free area reached 68,306 ha and the remaining inundated area was 0.040 ha..

A volume-conservation particle shifting scheme for moving particle method simulating free-surface flow

AbstractIn this study, a novel particle shifting scheme for the moving particle method simulating free surface flow is developed. The overall method is based on the framework of least square moving particle semi-implicit (LSMPS) method, enabling accurate and stable treatment of wall boundary without configuration of dummy or virtual wall particles. To avoid volume expansion, a volume-conservation particle shifting (VCPS) model is developed. An additional term considering the variation of particle numerical density is incorporated into the VCPS model to avoid volume expansion. Several numerical simulations are calculated to validate the effectiveness of the VCPS. It is demonstrated that LSMPS incorporating with VCPS shows satisfactory accuracy and superior capability to conserve volume.

Confinement-induced drift in Marangoni-driven transport of surfactant: a Lagrangian perspective

Successive drops of coloured ink mixed with surfactant are deposited onto a thin film of water to create marbling patterns in the Japanese art technique of Suminagashi. To understand the physics behind this and other applications where surfactant transports adsorbed passive matter at gas–liquid interfaces, we investigate the Lagrangian trajectories of material particles on the surface of a thin film of a confined viscous liquid under Marangoni-driven spreading by an insoluble surfactant. We study a model problem in which several deposits of exogenous surfactant simultaneously spread on a bounded rectangular surface containing a pre-existing endogenous surfactant. We derive Eulerian and Lagrangian formulations of the equations governing the Marangoni-driven surface flow. Both descriptions show how confinement can induce drift and flow reversal during spreading. The Lagrangian formulation captures trajectories without the need to calculate surfactant concentrations; however, concentrations can still be inferred from the Jacobian of the map from initial to current particle position. We explore a link between thin-film surfactant dynamics and optimal transport theory to find the approximate equilibrium locations of material particles for any given initial condition by solving a Monge–Ampère equation. We find that as the endogenous surfactant concentration $\delta$ vanishes, the equilibrium shapes of deposits using the Monge–Ampère approximation approach polygons with corners curving in a self-similar manner over lengths scaling as $\delta ^{1/2}$ . We explore how Suminagashi patterns may be produced by using computationally efficient successive solutions of the Monge–Ampère equation.

Quantifying of Radioactive Elements in Soil, Water and Plant Samples using Laser Induced Breakdown Spectroscopy (LIBS) Technique

In this study, the laser-induced breakdown spectroscopy (LIBS) technique was applied to detect radioactive elements on surface soil, plant leaves (Tamarindus indica), and flow water collected from Bala, Mayo-Kebbi Quest, Chad. The survey process indicates that the samples contain radioactive elements U, Fr, Ce, Ac, Cm, Tb, and Pm. The concentrations of the detected elements in samples were determined using the calibration curve method, while the plasma temperatures (Texc) and electron density (Ne) of the detected radioactive elements were calculated from Boltzmann linear plots. Except for uranium, all elements in the soil sample had amounts below the acceptable limit. The concentration of uranium in soil samples ranges from 0.024236 to 0.23439 ppm, with a safety limit of 0.03 ppm. The average plasma temperature Texc for uranium (U) is 700, while the electron density Ne is 500 cm-3 .

Unfolding via Mesh Approximation using Surface Flows

AbstractManufacturing a 3D object by folding from a 2D material is typically done in four steps: 3D surface approximation, unfolding the surface into a plane, printing and cutting the outline of the unfolded shape, and refolding it to a 3D object. Usually, these steps are treated separately from each other. In this work we jointly address the first two pipeline steps by allowing the 3D representation to smoothly change while unfolding. This way, we increase the chances to overcome possible ununfoldability issues. To join the two pipeline steps, our work proposes and combines different surface flows with a Tabu Unfolder. We empirically investigate the effects that different surface flows have on the performance as well as on the quality of the unfoldings. Additionally, we demonstrate the ability to solve cases by approximation which comparable algorithms either have to segment or can not solve at all.