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$$ . The velocity equation in this system is the 3D Navier–Stokes equation with dissipation only in the$$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)$$ . In addition, explicit decay rates in$$H^2({\mathbb {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$$ 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
Higher-order error estimates for physics-informed neural networks approximating the primitive equations
Abstract Large-scale dynamics of the oceans and the atmosphere are governed by primitive equations (PEs). Due to the nonlinearity and nonlocality, the numerical study of the PEs is generally challenging. Neural networks have been shown to be a promising machine learning tool to tackle this challenge. In this work, we employ physics-informed neural networks (PINNs) to approximate the solutions to the PEs and study the error estimates. We first establish the higher-order regularity for the global solutions to the PEs with either full viscosity and diffusivity, or with only the horizontal ones. Such a result for the case with only the horizontal ones is new and required in the analysis under the PINNs framework. Then we prove the existence of two-layer tanh PINNs of which the corresponding training error can be arbitrarily small by taking the width of PINNs to be sufficiently wide, and the error between the true solution and its approximation can be arbitrarily small provided that the training error is small enough and the sample set is large enough. In particular, all the estimates area priori, and our analysis includes higher-order (in spatial Sobolev norm) error estimates. Numerical results on prototype systems are presented to further illustrate the advantage of using the$$H^s$$ norm during the training.
more »
« less
- PAR ID:
- 10433266
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- Partial Differential Equations and Applications
- Volume:
- 4
- Issue:
- 4
- ISSN:
- 2662-2963
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We study the sparsity of the solutions to systems of linear Diophantine equations with and without non-negativity constraints. The sparsity of a solution vector is the number of its nonzero entries, which is referred to as the$$\ell _0$$ -norm of the vector. Our main results are new improved bounds on the minimal$$\ell _0$$ -norm of solutions to systems$$A\varvec{x}=\varvec{b}$$ , where$$A\in \mathbb {Z}^{m\times n}$$ ,$${\varvec{b}}\in \mathbb {Z}^m$$ and$$\varvec{x}$$ is either a general integer vector (lattice case) or a non-negative integer vector (semigroup case). In certain cases, we give polynomial time algorithms for computing solutions with$$\ell _0$$ -norm satisfying the obtained bounds. We show that our bounds are tight. Our bounds can be seen as functions naturally generalizing the rank of a matrix over$$\mathbb {R}$$ , to other subdomains such as$$\mathbb {Z}$$ . We show that these new rank-like functions are all NP-hard to compute in general, but polynomial-time computable for fixed number of variables.more » « less
-
Abstract Energy-conserving Hermite methods for solving Maxwell’s equations in dielectric and dispersive media are described and analyzed. In three space dimensions, methods of order 2mto$$2m+2$$ require$$(m+1)^3$$ degrees-of-freedom per node for each field variable and can be explicitly marched in time with steps independent ofm. We prove the stability for time steps limited only by domain-of-dependence requirements along with error estimates in a special semi-norm associated with the interpolation process. Numerical experiments are presented which demonstrate that Hermite methods of very high order enable the efficient simulation of the electromagnetic wave propagation over thousands of wavelengths.more » « less
-
Abstract We study commutators with the Riesz transforms on the Heisenberg group$${\mathbb {H}}^{n}$$ . The Schatten norm of these commutators is characterized in terms of Besov norms of the symbol. This generalizes the classical Euclidean results of Peller, Janson–Wolff and Rochberg–Semmes. The method in proof bypasses the use of Fourier analysis, allowing us to address not just the Riesz transforms, but also the Cauchy–Szegő projection and second order Riesz transforms on$${\mathbb {H}}^{n}$$ among other settings.more » « less
-
Abstract We find examples of cohomogeneity one metrics on$$S^4$$ and$$\mathbb {C}P^2$$ with positive sectional curvature that lose this property when evolved via Ricci flow. These metrics are arbitrarily small perturbations of Grove–Ziller metrics with flat planes that become instantly negatively curved under Ricci flow.more » « less
An official website of the United States government
