Abstract Imagebased computational fluid dynamics (CFD) has become a new capability for determining wall stresses of pulsatile flows. However, a computational platform that directly connects image information to pulsatile wall stresses is lacking. Prevailing methods rely on manual crafting of a hodgepodge of multidisciplinary software packages, which is usually laborious and errorprone. We present a new computational platform, to compute wall stresses in imagebased pulsatile flows using the volumetric lattice Boltzmann method (VLBM). The novelty includes: (1) a unique image processing to extract flow domain and local wall normality, (2) a seamless connection between image extraction and VLBM, (3) an enroute calculation of strainrate tensor, and (4) GPU acceleration (not included here). We first generalize the streaming operation in the VLBM and then conduct application studies to demonstrate its reliability and applicability. A benchmark study is for laminar and turbulent pulsatile flows in an imagebased pipe (Reynolds number: 10 to 5000). The computed pulsatile velocity and shear stress are in good agreements with Womersley's analytical solutions for laminar pulsatile flows and concurrent laboratory measurements for turbulent pulsatile flows. An application study is to quantify the pulsatile hemodynamics in imagebased human vertebral and carotid arteries including velocity vector, pressure, and wallshearmore »
This content will become publicly available on November 25, 2023
Incompressible active phases at an interface. Part 1. Formulation and axisymmetric odd flows
Inspired by the recent realization of a twodimensional (2D) chiral fluid as an active monolayer droplet moving atop a 3D Stokesian fluid, we formulate mathematically its freeboundary dynamics. The surface droplet is described as a general 2D linear, incompressible and isotropic fluid, having a viscous shear stress, an active chiral driving stress and a Hall stress allowed by the lack of timereversal symmetry. The droplet interacts with itself through its driven internal mechanics and by driving flows in the underlying 3D Stokes phase. We pose the dynamics as the solution to a singular integral–differential equation, over the droplet surface, using the mapping from surface stress to surface velocity for the 3D Stokes equations. Specializing to the case of axisymmetric droplets, exact representations for the chiral surface flow are given in terms of solutions to a singular integral equation, solved using both analytical and numerical techniques. For a discshaped monolayer, we additionally employ a semianalytical solution that hinges on an orthogonal basis of Bessel functions and allows for efficient computation of the monolayer velocity field, which ranges from a nearly solidbody rotation to a unidirectional edge current, depending on the subphase depth and the Saffman–Delbrück length. Except in the nearwall limit, more »
 Award ID(s):
 2011854
 Publication Date:
 NSFPAR ID:
 10400882
 Journal Name:
 Journal of Fluid Mechanics
 Volume:
 951
 ISSN:
 00221120
 Sponsoring Org:
 National Science Foundation
More Like this


The relative velocities and positions of monodisperse highinertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al. , J. Fluid Mech. , vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian twotime correlation of fluid relative velocities seen by pairs that are nearly stationary. The twotime correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by largescale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integralscale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS.more »

We analyse the pressuredriven flow of the OldroydB fluid in slowly varying arbitrarily shaped, narrow channels and present a theoretical framework for calculating the relationship between the flow rate $q$ and pressure drop $\Delta p$ . We first identify the characteristic scales and dimensionless parameters governing the flow in the lubrication limit. Employing a perturbation expansion in powers of the Deborah number ( $De$ ), we provide analytical expressions for the velocity, stress and the $q$ – $\Delta p$ relation in the weakly viscoelastic limit up to $O(De^2)$ . Furthermore, we exploit the reciprocal theorem derived by Boyko $\&$ Stone ( Phys. Rev. Fluids , vol. 6, 2021, L081301) to obtain the $q$ – $\Delta p$ relation at the next order, $O(De^3)$ , using only the velocity and stress fields at the previous orders. We validate our analytical results with twodimensional numerical simulations in the case of a hyperbolic, symmetric contracting channel and find excellent agreement. While the velocity remains approximately Newtonian in the weakly viscoelastic limit (i.e. the theorem of Tanner and Pipkin), we reveal that the pressure drop strongly depends on the viscoelastic effects and decreases with $De$ . We elucidate the relative importance of different terms inmore »

Abstract This study investigates the dynamics of velocity shear and Reynolds stress in the ocean surface boundary layer for idealized misaligned wind and wave fields using a largeeddy simulation (LES) model based on the Craik–Leibovich equations, which captures Langmuir turbulence (LT). To focus on the role of LT, the LES experiments omit the Coriolis force, which obscures a stress–currentrelation analysis. Furthermore, a vertically uniform body force is imposed so that the volumeaveraged Eulerian flow does not accelerate but is steady. All simulations are first spunup without windwave misalignment to reach a fully developed stationary turbulent state. Then, a crosswind Stokes drift profile is abruptly imposed, which drives crosswind stresses and associated crosswind currents without generating volumeaveraged crosswind currents. The flow evolves to a new stationary state, in which the crosswind Reynolds stress vanishes while the crosswind Eulerian shear and Stokes drift shear are still present, yielding a misalignment between Reynolds stress and Lagrangian shear (sum of Eulerian current and Stokes drift). A Reynolds stress budgets analysis reveals a balance between stress production and velocity–pressure gradient terms (VPG) that encloses crosswind Eulerian shear, demonstrating a complex relation between shear and stress. In addition, the misalignment between Reynolds stress and Eulerian shearmore »

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinchoff happens. A onedimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a onedimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the selfsimilar nature of the solution. Numerical models for the onedimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a onedimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of onedimensional and twodimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for onedimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previousmore »