An official website of the United States government
We develop a theory of fluid--structure interaction (FSI) between an oscillatory Newtonian fluid flow and a compliant conduit. We consider the canonical geometries of a 2D channel with a deformable top wall and an axisymmetric deformable tube. Focusing on the hydrodynamics, we employ a linear relationship between wall displacement and hydrodynamic pressure, which has been shown to be suitable for a leading-order-in-slenderness theory. The slenderness assumption also allows the use of lubrication theory, and the flow rate is related to the pressure gradient (and the tube/wall deformation) via the classical solutions for oscillatory flow in a channel and in a tube (attributed to Womersley). Then, by two-way coupling the oscillatory flow and the wall deformation via the continuity equation, a one-dimensional nonlinear partial differential equation (PDE) governing the instantaneous pressure distribution along the conduit is obtained, without \textit{a priori} assumptions on the magnitude of the oscillation frequency (\textit{i.e.}, at arbitrary Womersley number). We find that the cycle-averaged pressure (for harmonic pressure-controlled conditions) deviates from the expected steady pressure distribution, suggesting the presence of a streaming flow. An analytical perturbative solution for a weakly deformable conduit is obtained to rationalize how FSI induces such streaming. In the case of a compliant tube, the results obtained from the proposed reduced-order PDE and its perturbative solutions are validated against three-dimensional, two-way-coupled direct numerical simulations. We find good agreement between theory and simulations for a range of dimensionless parameters characterizing the oscillatory flow and the FSI, demonstrating the validity of the proposed theory of oscillatory flows in compliant conduits at arbitrary Womersley number.
more »
« less
Elasto-inertial rectification of oscillatory flow in an elastic tube
The interaction between deformable surfaces and oscillatory driving is known to produce complex secondary time-averaged flows due to inertial and elastic nonlinearities. Here, we revisit the problem of oscillatory flow in a cylindrical tube with a deformable wall, and analyse it under a long-wave theory for small deformations, but for arbitrary Womersley numbers. We find that the oscillatory pressure does not vary linearly along the length of a deformable channel, but instead decays exponentially with spatial oscillations. We show that this decay occurs over an elasto-visco-inertial length scale that depends on the material properties of the fluid and the elastic walls, the geometry of the system, and the frequency of the oscillatory flow, but is independent of the amplitude of deformation. Inertial and geometric nonlinearities associated with the elastic deformation of the channel drive a time-averaged secondary flow. We quantify the flow using numerical solutions of the perturbation theory, and gain insight by developing analytic approximations. The theory identifies a complex non-monotonic dependence of the time-averaged flux on the elastic compliance and inertia, including a reversal of the flow. Finally, we show that our analytic theory is in excellent quantitative agreement with the three-dimensional direct numerical simulations of Pandeet al.(Phys. Rev. Fluids, vol. 8, no. 12, 2023, 124102).
more »
« less
- PAR ID:
- 10549853
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 996
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present a new computational fluid dynamics approach for simulating two‐phase flow in hybrid systems containing solid‐free regions and deformable porous matrices. Our approach is based on the derivation of a unique set of volume‐averaged partial differential equations that asymptotically approach the Navier‐Stokes Volume‐of‐Fluid equations in solid‐free regions and multiphase Biot Theory in porous regions. The resulting equations extend our recently developed Darcy‐Brinkman‐Biot framework to multiphase flow. Through careful consideration of interfacial dynamics (relative permeability and capillary effects) and extensive benchmarking, we show that the resulting model accurately captures the strong two‐way coupling that is often exhibited between multiple fluids and deformable porous media. Thus, it can be used to represent flow‐induced material deformation (swelling, compression) and failure (cracking, fracturing). The model's open‐source numerical implementation,hybridBiotInterFoam, effectively marks the extension of computational fluid mechanics into modeling multiscale multiphase flow in deformable porous systems. The versatility of the solver is illustrated through applications related to material failure in poroelastic coastal barriers and surface deformation due to fluid injection in poro‐visco‐plastic systems.more » « less
-
In this work, full 3-D numerical simulations are performed to study the combined effects of elastic and inertial forces along the Y and Z-midline of the channel. Ultimately, simulation results are compared and matched with experimental fluorescent streak images of the focusing of particles under the same parametric conditions. We reported that shear-gradient (FSG), N2-induced secondary flow transversal drag (FSF), and elastic (FEL) lift are the main forces responsible for the focusing of particles in the elasto-inertial regime.more » « less
-
Driving oscillatory flow around an obstacle generates, due to inertial rectification, a steady ‘streaming’ flow that is useful in a host of microfluidic applications. While theory has focused largely on two-dimensional flows, streaming in many practical microfluidic devices is three-dimensional due to confinement. We develop a three-dimensional streaming theory around an obstacle in a microchannel with a Hele-Shaw-like geometry, where one dimension (depth) is much shorter than the other two dimensions. Utilizing inertial lubrication theory, we demonstrate that the time-averaged streaming flow has a three-dimensional structure. Notably, the flow reverses direction across the depth of the channel, which is a feature not observed in less confined streaming set-ups. This feature is confirmed by our experiments of streaming around a cylinder sandwiched in a microchannel. Our theory also predicts that the streaming velocity decays as the inverse cube of the distance from the cylinder, faster than that expected from previous two-dimensional approaches. We verify this velocity decay quantitatively using particle tracking measurements from experiments of streaming around cylinders with different aspect ratios at different driving frequencies.more » « less
-
This paper investigates flow and transport in a slender wavy-walled vertical channel subject to a prescribed oscillatory pressure difference between its ends. When the ratio of the stroke length of the pulsatile flow to the channel wavelength is small, the resulting flow velocity is known to include a slow steady-streaming component resulting from the effect of the convective acceleration. Our study considers the additional effect of gravitational forces in configurations with a non-uniform density distribution. Specific attention is given to the slowly evolving buoyancy-modulated flow emerging after the deposition of a finite amount of solute whose density is different from that of the fluid contained in the channel, a relevant problem in connection with drug dispersion in intrathecal drug delivery (ITDD) processes, involving the injection of the drug into the cerebrospinal fluid that fills the spinal canal. It is shown that when the Richardson number is of order unity, the relevant limit in ITDD applications, the resulting buoyancy-induced velocities are comparable to those of steady streaming. As a consequence, the slow time-averaged Lagrangian motion of the fluid, involving the sum of the Stokes drift and the time-averaged Eulerian velocity, is intimately coupled with the transport of the solute, resulting in a slowly evolving problem that can be treated with two-time-scale methods. The asymptotic development leads to a time-averaged, nonlinear integro-differential transport equation that describes the slow dispersion of the solute, thereby circumventing the need to describe the small concentration fluctuations associated with the fast oscillatory motion. The ideas presented here can find application in developing reduced models for future quantitative analyses of drug dispersion in the spinal canal.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2024.612
