Source author record

Andrew R. Winters

Andrew R. Winters appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

A new class of entropy stable fluctuations for the discontinuous Galerkin method with application to the Saint-Venant-Exner model

In this work we consider entropy stable discontinuous Galerkin methods applied to nonconservative hyperbolic systems. We introduce a new class of entropy conservative fluctuations that allow us to construct entropy conservative schemes without any system-specific derivations. We demonstrate that a loss of entropy symmetrization for nonconservative systems restricts the design of entropy stable fluctuations and propose a novel blending procedure to construct entropy stable dissipation terms from general numerical viscosity matrices. The resulting methodology is applied to develop a high-order, entropy stable, and well-balanced approximation for the Saint-Venant-Exner system. Numerical tests are presented to verify the theoretical findings and demonstrate the performance and robustness of the proposed scheme.

preprint2022arXiv

A linear and nonlinear analysis of the shallow water equations and its impact on boundary conditions

We derive boundary conditions and estimates based on the energy and entropy analysis of systems of the nonlinear shallow water equations in two spatial dimensions. It is shown that the energy method provides more details, but is fully consistent with the entropy analysis. The details brought forward by the nonlinear energy analysis allow us to pinpoint where the difference between the linear and nonlinear analysis originate. We find that the result from the linear analysis does not necessarily hold in the nonlinear case. The nonlinear analysis leads in general to a different minimal number of boundary conditions compared with the linear analysis. In particular, and contrary to the linear case, the magnitude of the flow does not influence the number of required boundary conditions.

preprint2022arXiv

Adaptive numerical simulations with Trixi.jl: A case study of Julia for scientific computing

We present Trixi.jl, a Julia package for adaptive high-order numerical simulations of hyperbolic partial differential equations. Utilizing Julia's strengths, Trixi.jl is extensible, easy to use, and fast. We describe the main design choices that enable these features and compare Trixi.jl with a mature open source Fortran code that uses the same numerical methods. We conclude with an assessment of Julia for simulation-focused scientific computing, an area that is still dominated by traditional high-performance computing languages such as C, C++, and Fortran.

preprint2020arXiv

Construction of Modern Robust Nodal Discontinuous Galerkin Spectral Element Methods for the Compressible Navier-Stokes Equations

Discontinuous Galerkin (DG) methods have a long history in computational physics and engineering to approximate solutions of partial differential equations due to their high-order accuracy and geometric flexibility. However, DG is not perfect and there remain some issues. Concerning robustness, DG has undergone an extensive transformation over the past seven years into its modern form that provides statements on solution boundedness for linear and nonlinear problems. This chapter takes a constructive approach to introduce a modern incarnation of the DG spectral element method for the compressible Navier-Stokes equations in a three-dimensional curvilinear context. The groundwork of the numerical scheme comes from classic principles of spectral methods including polynomial approximations and Gauss-type quadratures. We identify aliasing as one underlying cause of the robustness issues for classical DG spectral methods. Removing said aliasing errors requires a particular differentiation matrix and careful discretization of the advective flux terms in the governing equations.

preprint2020arXiv

Stable filtering procedures for nodal discontinuous Galerkin methods

We prove that the most common filtering procedure for nodal discontinuous Galerkin (DG) methods is stable. The proof exploits that the DG approximation is constructed from polynomial basis functions and that integrals are approximated with high-order accurate Legendre-Gauss-Lobatto quadrature. The theoretical discussion serves to re-contextualize stable filtering results for finite difference methods into the DG setting. It is shown that the stability of the filtering is equivalent to a particular contractivity condition borrowed from the analysis of so-called transmission problems. As such, the temporal stability proof relies on the fact that the underlying spatial discretization of the problem possesses a semi-discrete bound on the solution. Numerical tests are provided to verify and validate the underlying theoretical results.

preprint2016arXiv

A Novel High-Order, Entropy Stable, 3D AMR MHD Solver with Guaranteed Positive Pressure

We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum, and energy and is entropy stable. The method is designed to treat the divergence-free constraint on the magnetic field in a similar fashion to a hyperbolic divergence cleaning technique. The solver described herein is especially well-suited for flows involving strong discontinuities. Furthermore, we present a new formulation to guarantee positivity of the pressure. We present the underlying theory and implementation of the new solver into the multi-physics, multi-scale adaptive mesh refinement (AMR) simulation code $\texttt{FLASH}$ (http://flash.uchicago.edu). The accuracy, robustness and computational efficiency is demonstrated with a number of tests, including comparisons to available MHD implementations in $\texttt{FLASH}$.

preprint2016arXiv

An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry

We design an arbitrary high-order accurate nodal discontinuous Galerkin spectral element approximation for the nonlinear two dimensional shallow water equations with non-constant, possibly discontinuous, bathymetry on unstructured, possibly curved, quadrilateral meshes. The scheme is derived from an equivalent flux differencing formulation of the split form of the equations. We prove that this discretisation exactly preserves the local mass and momentum. Furthermore, combined with a special numerical interface flux function, the method exactly preserves the mathematical entropy, which is the total energy for the shallow water equations. By adding a specific form of interface dissipation to the baseline entropy conserving scheme we create a provably entropy stable scheme. That is, the numerical scheme discretely satisfies the second law of thermodynamics. Finally, with a particular discretisation of the bathymetry source term we prove that the numerical approximation is well-balanced. We provide numerical examples that verify the theoretical findings and furthermore provide an application of the scheme for a partial break of a curved dam test problem.

preprint2016arXiv

Hybrid Entropy Stable HLL-Type Riemann Solvers for Hyperbolic Conservation Laws

It is known that HLL-type schemes are more dissipative than schemes based on characteristic decompositions. However, HLL-type methods offer greater flexibility to large systems of hyperbolic conservation laws because the eigenstructure of the flux Jacobian is not needed. We demonstrate in the present work that several HLL-type Riemann solvers are provably entropy stable. Further, we provide convex combinations of standard dissipation terms to create hybrid HLL-type methods that have less dissipation while retaining entropy stability. The decrease in dissipation is demonstrated for the ideal MHD equations with a numerical example.

preprint2016arXiv

Split Form Nodal Discontinuous Galerkin Schemes with Summation-By-Parts Property for the Compressible Euler Equations

Fisher and Carpenter (\textit{High-order entropy stable finite difference schemes for non-linear conservation laws: Finite domains, Journal of Computational Physics, 252:518--557, 2013}) found a remarkable equivalence of general diagonal norm high-order summation-by-parts operators to a subcell based high-order finite volume formulation. This equivalence enables the construction of provably entropy stable schemes by a specific choice of the subcell finite volume flux. We show that besides the construction of entropy stable high order schemes, a careful choice of subcell finite volume fluxes generates split formulations of quadratic or cubic terms. Thus, by changing the subcell finite volume flux to a specific choice, we are able to generate, in a systematic way, all common split forms of the compressible Euler advection terms, such as the Ducros splitting and the Kennedy and Gruber splitting. Although these split forms are not entropy stable, we present a systematic way to prove which of those split forms are at least kinetic energy preserving. With this, we show we construct a unified high-order split form DG framework. We investigate with three dimensional numerical simulations of the inviscid Taylor-Green vortex and show that the new split forms enhance the robustness of high order simulations in comparison to the standard scheme when solving turbulent vortex dominated flows. In fact, we show that for certain test cases, the novel split form discontinuous Galerkin schemes are more robust than the discontinuous Galerkin scheme with over-integration.

preprint2015arXiv

A Provably Stable Discontinuous Galerkin Spectral Element Approximation for Moving Hexahedral Meshes

We design a novel provably stable discontinuous Galerkin spectral element (DGSEM) approximation to solve systems of conservation laws on moving domains. To incorporate the motion of the domain, we use an arbitrary Lagrangian-Eulerian formulation to map the governing equations to a fixed reference domain. The approximation is made stable by a discretization of a skew-symmetric formulation of the problem. We prove that the discrete approximation is stable, conservative and, for constant coefficient problems, maintains the free-stream preservation property. We also provide details on how to add the new skew-symmetric ALE approximation to an existing discontinuous Galerkin spectral element code. Lastly, we provide numerical support of the theoretical results.

preprint2015arXiv

Affordable, Entropy Conserving and Entropy Stable Flux Functions for the Ideal MHD Equations

In this work, we design an entropy stable, finite volume approximation for the ideal magnetohydrodynamics (MHD) equations. The method is novel as we design an affordable analytical expression of the numerical interface flux function that discretely preserves the entropy of the system. To guarantee the discrete conservation of entropy requires the addition of a particular source term to the ideal MHD system. Exact entropy conserving schemes cannot dissipate energy at shocks, thus to compute accurate solutions to problems that may develop shocks, we determine a dissipation term to guarantee entropy stability for the numerical scheme. Numerical tests are performed to demonstrate the theoretical findings of entropy conservation and robustness.

preprint2015arXiv

An Entropy Stable Finite Volume Scheme for the Equations of Shallow Water Magnetohydrodynamics

In this work, we design an entropy stable, finite volume approximation for the shallow water magnetohydrodynamics (SWMHD) equations. The method is novel as we design an affordable analytical expression of the numerical interface flux function that exactly preserves the entropy, which is also the total energy for the SWMHD equations. To guarantee the discrete conservation of entropy requires a special treatment of a consistent source term for the SWMHD equations. With the goal of solving problems that may develop shocks, we determine a dissipation term to guarantee entropy stability for the numerical scheme. Numerical tests are performed to demonstrate the theoretical findings of entropy conservation and robustness.