Information geometric regularization of the barotropic Euler equation
Fuente:
arXiv
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2023
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866915824791453696 |
|---|---|
| author | Cao, Ruijia Schäfer, Florian |
| author_facet | Cao, Ruijia Schäfer, Florian |
| contents | Shock waves in gas dynamics feature jump discontinuities that hinder numerical simulations. Viscous regularizations are prone to excessive dissipation of fine-scale structures. In this work, we propose the first inviscid regularization of the multidimensional Euler equation based on ideas from semidefinite programming, information geometry, geometric hydrodynamics, and nonlinear elasticity. The Lagrangian flow maps of Euler solutions are a dynamical system on the manifold of diffeomorphisms. We observe that shock formation arises from the manifold's geodesic incompleteness. Our regularization embeds it into an ambient space equipped with the information geometry of the logarithmic barrier function. Thus, the diffeomorphism manifold inherits a geodesically complete geometry. The resulting regularized conservation law replaces shocks with smooth profiles without affecting oscillatory structures. One and two-dimensional numerical experiments show its practical potential to enable higher-order methods without explicit shock capturing. While we focus on the barotropic Euler equations for concreteness and simplicity of exposition, our regularization easily extends to more general Euler and Navier-Stokes-type equations. Our approach regularizes the Wasserstein geometry of the mass density with its information geometry. The former captures the natural trajectories of physical particles and the latter that of statistical estimators. Information geometric regularization accounts for the mass density's dual nature as a statistical/computational tool summarizing the motion of physical particles. Thus, our work is a starting point for information geometric mechanics that views solutions of continuum mechanical PDEs as parameters of statistical models for unresolved scales and uses their information geometry to evolve them in time. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2308_14127 |
| institution | arXiv |
| publishDate | 2023 |
| record_format | arxiv |
| spellingShingle | Information geometric regularization of the barotropic Euler equation Cao, Ruijia Schäfer, Florian Numerical Analysis Analysis of PDEs 35L65, 76L05, 65M25, 76J20, 58B20 Shock waves in gas dynamics feature jump discontinuities that hinder numerical simulations. Viscous regularizations are prone to excessive dissipation of fine-scale structures. In this work, we propose the first inviscid regularization of the multidimensional Euler equation based on ideas from semidefinite programming, information geometry, geometric hydrodynamics, and nonlinear elasticity. The Lagrangian flow maps of Euler solutions are a dynamical system on the manifold of diffeomorphisms. We observe that shock formation arises from the manifold's geodesic incompleteness. Our regularization embeds it into an ambient space equipped with the information geometry of the logarithmic barrier function. Thus, the diffeomorphism manifold inherits a geodesically complete geometry. The resulting regularized conservation law replaces shocks with smooth profiles without affecting oscillatory structures. One and two-dimensional numerical experiments show its practical potential to enable higher-order methods without explicit shock capturing. While we focus on the barotropic Euler equations for concreteness and simplicity of exposition, our regularization easily extends to more general Euler and Navier-Stokes-type equations. Our approach regularizes the Wasserstein geometry of the mass density with its information geometry. The former captures the natural trajectories of physical particles and the latter that of statistical estimators. Information geometric regularization accounts for the mass density's dual nature as a statistical/computational tool summarizing the motion of physical particles. Thus, our work is a starting point for information geometric mechanics that views solutions of continuum mechanical PDEs as parameters of statistical models for unresolved scales and uses their information geometry to evolve them in time. |
| title | Information geometric regularization of the barotropic Euler equation |
| topic | Numerical Analysis Analysis of PDEs 35L65, 76L05, 65M25, 76J20, 58B20 |
| url | https://arxiv.org/abs/2308.14127 |