The geometry of r-adaptive meshes generated using optimal transport methods

JOURNAL OF COMPUTATIONAL PHYSICS

卷:282

Article

Budd, C. J.¹  Russell, R. D.²  Walsh, E.² 

[1] Univ Bath, Bath BA2 7AY, Avon, England

[2] Simon Fraser Univ, Burnaby, BC V5N 1S6, Canada

关键词: Alignment;    Anisotropy;    Mesh adaptation;    Metric tensor;    Monge-Ampere;   

DOI  :  10.1016/j.jcp.2014.11.007

来源: Elsevier

【 摘 要 】

The principles of mesh equidistribution and alignment play a fundamental role in the design of adaptive methods, and a metric tensor and mesh metric are useful theoretical tools for understanding a method's level of mesh alignment, or anisotropy. We consider a mesh redistribution method based on the Monge-Ampere equation which combines equidistribution of a given scalar density function with optimal transport. It does not involve explicit use of a metric tensor, although such a tensor must exist for the method, and an interesting question to ask is whether or not the alignment produced by the metric gives an anisotropic mesh. For model problems with a linear feature and with a radially symmetric feature, we derive the exact form of the metric, which involves expressions for its eigenvalues and eigenvectors. The eigenvectors are shown to be orthogonal and tangential to the feature, and the ratio of the eigenvalues (corresponding to the level of anisotropy) is shown to depend, both locally and globally, on the value of the density function and the amount of curvature. We thereby demonstrate how the optimal transport method produces an anisotropic mesh along a given feature while equidistributing a suitably chosen scalar density function. Numerical results are given to verify these results and to demonstrate how the analysis is useful for problems involving more complex features, including for a non-trivial time dependant nonlinear PDE which evolves narrow and curved reaction fronts. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free