11,125 research outputs found
Automated code generation for discontinuous Galerkin methods
A compiler approach for generating low-level computer code from high-level
input for discontinuous Galerkin finite element forms is presented. The input
language mirrors conventional mathematical notation, and the compiler generates
efficient code in a standard programming language. This facilitates the rapid
generation of efficient code for general equations in varying spatial
dimensions. Key concepts underlying the compiler approach and the automated
generation of computer code are elaborated. The approach is demonstrated for a
range of common problems, including the Poisson, biharmonic,
advection--diffusion and Stokes equations
Efficient Localization of Discontinuities in Complex Computational Simulations
Surrogate models for computational simulations are input-output
approximations that allow computationally intensive analyses, such as
uncertainty propagation and inference, to be performed efficiently. When a
simulation output does not depend smoothly on its inputs, the error and
convergence rate of many approximation methods deteriorate substantially. This
paper details a method for efficiently localizing discontinuities in the input
parameter domain, so that the model output can be approximated as a piecewise
smooth function. The approach comprises an initialization phase, which uses
polynomial annihilation to assign function values to different regions and thus
seed an automated labeling procedure, followed by a refinement phase that
adaptively updates a kernel support vector machine representation of the
separating surface via active learning. The overall approach avoids structured
grids and exploits any available simplicity in the geometry of the separating
surface, thus reducing the number of model evaluations required to localize the
discontinuity. The method is illustrated on examples of up to eleven
dimensions, including algebraic models and ODE/PDE systems, and demonstrates
improved scaling and efficiency over other discontinuity localization
approaches
A high-order semi-explicit discontinuous Galerkin solver for 3D incompressible flow with application to DNS and LES of turbulent channel flow
We present an efficient discontinuous Galerkin scheme for simulation of the
incompressible Navier-Stokes equations including laminar and turbulent flow. We
consider a semi-explicit high-order velocity-correction method for time
integration as well as nodal equal-order discretizations for velocity and
pressure. The non-linear convective term is treated explicitly while a linear
system is solved for the pressure Poisson equation and the viscous term. The
key feature of our solver is a consistent penalty term reducing the local
divergence error in order to overcome recently reported instabilities in
spatially under-resolved high-Reynolds-number flows as well as small time
steps. This penalty method is similar to the grad-div stabilization widely used
in continuous finite elements. We further review and compare our method to
several other techniques recently proposed in literature to stabilize the
method for such flow configurations. The solver is specifically designed for
large-scale computations through matrix-free linear solvers including efficient
preconditioning strategies and tensor-product elements, which have allowed us
to scale this code up to 34.4 billion degrees of freedom and 147,456 CPU cores.
We validate our code and demonstrate optimal convergence rates with laminar
flows present in a vortex problem and flow past a cylinder and show
applicability of our solver to direct numerical simulation as well as implicit
large-eddy simulation of turbulent channel flow at as well as
.Comment: 28 pages, in preparation for submission to Journal of Computational
Physic
GPU-accelerated discontinuous Galerkin methods on hybrid meshes
We present a time-explicit discontinuous Galerkin (DG) solver for the
time-domain acoustic wave equation on hybrid meshes containing vertex-mapped
hexahedral, wedge, pyramidal and tetrahedral elements. Discretely energy-stable
formulations are presented for both Gauss-Legendre and Gauss-Legendre-Lobatto
(Spectral Element) nodal bases for the hexahedron. Stable timestep restrictions
for hybrid meshes are derived by bounding the spectral radius of the DG
operator using order-dependent constants in trace and Markov inequalities.
Computational efficiency is achieved under a combination of element-specific
kernels (including new quadrature-free operators for the pyramid), multi-rate
timestepping, and acceleration using Graphics Processing Units.Comment: Submitted to CMAM
Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
We present a class of spline finite element methods for time-domain wave
propagation which are particularly amenable to explicit time-stepping. The
proposed methods utilize a discontinuous Galerkin discretization to enforce
continuity of the solution field across geometric patches in a multi-patch
setting, which yields a mass matrix with convenient block diagonal structure.
Over each patch, we show how to accurately and efficiently invert mass matrices
in the presence of curved geometries by using a weight-adjusted approximation
of the mass matrix inverse. This approximation restores a tensor product
structure while retaining provable high order accuracy and semi-discrete energy
stability. We also estimate the maximum stable timestep for spline-based finite
elements and show that the use of spline spaces result in less stringent CFL
restrictions than equivalent piecewise continuous or discontinuous finite
element spaces. Finally, we explore the use of optimal knot vectors based on L2
n-widths. We show how the use of optimal knot vectors can improve both
approximation properties and the maximum stable timestep, and present a simple
heuristic method for approximating optimal knot positions. Numerical
experiments confirm the accuracy and stability of the proposed methods
Higher-order finite element methods for elliptic problems with interfaces
We present higher-order piecewise continuous finite element methods for
solving a class of interface problems in two dimensions. The method is based on
correction terms added to the right-hand side in the standard variational
formulation of the problem. We prove optimal error estimates of the methods on
general quasi-uniform and shape regular meshes in maximum norms. In addition,
we apply the method to a Stokes interface problem, adding correction terms for
the velocity and the pressure, obtaining optimal convergence results.Comment: 26 pages, 6 figures. An earlier version of this paper appeared on
November 13, 2014 in
http://www.brown.edu/research/projects/scientific-computing/reports/201
A free-energy stable nodal discontinuous Galerkin approximation with summation-by-parts property for the Cahn-Hilliard equation
We present a nodal Discontinuous Galerkin (DG) scheme for the Cahn-Hilliard
equation that satisfies the summation-by-parts simultaneous-approximation-term
(SBP-SAT) property. The latter permits us to show that the discrete free-energy
is bounded, and as a result, the scheme is provably stable. The scheme and the
stability proof are presented for general curvilinear three-dimensional
hexahedral meshes. We use the Bassi-Rebay 1 (BR1) scheme to compute interface
fluxes, and an IMplicit-EXplicit (IMEX) scheme to integrate in time. Lastly, we
test the theoretical findings numerically and present examples for two and
three-dimensional problems
- …