Optimal-transport-based mesh adaptivity on the plane and sphere using
finite elements
finite elements
File(s)paper-3.pdf (1.88 MB)
Accepted version
Author(s)
McRae, Andrew TT
Cotter, Colin J
Budd, Chris J
Type
Journal Article
Abstract
In moving mesh methods, the underlying mesh is dynamically adapted without
changing the connectivity of the mesh. We specifically consider the generation
of meshes which are adapted to a scalar monitor function through
equidistribution. Together with an optimal transport condition, this leads to a
Monge-Amp\`ere equation for a scalar mesh potential. We adapt an existing
finite element scheme for the standard Monge-Amp\`ere equation to this mesh
generation problem; this is a mixed finite element scheme, in which an extra
discrete variable is introduced to represent the Hessian matrix of second
derivatives. The problem we consider has additional nonlinearities over the
basic Monge-Amp\`ere equation due to the implicit dependence of the monitor
function on the resulting mesh. We also derive the equivalent
Monge-Amp\`ere-like equation for generating meshes on the sphere. The finite
element scheme is extended to the sphere, and we provide numerical examples.
All numerical experiments are performed using the open-source finite element
framework Firedrake.
changing the connectivity of the mesh. We specifically consider the generation
of meshes which are adapted to a scalar monitor function through
equidistribution. Together with an optimal transport condition, this leads to a
Monge-Amp\`ere equation for a scalar mesh potential. We adapt an existing
finite element scheme for the standard Monge-Amp\`ere equation to this mesh
generation problem; this is a mixed finite element scheme, in which an extra
discrete variable is introduced to represent the Hessian matrix of second
derivatives. The problem we consider has additional nonlinearities over the
basic Monge-Amp\`ere equation due to the implicit dependence of the monitor
function on the resulting mesh. We also derive the equivalent
Monge-Amp\`ere-like equation for generating meshes on the sphere. The finite
element scheme is extended to the sphere, and we provide numerical examples.
All numerical experiments are performed using the open-source finite element
framework Firedrake.
Date Issued
2018-03-01
Date Acceptance
2018-01-16
Citation
SIAM Journal on Scientific Computing, 2018, 40, pp.A1121-A1148
Start Page
A1121
End Page
A1148
Journal / Book Title
SIAM Journal on Scientific Computing
Volume
40
Copyright Statement
© The Authors
Sponsor
Natural Environment Research Council (NERC)
Identifier
http://arxiv.org/abs/1612.08077v3
Grant Number
NE/M013634/1
Subjects
math.NA
math.NA
Notes
Updated following reviews, 36 pages
Date Publish Online
2018-04-24