An optimal robust equidistribution method for two-dimensional grid generation based on Monge-Kantorovich optimization

Speaker: Dr. John M. Finn Lawrence Livermore National laboratory

Abstract: I will present a new cell-area equidistribution method for grid adaptation, based on Monge-Kantorovich optimization (or Monge-Kantorovich optimal transport). The method is based on a rigorous variational principle, in which the L_2 norm of the grid displacement is minimized, constrained locally to produce a prescribed positive-definite cell volume distribution. The procedure involves solving the Monge-Ampere equation,a single, nonlinear, elliptic scalar equation with no free parameters, and with proved solution existence and uniqueness theorems. We show that, for sufficiently small grid displacement, this method also minimizes the mean grid-cell distortion, measured by the trace of the covariant metric tensor. We solve the Monge-Ampere equation numerically with a Jacobian-Free Newton-Krylov method. The ellipticity property of the Monge-Ampere equation allows multigrid preconditioning techniques to be used effectively, delivering a scalable algorithm under grid refinement. Several challenging test cases demonstrate that this method produces optimal grids in which the constraint is satisfied numerically to truncation error. We also compare this method to the well known deformation method [G. Liao and D. Anderson, Appl. Anal., v.44, p.285 (1992)]. We show that the new method achieves the desired equidistributed grid using comparable computational time, but with considerably better grid quality than the deformation method. I will present recent work with more general boundaries in 2D and in a cube in 3D.

Last modified: Wed Jul 2 11:28:43 EDT 2008