1,114 research outputs found
Dynamic p-enrichment schemes for multicomponent reactive flows
We present a family of p-enrichment schemes. These schemes may be separated
into two basic classes: the first, called \emph{fixed tolerance schemes}, rely
on setting global scalar tolerances on the local regularity of the solution,
and the second, called \emph{dioristic schemes}, rely on time-evolving bounds
on the local variation in the solution. Each class of -enrichment scheme is
further divided into two basic types. The first type (the Type I schemes)
enrich along lines of maximal variation, striving to enhance stable solutions
in "areas of highest interest." The second type (the Type II schemes) enrich
along lines of maximal regularity in order to maximize the stability of the
enrichment process. Each of these schemes are tested over a pair of model
problems arising in coastal hydrology. The first is a contaminant transport
model, which addresses a declinature problem for a contaminant plume with
respect to a bay inlet setting. The second is a multicomponent chemically
reactive flow model of estuary eutrophication arising in the Gulf of Mexico.Comment: 29 pages, 7 figures, 3 table
An entropy stable discontinuous Galerkin method for the shallow water equations on curvilinear meshes with wet/dry fronts accelerated by GPUs
We extend the entropy stable high order nodal discontinuous Galerkin spectral
element approximation for the non-linear two dimensional shallow water
equations presented by Wintermeyer et al. [N. Wintermeyer, A. R. Winters, G. J.
Gassner, and D. A. Kopriva. An entropy stable nodal discontinuous Galerkin
method for the two dimensional shallow water equations on unstructured
curvilinear meshes with discontinuous bathymetry. Journal of Computational
Physics, 340:200-242, 2017] with a shock capturing technique and a positivity
preservation capability to handle dry areas. The scheme preserves the entropy
inequality, is well-balanced and works on unstructured, possibly curved,
quadrilateral meshes. For the shock capturing, we introduce an artificial
viscosity to the equations and prove that the numerical scheme remains entropy
stable. We add a positivity preserving limiter to guarantee non-negative water
heights as long as the mean water height is non-negative. We prove that
non-negative mean water heights are guaranteed under a certain additional time
step restriction for the entropy stable numerical interface flux. We implement
the method on GPU architectures using the abstract language OCCA, a unified
approach to multi-threading languages. We show that the entropy stable scheme
is well suited to GPUs as the necessary extra calculations do not negatively
impact the runtime up to reasonably high polynomial degrees (around ). We
provide numerical examples that challenge the shock capturing and positivity
properties of our scheme to verify our theoretical findings
Space-time discontinuous Galerkin discretization of rotating shallow water equations on moving grids
A space-time discontinuous Galerkin (DG) discretization is presented for the (rotating) shallow water equations over varying topography. We formulate the space-time DG finite element discretization in an efficient and conservative discretization. The HLLC flux is used as numerical flux through the finite element boundaries. When discontinuities are present, we locally apply dissipation around these discontinuities with the help of Krivodonova's discontinuity indicator such that spurious oscillations are suppressed. The non-linear algebraic system resulting from the discretization is solved using a pseudo-time integration with a second-order five-stage Runge-Kutta method. A thorough verification of the space-time DG finite element method is undertaken by comparing numerical and exact solutions. We also carry out a discrete Fourier analysis of the one dimensional linear rotating shallow water equations to show that the method is unconditionally stable with minimal dispersion and dissipation error. The numerical scheme is validated in a novel way by considering various simulations of bore-vortex interactions in combination with a qualitative analysis of PV generation by non-uniform bores. Finally, the space-time DG method is particularly suited for problems where dynamic grid motion is required. To demonstrate this we simulate waves generated by a wave maker and verify these for low amplitude waves where linear theory is approximately valid
Well-balanced -adaptive and moving mesh space-time discontinuous Galerkin method for the shallow water equations
In this article we introduce a well-balanced discontinuous Galerkin method for the shallow water equations on moving meshes. Particular emphasis will be given on -adaptation in which mesh points of an initially uniform mesh move to concentrate in regions where interesting behaviour of the solution is observed. Obtaining well-balanced numerical schemes for the shallow water equations on fixed meshes is nontrivial and has been a topic of much research. In [S. Rhebergen, O. Bokhove, J.J.W. van der Vegt, Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations, J. Comput. Phys. 227 (2008) 1887–1922] we introduced a well-balanced discontinuous Galerkin method using the theory of weak solutions for nonconservative products introduced in [G. Dal Maso, P.G. LeFloch, F. Murat, Definition and weak stability of nonconservative products, J. Math. Pures Appl. 74 (1995) 483–548]. In this article we continue this approach and prove well-balancedness of a discontinuous Galerkin method for the shallow water equations on moving meshes. Numerical simulations are then performed to verify the -adaptive method in combination with the space-time discontinuous Galerkin method against analytical solutions and showing its robustness on more complex problems
- …