skip to main content

Title: Revisiting steady viscous flow of a generalized Newtonian fluid through a slender elastic tube using shell theory
A flow vessel with an elastic wall can deform significantly due to viscous fluid flow within it, even at vanishing Reynolds number (no fluid inertia). Deformation leads to an enhancement of throughput due to the change in cross‐sectional area. The latter gives rise to a non‐constant pressure gradient in the flow‐wise direction and, hence, to a nonlinear flow rate–pressure drop relation (unlike the Hagen–Poiseuille law for a rigid tube). Many biofluids are non‐Newtonian, and are well approximated by generalized Newtonian (say, power‐law) rheological models. Consequently, we analyze the problem of steady low Reynolds number flow of a generalized Newtonian fluid through a slender elastic tube by coupling fluid lubrication theory to a structural problem posed in terms of Donnell shell theory. A perturbative approach (in the slenderness parameter) yields analytical solutions for both the flow and the deformation. Using matched asymptotics, we obtain a uniformly valid solution for the tube's radial displacement, which features both a boundary layer and a corner layer caused by localized bending near the clamped ends. In doing so, we obtain a “generalized Hagen–Poiseuille law” for soft microtubes. We benchmark the mathematical predictions against three‐dimensional two‐way coupled direct numerical simulations (DNS) of flow and deformation performed more » using the commercial computational engineering platform by ANSYS. The simulations show good agreement and establish the range of validity of the theory. Finally, we discuss the implications of the theory on the problem of the flow‐induced deformation of a blood vessel, which is featured in some textbooks. « less
Award ID(s):
Publication Date:
Journal Name:
ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Characterizing the elastic properties of soft materials through bulge testing relies on accurate measurement of deformation, which is experimentally challenging. To avoid measuring deformation, we propose a hydrodynamic bulge test for characterizing the material properties of thick, pre-stressed elastic sheets via their fluid–structure interaction with a steady viscous fluid flow. Specifically, the hydrodynamic bulge test relies on a pressure drop measurement across a rectangular microchannel with a deformable top wall. We develop a mathematical model using first-order shear deformation theory of plates with stretching and the lubrication approximation for the Newtonian fluid flow. Specifically, a relationship is derived between the imposed flowrate and the total pressure drop. Then, this relationship is inverted numerically to yield estimates of the Young’s modulus (given the Poisson ratio) if the pressure drop is measured (given the steady flowrate). Direct numerical simulations of two-way-coupled fluid–structure interaction are carried out in ansys to determine the cross-sectional membrane deformation and the hydrodynamic pressure distribution. Taking the simulations as “ground truth,” a hydrodynamic bulge test is performed using the simulation data to ascertain the accuracy and the validity of the proposed methodology for estimating material properties. An error propagation analysis is performed via Monte Carlo simulation tomore »characterize the susceptibility of the hydrodynamic bulge test estimates to noise. We find that, while a hydrodynamic bulge test is less accurate in characterizing material properties, it is less susceptible to noise, in the input (measured) variable, than a hydrostatic bulge test.« less
  2. Non-Newtonian fluids are characterized by complex rheological behaviour that affects the hydrodynamic features, such as the flow rate–pressure drop relation. While flow rate–pressure drop measurements of such fluids are common in the literature, a comparison of experimental data with theory is rare, even for shear-thinning fluids at low Reynolds number, presumably due to the lack of analytical expressions for the flow rate–pressure drop relation covering the entire range of pressures and flow rates. Such a comparison, however, is of fundamental importance as it may provide insight into the adequacy of the constitutive model that was used and the values of the rheological parameters. In this work, we present a theoretical approach to calculating the flow rate–pressure drop relation of shear-thinning fluids in long, narrow channels that can be used for comparison with experimental measurements. We utilize the Carreau constitutive model and provide a semi-analytical expression for the flow rate–pressure drop relation. In particular, we derive three asymptotic solutions for small, intermediate and large values of the dimensionless pressures or flow rates, which agree with distinct limits previously known and allow us to approximate analytically the entire flow rate–pressure drop curve. We compare our semi-analytical and asymptotic results with the experimentalmore »measurements of Pipe et al. ( Rheol. Acta , vol. 47, 2008, pp. 621–642) and find excellent agreement. Our results rationalize the change in the slope of the flow rate–pressure drop data, when reported in log–log coordinates, at high flow rates, which cannot be explained using a simple power-law model.« less
  3. A phenomenological model is proposed to describe the deformation and orientation dynamics of finite-sized bubbles in both quiescent and turbulent aqueous media. This model extends and generalizes a previous work that is limited to only the viscous deformation of neutrally buoyant droplets, conducted by Maffettone & Minale ( J. Non-Newtonian Fluid Mech. , vol. 78, 1998, pp. 227–241), into a high Reynolds number regime where the bubble deformation is dominated by flow inertia. By deliberately dividing flow inertia into contributions from the slip velocity and velocity gradients, a new formulation for bubble deformation is constructed and validated against two experiments designed to capture the deformation and orientation dynamics of bubbles simultaneously with two types of surrounding flows. The relative importance of each deformation mechanism is measured by its respective dimensionless coefficient, which can be isolated and evaluated independently through several experimental constraints without multi-variable fitting, and the results agree with the model predictions well. The acquired coefficients imply that bubbles reorient through body rotation as they rise in water at rest but through deformation along a different direction in turbulence. Finally, we provide suggestions on how to implement the proposed framework for characterizing the dynamics of deformable bubbles/drops in simulations.
  4. Flight vehicles that operate in the supersonic regime can be subject to adverse fluid-structure interactions due to their lightweight design. The presence of geometric obstructions such as control surfaces or fins can induce compression shocks that can interact with the boundary layer, leading to flow separation. The interaction of flow, compression shock and structural dynamics is very difficult to model and currently only poorly understood. This work investigates experimentally the interaction between a compliant panel in a Mach 2 flow under a ramp-induced shock-wave/boundary layer interaction (SWBLI). Brass panels of length 4.8" and width 2.5" and different thicknesses (h=0.020", 0.016", 0.012" and 0.010") are investigated. Tests are performed both with and without a compression ramp installed. This direct comparison allows characterization of the effect of the SWBLI on the system dynamics. High-speed stereoscopic digital image correlation (DIC) and fast-response pressure sensitive paint (PSP) are used to obtain simultaneous full field deformation and surface pressure of the panels. The results show that the shock induced by the 20compression ramp leads to separation of the turbulent boundary layer close to the ramp starting at about 80% of the panel length. This results in a region of large pressure fluctuations which primarily increasemore »the vibration amplitude of the second panel mode. Analysis of the fundamental mode, which contains most of the vibration energy of the panel, shows that the SWBLI does not lead to changes of this mode, neither in frequency, amplitude or mode shape. On the other hand, analysis of the shock foot motion shows that the shock primarily oscillates at the fundamental frequency of the panel. This means that while the shock and panel oscillate at the same frequency, it is not two-way coupling. The panel vibration dictates the motion of the shock, but the shock (or rather the SWBLI) does not modify the fundamental panel vibration beyond the forcing provided by the turbulent boundary layer. Full field surface pressure predictions are made using linearized potential flow theory, which relates the local slope of the panel to the surface pressure. Results are found to be in good agreement in the region of attached flow.« less
  5. For plane channel flow, thermal stratification resulting from a wall-normal temperature gradient together with an opposing gravitational field can lead to buoyancy-driven instability of three-dimensional waves. Moreover, viscosity-driven instability can lead to the amplification of two-dimensional Tollmien-Schlichting waves. Temporal stability simulations considering different combinations of Reynolds number and Rayleigh number were performed to investigate both the buoyancy and viscosity-driven instability of Rayleigh-Benard-Poiseuille flow. The investigated cases are either (1) stable, (2) unstable with respect to three-dimensional waves (buoyancy-driven instability), or (3) unstable with respect to two-dimensional waves (viscosity-driven instability). Two new and highly accurate computational fluid dynamics codes have been developed for solving the full and linearized unsteady compressible Navier-Stokes equations in Cartesian coordinates. The codes employ fifth-order-accurate upwind-biased compact finite differences for the convective terms and fourth-order-accurate compact finite differences for the viscous terms. For the case with buoyancy-driven instability, strong linear growth is observed for a broad range of spanwise wavenumbers and the wavelength of the spanwise mode with the strongest non-linear growth is gradually decreasing in time. For the case with viscosity-driven instability, the linear growth rates are lower and the first mode to experience non-linear growth is a higher harmonic with half the wavelength of themore »primary wave. The present results are consistent with the neutral curves from the linear stability theory analysis by Gage and Reid.« less