Hamiltonian-Preserving Discontinuous Galerkin Methods for the Liouville Equation With Discontinuous Potential
- Award ID(s):
- 1753581
- PAR ID:
- 10416402
- Date Published:
- Journal Name:
- SIAM Journal on Scientific Computing
- Volume:
- 44
- Issue:
- 5
- ISSN:
- 1064-8275
- Page Range / eLocation ID:
- A3317 to A3340
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
The discontinuous Petrov–Galerkin (DPG) method is a Petrov–Galerkin finite element method with test functions designed for obtaining stability. These test functions are computable locally, element by element, and are motivated by optimal test functions which attain the supremum in an inf-sup condition. A profound consequence of the use of nearly optimal test functions is that the DPG method can inherit the stability of the (undiscretized) variational formulation, be it coercive or not. This paper combines a presentation of the fundamentals of the DPG ideas with a review of the ongoing research on theory and applications of the DPG methodology. The scope of the presented theory is restricted to linear problems on Hilbert spaces, but pointers to extensions are provided. Multiple viewpoints to the basic theory are provided. They show that the DPG method is equivalent to a method which minimizes a residual in a dual norm, as well as to a mixed method where one solution component is an approximate error representation function. Being a residual minimization method, the DPG method yields Hermitian positive definite stiffness matrix systems even for non-self-adjoint boundary value problems. Having a built-in error representation, the method has the out-of-the-box feature that it can immediately be used in automatic adaptive algorithms. Contrary to standard Galerkin methods, which are uninformed about test and trial norms, the DPG method must be equipped with a concrete test norm which enters the computations. Of particular interest are variational formulations in which one can tailor the norm to obtain robust stability. Key techniques to rigorously prove convergence of DPG schemes, including construction of Fortin operators, which in the DPG case can be done element by element, are discussed in detail. Pointers to open frontiers are presented.more » « less
-
Abstract Barrier islands and spits tend to migrate landward in response to sea‐level rise through the storm‐driven process of overwash, but overwash flux depends on the height of the frontal dunes. Here, we explore this fundamental linkage between dune dynamics and barrier migration using the new model Barrier3D. Our experiments demonstrate that discontinuous barrier retreat is a prevalent behavior that can arise directly from the bistability of foredune height, occurring most likely when the storm return period and characteristic time scale of dune growth are of similar magnitudes. Under conditions of greater storm intensity, discontinuous retreat becomes the dominant behavior of barriers that were previously stable. Alternatively, higher rates of sea‐level rise decrease the overall likelihood of discontinuous retreat in favor of continuous transgression. We find that internal dune dynamics, while previously neglected in exploratory barrier modeling, are an essential component of barrier evolution on time scales relevant to coastal management.more » « less
An official website of the United States government

