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.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.


Title: An inf-sup approach to $$C_0$$-semigroup generation for an interactive composite structure-Stokes PDE dynamics
Abstract In this work, we investigate the existence and uniqueness properties of a composite structure (multilayered)–fluid interaction PDE system which arises in multi-physics problems, and particularly in biofluidic applications related to the mammalian blood transportation process. The PDE system under consideration consists of the interactive coupling of 3D Stokes flow and 3D elastic dynamics which gives rise to an additional 2D elastic equation on the boundary interface between these 3D PDE systems. By means of a nonstandard mixed variational formulation, we show that the PDE system generates a$$C_0$$ C 0 -semigroup on the associated finite energy space of data. In this work, the presence of the pressure term in the 3D Stokes equation adds a great challenge to our analysis. To overcome this difficulty, we follow a methodology which is based on the necessarily non-Leray-based elimination of the associated pressure term, via appropriate nonlocal operators. Moreover, while we express the fluid solution variable via decoupling of the Stokes equation, we construct the elastic solution variables by solving a mixed variational formulation via a Babuska–Brezzi approach.  more » « less
Award ID(s):
2348312
PAR ID:
10533390
Author(s) / Creator(s):
Publisher / Repository:
Journal of Evolution Equations
Date Published:
Journal Name:
Journal of Evolution Equations
Volume:
24
Issue:
3
ISSN:
1424-3199
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Physical experiments and numerical simulations have observed a remarkable stabilizing phenomenon: a background magnetic field stabilizes and dampens electrically conducting fluids. This paper intends to establish this phenomenon as a mathematically rigorous fact on a magnetohydrodynamic (MHD) system with anisotropic dissipation in$$\mathbb R^3$$ R 3 . The velocity equation in this system is the 3D Navier–Stokes equation with dissipation only in the$$x_1$$ x 1 -direction, while the magnetic field obeys the induction equation with magnetic diffusion in two horizontal directions. We establish that any perturbation near the background magnetic field (0, 1, 0) is globally stable in the Sobolev setting$$H^3({\mathbb {R}}^3)$$ H 3 ( R 3 ) . In addition, explicit decay rates in$$H^2({\mathbb {R}}^3)$$ H 2 ( R 3 ) are also obtained. For when there is no presence of a magnetic field, the 3D anisotropic Navier–Stokes equation is not well understood and the small data global well-posedness in$$\mathbb R^3$$ R 3 remains an intriguing open problem. This paper reveals the mechanism of how the magnetic field generates enhanced dissipation and helps to stabilize the fluid. 
    more » « less
  2. Abstract This paper proposes and investigates the two-grid stabilized lowest equal-order finite element method for the time-independent dual-permeability-Stokes model with the Beavers-Joseph-Saffman-Jones interface conditions. This method is mainly based on the idea of combining the two-grid and the two local Gauss integrals for the dual-permeability-Stokes system. In this technique, we use a difference between a consistent mass matrix and an under-integrated mass matrix for the pressure variable of the dual-permeability-Stokes model using the lowest equal-order finite element quadruples. In the two-grid scheme, the global problem is solved using the standard$$ P_1-P_1-P_1-P_1 $$ P 1 - P 1 - P 1 - P 1 finite element approximations only on a coarse grid with grid sizeH. Then, a coarse grid solution is applied on a fine grid of sizehto decouple the interface terms and the mass exchange terms for solving the three independent subproblems such as the Stokes equations, microfracture equations, and the matrix equations on the fine grid. On the other hand, microfracture and matrix equations are decoupled through the mass exchange terms. The weak formulation is reported, and the optimal error estimate is derived for the two-grid schemes. Furthermore, the numerical results validate that the two-grid stabilized lowest equal-order finite element method is effective and has the same accuracy as the coupling scheme when we choose$$ h=H^2 $$ h = H 2
    more » « less
  3. Martin Blunt (Ed.)
    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
  4. Abstract In the (special) smoothing spline problem one considers a variational problem with a quadratic data fidelity penalty and Laplacian regularization. Higher order regularity can be obtained via replacing the Laplacian regulariser with a poly-Laplacian regulariser. The methodology is readily adapted to graphs and here we consider graph poly-Laplacian regularization in a fully supervised, non-parametric, noise corrupted, regression problem. In particular, given a dataset$$\{x_i\}_{i=1}^n$$ { x i } i = 1 n and a set of noisy labels$$\{y_i\}_{i=1}^n\subset \mathbb {R}$$ { y i } i = 1 n R we let$$u_n{:}\{x_i\}_{i=1}^n\rightarrow \mathbb {R}$$ u n : { x i } i = 1 n R be the minimizer of an energy which consists of a data fidelity term and an appropriately scaled graph poly-Laplacian term. When$$y_i = g(x_i)+\xi _i$$ y i = g ( x i ) + ξ i , for iid noise$$\xi _i$$ ξ i , and using the geometric random graph, we identify (with high probability) the rate of convergence of$$u_n$$ u n togin the large data limit$$n\rightarrow \infty $$ n . Furthermore, our rate is close to the known rate of convergence in the usual smoothing spline model. 
    more » « less
  5. A<sc>bstract</sc> A number of supersymmetric Jackiw-Teitelboim (JT) gravity theories are known to be described (in the Euclidean path integral formulation) by double-scaled random matrix models. Such matrix models can be characterized using a certain “string equation”. It was shown recently that in extended supergravity, when the number of BPS states scales as$$ {\textrm{e}}^{S_0} $$ e S 0 , whereS0is the extremal entropy, a special ansatz for the leading order solution of the string equation yields the supergravity spectrum. Somewhat miraculously, the construction showed that the functional form of the non-BPS (continuum) sector predicts the precise form of the BPS sector, showing the robustness of the supergravity/matrix-model correspondence. In this paper, we refine the analysis and show that the string equation, combined with some simple requirements on solutions, are powerful tools for constraining the spectrum of extended JT supergravity theories. We re-explore the cases of$$ \mathcal{N}=2 $$ N = 2 and (small)$$ \mathcal{N}=4 $$ N = 4 JT supergravity, and then explore the new cases of spectra from$$ \mathcal{N}=3 $$ N = 3 and large$$ \mathcal{N}=4 $$ N = 4 JT supergravity (recently derived by Heydeman, Shi, and Turiaci) showing that our approach also works naturally for (nearly) all the models. Based on this success, we conjecture that these new supergravity models also have matrix model descriptions. 
    more » « less