Adaptive discrete harmonic grid generation (Q1005668)

The grid generation problem can be described as the construction of a continuous parametrization from the unit square of a simple connected domain, \(\Omega\), in the plane in such a way that the boundary of the unit square parametrizes the boundary of the region \(\Omega\). The authors present a new adaptive-harmonic structured grid generation method. It is based on a functional that shares a common set of minimizers with Ivanenko's harmonic functional [cf. \textit{A. A. Charakhch'yan} and \textit{S. A. Ivanenko}, U.S.S.R. Comput. Math. Math. Phys. 28, No.~2, 126--133 (1988; Zbl 0671.65095)]. A discretization of the variational principle is used instead of the Euler-Lagrange equations. This allows to control directly the interior nodes of the grid. An unconstrained optimization process related to a continuation parameter is used to guarantee the convexity of the grid cells. Several numerical examples of grids generated on highly irregular 2D regions, show the effective performance of the proposed method and its robustness when preset monitor functions are given.