skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Pressure-driven flow across a hyperelastic porous membrane
We report an experimental investigation of pressure-driven flow of a viscous liquid across thin polydimethylsiloxane (PDMS) membranes. Our experiments revealed a nonlinear relation between the flow rate $$Q$$ and the applied pressure drop $$\unicode[STIX]{x0394}p$$ , in apparent disagreement with Darcy’s law, which dictates a linear relationship between flow rate, or average velocity, and pressure drop. These observations suggest that the effective permeability of the membrane decreases with pressure due to deformation of the nanochannels in the PDMS polymeric network. We propose a model that incorporates the effects of pressure-induced deformation of the hyperelastic porous membrane at three distinct scales: the membrane surface area, which increases with pressure, the membrane thickness, which decreases with pressure, and the structure of the porous material, which is deformed at the nanoscale. With this model, we are able to rationalize the deviation between Darcy’s law and the data. Our result represents a novel case in which macroscopic deformations can impact the microstructure and transport properties of soft materials.  more » « less
Award ID(s):
1661672
PAR ID:
10140511
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
871
ISSN:
0022-1120
Page Range / eLocation ID:
742 to 754
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 experimental 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. 
    more » « less
  2. Abstract When a water drop is placed onto a soft polymer network, a wetting ridge develops at the drop periphery. The height of this wetting ridge is typically governed by the drop surface tension balanced by elastic restoring forces of the polymer network. However, the situation is more complex when the network is swollen with fluid, because the fluid may separate from the network at the contact line. Here we study the fluid separation and network deformation at the contact line of a soft polydimethylsiloxane (PDMS) network, swollen with silicone oil. By controlling both the degrees of crosslinking and swelling, we find that more fluid separates from the network with increasing swelling. Above a certain swelling, network deformation decreases while fluid separation increases, demonstrating synergy between network deformation and fluid separation. When the PDMS network is swollen with a fluid having a negative spreading parameter, such as hexadecane, no fluid separation is observed. A simple balance of interfacial, elastic, and mixing energies can describe this fluid separation behavior. Our results reveal that a swelling fluid, commonly found in soft networks, plays a critical role in a wetting ridge. 
    more » « less
  3. Abstract 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 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. 
    more » « less
  4. Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O(nlog 𝑛) instead of the typical O(𝑛^3 ) rate in the number of samples 𝑛 used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients. 
    more » « less
  5. Abstract Using extensive numerical simulation of the Navier–Stokes equations, we study the transition from the Darcy’s law for slow flow of fluids through a disordered porous medium to the nonlinear flow regime in which the effect of inertia cannot be neglected. The porous medium is represented by two-dimensional slices of a three-dimensional image of a sandstone. We study the problem over wide ranges of porosity and the Reynolds number, as well as two types of boundary conditions, and compute essential features of fluid flow, namely, the strength of the vorticity, the effective permeability of the pore space, the frictional drag, and the relationship between the macroscopic pressure gradient$${\varvec{\nabla }}P$$ P and the fluid velocityv. The results indicate that when the Reynolds number Re is low enough that the Darcy’s law holds, the magnitude$$\omega _z$$ ω z of the vorticity is nearly zero. As Re increases, however, so also does$$\omega _z$$ ω z , and its rise from nearly zero begins at the same Re at which the Darcy’s law breaks down. We also show that a nonlinear relation between the macroscopic pressure gradient and the fluid velocityv, given by,$$-{\varvec{\nabla }}P=(\mu /K_e)\textbf{v}+\beta _n\rho |\textbf{v}|^2\textbf{v}$$ - P = ( μ / K e ) v + β n ρ | v | 2 v , provides accurate representation of the numerical data, where$$\mu$$ μ and$$\rho$$ ρ are the fluid’s viscosity and density,$$K_e$$ K e is the effective Darcy permeability in the linear regime, and$$\beta _n$$ β n is a generalized nonlinear resistance. Theoretical justification for the relation is presented, and its predictions are also compared with those of the Forchheimer’s equation. 
    more » « less