5,899 research outputs found
Graph Element Networks: adaptive, structured computation and memory
We explore the use of graph neural networks (GNNs) to model spatial processes
in which there is no a priori graphical structure. Similar to finite element
analysis, we assign nodes of a GNN to spatial locations and use a computational
process defined on the graph to model the relationship between an initial
function defined over a space and a resulting function in the same space. We
use GNNs as a computational substrate, and show that the locations of the nodes
in space as well as their connectivity can be optimized to focus on the most
complex parts of the space. Moreover, this representational strategy allows the
learned input-output relationship to generalize over the size of the underlying
space and run the same model at different levels of precision, trading
computation for accuracy. We demonstrate this method on a traditional PDE
problem, a physical prediction problem from robotics, and learning to predict
scene images from novel viewpoints.Comment: Accepted to ICML 201
Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes
We present a new family of high order accurate fully discrete one-step
Discontinuous Galerkin (DG) finite element schemes on moving unstructured
meshes for the solution of nonlinear hyperbolic PDE in multiple space
dimensions, which may also include parabolic terms in order to model
dissipative transport processes. High order piecewise polynomials are adopted
to represent the discrete solution at each time level and within each spatial
control volume of the computational grid, while high order of accuracy in time
is achieved by the ADER approach. In our algorithm the spatial mesh
configuration can be defined in two different ways: either by an isoparametric
approach that generates curved control volumes, or by a piecewise linear
decomposition of each spatial control volume into simplex sub-elements. Our
numerical method belongs to the category of direct
Arbitrary-Lagrangian-Eulerian (ALE) schemes, where a space-time conservation
formulation of the governing PDE system is considered and which already takes
into account the new grid geometry directly during the computation of the
numerical fluxes. Our new Lagrangian-type DG scheme adopts the novel a
posteriori sub-cell finite volume limiter method, in which the validity of the
candidate solution produced in each cell by an unlimited ADER-DG scheme is
verified against a set of physical and numerical detection criteria. Those
cells which do not satisfy all of the above criteria are flagged as troubled
cells and are recomputed with a second order TVD finite volume scheme. The
numerical convergence rates of the new ALE ADER-DG schemes are studied up to
fourth order in space and time and several test problems are simulated.
Finally, an application inspired by Inertial Confinement Fusion (ICF) type
flows is considered by solving the Euler equations and the PDE of viscous and
resistive magnetohydrodynamics (VRMHD).Comment: 39 pages, 21 figure
On the mesh nonsingularity of the moving mesh PDE method
The moving mesh PDE (MMPDE) method for variational mesh generation and
adaptation is studied theoretically at the discrete level, in particular the
nonsingularity of the obtained meshes. Meshing functionals are discretized
geometrically and the MMPDE is formulated as a modified gradient system of the
corresponding discrete functionals for the location of mesh vertices. It is
shown that if the meshing functional satisfies a coercivity condition, then the
mesh of the semi-discrete MMPDE is nonsingular for all time if it is
nonsingular initially. Moreover, the altitudes and volumes of its elements are
bounded below by positive numbers depending only on the number of elements, the
metric tensor, and the initial mesh. Furthermore, the value of the discrete
meshing functional is convergent as time increases, which can be used as a
stopping criterion in computation. Finally, the mesh trajectory has limiting
meshes which are critical points of the discrete functional. The convergence of
the mesh trajectory can be guaranteed when a stronger condition is placed on
the meshing functional. Two meshing functionals based on alignment and
equidistribution are known to satisfy the coercivity condition. The results
also hold for fully discrete systems of the MMPDE provided that the time step
is sufficiently small and a numerical scheme preserving the property of
monotonically decreasing energy is used for the temporal discretization of the
semi-discrete MMPDE. Numerical examples are presented.Comment: Revised and improved version of the WIAS preprin
- …