Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Thursday, July 16 until 12:00 AM ET on Friday 17 due to maintenance. We apologize for the inconvenience.


Title: Code and data for "Theory and simulation of elastoinertial rectification of oscillatory flows in two-dimensional deformable rectangular channels"
Oscillatory flows between deformable boundaries frequently occur in various processes within engineered and biological systems. The fluid–structure interaction problem, and specifically the two-way pressure-deformation coupling that arises due to the flow in a compliant channel, is central to applications. A complete understanding of oscillatory flow in a compliant 2D channel (a canonical configuration) has not been achieved; specifically, developing the elastoinertial theory and benchmarking its predictions against detailed simulation data remains a knowledge gap. The goal of this work is to address these knowledge gaps.  more » « less
Award ID(s):
2245343
PAR ID:
10675534
Author(s) / Creator(s):
; ;
Publisher / Repository:
Purdue University Research Repository
Date Published:
Edition / Version:
1.0
Subject(s) / Keyword(s):
CFD Mechanical Engineering FOS: Mechanical engineering fluid-structure interactions microfluidics lubrication theory viscous fluid flow soft hydraulics elastic deformation mathematical modeling oscillatory flow Finite-element method
Format(s):
Medium: X Size: 24 files; 467 MB Other: text/plain; text/x-script.python; application/zip; image/png
Size(s):
24 files 467 MB
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. This paper presents a theoretical and experimental study of oscillatory flow through a finite compliant tube, aiming to investigate fluid–structure interactions in structures like blood vessels. The setup consists of a horizontally submerged polydimethylsiloxane tube mounted between two fixed ends, each instrumented with a pressure sensor. The device that provides the flow was custom-designed and manufactured by our team. It consists of two rigidly connected piston–cylinder assemblies that operate 180° out of phase, ensuring that the amount of fluid pushed into the tube from one end is identical to the amount pulled out from the other end, thereby making the flow purely oscillatory. Visualizing the flow inside a narrow tube with a constantly changing shape can be quite cumbersome, making corrections for light refraction distinctly challenging. Therefore, the tube's undulating shape is monitored using a high-speed camera. The instantaneous tube profile is obtained by processing the video offline in MATLAB. A theoretical model was developed to describe the tube's wall radial motion. The resulting equation was solved analytically for small deformations to obtain the predicted local deformation history. For radial deformations within ∼7%, the analytical solution agreed well with observation. Driven by this validation, a complete solution for the flow field is proposed within the small deformation limit. This study showcases a setup where a compliant tube's flow conditions are measured in tandem with its deformation, offering a unique avenue through which the fluid–structure interaction models in compliant tubes can be tested and refined. 
    more » « less
  3. We study the steady-flow-induced deformation between an incompressible non-Newtonian fluid and a three-dimensional deformable channel. Specifically, we provide a comprehensive experimental-theoretical framework for such flows of constant-viscosity viscoelastic (i.e., Boger) fluids, which allows us to quantify the influence of the fluid’s viscoelasticity on the flow rate–pressure drop (𝑞–Δ⁢𝑝) relation. For a flow-rate-controlled regime and weakly viscoelastic flow, we find excellent agreement between the theoretical prediction and experimental results, specifically providing an experimental demonstration of how both wall compliance and fluid viscoelasticity decrease the pressure drop. To definitively demonstrate the coupled effects of wall compliance and fluid viscoelasticity (change of cross-sectional area and curving of streamlines), we perform experiments in an equivalent rigid channel, showing that both the Boger fluid and a viscosity-matched Newtonian fluid have the same pressure drop. Thus, the experimentally verified theory for the 𝑞–Δ⁢𝑝 relation of weakly viscoelastic flows of Boger fluids in compliant channels can now be used for the design of microfluidics and soft hydraulic systems operating with complex fluids. 
    more » « less
  4. 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
  5. Experiments have shown that flow in compliant microchannels can become unstable at a much lower Reynolds number than the corresponding flow in a rigid conduit. Therefore, it has been suggested that the wall's elastic compliance can be exploited towards new modalities of microscale mixing. While previous studies mainly focused on the local instability induced by the fluid–structure interactions (FSIs) in the system, we derive a one-dimensional (1-D) model to study the FSI's effect on the global instability. The proposed 1-D FSI model is tailored to long, shallow rectangular microchannels with a deformable top wall, similar to the experiments. Going beyond the usual lubrication flows analysed in these geometries, we include finite fluid inertia and couple the reduced flow equations to a novel reduced 1-D wall deformation equation. Although a quantitative comparison with previous experiments is difficult, the behaviours of the proposed model show, qualitatively, agreement with the experimental observations, and capture several key effects. Specifically, we find the critical conditions under which the inflated base state of the 1-D FSI model is linearly unstable to infinitesimal perturbations. The critical Reynolds numbers predicted are in agreement with experimental observations. The unstable modes are highly oscillatory, with frequencies close to the natural frequency of the wall, suggesting that the observed instabilities are resonance phenomena. Furthermore, during the start-up from an undeformed initial state, self-sustained oscillations can be triggered by FSI. Our modelling framework can be applied to other microfluidic systems with similar geometric scale separation under different operating conditions. 
    more » « less