skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Robust, multidimensional mesh motion based on Monge-Kantorovich equidistribution

Abstract

Mesh-motion (r-refinement) grid adaptivity schemes are attractive due to their potential to minimize the numerical error for a prescribed number of degrees of freedom. However, a key roadblock to a widespread deployment of this class of techniques has been the formulation of robust, reliable mesh-motion governing principles, which (1) guarantee a solution in multiple dimensions (2D and 3D), (2) avoid grid tangling (or folding of the mesh, whereby edges of a grid cell cross somewhere in the domain), and (3) can be solved effectively and efficiently. In this study, we formulate such a mesh-motion governing principle, based on volume equidistribution via Monge-Kantorovich optimization (MK). In earlier publications [1] and [2], the advantages of this approach with regard to these points have been demonstrated for the time-independent case. In this study, we demonstrate that Monge-Kantorovich equidistribution can in fact be used effectively in a time-stepping context, and delivers an elegant solution to the otherwise pervasive problem of grid tangling in mesh-motion approaches, without resorting to ad hoc time-dependent terms (as in moving-mesh PDEs, or MMPDEs [3] and [4]). We explore two distinct r-refinement implementations of MK: the direct method, where the current mesh relates to an initial, unchanging mesh, and themore » sequential method, where the current mesh is related to the previous one in time. We demonstrate that the direct approach is superior with regard to mesh distortion and robustness. The properties of the approach are illustrated with a hyperbolic PDE, the advection of a passive scalar, in 2D and 3D. Velocity flow fields with and without flow shear are considered. Three-dimensional grid, time-step, and nonlinear tolerance convergence studies are presented which demonstrate the optimality of the approach.« less

Authors:
 [1];  [1]
  1. Los Alamos National Laboratory (LANL)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
991694
DOE Contract Number:  
DE-AC05-00OR22725
Resource Type:
Journal Article
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 230; Journal Issue: 1; Journal ID: ISSN 0021-9991
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ADVECTION; CONVERGENCE; DEGREES OF FREEDOM; DIMENSIONS; OPTIMIZATION; SHEAR; TOLERANCE; VELOCITY

Citation Formats

Chacon De La Rosa, Luis, Delzanno, Gian Luca, and Finn, John M. Robust, multidimensional mesh motion based on Monge-Kantorovich equidistribution. United States: N. p., 2011. Web. doi:10.1016/j.jcp.2010.09.013.
Chacon De La Rosa, Luis, Delzanno, Gian Luca, & Finn, John M. Robust, multidimensional mesh motion based on Monge-Kantorovich equidistribution. United States. https://doi.org/10.1016/j.jcp.2010.09.013
Chacon De La Rosa, Luis, Delzanno, Gian Luca, and Finn, John M. 2011. "Robust, multidimensional mesh motion based on Monge-Kantorovich equidistribution". United States. https://doi.org/10.1016/j.jcp.2010.09.013.
@article{osti_991694,
title = {Robust, multidimensional mesh motion based on Monge-Kantorovich equidistribution},
author = {Chacon De La Rosa, Luis and Delzanno, Gian Luca and Finn, John M},
abstractNote = {Mesh-motion (r-refinement) grid adaptivity schemes are attractive due to their potential to minimize the numerical error for a prescribed number of degrees of freedom. However, a key roadblock to a widespread deployment of this class of techniques has been the formulation of robust, reliable mesh-motion governing principles, which (1) guarantee a solution in multiple dimensions (2D and 3D), (2) avoid grid tangling (or folding of the mesh, whereby edges of a grid cell cross somewhere in the domain), and (3) can be solved effectively and efficiently. In this study, we formulate such a mesh-motion governing principle, based on volume equidistribution via Monge-Kantorovich optimization (MK). In earlier publications [1] and [2], the advantages of this approach with regard to these points have been demonstrated for the time-independent case. In this study, we demonstrate that Monge-Kantorovich equidistribution can in fact be used effectively in a time-stepping context, and delivers an elegant solution to the otherwise pervasive problem of grid tangling in mesh-motion approaches, without resorting to ad hoc time-dependent terms (as in moving-mesh PDEs, or MMPDEs [3] and [4]). We explore two distinct r-refinement implementations of MK: the direct method, where the current mesh relates to an initial, unchanging mesh, and the sequential method, where the current mesh is related to the previous one in time. We demonstrate that the direct approach is superior with regard to mesh distortion and robustness. The properties of the approach are illustrated with a hyperbolic PDE, the advection of a passive scalar, in 2D and 3D. Velocity flow fields with and without flow shear are considered. Three-dimensional grid, time-step, and nonlinear tolerance convergence studies are presented which demonstrate the optimality of the approach.},
doi = {10.1016/j.jcp.2010.09.013},
url = {https://www.osti.gov/biblio/991694}, journal = {Journal of Computational Physics},
issn = {0021-9991},
number = 1,
volume = 230,
place = {United States},
year = {2011},
month = {1}
}