Optimization-based mesh correction with volume and convexity constraints
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
In this study, we consider the problem of finding a mesh such that 1) it is the closest, with respect to a suitable metric, to a given source mesh having the same connectivity, and 2) the volumes of its cells match a set of prescribed positive values that are not necessarily equal to the cell volumes in the source mesh. This volume correction problem arises in important simulation contexts, such as satisfying a discrete geometric conservation law and solving transport equations by incremental remapping or similar semi-Lagrangian transport schemes. In this paper we formulate volume correction as a constrained optimization problem in which the distance to the source mesh defines an optimization objective, while the prescribed cell volumes, mesh validity and/or cell convexity specify the constraints. We solve this problem numerically using a sequential quadratic programming (SQP) method whose performance scales with the mesh size. To achieve scalable performance we develop a specialized multigrid-based preconditioner for optimality systems that arise in the application of the SQP method to the volume correction problem. Numerical examples illustrate the importance of volume correction, and showcase the accuracy, robustness and scalability of our approach.
- Research Organization:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE National Nuclear Security Administration (NNSA)
- Grant/Contract Number:
- 14-017511; AC52-06NA25396; AC04-94AL85000
- OSTI ID:
- 1240603
- Alternate ID(s):
- OSTI ID: 1254321; OSTI ID: 1348256
- Report Number(s):
- LA-UR-15-22070; SAND-2015-1521J; PII: S0021999116001224
- Journal Information:
- Journal of Computational Physics, Vol. 313, Issue C; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
2D Mesh smoothing based on Markov chain method
|
journal | May 2019 |
Mixing Local and Nonlocal Evolution Equations
|
journal | January 2023 |
Homogenization for Nonlocal Evolution Problems with Three Different Smooth Kernels
|
journal | February 2023 |
An efficient implementation of mass conserving characteristic-based schemes in 2D and 3D | preprint | January 2019 |
Homogenization for nonlocal problems with smooth kernels | preprint | January 2020 |
Similar Records
Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid–structure interaction, and free surface flow: Part II
Towards Adaptive Mesh Refinement for the Spectral Element Solver Nek5000