Generalized Monge-Kantorovich Optimization for Grid Generation and Adaptation in Lp

The Monge-Kantorovich grid generation and adaptation scheme of is generalized from a variational principle based on L{sub 2} to a variational principle based on L{sub p}. A generalized Monge-Ampere (MA) equation is derived and its properties are discussed. Results for p > 1 are obtained and compared in terms of the quality of the resulting grid. We conclude that for the grid generation application, the formulation based on L{sub p} for p close to unity leads to serious problems associated with the boundary. Results for 1.5 {approx}< p {approx}< 2.5 are quite good, but there is a fairly narrow range around p = 2 where the results are close to optimal with respect to grid distortion. Furthermore, the Newton-Krylov methods used to solve the generalized MA equation perform best for p = 2.

[1]  J. Scanlon.  The Dirichlet problem for nonlinear elliptic equations. , 1955 .

[2]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

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

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

[5]  G. M.,et al.  Partial Differential Equations I , 2023, Applied Mathematical Sciences.

[6]  Giovanni Lapenta.  Variational grid adaptation based on the minimization of local truncation error: time-independent problems , 2004 .

[7]  J. U. Brackbill,et al.  FLIP MHD: a particle-cell method for magnetohydrodynamics , 1991 .

[8]  HakerSteven,et al.  Optimal Mass Transport for Registration and Warping , 2004 .

[9]  Lei Zhu,et al.  Optimal Mass Transport for Registration and Warping , 2004, International Journal of Computer Vision.

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

[11]  D. A. Knoll,et al.  An Implicit Energy-Conservative 2D Fokker—Planck Algorithm , 2000 .

[12]  William J. Rider,et al.  A Multigrid Preconditioned Newton-Krylov Method , 1999, SIAM J. Sci. Comput..

[13]  J. Brackbill,et al.  FLIP: A method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions , 1986 .

[14]  Gian Luca Delzanno,et al.  The fluid dynamic approach to equidistribution methods for grid adaptation , 2011, Comput. Phys. Commun..

[15]  Lei Zhu,et al.  An Image Morphing Technique Based on Optimal Mass Preserving Mapping , 2007, IEEE Transactions on Image Processing.

[16]  Gian Luca Delzanno,et al.  Robust, multidimensional mesh-motion based on Monge-Kantorovich equidistribution , 2011, J. Comput. Phys..

[17]  L. Evans,et al.  Differential equations methods for the Monge-Kantorovich mass transfer problem , 1999 .

[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]  J. Brackbill,et al.  A numerical method for suspension flow , 1991 .

[20]  Steven Haker,et al.  Minimizing Flows for the Monge-Kantorovich Problem , 2003, SIAM J. Math. Anal..