Flow Dynamics Research Articles

Investigating Pressure Gradient Dynamics in Two-phase Fluid Flow through Porous Media: An Experimental and Numerical Analysis

Investigating Pressure Gradient Dynamics in Two-phase Fluid Flow through Porous Media: An Experimental and Numerical Analysis

Wastewater discharges and urban land cover dominate urban hydrology signals across England and Wales

Abstract Urbanisation is an important driver of changes in streamflow. These changes are not uniform across catchments due to the diverse nature of water sources, storage, and pathways in urban river systems. While land cover data are typically used in urban hydrology analyses, other characteristics of urban systems (such as water management practices) are poorly quantified which means that urbanisation impacts on streamflow are often difficult to detect and quantify. Here, we assess urban impacts on streamflow dynamics for 711 catchments across England and Wales. We use the CAMELS-GB dataset, which is a large-sample hydrology dataset containing hydro-meteorological timeseries and catchment attributes characterising climate, geology, water management practices and land cover. We quantify urban impacts on a wide range of streamflow dynamics (flow magnitudes, variability, frequency, and duration) using random forest models. We demonstrate that wastewater discharges from sewage treatment plants and urban land cover dominate urban hydrology signals across England and Wales. Wastewater discharges increase low flows and reduce flashiness in urban catchments. In contrast, urban land cover increases flashiness and frequency of medium and high flow events. We highlight the need to move beyond land cover metrics and include other features of urban river systems in hydrological analyses to quantify current and future drivers of urban streamflow.

Porous media flow of polymeric nanofluids in multilayered formations: a simulation study

ABSTRACT The injection of viscoelastic fluids into the oil-bearing formations to push the bypassed and entrapped oil toward the producers is a widely accepted technique known for creating a stable displacement front and piston-like displacement, leading to Enhanced Oil Recovery (EOR). However, the sweep efficiency of such complex fluids in stratified formations is still a topic of debate in the literature and is problematic, particularly in harsh reservoir conditions. This research concerns how the inclusion of nanoparticles into the polymer solution could affect the porous media flow of polymeric fluids in multilayered reservoirs. Accordingly, validated core-scale numerical models for polymer flood and polymer nanohybrid flood were upscaled to develop a reservoir-scale model. Through numerical simulations using a compositional reservoir simulator, the simultaneous flow of crude oil and polymer nanohybrids in both low and high-permeable layers of the reservoir were scrutinized. The numerical outputs showed that, compared with polymer flood, the flow dynamics of the polymer nanohybrids in different layers of the reservoir are more favorable, which resulted in improved oil recovery of over 20% of Residual Oil In Place (ROIP). Additionally, our numerical findings indicated that the polymer nanohybrids could exhibit lower polymer adsorption and acceptable salt tolerance in porous media.

Nonlinear radiative thermal analysis of magneto-micropolar fluid over a bidirectionally extending sheet with varied thermal conditions

Thermal fluid flow analysis occasioned by stretching surfaces offers practical industrial and engineering applications, such as extrusion and drawing processes and materials experiencing heating and cooling conditions. Thus, this study explores magnetohydrodynamic micropolar fluid flow under varied thermal conditions and nonlinear thermal radiation over a dually stretched sheet subject to surface mass flux and an internal heat source. The study thus provides insight into the dynamics of heat transfer mechanisms for improving material processing techniques in many engineering processes. The flow dynamics with the thermal analysis are evaluated by applying two thermal conditions in the heat equation: the prescribed surface temperature (PST) and the prescribed heat flux (PHF). The formulated partial derivative equations that describe the current problem are changed into ordinary derivatives using some similarity quantities, while the converted equations are tackled with the combined techniques of shooting and Runge–Kutta Fehlberg. Consequently, diverse tables and figures are displayed to deliberate on the impacts of the physical terms on the dynamics of flow and heat transmission processes. The consequence of the analysis reveals a higher heat gradient in the prescribed surface temperature case than the uniform temperature distribution, as the Prandtl number magnifies. The micropolar material term reduces the surface frictional factor and depletes heat transmission, but the magnetic field term raises the skin friction coefficient.

Investigation and optimization into flow dynamics for an axial flow pump as turbine (PAT) with ultra-low water head

When the water downstream is abundant, axial flow pumps are often used as turbines (PAT) for reverse power generation to obtain additional benefits. However, due to the actual water level difference, PAT usually operates at off-design conditions. The CFD technology and entropy production theory were used to study the working performance, flow characteristics, and energy loss mechanism of the blade adjustable axial-flow PAT under a low water head. The results indicate that the strength of the tubular vortex at the pressure surface and the separation position of the tip-leakage-vortex at the suction surface depend on the inlet angle of the blade. The energy loss in the draft tube is influenced by the velocity circulation at the runner’s outlet, which depends on the blade’s outlet angle. Subsequently, the blade inlet and outlet angles in the runner were optimized. In turbine mode, the overall efficiency of the optimized runner increased by an average of 0.99 %. The blade pressure load distribution of the optimized runner has improved, and the pressure fluctuation in the draft tube is also under control. In pump mode, due to the increasing of the inlet angle within the same demand pump head range, the optimized runner’s pumping flow increases by an average of 2 m3/s.

Steady-state relative permeability measurements in rough-walled fractures: The effects of wettability and aperture

Tight subsurface geosystems provide significant prospects for the recovery and storage of mass and energy. Because of the tightness of these systems, natural and engineered fractures are the main conveyance conduits transporting the fluids from the tight rock matrix to the wellbore or vice versa. During the fracturing process, most of the newly induced fractures are considered water-wet, however during injection/production, the wettability of the fractures gradually changes. These changes are mainly due to the interactions between the resident and injected fluids, or the exposure of the virgin fracture surfaces to the resident fluids during production. This alteration in wettability can affect the trapping and flow dynamics, e.g., relative permeability, of the fluids flowing in the fracture. Furthermore, fluid injection and production into and from these systems impact the fractures pore pressure and hence the effective stress they are exposed to, resulting in variations in the fracture aperture. In this study, we present a comprehensive macro-scale experimental investigation focusing on the impact of wettability and aperture on steady-state two-phase relative permeability in rough-walled fractures. Using the stead-state measurement technique, we meticulously characterized oil-brine relative permeabilities and residual trapping in water-wet, oil-wet, and mixed-wet rough-walled fractures induced in Eagle Ford shale rock samples. Additionally, we assessed the influence of fracture aperture size on these properties during both drainage and imbibition flow processes. Our results indicate significant phase interference in oil-brine flow in non-water-wet rough-walled fractures, which renders the commonly used x-curve and Corey models inadequate to represent the steady-state oil-brine relative permeabilities measured in our study, leading to substantial overestimations. Mixed-wet fractures exhibited reduced relative permeabilities compared to oil-wet fractures, which in turn demonstrated lower relative permeabilities than water-wet counterparts. The observed trends stem from evolving fluid configuration and transport dynamics within the fractures as surface wettability transitions. Furthermore, an increase in fracture aperture corresponded with an increase in relative permeability. This was attributed to a broader flow area, diminished capillarity, attenuated impact of fracture wall roughness and tortuosity effects, and a reduction in bypassing and trapping events, all synergistically contributing to the mitigation of phase interference. Similarly, an increase in fracture aperture led to a marked decrease in post-waterflooding endpoint oil saturation and post-oil flooding endpoint brine saturation. Finally, we provided generalized correlation models based on experimental results from fractures with varying conductivities and wettability, yielding a more accurate representation of fracture relative permeabilities compared to traditional models. The data and correlation models produced in this study improve our understanding and prediction of multiphase flow in subsurface fractures and offer insights into rapid production decline in unconventional shale reservoirs.

Introducing a Novel Innovative Technique for the Recording and Interpretation of Dynamic Coronary Angiography.

In the study of coronary artery disease (CAD), the mechanism of plaque formation and development is still an important subject for investigation. A limitation of current coronary angiography (CAG) is that it can only show static images of the narrowing of arterial channels without identifying the mechanism of the disease or predicting its progression or regression. To address this limitation, the CAG technique has been modified. The new approach emphasizes identifying and analyzing blood flow patterns, employing methodologies akin to those used by hydraulic engineers for fluid or gas movement through domestic or industrial pipes and pumps. With the new technique, various flow patterns and arterial phenomena-such as laminar, turbulent, antegrade, retrograde, and recirculating flow and potentially water hammer shock and vortex formation-are identified, recorded, and classified. These phenomena are then correlated with the presence of lesions at different locations within the coronary vasculature. The formation and growth of these lesions are explained from the perspective of fluid mechanics. As the pathophysiology of CAD and other cardiovascular conditions becomes clearer, new medical, surgical, and interventional treatments could be developed to reverse abnormal coronary flow dynamics and restore laminar flow, leading to improved clinical outcomes.

INTELLIGENT NUMERICAL METHOD FOR STUDYING MAXWELL WILLIAMSON NANOFLUID FLOW WITH ACTIVATION ENERGY

The use of artificial intelligence and its techniques has become increasingly widespread in recent times. It is being used to solve stiff non-linear equations. Additionally, nanofluids play a pivotal role in studying heat transfer. All of this was the motivation for doing this work. This work investigates a two-dimensional magnetohydrodynamic stretched flow (2D-MHDSF) of Maxwell Williamson nanofluid (MWNF) affected by bioconvection and activation energy numerically through Levenberg-Marquardt backpropagation method (LMBM)-based artificial neural network approach. The mathematical formulation for the problem was obtained through non-linear partial differential equations (PDEs). The leading PDEs were transmitted into non-linear ordinary differential equations by similarity transformation variables. The reference results for the 2D-MHDSF-MWNF model are produced by the Lobatto IIIA method through different scenarios of specific parameters for the flow velocity, fluid temperature, nanoparticle concentration, and motile density profiles. Using obtained results as a dataset to apply the testing, training, and validation steps of the suggested LMBM for the 2D-MHDSF-MWNF model. The mean squared error, analysis of regression, and error histograms are presented to prove the efficiency and precision of the proposed method. The numerical results of LMBM are displayed as a study of the effects of different physical factors on flow dynamics for 2D-MHDSF-MWNF.

Direct numerical simulation of the 7-7-7 film cooling at a range of compound angles

Film cooling is essential for protecting substrate materials from the hot gases in turbomachines. Recent experiments have explored the effects of compound angles for the 7-7-7 film cooling hole, a baseline shaped film cooling hole developed at the START lab. However, differences in facility conditions, manufacturing process uncertainties, uncertainties in flow conditions, and the absence of internal cooling hole measurements necessitate further investigation to fully understand the impact of compound angles and the associated flow physics. This study presents direct numerical simulations (DNS) for the baseline 7-7-7 shaped film cooling hole, examining a range of compound angles and flow conditions paralleling those in prior lab experiments. The data gleaned offers insights into the flow dynamics within and in the vicinity of the cooling hole, particularly in areas proximate to the wall. Our findings indicate that increasing the compound angle adversely affects adiabatic cooling efficiency. A detailed analysis reveals that this decline in efficiency is primarily due to hot gas intrusion into the cooling hole, the elevation of the cooling jet from the wall, and intensified secondary flows at higher compound angles. Additionally, the DNS data reveals a logarithmic scaling in the laterally averaged mean velocity profile and a self-similarity in temperature profiles beyond a certain distance from the cooling hole.

Experimental and numerical study of combined impinging jets at low nozzle-plate distance for confined and unconfined impingement plate

This study investigates the heat transfer and flow dynamics of combined impinging jets at low nozzle-plate distances. The combined jets comprise (i) circular and quadruple helical jets and (ii) circular and double helical jets. A conventional jet is generated from the only circular part of the combined jet. Experiments are conducted with a constant total flow rate (Qt = 85 LPM) and low nozzle-plate distances (H/D = 0.1, 0.2, 0.3, 0.4) for confined and unconfined jet cases. Flow rate ratios (Q*=0.176, 0.235, 0.352, 0.588, 0.823) are used for different swirl intensities. A numerical study was performed using ANSYS Fluent software. The results show that the numerical analysis is consistent with the experimental research. The combined impingement jet shows a potential to improve the low heat transfer occurring at the stagnation zone at low nozzle-plate distance in the swirl jet, leading to non-uniform heat transfer. The confined jet case at the low nozzle-plate distance reduces the heat transfer at the stagnation point and over the entire plate as regards the unconfined jet case. Increasing the H/D distance in the combined impinging jet improves the heat transfer in the confined jet case but reduces it in the unconfined case. Correlations are given for heat transfer prediction depending on H/D and Q* for unconfined and confined cases of the combined jet. Increasing the flow rate ratio and the nozzle-plate distance increases the heat transfer uniformity for the confined and unconfined jet cases. In contrast, the unconfined jet case is better regarding heat transfer uniformity. The confinement plate creates a negative pressure on the impingement plate in all conditions. In contrast, the unconfined jet case is precisely the opposite.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage.

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.

The characteristics of flows past a circular cylinder with heat transfer in the proximity to the wall at different gap ratios

In this study, two-dimensional laminar flows past a circular cylinder under isothermal and non-isothermal conditions over a flat-plate were investigated numerically at Re = 200. The gap ratio G / D ranges from 0.5 to 3.0, with G representing the distance between the lower surface of the cylinder and the wall. The effects of G / D on the flow dynamics were explored. For G / D = 0.5 , the vortex shedding of the lower cylinder surface is suppressed and there is a large recirculation region on the wall. For G / D = 1.0 , the recirculation region is small and the lower vortex begins to shed. For G / D ⩾ 1.5 , the flow exhibits similarities to that of an unbounded cylinder. The heat transfer characteristics was evaluated. When the cylinder temperature T c is higher than the air temperature T a , the mean drag coefficient C D ¯ is reduced and the root-mean-square lift coefficient C L ′ is increased. The temperature field is reconstructed using the proper orthogonal decomposition method, which is influenced by the large-scale vortical structures.

Numerical investigation of mesoscale multiphase mass transport mechanism in fibrous porous media

At present, the proton exchange membrane fuel cell (PEMFC) is one of the most promising new energy solutions. The structure and material properties of the gas diffusion layer (GDL) are important factors that restrict the efficient water management and structural stability of PEMFCs. The complex mesoscopic fibrous porous structure of the GDL results in complex multiphase flow dynamics problems that have highly nonlinear characteristics. It is difficult to describe the details of mesoscopic multiphase coupled transport dynamics and numerically solve problems in which there is a two-phase flow with a high-density ratio. To solve these problems, a mesoscopic multiphase coupled transport model based on a lattice Boltzmann method and volume-of-fluid (LBM-VOF) model is proposed. Moreover, we also discuss the mechanism of the interaction between the fibrous porous structure and internal multiphase flow. The results obtained illustrate that the proposed method can obtain dynamic details of the flow field for fibrous porous media. The surface wettability of fibres strongly influences both the distribution and stability of liquid water clusters within porous materials. Concurrently, the fibre diameter assumes a pivotal role in governing lateral water diffusion and exclusion efficiency. This work lays a foundation for efficient water and thermal management of PEMFCs.

Insights from Murine Studies on the Site Specificity of Atherosclerosis.

Atherosclerosis is an inflammatory reaction that develops at specific regions within the artery wall and at specific sites of the arterial tree over a varying time frame in response to a variety of risk factors. The mechanisms that account for the interaction of systemic factors and atherosclerosis-susceptible regions of the arterial tree to mediate this site-specific development of atherosclerosis are not clear. The dynamics of blood flow has a major influence on where in the arterial tree atherosclerosis develops, priming the site for interactions with atherosclerotic risk factors and inducing cellular and molecular participants in atherogenesis. But how this accounts for lesion development at various locations along the vascular tree across differing time frames still requires additional study. Currently, murine models are favored for the experimental study of atherogenesis and provide the most insight into the mechanisms that may contribute to the development of atherosclerosis. Based largely on these studies, in this review, we discuss the role of hemodynamic shear stress, SR-B1, and other factors that may contribute to the site-specific development of atherosclerosis.

Control of plug flow dynamics in microfluidic T-junction using pulsations of dispersed phase flow rate

Active control methods are commonly employed in droplet-based microfluidics. While passive methods of droplet generation face limitations, active methods, including flow rate pulsations provide an avenue for precise control. We investigate the influence of sinusoidal pulsations of dispersed phase flow rate on the plug flow dynamics in a T-junction microchannel. The study covers a range of fluid sets, including gas-liquid and liquid-liquid systems. In the absence of pulsations, distinct regimes of plug formation were identified, and plug length, velocity, and natural frequency of formation were determined. Introducing flow rate pulsations at frequencies proportional to the natural ones f = fpulse/fplug, we observed distinct plug length distributions, including double- or triple-mode, drop-on-demand mode, and multi-mode distribution. The influence of pulsations on plug length was particularly notable at low dimensionless frequencies (f < 1), especially for fluid combinations with higher viscosity ratios. The model was proposed to explain plug length distribution patterns under pulsatile conditions. Micro-PIV measurements provided insights into the velocity fields within plugs, revealing the stretching of streamlines without altering plug length. Our findings contribute to the understanding of active control methods in droplet-based microfluidics and suggest potential applications for precise control over microfluidic processes.

Dynamics and Chaos of Convective Fluid Flow

In this paper, a mathematical model of fluid flows in a convective thermal system is developed, and a five-dimensional dynamical system is developed for the investigation of the convective fluid dynamics. The analytical solutions of periodic motions to chaos of the convective fluid flows are developed for steady-state vortex flows, and the corresponding stability and bifurcations of periodic motions in the five-dimensional dynamical system are studied. The harmonic frequency-amplitude characteristics for periodic flows are obtained, which provide energy distribution in the parameter space. Analytical homoclinic orbits for the convective fluid flow systems are developed for the asymptotic convection through the infinite-many homoclinic orbits in the five-dimensional dynamical system. The dynamics of fluid flows in the convective thermal systems are revealed, and one can use such methodology to predict atmospheric and oceanic phenomena through thermal convections.

Dynamic simulation of unsteady Casson–Carreau fluid flow considering temperature-dependent variable-properties and thermo-solutal mixed convection over flat and slendering sheets

The present investigation delves into the intricate interplay between surface effects and fluid dynamics, focusing on the flow characteristics of Casson and Carreau fluids endowed with temperature-dependent thermophysical properties. Specifically, the study explores the ramifications of radiative, mixed convection, and magnetized processes on heat and mass transfer phenomena. A comprehensive mathematical model is developed to analyze the behavior of Casson–Carreau fluids in conjunction with microscopic particles and gyrotactic microorganisms over both flat and slender surfaces. Leveraging similarity transformations, the governing equations are transformed into dimensionless form, facilitating numerical solutions using the bvp5c solver in MATLAB. The ensuing numerical and graphical analyses elucidate key physical quantities, including local skin friction coefficients, gyrotactic microorganism density, Nusselt number, and Sherwood number, across a spectrum of physical parameters. Notably, the investigation underscores the profound influence of flat surfaces on flow patterns and highlights the heightened consumption of fluid particles with increasing variable diffusivity and thermal conductivity. Moreover, enhanced values of parameters such as variable motile microorganism diffusion coefficient, Brownian motion parameter, and temperature-dependent thermal conductivity parameter are found to augment the Nusselt number, indicative of intensified heat transfer rates. This research underscores the imperative of comprehending the intricate dynamics of fluid flow and heat transfer, offering actionable insights for enhancing engineering designs and addressing multifaceted challenges across domains encompassing chemical engineering, environmental science, and biomedical applications.

Time-Resolved Dynamic Optical Coherence Tomography for Retinal Blood Flow Analysis.

Optical coherence tomography (OCT) representations in clinical practice are static and do not allow for a dynamic visualization and quantification of blood flow. This study aims to present a method to analyze retinal blood flow dynamics using time-resolved structural OCT. We developed novel imaging protocols to acquire video-rate time-resolved OCT B-scans (1024 × 496 pixels, 10 degrees field of view) at four different sensor integration times (integration time of 44.8 µs at a nominal A-scan rate of 20 kHz, 22.4 µs at 40 kHz, 11.2 µs at 85 kHz, and 7.24 µs at 125 kHz). The vessel centers were manually annotated for each B-scan and surrounding subvolumes were extracted. We used a velocity model based on signal-to-noise ratio (SNR) drops due to fringe washout to calculate blood flow velocity profiles in vessels within five optic disc diameters of the optic disc rim. Time-resolved dynamic structural OCT revealed pulsatile SNR changes in the analyzed vessels and allowed the calculation of potential blood flow velocities at all integration times. Fringe washout was stronger in acquisitions with longer integration times; however, the ratio of the average SNR to the peak SNR inside the vessel was similar across all integration times. We demonstrated the feasibility of estimating blood flow profiles based on fringe washout analysis, showing pulsatile dynamics in vessels close to the optic nerve head using structural OCT. Time-resolved dynamic OCT has the potential to uncover valuable blood flow information in clinical settings.

On the prediction of the turbulent flow behind cylinder arrays via Echo State Networks

Abstract This study aims at the prediction of the turbulent flow behind cylinder arrays by the application of Echo State Networks (ESN). Three different arrangements of arrays of seven cylinders are chosen for the current study. These represent different flow regimes: single bluff body flow, transient flow, and co-shedding flow. This allows the investigation of turbulent flows that fundamentally originate from wake flows yet exhibit highly diverse dynamics. The data is reduced by Proper Orthogonal Decomposition (POD) which is optimal in terms of kinetic energy. The Time Coefficients of the POD Modes (TCPM) are predicted by the ESN. The network architecture is optimized with respect to its three main hyperparameters, Input Scaling (INS), Spectral Radius (SR), and Leaking Rate (LR), in order to produce the best predictions in terms of Weighted Prediction Score (WPS), a metric leveling statistic and deterministic prediction. In general, the ESN is capable of imitating the complex dynamics of turbulent flows even for longer periods of several vortex shedding cycles. Furthermore, the mutual interdependencies of the TCPM are well preserved. However, optimal hyperparameters depend strongly on the flow characteristics. Generally, as flow dynamics become faster and more intermittent, larger LR and INS values result in better predictions, whereas less clear trends for SR are observable.

Spatiotemporal dynamics of cortical somatosensory network in typically developing children.

Sense of touch is essential for our interactions with external objects and fine control of hand actions. Despite extensive research on human somatosensory processing, it is still elusive how involved brain regions interact as a dynamic network in processing tactile information. Few studies probed temporal dynamics of somatosensory information flow and reported inconsistent results. Here, we examined cortical somatosensory processing through magnetic source imaging and cortico-cortical coupling dynamics. We recorded magnetoencephalography signals from typically developing children during unilateral pneumatic stimulation. Neural activities underlying somatosensory evoked fields were mapped with dynamic statistical parametric mapping, assessed with spatiotemporal activation analysis, and modeled by Granger causality. Unilateral pneumatic stimulation evoked prominent and consistent activations in the contralateral primary and secondary somatosensory areas but weaker and less consistent activations in the ipsilateral primary and secondary somatosensory areas. Activations in the contralateral primary motor cortex and supramarginal gyrus were also consistently observed. Spatiotemporal activation and Granger causality analysis revealed initial serial information flow from contralateral primary to supramarginal gyrus, contralateral primary motor cortex, and contralateral secondary and later dynamic and parallel information flows between the consistently activated contralateral cortical areas. Our study reveals the spatiotemporal dynamics of cortical somatosensory processing in the normal developing brain.