skip to main content

This content will become publicly available on October 1, 2022

Title: General solutions of poroelastic equations with viscous stress
Mechanical properties of cellular structures, including the cell cytoskeleton, are increasingly used as biomarkers for disease diagnosis and fundamental studies in cell biology. Recent experiments suggest that the cell cytoskeleton and its permeating cytosol, can be described as a poroelastic (PE) material. Biot theory is the standard model used to describe PE materials. Yet, this theory does not account for the fluid viscous stress, which can lead to inaccurate predictions of the mechanics in the dilute filamentous network of the cytoskeleton. Here, we adopt a two-phase model that extends Biot theory by including the fluid viscous stresses in the fluid's momentum equation. We use generalized linear viscoelastic (VE) constitutive equations to describe the permeating fluid and the network stresses and assume a constant friction coefficient that couples the fluid and network displacement fields. As the first step in developing a computational framework for solving the resulting equations, we derive closed-form general solutions of the fluid and network displacement fields in spherical coordinates. To demonstrate the applicability of our results, we study the motion of a rigid sphere moving under a constant force inside a PE medium, composed of a linear elastic network and a Newtonian fluid. We find that the more » network compressibility introduces a slow relaxation of the sphere and a non-monotonic network displacements with time along the direction of the applied force. These novel features cannot be predicted if VE constitutive equation is used for the medium. We show that our results can be applied to particle-tracking microrheology to differentiate between PE and VE materials and to independently measure the permeability and VE properties of the fluid and the network phases. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Sponsoring Org:
National Science Foundation
More Like this
  1. This study explores thermal convection in suspensions of neutrally buoyant, non-colloidal suspensions confined between horizontal plates. A constitutive diffusion equation is used to model the dynamics of the particles suspended in a viscous fluid and it is coupled with the flow equations. We employ a simple model that was proposed by Metzger, Rahli & Yin ( J. Fluid Mech. , vol. 724, 2013, pp. 527–552) for the effective thermal diffusivity of suspensions. This model considers the effect of shear-induced diffusion and gives the thermal diffusivity increasing linearly with the thermal Péclet number ( Pe ) and the particle volume fractionmore »( ϕ ). Both linear stability analysis and numerical simulation based on the mathematical models are performed for various bulk particle volume fractions $({\phi _b})$ ranging from 0 to 0.3. The critical Rayleigh number $(R{a_c})$ grows gradually by increasing ${\phi _b}$ from the critical value $(R{a_c} = 1708)$ for a pure Newtonian fluid, while the critical wavenumber $({k_c})$ remains constant at 3.12. The transition from the conduction state of suspensions is subcritical, whereas it is supercritical for the convection in a pure Newtonian fluid $({\phi _b} = 0)$ . The heat transfer in moderately dense suspensions $({\phi _b} = 0.2\text{--}0.3)$ is significantly enhanced by convection rolls for small Rayleigh number ( Ra ) close to $R{a_c}$ . We also found a power-law increase of the Nusselt number ( Nu ) with Ra , namely, $Nu\sim R{a^b}$ for relatively large values of Ra where the scaling exponent b decreases with ${\phi _b}$ . Finally, it turns out that the shear-induced migration of particles can modify the heat transfer.« less
  2. Electrohydrodynamics of drops is a classic fluid mechanical problem where deformations and microscale flows are generated by application of an external electric field. In weak fields, electric stresses acting on the drop surface drive quadrupolar flows inside and outside and cause the drop to adopt a steady axisymmetric shape. This phenomenon is best explained by the leaky-dielectric model under the premise that a net surface charge is present at the interface while the bulk fluids are electroneutral. In the case of dielectric drops, increasing the electric field beyond a critical value can cause the drop to start rotating spontaneously andmore »assume a steady tilted shape. This symmetry-breaking phenomenon, called Quincke rotation, arises due to the action of the interfacial electric torque countering the viscous torque on the drop, giving rise to steady rotation in sufficiently strong fields. Here, we present a small-deformation theory for the electrohydrodynamics of dielectric drops for the complete Melcher–Taylor leaky-dielectric model in three dimensions. Our theory is valid in the limits of strong capillary forces and highly viscous drops and is able to capture the transition to Quincke rotation. A coupled set of nonlinear ordinary differential equations for the induced dipole moments and shape functions are derived whose solution matches well with experimental results in the appropriate small-deformation regime. Retention of both the straining and rotational components of the flow in the governing equation for charge transport enables us to perform a linear stability analysis and derive a criterion for the applied electric field strength that must be overcome for the onset of Quincke rotation of a viscous drop.« less
  3. The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. Themore »drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Potential applications and caveats to the usage of LD are included as a summary.« less
  4. Hydrogels are a class of soft, highly deformable materials formed by swelling a network of polymer chains in water. With mechanical properties that mimic biological materials, hydrogels are often proposed for load bearing biomedical or other applications in which their deformation and failure properties will be important. To study the failure of such materials a means for the measurement of deformation fields beyond simple uniaxial tension tests is required. As a non-contact, full-field deformation measurement method, Digital Image Correlation (DIC) is a good candidate for such studies. The application of DIC to hydrogels is studied here with the goal ofmore »establishing the accuracy of DIC when applied to hydrogels in the presence of large strains and large strain gradients. Experimental details such as how to form a durable speckle pattern on a material that is 90% water are discussed. DIC is used to measure the strain field in tension loaded samples containing a central hole, a circular edge notch and a sharp crack. Using a nonlinear, large deformation constitutive model, these experiments are modeled using the finite element method (FEM). Excellent agreement between FEM and DIC results for all three geometries shows that the DIC measurements are accurate up to strains of over 10, even in the presence of very high strain gradients near a crack tip. The method is then applied to verify a theoretical prediction that the deformation field in a cracked sample under relaxation loading, i.e. constant applied boundary displacement, is stationary in time even as the stress relaxes by a factor of three.« less
  5. The laminar boundary layer of a viscous incompressible fluid subject to a two-dimensional wall curvature is evaluated. It is well known that a curved surface induces streamwise pressure gradient as well as wall curvature driven pressure gradient. Under certain assumptions, a family of similarity solutions can be obtained under the influence of flow acceleration/deceleration, which is known as the Falkner-Skan similarity solutions. In this study, the effect of the wall normal pressure gradient is taken into consideration, and the freestream flow parameters are adjusted for flow over a curved surface. Present results are obtained by numerical solution of a generalizedmore »Falkner-Skan equation governing similar solutions for flows over curved surfaces. The Falkner-Skan equations are solved by an RK4 shooting algorithm. Additionally, the transport of a passive scalar is incorporated in the present analysis at different Prandtl numbers. The objective of this paper is to use the curvilinear or axisymmetric boundary layer and energy equations to assess the effect of Favorable, Adverse and Zero pressure gradient on the laminar momentum and thermal boundary layer development. Major conclusions are summarized as follows: (i) as the pressure gradient β increases from negative values (APG) towards positive (FPG) values, the displacement (Δ∗) and momentum (θ∗) thickness tend to decrease no matter the curvature type, and, (ii) the normalized wall shear stress (i.e., f′′) exhibits a linear decreasing behavior as the wall curvature switches from concave (negative) to convex (positive) at a constant pressure gradient.« less