skip to main content


Title: Weak solutions for a poro-elastic plate system
We consider a recent plate model obtained as a scaled limit of the three-dimensional Biot system of poro-elasticity. The result is a ‘2.5’-dimensional linear system that couples traditional Euler–Bernoulli plate dynamics to a pressure equation in three dimensions, where diffusion acts only transversely. We allow the permeability function to be time dependent, making the problem non-autonomous and disqualifying much of the standard abstract theory. Weak solutions are defined in the so-called quasi-static case, and the problem is framed abstractly as an implicit, degenerate evolution problem. Utilizing the theory for weak solutions for implicit evolution equations, we obtain existence of solutions. Uniqueness is obtained under additional hypotheses on the regularity of the permeability function. We address the inertial case in an appendix, by way of semigroup theory. The work here provides a baseline theory of weak solutions for the poro-elastic plate and exposits a variety of interesting related models and associated analytical investigations.  more » « less
Award ID(s):
1907620
NSF-PAR ID:
10285198
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Applicable Analysis
ISSN:
0003-6811
Page Range / eLocation ID:
1 to 20
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The dynamics of initially truncated and bent line solitons for the Kadomtsev–Petviashvili (KPII) equation modelling internal and surface gravity waves is analysed using modulation theory. In contrast to previous studies on obliquely interacting solitons that develop from acute incidence angles, this work focuses on initial value problems for the obtuse incidence of two or three partial line solitons, which propagate away from one another. Despite counterpropagation, significant residual soliton interactions are observed with novel physical consequences. The initial value problem for a truncated line soliton – describing the emergence of a quasi-one-dimensional soliton from a wide channel – is shown to be related to the interaction of oblique solitons. Analytical descriptions for the development of weak and strong interactions are obtained in terms of interacting simple wave solutions of modulation equations for the local soliton amplitude and slope. In the weak interaction case, the long-time evolution of truncated and large obtuse angle solitons exhibits a decaying, parabolic wave profile with temporally increasing focal length that asymptotes to a cylindrical Korteweg–de Vries soliton. In contrast, the strong interaction case of slightly obtuse interacting solitons evolves into a steady, one-dimensional line soliton with amplitude reduced by an amount proportional to the incidence slope. This strong interaction is identified with the ‘Mach expansion’ of a soliton with an expansive corner, contrasting with the well-known Mach reflection of a soliton with a compressive corner. Interestingly, the critical angles for Mach expansion and reflection are the same. Numerical simulations of the KPII equation quantitatively support the analytical findings. 
    more » « less
  2. Abstract

    Early research in aerodynamics and biological propulsion was dramatically advanced by the analytical solutions of Theodorsen, von Kármán, Wu and others. While these classical solutions apply only to isolated swimmers, the flow interactions between multiple swimmers are relevant to many practical applications, including the schooling and flocking of animal collectives. In this work, we derive a class of solutions that describe the hydrodynamic interactions between an arbitrary number of swimmers in a two‐dimensional inviscid fluid. Our approach is rooted in multiply‐connected complex analysis and exploits several recent results. Specifically, the transcendental (Schottky–Klein) prime function serves as the basic building block to construct the appropriate conformal maps and leading‐edge‐suction functions, which allows us to solve the modified Schwarz problem that arises. As such, our solutions generalize classical thin aerofoil theory, specifically Wu's waving‐plate analysis, to the case of multiple swimmers. For the case of a pair of interacting swimmers, we develop an efficient numerical implementation that allows rapid computations of the forces on each swimmer. We investigate flow‐mediated equilibria and find excellent agreement between our new solutions and previously reported experimental results. Our solutions recover and unify disparate results in the literature, thereby opening the door for future studies into the interactions between multiple swimmers.

     
    more » « less
  3. Abstract

    Efficient and accurate modeling of the coupled thermal‐hydraulic‐mechanical‐chemical (THMC) processes in various rock formations is indispensable for designing energy geo‐structures such as underground repositories for high‐level nuclear wastes. This work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user‐defined element (UEL) subroutine which is then combined with the general‐purpose FE software Abaqus Standard to solve initial‐boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic‐mechanical (HM) and thermal‐hydraulic‐mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed‐form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three‐dimensional case study is performed with reference to the in‐situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical‐concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

     
    more » « less
  4. The cell cytoskeleton is a dynamic assembly of semi-flexible filaments and motor proteins. The cytoskeleton mechanics is a determining factor in many cellular processes, including cell division, cell motility and migration, mechanotransduction and intracellular transport. Mechanical properties of the cell, which are determined partly by its cytoskeleton, are also used as biomarkers for disease diagnosis and cell sorting. Experimental studies suggest that in whole cell scale, the cell cytoskeleton and its permeating cytosol may be modelled as a two-phase poro-viscoelastic (PVE) material composed of a viscoelastic (VE) network permeated by a viscous cytosol. We present the first general solution to this two-phase system in spherical coordinates, where we assume that both the fluid and network phases are in their linear response regime. Specifically, we use generalized linear incompressible and compressible VE constitutive equations to describe the stress in the fluid and network phases, respectively. We assume a constant permeability that couples the fluid and network displacements. We use these general solutions to study the motion of a rigid sphere moving under a constant force inside a two-phase system, composed of a linear elastic network and a Newtonian fluid. It is shown that the network compressibility introduces a slow relaxation of the sphere and non-monotonic network displacements with time along the direction of the applied force. Our results can be applied to particle-tracking microrheology to differentiate between PVE and VE materials, and to measure the fluid permeability as well as VE properties of the fluid and the network phases. 
    more » « less
  5. Abstract

    Conduit models of volcanic eruptions simulate magma evolution through phase transitions and material changes during ascent. We present a time‐dependent one‐dimensional model of a chamber‐conduit system to examine the temporal evolution of dome‐forming eruptions. As magma ascends, volatiles exsolve and may escape vertically through the column or laterally through the conduit walls. Magma solidifies which increases viscosity, leading to a natural transition from viscous flow at depth to frictional sliding along the conduit walls near the surface, resulting in the extrusion of a semisolid plug. The model evaluates time‐ and depth‐dependent pressure, velocity, porosity, and relative amounts of exsolved water to carbon dioxide. Transient effects arise when magma outflux from the chamber appreciably decreases pressure over the magma ascent timescale. For low magma permeability, transient effects increase porosity and velocity relative to steady‐state solutions. For high magma permeability, efficient vertical and lateral gas escape depresses porosity and velocity at later times. We use the model to predict three time series data sets from the 2004–2008 eruption of Mount St. Helens: extruded volume, ground deformation, and carbon dioxide emissions. We quantify sensitivity of model predictions to input parameters using the distance‐based generalized sensitivity analysis. Chamber volatile content, volume, and excess pressure influence the amplitude of observables, while conduit radius, frictional rate dependence and magma permeability influence temporal evolution. High magma permeability can cause marked departures from exponentially decaying flux and may explain the unique temporal evolution of deformation observed at the only nearby continuous GPS station in operation at the eruption onset.

     
    more » « less