Hydrodynamic Fluid Research Articles

A GPU-based hydrodynamic simulator with boid interactions

We present a hydrodynamic simulation system using the GPU compute shaders of DirectX for simulating virtual agent behaviors and navigation inside a smoothed particle hydrodynamical (SPH) fluid environment with real-time water mesh surface reconstruction. The current SPH literature includes interactions between SPH and heterogeneous meshes but seldom involves interactions between SPH and virtual boid agents. The contribution of the system lies in the combination of the parallel smoothed particle hydrodynamics model with the distributed boid model of virtual agents to enable agents to interact with fluids. The agents based on the boid algorithm influence the motion of SPH fluid particles, and the forces from the SPH algorithm affect the movement of the boids. To enable realistic fluid rendering and simulation in a particle-based system, it is essential to construct a mesh from the particle attributes. Our system also contributes to the surface reconstruction aspect of the pipeline, in which we performed a set of experiments with the parallel marching cubes algorithm per frame for constructing the mesh from the fluid particles in a real-time compute and memory-intensive application, producing a wide range of triangle configurations. We also demonstrate that our system is versatile enough for reinforced robotic agents instead of boid agents to interact with the fluid environment for underwater navigation and remote control engineering purposes.

Flow around an Aircraft Model—Comparison between Hydrodynamic Tunnel Tests and Computational Fluid Dynamics Simulations

This paper presents a comparative analysis of the results obtained using the computational fluid dynamic method and a hydrodynamic tunnel. The aerodynamic properties of the General Dynamics F-111 Aardvark aircraft were subjected to both qualitative and quantitative assessment. In both cases, the same geometric model of the aircraft was used. For the hydrodynamic tunnel tests, the model was 3D printed, while for the CFD simulations, the geometry was imported into ANSYS Fluent 2021 R2 software. During quantitative studies, aerodynamic characteristics, including the coefficients of lift, drag, and pitching moment for the aircraft, were determined over a wide range of angles of attack. In turn, flow field images and pathline visualization of the vortex were obtained following qualitative studies. The comparative analysis showed that high consistency between the simulation and experimental results was achieved.

Two-dimensional electron gases as non-Newtonian fluids

Two-dimensional electron systems offer an appealing platform to explore long-lived excitations arising due to collinear carrier scattering enabled by phase-space constraints at the Fermi surface. Recently it was found that these effects can boost excitation lifetimes over the fundamental bound set by Landau’s Fermi-liquid theory by a factor as large as (TF/T)α with α≈2. Long-lived degrees of freedom possess the capability to amplify the response to weak perturbations, producing lasting collective memory effects. This leads to non-Newtonian hydrodynamics in 2D electron fluids driven by multiple viscous modes with scale-dependent viscosity. We describe these modes as Fermi surface modulations of odd parity evolving in space and time, and discuss their implications for experimental studies of electron hydrodynamics.

Smoothed particle hydrodynamics simulations of microstructure induced stress overshoot in structured fluids

Stress overshoot/undershoot is an important phenomenon in structured fluids undergoing dynamic and transient flow. To accurately capture the flow process, it is important to have a better understanding of and include in the numerical modeling the microstructure evolution that leads to the stress overshoot/undershoot phenomenon. We present a procedure for incorporating a microstructure model into a Lagrangian framework based on the smoothed particle hydrodynamics fluid solver. The numerical simulation is performed for a typical structured fluid under an applied strain rate history flow. Good agreement between the numerical results and the experimental data lends credence to and validates the proposed procedure for simulations of complex mixture flows. Additionally, the interaction between a flow of structured fluid and a circular cylinder placed in a channel is investigated. The viscous force is found to overshoot together with the applied gradient pressure and decrease over time as the fluid approaches the equilibrium state.

Pore-scale hydrodynamics of non-Newtonian power-law fluids across a partially blocked porous medium in a confined channel

Transport of non-Newtonian fluids in porous media is pervasive in many natural and industrial applications. However, capturing the rheological behaviors of fluids by direct experimental techniques is challenging at the pore-scale. This paper outlines the pore-scale hydrodynamic interactions of non-Newtonian power-law fluids across a partially blocked porous medium in the laminar flow regime by computational fluid dynamics. The porous medium consists of an array of uniformly arranged square pillars. We explore the complex interplay of power-law rheology and Reynolds number on the microscopic flow field at the pore-scale. We capture the momentum transfer at the permeable interface between the porous and non-porous regions through stream-wise and span-wise velocity components and average volumetric flow rates at each pore-throat. Our results unveil a significant augmentation in stream-wise momentum by shear-thinning behavior and a diminution in momentum by shear-thickening behavior of the fluid through the porous medium. Further, the flow-leakage at the top interface purely depends on the combined effects of Reynolds number and power-law index. The channel pressure drop between the windward and leeward faces of the porous medium increases with the power-law index at low Reynolds number, while it decreases at high Reynolds number. Moreover, we provide a simple numerical framework to comprehend how the power-law behavior of the fluid dynamically regulates the flow field at the pore-scale.

Impaired hemodynamic renal reserve response following recovery from established acute kidney injury and improvement by hydrodynamic isotonic fluid delivery.

Renal reserve capacity may be compromised following recovery from acute kidney injury (AKI) and could be used to identify impaired renal function in the face of restored glomerular filtration rate (GFR) or plasma creatinine. To investigate the loss of hemodynamic renal reserve responses following recovery in a model of AKI, rats were subjected to left unilateral renal ischemia-reperfusion (I/R) injury and contralateral nephrectomy and allowed to recover for 5 wk. Some rats were treated 24 h post-I/R by hydrodynamic isotonic fluid delivery (AKI-HIFD) of saline through the renal vein, previously shown to improve recovery and inflammation relative to control rats that received saline through the vena cava (AKI-VC). At 5 wk after surgery, plasma creatinine and GFR recovered to levels observed in uninephrectomized sham controls. Baseline renal blood flow (RBF) was not different between AKI or sham groups, but infusion of l-arginine (7.5 mg/kg/min) significantly increased RBF in sham controls, whereas the RBF response to l-arginine was significantly reduced in AKI-VC rats relative to sham rats (22.6 ± 2.2% vs. 13.8 ± 1.8%, P < 0.05). RBF responses were partially protected in AKI-HIFD rats relative to AKI-VC rats (17.0 ± 2.2%) and were not significantly different from sham rats. Capillary rarefaction observed in AKI-VC rats was significantly protected in AKI-HIFD rats. There was also a significant increase in T helper 17 cell infiltration and interstitial fibrosis in AKI-VC rats versus sham rats, which was not present in AKI-HIFD rats. These data suggest that recovery from AKI results in impaired hemodynamic reserve and that associated CKD progression may be mitigated by HIFD in the early post-AKI period.NEW & NOTEWORTHY Despite the apparent recovery of renal filtration function following acute kidney injury (AKI) in rats, the renal hemodynamic reserve response is significantly attenuated, suggesting that clinical evaluation of this parameter may provide information on the potential development of chronic kidney disease. Treatments such as hydrodynamic isotonic fluid delivery, or other treatments in the early post-AKI period, could minimize chronic inflammation or loss of microvessels with the potential to promote a more favorable outcome on long-term function.

Jeans instability of predictive model of baryonic and fermionic dark matter with the impact of dissipative force

This study examines the Jeans instability and growth rate of a self-gravitating system made up of baryonic particles and dark matter with a source of dissipation. While treating dark matter using a quantum damped-hydrodynamic technique, we employed the classical hydrodynamic fluid approach for baryonic particles. We get the coupled system's dispersion relation by evaluating the equations. We examine the dispersion relation and the Jeans condition in various limiting circumstances and demonstrate how parameters influence growth rate with the use of graphical instances. We discovered that, irrespective of the jeans wave number, the dissipation source damps the growth rate. The study highlights how the introduction of dissipation influences the growth rate and damping effect observed in the system. Overall, the findings of this study provide valuable insights into the dynamics of self-gravitating systems and their evolution over time.

Transformation of a fully nonlinear breather-like package of internal waves over a bottom step in a layered fluid

In this paper, we study the process of transformation of a localized wave packet over a bottom step in a three-layer fluid, in which the height of the step is equal to or exceeds the thickness of the lower layer; therefore, density stratification becomes two-layer in the shallow water zone. In numerical experiments, both the height of the step and the width of the step were varied. The problem is solved in the framework of a fully nonlinear model of hydrodynamics of an inviscid incompressible stratified fluid. The primary analysis consisted in estimating the values of dimensionless parameters used, as a rule, in runup problems: the Froude and Iribarren numbers, the ratio of the characteristic wavelength to the characteristic slope width, the ratio of the topographic slope to the characteristic wave beam angle. Since the “cutoff” line for the lower pycnocline is partially or completely located on a step, one could expect the effects of run-up, breaking or reflection of waves propagating along the lower pycnocline, but this doesn’t happen. It is shown that the reflection of the wave packet from the step is minimal in all cases considered, a strong steepening of the wave is observed, but no breaking occurs in this case — the wave then just quickly decays on the lower pycnocline. An analysis of the spectral amplitudes and energy fields allows us to conclude that there is a transfer of energy from the lower pycnocline to the upper one. The breather in a two-layer fluid cannot exist, but the wave packet formed in the upper pycnocline after its destruction has much higher energy than it has before the step.

Radiation Absorption and Diffusion Thermo Effects on Unsteady MHD Kuvshinski Fluid Flow Past an Inclined Porous Plate in the Presence of Thermal Radiation and Chemical Reaction

The aim of the present paper is to investigate investigated the Radiation absorption and diffusion thermo on unsteady magneto hydrodynamic Kuvshinski fluid flow past an infinite vertical permeable moving plate in Presence of thermal radiation, heat absorption and homogenous chemical reaction, subjected to variable suction. The plate is assumed to be embedded in a uniform porous medium and Moves with a constant velocity in the flow direction in the presence of a transverse magnetic Field. The equations governing the flow are transformed into a system of nonlinear ordinary differential equations by using perturbation technique. Graphical results for the velocity distribution, temperature distribution and concentration distribution based on the numerical solutions are presented and discussed. Also discuss the effects of various parameters on the skin-friction coefficient and the rate of heat transfer in the form of Nusselt number and rate of mass transfer in the form of Sherwood number at the surface. Velocity and temperature distribution is observed to increase with an increase in Radiation absorption and Dufour parameter, whereas diminishes with enhances of Magnetic field, heat absorption coefficient, radiation parameter and chemical reaction parameter in respective velocity, temperature and concentration distributions.

A calibration-free model of micropipette aspiration for measuring properties of protein condensates

There is growing evidence that biological condensates, which are also referred to as membraneless organelles, and liquid-liquid phase separation play critical roles regulating many important cellular processes. Understanding the roles these condensates play in biology is predicated on understanding the material properties of these complex substances. Recently, micropipette aspiration (MPA) has been proposed as a tool to assay the viscosity and surface tension of condensates. This tool allows the measurement of both material properties in one relatively simple experiment, in contrast to many other techniques that only provide one or a ratio of parameters. While this technique has been commonly used in the literature to determine the material properties of membrane-bound objects dating back decades, the model describing the dynamics of MPA for objects with an external membrane does not correctly capture the hydrodynamics of unbounded fluids, leading to a calibration parameter several orders of magnitude larger than predicted. In this work we derive a new model for MPA of biological condensates that does not require any calibration and is consistent with the hydrodynamics of the MPA geometry. We validate the predictions of this model by conducting MPA experiments on a standard silicone oil of known material properties and are able to predict the viscosity and surface tension using MPA. Finally, we reanalyze with this new model the MPA data presented in previous works for condensates formed from LAF-1 RGG domains.

Korteweg-Type Fluids and Thermodynamic Modelling via Higher-Order Gradients

This paper investigates the modelling of Korteweg-type fluids and hence the dependence of the stress tensor on gradients of mass density. This topic, originating from the need for describing capillarity effects, is mainly of interest in connection with nanosystems where the mean free path may be comparable with the geometric dimensions of the system. In addition to the Korteweg fluid model, the paper gives a review of the stress tensor function arising in quantum fluid hydrodynamics. Next, thermodynamic consistency is established for a fluid involving first- and second-order density gradients. The modelling investigated is a generalization of the classical Korteweg fluid and allows a better understanding of previous thermodynamic restrictions. The restrictions determined for the general scheme with second-order gradients are applied to the particular cases of the Korteweg fluid and the quantum fluid. Further, to allow for discontinuity wave solutions with finite speed of propagation, a model is established which involves higher-order derivatives and reduces to the Korteweg fluid in stationary conditions.

Hydrodynamics of shear thinning fluid in a square microchannel: a numerical approach

Abstract In this present work, a numerical study was conducted for the formation of a slug bubble for shear thinning non-Newtonian fluid in a cross-junction 2-D square horizontal microchannel. Carboxymethyl cellulose (CMC) of concentration 0.2 (w/w%) percent was used as a continuous phase that shows the shear thinning behavior of non-Newtonian fluid and Nitrogen (N2) was used as the discrete phase. The pressure-based double precision solver was used in ANSYS FLUENT 2021 R2 with the volume of fluid (VOF) method. The finite volume method is applied for the discretization of the continuity and momentum equation. This article also focuses on the fluctuation of static pressure, mechanism of slug, annular, and churn annular flow i.e., obtained by the variation in the inlet velocities. On the other hand, a concept that was applied in this work was also validated with the prior literature data.

Natural Convection Couette Flow in the Presence of Magnetic Field and Thermal Property

This article examines the natural convection Couette flow in a hydrodynamic viscous fluid that is electrically conductive due to the thermal radiation effect. The governing flow in this study is modelled in the form of partial differential equations (PDEs) in dimensional form with initial and boundary conditions and the Couette fluid model is also be used to characterize the fluid behavior. Then, using suitable non-dimensional quantities, the governing non-linear PDEs are transformed. Since the flow governing equations of the problem under study are extremely complex and complicated, techniques that complement experimental and theoretical fluid dynamics by providing alternative potentially cheaper means of testing fluid flow systems is employed. Therefore, the Finite Element Method (FEM) is used after discretization of the PDEs. With Graphs and tables, the effects of embedded thermo physical parameters of engineering interests associated with the flow quantities viz. velocity, temperature, concentration of the fluid were examined through series of numerical experiments and discussed. This research also analyzes and compares the results obtained by Omokhuale and Jabaka (2022b). It is interesting to report that an excellent agreement was established, thereby authenticating and validating the accuracy of FEM as a strong tool. According to the results of this study, the actions of thermal radiation on the thermal and momentum boundary layers for increasing values are significant, also, increasing the magnetic field parameter impedes the fluid movement due to the Lorentz force action.

A fully resolved smoothed particle hydrodynamics-discrete element method study of the rheology of suspensions: The role of inertia and grain shape

This paper presents a numerical study on suspensions of monodisperse non-Brownian grains in a Couette flow. The fully resolved coupled smoothed particle hydrodynamics and discrete element method is employed to model the motion of arbitrarily shaped grains in a viscous fluid. The numerical method is benchmarked against its capability in accurately handling grain–fluid hydrodynamics and inter-grain collisions. It is then used to simulate suspension flows of varying particle Reynolds and Bagnold numbers subjected to different shear rates, solid concentrations, and solid-to-fluid density ratios. A special focus is placed on the effect of grain shape with different aspect ratios and convexities on the flow behavior. Both the inertia and the grain shape are found to affect the grain–fluid and inter-grain interactions and uniquely contribute to the overall shear stress and the rheology of the suspension. The local profiles of solid concentration suggest the presence of grain layering near the boundary walls, which becomes more pronounced with higher solid concentration and inertia, and increased non-circularity in grain shape. A further examination of the pair distribution function and average particle rotation reveals a strong correlation between suspension viscosity and grain microstructure and kinematics.

Homotopy time-fractional enhancement soliton solution for six component complex cometary plasma with thermal heavy ions

A dusty plasma system composed of three inertial species and three Kappa-distributed inertial-less species is considered. The hydrodynamic fluid equation used to describe investigated plasma system. The plasma system under consideration is a very rich system, so in this paper the effective Debye length and effective plasma frequency are considered to normalize all physical quantities. The Kortweg de-Veris equation is derived using the reductive perturbation method. Two envelope solitons (bright and dark) are propagating, and their appearance depends mainly on the values thermal parameter of heavy ions (oxygen ions) as well as the mass ratio for dust plasma. The homotopy perturbed transform method (HPTM) is used to solve the reduced time-fractional Korteweg-de Vries (TFKdV) equation. The effects of fractional parameters on the solution solutions are examined. Our results show that the shape of the solitary pulse is highly controlled by the fractional order, as well as the plasma parameters in addition to that the HPTM is very useful method, where few terms approximately three terms are enough to get successful approximate solution for fractional type equations.

Fluid hydrodynamics and thermal transports in nanofluids pulsating heat pipes applied for building energy Exploitations: Experimental investigations and full numerical simulations

Fluid hydrodynamics and thermal transports in nanofluids pulsating heat pipes applied for building energy Exploitations: Experimental investigations and full numerical simulations

Development of a scale-up strategy for an aerated coaxial mixer containing a non-Newtonian fluid: A mass transfer approach

Coaxial mixers have been shown to be effective in enhancing the hydrodynamic stress and shear environment inside the aerated systems. However, the scale-up study of the aerated coaxial mixing reactors based on a constant mass transfer coefficient has never been reported in the literature. In this study, for the first time, a practical technique is suggested to evaluate the scalability of these systems in terms of a constant mass transfer coefficient. The effects of impeller speed, impeller type, aeration rate, and pumping direction on the mass transfer, power consumption, gas holdup profile, fluid hydrodynamics, and energy dissipation rate were explored for gas dispersion in non-Newtonian fluids inside coaxial mixers through tomography, dynamic gassing-in, and computational fluid dynamics. It was found that a practical approach to preserve the mass transfer coefficient of the large-scale coaxial mixer the same as its small-scale counterpart was to maintain the volumetric aeration rate per working fluid volume constant.

On the statistical theory of self-gravitating collisionless dark matter flow

Dark matter, if it exists, accounts for five times as much as the ordinary baryonic matter. Compared to hydrodynamic turbulence, the flow of dark matter might possess the widest presence in our universe. This paper presents a statistical theory for the flow of dark matter that is compared with N-body simulations. By contrast to hydrodynamics of normal fluids, dark matter flow is self-gravitating, long-range, and collisionless with a scale-dependent flow behavior. The peculiar velocity field is of constant divergence nature on small scale and irrotational on large scale. The statistical measures, i.e., correlation, structure, dispersion, and spectrum functions, are modeled on both small and large scales, respectively. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. Incompressible and constant divergence flow share the same kinematic relations for even order correlations. The limiting correlation of velocity ρL=1/2 on the smallest scale (r = 0) is a unique feature of collisionless flow (ρL=1 for incompressible flow). On large scale, transverse velocity correlation has an exponential form T2∝e−r/r2 with a constant comoving scale r2=21.3 Mpc/h that may be related to the horizon size at matter–radiation equality. All other correlation, structure, dispersion, and spectrum functions for velocity, density, and potential fields are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of S2l∝r1/4. All other statistical measures can be obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for scale r between 1 and 7 Mpc/h. Divergence is negatively correlated for r &amp;gt; 30 Mpc/h that leads to a negative density correlation.

Vortex shedding characteristics and hydrodynamic forces of stationary and elastically mounted side-by-side cylinders fitted with small diameter control rods

Vortex shedding behind bluff bodies can induce noise and vibration problems when it engages in a feedback cycle of oscillations with an acoustic or structural mode. Therefore, several passive and active flow control techniques have been developed to suppress vortex shedding behind a single bluff body. However, less attention was given to the case of multiple bluff bodies in cross-flow. Therefore, in this paper, a numerical investigation is performed to study the effect of adding small diameter control rods in the vicinity of two side-by-side cylinders on the vortex shedding characteristics, the hydrodynamic fluid forces, as well as the suppression of vortex-induced vibration. The simulations are conducted at a Reynolds number of Re = 200 and the side-by-side cylinders have a spacing ratio of T∗=2.5. The number of control rods and their angular positions around the main cylinders are varied. Fitting the side-by-side cylinders with control rods has resulted in three different flow regimes compared to the base case with no control rods. These flow regimes depend on the angular positions of the control rods and they are classified as; bistable flow, stable biased flow, and merged wake flow regimes. The stable biased flow regime has resulted in the maximum reduction of the fluctuating lift coefficient of the side-by-side cylinders, especially when the control rods are placed at θ=135° from the flow stagnation point. Increasing the number of control rods around the cylinders has further reduced the fluctuating lift coefficient and produced a more symmetric wake at the rod angles which produce the stable biased flow, and the merged wake flow regimes. The vortex-induced vibration of the side-by-side cylinders is investigated when the control rods are placed at the angular position that results in a stable biased flow regime. It is revealed that fitting the cylinders with control rods at an angular position of θ=135° has resulted in reducing the vibration amplitude over the whole range of the simulated reduced flow velocities, which makes this technique a viable one to control the vortex shedding and the vortex-induced vibration of side-by-side cylinders.

Subject-specific tribo-contact conditions in total knee replacements: a simulation framework across scales.

Fundamental knowledge about invivo kinematics and contact conditions at the articulating interfaces of total knee replacements are essential for predicting and optimizing their behavior and durability. However, the prevailing motions and contact stresses in total knee replacements cannot be precisely determined using conventional invivo measurement methods. Insilico modeling, in turn, allows for a prediction of the loads, velocities, deformations, stress, and lubrication conditions across the scales during gait. Within the scope of this paper, we therefore combine musculoskeletal modeling with tribo-contact modeling. In the first step, we compute contact forces and sliding velocities by means of inverse dynamics approach and force-dependent kinematic solver based upon experimental gait data, revealing contact forces during healthy/physiological gait of young subjects. In a second step, the derived data are employed as input data for an elastohydrodynamic model based upon the finite element method full-system approach taking into account elastic deformation, the synovial fluid's hydrodynamics as well as mixed lubrication to predict and discuss the subject-specific pressure and lubrication conditions.