An official website of the United States government
We propose an operator-splitting scheme to approximate scalar conservation equations with stiff source terms having multiple (at least two) stable equilibrium points. The scheme com- bines a (reaction-free) transport substep followed by a (transport-free) reaction substep. The transport substep is approximated using the forward Euler method with continuous finite elements and graph viscosity. The reaction substep is approximated using an exponential integrator. The crucial idea of the paper is to use a mesh-dependent cutoff of the reaction time-scale in the reaction substep. We establish a bound on the entropy residual motivating the design of the scheme. We show that the proposed scheme is invariant-domain preserv- ing under the same CFL restriction on the time step as in the nonreactive case. Numerical experiments in one and two space dimensions using linear, convex, and nonconvex fluxes with smooth and nonsmooth initial data in various regimes show that the proposed scheme is asymptotic preserving.
more »
« less
First-Order Greedy Invariant-Domain Preserving Approximation for Hyperbolic Problems: Scalar Conservation Laws, and p-System
The paper focuses on first-order invariant-domain preserving approximations of hyperbolic systems. We propose a new way to estimate the artificial viscosity that has to be added to make explicit, conservative, consistent numerical methods invariant-domain preserving and entropy inequality compliant. Instead of computing an upper bound on the maximum wave speed in Riemann problems, we estimate a minimum wave speed in the said Riemann problems such that the approximation satisfies predefined invariant-domain properties and predefined entropy inequalities. This technique eliminates non-essential fast waves from the construction of the artificial viscosity, while preserving pre-assigned invariant-domain properties and entropy inequalities.
more »
« less
- PAR ID:
- 10525706
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Journal of Scientific Computing
- Volume:
- 100
- Issue:
- 2
- ISSN:
- 0885-7474
- Subject(s) / Keyword(s):
- conservation equations hyperbolic systems invariant domains convex limiting finite element method
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Kolmogorov–Sinai entropy is an invariant of measure-preserving actions of the group of integers that is central to classification theory. There are two recently developed invariants, sofic entropy and Rokhlin entropy, that generalize classical entropy to actions of countable groups. These new theories have counterintuitive properties such as factor maps that increase entropy. This survey article focusses on examples, many of which have not appeared before, that highlight the differences and similarities with classical theory.more » « less
-
Abstract High order schemes are known to be unstable in the presence of shock discontinuities or under‐resolved solution features, and have traditionally required additional filtering, limiting, or artificial viscosity to avoid solution blow up. Entropy stable schemes address this instability by ensuring that physically relevant solutions satisfy a semi‐discrete entropy inequality independently of discretization parameters. However, additional measures must be taken to ensure that solutions satisfy physical constraints such as positivity. In this work, we present a high order entropy stable discontinuous Galerkin (ESDG) method for the nonlinear shallow water equations (SWE) on two‐dimensional (2D) triangular meshes which preserves the positivity of the water heights. The scheme combines a low order positivity preserving method with a high order entropy stable method using convex limiting. This method is entropy stable and well‐balanced for fitted meshes with continuous bathymetry profiles.more » « less
-
In this paper we study the $$\mathbb {C}^*$$ -fixed points in moduli spaces of Higgs bundles over a compact Riemann surface for a complex semisimple Lie group and its real forms. These fixed points are called Hodge bundles and correspond to complex variations of Hodge structure. We introduce a topological invariant for Hodge bundles that generalizes the Toledo invariant appearing for Hermitian Lie groups. An important result of this paper is a bound on this invariant which generalizes the Milnor–Wood inequality for a Hodge bundle in the Hermitian case, and is analogous to the Arakelov inequalities of classical variations of Hodge structure. When the generalized Toledo invariant is maximal, we establish rigidity results for the associated variations of Hodge structure which generalize known rigidity results for maximal Higgs bundles and their associated maximal representations in the Hermitian case.more » « less
-
This paper proposes an invariant-domain preserving approximation technique for nonlinear conservation systems that is high-order accurate in space and time. The algorithm mixes a high order finite element method with an invariant-domain preserving low-order method that uses the closest neighbor stencil. The construction of the flux of the low-order method is based on an idea from Abgrall et al. (2017). The mass flux of the low-order and the high-order methods are identical on each finite element cell. This allows for mass preserving and invariant-domain preserving limiting.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s10915-024-02592-4
