Simple Fluid Research Articles

No yield stress required: Stress-activated flow in simple yield-stress fluids

An elastoviscoplastic constitutive equation is proposed to describe both the elastic and rate-dependent plastic deformation behavior of Carbopol® dispersions, commonly used to study yield-stress fluids. The model, a variant of the nonlinear Maxwell model with stress-dependent relaxation time, eliminates the need for a separate Herschel–Bulkley yield stress. The stress dependence of the viscosity was determined experimentally by evaluating the steady-state flow stress at a constant applied shear rate and by measuring the steady-state creep rate at constant applied shear stress. Experimentally, the viscosity’s stress-dependence was confirmed to follow the Ree–Eyring model. Furthermore, it is shown that the Carbopol® dispersions used here obey time-stress superposition, indicating that all relaxation times experience the same stress dependence. This was demonstrated by building a compliance mastercurve using horizontal shifting on a logarithmic time axis of creep curves measured at different stress levels and by constructing mastercurves of the storage- and loss-modulus curves determined independently by orthogonal superposition measurements at different applied constant shear stresses. Overall, the key feature of the proposed constitutive equation is its incorporation of a nonlinear stress-activated change in relaxation time, which enables a smooth transition from elastic to viscous behavior during start-up flow experiments. This approach bypasses the need for a distinct Herschel–Bulkley yield stress as a separate material characteristic. Additionally, the model successfully replicates the observed steady-state flow stress in transient-flow scenarios and the steady-state flow rate in creep experiments, underlining its effectiveness in capturing the material’s dynamic response. Finally, the one-dimensional description is readily extended to a full three-dimensional finite-strain elastoviscoplastic constitutive equation.

Tuning Wetting-Dewetting Thermomechanical Energy for Hydrophobic Nanopores via Preferential Intrusion.

Heat and the work of compression/decompression are among the basic properties of thermodynamic systems. Being relevant to many industrial and natural processes, this thermomechanical energy is challenging to tune due to fundamental boundaries for simple fluids. Here via direct experimental and atomistic observations, we demonstrate, for fluids consisting of nanoporous material and a liquid, one can overcome these limitations and noticeably affect both thermal and mechanical energies of compression/decompression exploiting preferential intrusion of water from aqueous solutions into subnanometer pores. We hypothesize that this effect is due to the enthalpy of dilution manifesting itself as the aqueous solution concentrates upon the preferential intrusion of pure water into pores. We suggest this genuinely subnanoscale phenomenon can be potentially a strategy for controlling the thermomechanical energy of microporous liquids and tuning the wetting/dewetting heat of nanopores relevant to a variety of natural and technological processes spanning from biomedical applications to oil-extraction and renewable energy.

Enhanced thermal study in hybrid nanofluid flow in a channel motivated by graphene/Fe3O4 and Newtonian heating

The significance of advanced hybrid type nanofluids cannot be avoided because of their multiple heat transport applications. These are not only in applied thermal, mechanical, chemical engineering but also use as coolant for nuclear reactors, paint, medical sciences and manufacturing of various home products. The dynamical study of fluid characteristics inside a channel with slowly expanding/contracting walls is a fascinating research topic. Therefore, the key aims of this research are to formulate a hybrid nanoliquid model using graphene and Fe3O4 nanoparticles. Further, physical effects of practical interest like Joule heating, magnetic field are merged in the model. After that, numerical tool used to solve it and then demonstrated the results via graphical aid. The results shown that the when the nanoparticles amount taken from 0.01 % to 0.06 % then the density changes from 1.00042 to 1.00252 (nanofluid), 1.00067 to 1.00277 (hybrid) and 1.00025 to 1.0015 (nano), 1.00125 to 1.0025 (hybrid). In similar way, thermal conductance boosted from 1.00025 to 1.0015 (nano) and 1.00085 to 1.0021 (hybrid). For elongating walls (α=2.0,4.0,6.0,8.0), the velocity increased and examined maximum at the center. The absorber walls (S=0.1,0.2,0.3,0.4) highly controlled the fluid motion and the particles at η=−1.0 and η=1.0 accelerates inconsequentially. The high energy dissipation, thermal radiations (Rd=0.1,0.3,0.5,0.7) and Re1=0.1,0.2,0.3,0.4 augmented the temperature mechanism in nano, hybrid and simple fluids.

On De Gennes narrowing of fluids confined at the molecular scale in nanoporous materials.

Beyond well-documented confinement and surface effects arising from the large internal surface and severely confining porosity of nanoporous hosts, the transport of nanoconfined fluids remains puzzling in many aspects. With striking examples such as memory, i.e., non-viscous effects, intermittent dynamics, and surface barriers, the dynamics of fluids in nanoconfinement challenge classical formalisms (e.g., random walk, viscous/advective transport)-especially for molecular pore sizes. In this context, while molecular frameworks such as intermittent Brownian motion, free volume theory, and surface diffusion are available to describe the self-diffusion of a molecularly confined fluid, a microscopic theory for collective diffusion (i.e., permeability), which characterizes the flow induced by a thermodynamic gradient, is lacking. Here, to fill this knowledge gap, we invoke the concept of "De Gennes narrowing," which relates the wavevector-dependent collective diffusivity D0(q) to the fluid structure factor S(q). First, using molecular simulation for a simple yet representative fluid confined in a prototypical solid (zeolite), we unravel an essential coupling between the wavevector-dependent collective diffusivity and the structural ordering imposed on the fluid by the crystalline nanoporous host. Second, despite this complex interplay with marked Bragg peaks in the fluid structure, the fluid collective dynamics is shown to be accurately described through De Gennes narrowing. Moreover, in contrast to the bulk fluid, the departure from De Gennes narrowing for the confined fluid in the macroscopic limit remains small as the fluid/solid interactions in severe confinement screen collective effects and, hence, weaken the wavevector dependence of collective transport.

An analytic model for the flow induced in syringomyelia cavities

A simple two-dimensional fluid–structure interaction problem, involving viscous oscillatory flow in a channel separated by an elastic membrane from a fluid-filled slender cavity, is analysed to shed light on the flow dynamics pertaining to syringomyelia, a neurological disorder characterized by the appearance of a large tubular cavity (syrinx) within the spinal cord. The focus is on configurations in which the velocity induced in the cavity, representing the syrinx, is comparable to that found in the channel, representing the subarachnoid space surrounding the spinal cord, both flows being coupled through a linear elastic equation describing the membrane deformation. An asymptotic analysis for small stroke lengths leads to closed-form expressions for the leading-order oscillatory flow, and also for the stationary flow associated with the first-order corrections, the latter involving a steady distribution of transmembrane pressure. The magnitude of the induced flow is found to depend strongly on the frequency, with the result that for channel flow rates of non-sinusoidal waveform, as those found in the spinal canal, higher harmonics can dominate the sloshing motion in the cavity, in agreement with previous in vivo observations. Under some conditions, the cycle-averaged transmembrane pressure, also showing a marked dependence on the frequency, changes sign on increasing the cavity transverse dimension (i.e. orthogonal to the cord axis), underscoring the importance of cavity size in connection with the underlying hydrodynamics. The analytic results presented here can be instrumental in guiding future numerical investigations, needed to clarify the pathogenesis of syringomyelia cavities.

Whole-heart modelling with valves in a fluid–structure interaction framework

Computational modelling of whole-heart function is a useful tool to study heart mechanics and haemodynamics. Many existing heart models focus on electromechanical aspect without considering physiological valves and use simplified fluid models instead. In this study we develop a four-chamber heart model featuring realistic chamber geometry, detailed valve modelling, hyperelasticity with fibre architecture and fluid–structure interaction analysis. Our model is used to investigate heart behaviours with different modelling assumptions including restricted/free valve annular dynamics, and with/without heart-pericardium interactions. Our simulation results capture the interactions between valve leaflet and surrounding flow, typical left ventricular flow vortices, typical venous and transvalvular flow waveform, and physiological heart deformations such as atrioventricular plane movement. The improvement of ventricular filling and atrial emptying at early diastole is evident with free annulus. In addition, we find that the added pericardial forces on the heart have a predominant effect on atrial wall deformation especially during atrial contraction, and further help with the atrial filling process. Most importantly, the current study provides a framework for comprehensive multi-physics whole-heart modelling considering all heart valves and fluid–structure interactions.

Detachment forces during parallel-plate gap separation mediated by a simple yield-stress fluid.

In this work we have monitored the multiple stages of the normal traction force response of a yield-stress fluid confined between two circular parallel plates that are separated at constant velocity. At narrow initial gaps, the air-fluid interface suffers from the Saffman-Taylor instability, confirmed by visual inspection of fingering patterns imprinted on the fluid. At larger initial gaps, the fluid preserves the initially imposed circular symmetry of the confining plates, indicating the absence of instability. Due to the system characteristics and experimental environment, the multiple traction force contributions occurred in cascade, permitting us to isolate the adhesion responses associated with viscosity, capillarity, and yield stress. Employing a standard Herschel-Bulkley model, we assessed the scaling of the traction force in multiple regimes-specifically, evaluating the dependencies of the fingering to yield-stress transitions.

Regulating droplet impact dynamics of nanoparticle suspension: Phenomena, mechanisms, and implications

Droplet impact on solid substrates is a ubiquitous phenomenon in nature, agriculture, and industrial processes, playing a crucial role in numerous applications including self-cleaning, pesticide utilization, and inkjet printing. As a promising technique, adding nanoparticles into simple fluids to form nanofluids can effectively manipulate droplet impact behaviors. However, a comprehensive understanding of how nanoparticles modify the droplet impact dynamics, especially on the nanoscale, is still far from being fully explored. Hence, in this work, through the combined effort of molecular dynamics simulations and theoretical analysis, we elaborate on the influences of nanoparticles on droplet impact process. Using simple droplets as a control, we summarize four typical droplet impact modes and reveal how nanoparticles alter the impact behaviors of droplets, taking into account the key parameters including substrate wettability, impact velocity, volume fraction, and mass fraction of nanoparticles. We also demonstrate that with appropriate modifications, the theoretical/empirical models to predict the maximum contact diameter and the occurrence of breakup for simple droplets still hold to predict those of nanofluid droplets. Our findings and results enhance the understanding of the impact of nanoparticles on the droplet impact dynamics, with promising possibilities for various applications where regulating droplet impact behaviors is desired.

Influence of Fluid Rheology on Blood Flow Haemodynamics in Patient-Specific Arterial Networks of Varied Complexity – In-Silico Studies

Abstract Results obtained with computational fluid dynamics (CFD) rely on assumptions made during a pre-processing stage, including a mathematical description of a fluid rheology. Up to this date there is no clear answer to several aspects, mainly related to the question of whether and under what conditions blood can be simplified to a Newtonian fluid during CFD analyses. Different research groups present contradictory results, leaving the question unanswered. Therefore, the objective of this research was to perform steady-state and pulsatile blood flow simulations using eight different rheological models in geometries of varying complexity. A qualitative comparison of shear- and viscosity-related parameters showed no meaningful discrepancies, but a quantitative analysis revealed significant differences, especially in the magnitudes of wall shear stress (WSS) and its gradient (WSSG). We suggest that for the large arteries blood should be modelled as a non-Newtonian fluid, whereas for the cerebral vasculature the assumption of blood as a simple Newtonian fluid can be treated as a valid simplification.

Biomagnetic flow over a flat plate embedded by a magnetic dipole in a porous media

The present study used a simplified axisymmetric biomagnetic fluid dynamics and porous mediamodel which includes FHD (Ferrohydrodynamics), porosity and inertia effects saturated by magnetic dipole to study the influences of the leading parameters on various flow variables along a flat plate.The governing equations are simplified and solved by finite difference approach. We clarify how the ferromagnetic interaction parameter, B and porosity,ε assumptions contribute in the bio-background of the problem of interest. Moreover, from the results of the flow profiles, accelerating and decelerating phenomena are noticed for B, and ε interaction.

Thermodynamic properties of fluids with Mie n − m potentials and application to tune effective Mie potentials for simple real fluids

ABSTRACT Monte Carlo computer simulations have been performed to obtain the equation of state and internal energy of fluids with Mie n − m potentials, with m = 6 and n = 9, 12 and 15, for different supercritical temperatures and densities covering most of the fluid phase. These data have been used to test the performance of a third-order perturbation theory with the first three perturbation terms determined from computer simulation, giving rise to a very good agreement between theory and simulation. Next, the theory has been used to tune effective two-body Mie n − 6 potentials for the real fluids Ne, Ar, Kr, Xe and C H 4 at supercritical temperatures for wide range of densities and pressures. It has been shown that the theory is able to accurately fit the experimental data for the energy and pressure of these fluids within the temperature and density ranges considered. This theory might be useful to obtain the equation of state and the internal energy of real molecular fluids within the context of the SAFT-VR Mie theory.

Analytic expressions for correlations in coarse-grained simple fluids

Coarse-graining of fluids is challenging because fluid particles are unbound and diffuse long distances in time. One approach creates coarse-grain variables that group all particles within a region centered on specific points in space and accounts for the movement of particles among such regions. In our previous work, we showed that in many cases, potential interactions for such a scheme adopted a generalized quadratic form, whose parameters depend on means, variances, and correlation coefficients among the coarse-grain variables. In this work, we use statistical mechanics to derive analytic expressions for these parameters, using properties of the fluid, including pair distribution functions. These expressions are compared against simulation-derived values and shown to be in good agreement. This approach can be used to calculate a priori the potential for any homogeneous, simple fluid, without the need for fitting procedures or matching, thus increasing the ease of use of this coarse-grain scheme and creating a foundation for large-scale bottom-up simulations. Furthermore, these expressions provide a quantitative way of studying the boundary between discrete (atomic) and continuum models of fluids.

Central vs peripheral vascular factors determining risk of aortic dissection

Hypertension is the commonest cause of “central” aortic dissection, despite being an “peripheral” arteriole disease. We hypothesised that known regional variations in peripheral organ perfusion, due to variation in arteriole tone, via its effect on differential head and neck vessel flows can alter the risk and location of central aortic dissection. This might help to explain why some hypertensive patients dissect their aorta and others don't despite having similar aortic geometry. We utilised a simple 2D computational fluid dynamic model to test our hypothesis. We hypothesise (1) that variations in “peripheral” organ flow can alter the site of “central” aortic dissection, despite the blood pressure and cardiac output being identical, (2) variations in cardiac output at rest may also help to explain the risk of dissection but not the site of the initial intimal tear, and (3) atrial fibrillation, a known factor for aortic dissection, increases the risk via stroke volume variation. MRI scanning, via regional aortic branch flow quantification, may help to identify individuals who have a family history of dissection/aneurysm, who are normotensive with normal aortic anatomy that have a high potential risk for dissection due to calculated high wall velocity/shear stresses. This may allow selective beta blocker prescriptions in an asymptomatic patient population that is at a predicted high risk of aortic dissection – a hither to unidentified patient subpopulation.

Comparison of dimensionality reduction techniques for multi-variable spatiotemporal flow fields

In the field of fluid mechanics, it is a potential consensus that nonlinear dimensionality reduction (DR) techniques outperform linear methods. However, this conclusion has been obtained based on simple fluid phenomena and an incomplete evaluation system of dimensionality reduction algorithms. In this study, we use an improved evaluation system of DR methods to compare and evaluate the performance of four DR methods, including two linear techniques: Principal Component Analysis (PCA) and Independent Component Analysis (ICA), and two non-linear techniques: Isometric Mapping (ISOMAP) and Locally Linear Embedding (LLE). The four methods are applied to analyze a complex hydrodynamic flow field with cavitation by considering the joint features of multiple variables. Results show that LLE can capture redundant features that do not contribute to understanding the characteristics of flow fields, while ISOMAP is more suitable for handling datasets with multiple scales than LLE. PCA and ISOMAP can successfully capture the characteristics and evolution of flow fields. In addition, 3D supplementary information can assist ICA in improving the problem of identifying unstable flow field states.

Statistical field theory of mechanical stresses in Coulomb fluids: general covariant approach vs Noether’s theorem

In this paper, we introduce a statistical field theory that describes the macroscopic mechanical forces in inhomogeneous Coulomb fluids. Our approach employs the generalization of Noether’s first theorem for the case of a fluctuating order parameter to calculate the stress tensor for Coulomb fluids. This tensor encompasses the mean-field stress tensor and fluctuation corrections derived through the one-loop approximation. The correction for fluctuations includes a term that accounts for the thermal fluctuations of the local electrostatic potential and field in the vicinity of the mean-field configuration. This correlation stress tensor determines how electrostatic correlation affects local stresses in a nonuniform Coulomb fluid. We also use a previously formulated general covariant methodology (Brandyshev and Budkov 2023 J. Chem. Phys. 158 174114) in conjunction with a functional Legendre transformation method and derive within it the same total stress tensor. We would like to emphasize that our general approaches are applicable not only to Coulomb fluids but also to nonionic simple or complex fluids, for which the field-theoretic Hamiltonian is known as a function of the relevant scalar order parameters.

Simplified models of aerosol collision and deposition for disease transmission

Fluid-mechanics research has focused primarily on droplets/aerosols being expelled from infected individuals and transmission of well-mixed aerosols indoors. However, aerosol collisions with susceptible hosts earlier in the spread, as well as aerosol deposition in the nasal cavity, have been relatively overlooked. In this paper, two simple fluid models are presented to gain a better understanding of the collision and deposition between a human and aerosols. The first model is based on the impact of turbulent diffusion coefficients and air flow in a room on the collisions between aerosols and humans. Infection rates can be determined based on factors such as air circulation and geometry as an infection zone expands from an infected host. The second model clarifies how aerosols of different sizes adhere to different parts of the respiratory tract. Based on the inhalation rate and the nasal cavity shape, the critical particle size and the deposition location can be determined. Our study offers simple fluid models to understand the effects of geometric factors and air flows on the aerosol transmission and deposition.

Effect of amphiphilic polymers on phase separating binary mixtures: A DPD simulation study.

We present the phase separation dynamics of a binary (AB), simple fluid (SF), and amphiphilic polymer (AP) mixture using dissipative particle dynamics simulation at d = 3. We study the effect of different AP topologies, including block copolymers, ring block copolymers (RCP), and miktoarm star polymers, on the evolution morphologies, dynamic scaling functions, and length scale of the AB mixture. Our results demonstrate that the presence of APs leads to significantly different evolution morphologies in SF. However, the deviation from dynamical scaling is prominent, mainly for RCP. Typically, the characteristic length scale for SF follows the power law R(t) ∼ tϕ, where ϕ is the growth exponent. In the presence of high AP, we observe diffusive growth (ϕ → 1/3) at early times, followed by saturation in length scale (ϕ → 0) at late times. The extent of saturation varies with constraints imposed on the APs, such as topology, composition ratio, chain length, and stiffness. At lower composition ratios, the system exhibits inertial hydrodynamic growth (ϕ → 2/3) at asymptotic times without clearly exhibiting the viscous hydrodynamic regime (ϕ → 1) at earlier times in our simulations. Our results firmly establish the existence of hydrodynamic growth regimes in low surfactant-influenced phase separation kinetics of binary fluids and settle the related ambiguity in d = 3 systems.

Survival probabilities and first-passage distributions of self-propelled particles in spherical cavities.

A model of self-propelled motion in a closed compartment containing simple or complex fluids is formulated in this paper in terms of the dynamics of a point particle moving in a spherical cavity under the action of random thermal forces and exponentially correlated noise. The particle's time evolution is governed by a generalized Langevin equation (GLE) in which the memory function, connected to the thermal forces by a fluctuation-dissipation relation, is described by Jeffrey's model of viscoelasticity (which reduces to a model of ordinary viscous dynamics in a suitable limit). The GLE is transformed exactly to a Fokker-Planck equation that in spherical polar coordinates is in turn found to admit of an exact solution for the particle's probability density function under absorbing boundary conditions at the surface of the sphere. The solution is used to derive an expression (that is also exact) for the survival probability of the particle in the sphere, starting from its center, which is then used to calculate the distribution of the particle's first-passage times to the boundary. The behavior of these quantities is investigated as a function of the Péclet number and the persistence time of the athermal forces, providing insight into the effects of nonequilibrium fluctuations on confined particle motion in three dimensions.

Measuring many-body distribution functions in fluids using test-particle insertion.

We derive a hierarchy of equations, which allow a general n-body distribution function to be measured by test-particle insertion of between 1 and n particles. We apply it to measure the pair and three-body distribution functions in a simple fluid using snapshots from Monte Carlo simulations in the grand canonical ensemble. The resulting distribution functions obtained from insertion methods are compared with the conventional distance-histogram method: the insertion approach is shown to overcome the drawbacks of the histogram method, offering enhanced structural resolution and a more straightforward normalization. At high particle densities, the insertion method starts breaking down, which can be delayed by utilizing the underlying hierarchical structure of the insertion method. Our method will be especially useful in characterizing the structure of inhomogeneous fluids and investigating closure approximations in liquid state theory.

Significant influence of fluid viscoelasticity on flow dynamics past an oscillating cylinder

This study focuses on numerically investigating the impact of fluid viscoelasticity on the flow dynamics around a transversely forced oscillating cylinder operating in the laminar vortex shedding regime at a fixed Reynolds number of $Re = 100$ . Specifically, we explore how fluid viscoelasticity affects the boundary between the lock-in and no lock-in regions and the corresponding wake characteristics compared with a simple Newtonian fluid. Our findings reveal that fluid viscoelasticity enables the synchronization of the vortex street with the cylinder motion at lower oscillation frequencies than those required for a Newtonian fluid. Consequently, the lock-in region boundary for a viscoelastic fluid differs from that of a Newtonian fluid and expands in the non-dimensional cylinder oscillation amplitude and frequency parameter space. In the primary synchronization region, the wake of a Newtonian fluid exhibits ‘2S’ (two single vortices) and ‘P+S’ (a pair of vortices and a single vortex) shedding modes. In contrast, a ‘2P’ (two pairs of vortices) vortex mode is observed for a viscoelastic fluid within the same region. To gain a deeper understanding of the differences in the coherent flow structures and their associated frequencies between the two fluids, we employ the data-driven reduced-order modelling technique, known as the dynamic mode decomposition (DMD) technique. Utilizing this technique, we successfully extract and visualize the two competing fundamental frequencies (cylinder oscillation and natural vortex shedding frequencies) and their associated flow structures in the case of the no lock-in state, whereas only the dominant cylinder oscillation frequency and associated flow structure in the case of the lock-in state. Furthermore, we propose that the presence of excess strain resulting from the stretching of polymer molecules in viscoelastic fluids leads to a distinct difference in the wake structure compared with Newtonian fluids. This observation aligns with the findings obtained from the $Q$ -criterion and vorticity transport analysis of the wake.