Optimal mass transport for higher dimensional adaptive grid generation

In this work, we describe an approach for higher dimensional adaptive grid generation based on solving the L^2 Monge-Kantorovich problem (MKP) which is a special case of the classical optimal mass transportation problem. Two methods are developed for computing the coordinate transformation used to define the grid adaptation. For the first method, the transformation is determined by solving a parabolic Monge-Ampere equation for a steady state solution. For the second method, the grid movement is determined from the velocity field obtained by solving a fluid dynamics formulation of the L^2 MKP. Several numerical experiments are presented to demonstrate the performance of the MKP methods and to compare them with some related adaptive grid methods. The experimental results demonstrate that the MKP methods show promise as effective and reliable methods for higher dimensional adaptive grid generation.

[1]  Gian Luca Delzanno,et al.  Grid Generation and Adaptation by Monge-Kantorovich Optimization in Two and Three Dimensions , 2008, IMR.

[2]  J. F. Williams,et al.  Parabolic Monge-Ampère methods for blow-up problems in several spatial dimensions , 2006 .

[3]  J. F. Williams,et al.  Moving Mesh Generation Using the Parabolic Monge--Amp[e-grave]re Equation , 2009, SIAM J. Sci. Comput..

[4]  WEIZHANG HUANG,et al.  A Moving Mesh Method Based on the Geometric Conservation Law , 2002, SIAM J. Sci. Comput..

[5]  Tao Tang,et al.  An adaptive mesh redistribution algorithm for convection-dominated problems , 2002 .

[6]  Robert D. Russell,et al.  Moving Mesh Methods for Problems with Blow-Up , 1996, SIAM J. Sci. Comput..

[7]  Weizhang Huang,et al.  A moving collocation method for solving time dependent partial differential equations , 1996 .

[8]  Allen R. Tannenbaum,et al.  An Efficient Numerical Method for the Solution of the L2 Optimal Mass Transfer Problem , 2010, SIAM J. Sci. Comput..

[9]  K. Guittet On the Time-Continuous Mass Transport Problem and Its Approximation by Augmented Lagrangian Techniques , 2003, SIAM J. Numer. Anal..

[10]  J. F. Williams,et al.  An efficient approach for the numerical solution of the Monge-Ampère equation , 2011 .

[11]  Mohamed H. Mahmoud Sulman Interpreting moving mesh methods and transformation techniques from various applications in terms of the mass transport problem , 2003 .

[12]  H. Schwetlick,et al.  Translating solutions for Gaußcurvature flows with Neumann boundary conditions , 2004 .

[13]  M. Knott,et al.  On the optimal mapping of distributions , 1984 .

[14]  Robert D. Russell,et al.  Moving Mesh Strategy Based on a Gradient Flow Equation for Two-Dimensional Problems , 1998, SIAM J. Sci. Comput..

[15]  Weizhang Huang,et al.  Moving Mesh Methods Based on Moving Mesh Partial Differential Equations , 1994 .

[16]  Guojun Liao,et al.  A new approach to grid generation , 1992 .

[17]  Weizhang Huang,et al.  Moving mesh partial differential equations (MMPDES) based on the equidistribution principle , 1994 .

[18]  Gian Luca Delzanno,et al.  An optimal robust equidistribution method for two-dimensional grid adaptation based on Monge-Kantorovich optimization , 2008, J. Comput. Phys..

[19]  P. Thomas,et al.  Geometric Conservation Law and Its Application to Flow Computations on Moving Grids , 1979 .

[20]  Thomas Y. Hou,et al.  An efficient dynamically adaptive mesh for potentially singular solutions , 2001 .

[21]  Yann Brenier,et al.  A computational fluid mechanics solution to the Monge-Kantorovich mass transfer problem , 2000, Numerische Mathematik.

[22]  Rao V. Garimella,et al.  Proceedings of the 17th International Meshing Roundtable , 2008 .

[23]  Knut Smoczyk,et al.  Neumann and second boundary value problems for Hessian and Gauß curvature flows , 2003 .

[24]  Oliver C. Schnürer,et al.  F Ur Mathematik in Den Naturwissenschaften Leipzig Translating Solutions to the Second Boundary Value Problem for Curvature Ows Translating Solutions to the Second Boundary Value Problem for Curvature Flows , 2022 .

[25]  Weizhang Huang,et al.  Analysis Of Moving Mesh Partial Differential Equations With Spatial Smoothing , 1997 .

[26]  Y. Brenier Polar Factorization and Monotone Rearrangement of Vector-Valued Functions , 1991 .

[27]  Weizhang Huang,et al.  Measuring Mesh Qualities and Application to Variational Mesh Adaptation , 2005, SIAM J. Sci. Comput..