On one problem of 2D regular grid generation based on mappings

B. N. Azarenok ; A. A. Charakhch'yan

Matematičeskoe modelirovanie, Tome 26 (2014) no. 12, pp. 48-64.

Résumé

To study the known problem of regular grid generation by the Winslow method in a rectangular domain with a boundary kink known as the backstep, a high-accuracy numerical method for the inverse harmonic mapping of the unit square onto the domain with a certain mapping of the domain boundaries is developed. Behavior of the level line of the mapping which enter the point of the boundary kink is studied. Near the kink, the angle between the boundary and the straight line connecting a point on the level line with the point of the boundary kink is found as the function of the coordinate of the point on the level line in the unit square. It is shown that the level line of the mapping is in tangent to the boundary at the kink point. Near the kink point the mapping is non-quasiisometric. The regular grid of the intersection points between the level lines connected by straight lines contains a self-intersecting cell that remains when the grid step along the boundary decreases. Basing on the universal elliptical equations reproducing any nondegenerate mapping of the parametric rectangle onto a given domain, it is suggested a simple two-parametric control of grid nodes in the backstep that allows one to control effectively the slope angle of the grid line entering the point of the kink, thereby removing escape of grid lines from the domain boundary. In the case of the small number of grid points $31\times 31$, the nondegenerate grid is generated by selection of a suitable value of one parameter. When increasing the number of grid points 8 times in both direction (the grid $241\times 241$), the non-convex cells appear within the domain which are easily removed by using the variational barrier method. Another possibility to avoid non-convex cells is to decrease the grid dimension along the second direction (the grid $241\times 121$).

Détail

Keywords: structured grids, harmonic mapping, control metric.

Bibliographie

[1] Azarenok B. N., “Postroenie setki s ortogonalizatsiei i zadannym sguscheniem uzlov okolo granitsy dvumernoi oblasti”, Matematicheskoe modelirovanie, 23:2 (2011), 62–74 | Zbl

[2] Winslow A. M., “Numerical solution of quasilinear Poisson equation in nonuniform triangle mesh”, J. Comput. Phys., 1 (1967), 149–172 | DOI

[3] Radó T., “Aufgabe 41”, Jahresber. Deutsche Math.-Verein., 35 (1926), 49

[4] Kneser H., “Lösung der Aufgabe 41”, Jahresber. Deutsche Math.-Verein., 35 (1926), 123–124

[5] Choquet G., “Sur un type de transformation analytique généralisant la représentation conforme et définie au moyen de fonctions harmoniques”, Bull. Sci. Math., 69 (1945), 156–165 | Zbl

[6] Godunov S. K., Zabrodin A. V., Ivanov M. Ya., Kraiko A. N., Prokopov G. P., Chislennoe reshenie mnogomernykh zadach gazovoi dinamiki, Nauka, M., 1976, 400 pp. | Zbl

[7] Ivanenko S. A., Charakhchyan A. A., “Algoritm postroeniya krivolineinykh setok iz vypuklykh chetyrekhugolnikov”, Dokl. AN SSSR, 295:2 (1987), 280–283 | Zbl

[8] Ivanenko S. A., “Variatsionnye metody postroeniya adaptivnykh setok”, Zh. vychisl. matem. i matem. fiz., 43:6 (2003), 830–844 | Zbl

[9] Azarenok B. N., “Variatsionnyi metod postroeniya geksaedralnykh setok s upravlyayuschei metrikoi”, Matem. modelirovanie, 20:9 (2008), 3–22

[10] Ivanenko S. A., “O suschestvovanii uravnenii dlya opisaniya klassov nevyrozhdennykh krivolineinykh koordinat v proizvolnoi oblasti”, ZhVMiMF, 42:1 (2002), 47–52 | Zbl

[11] Prokopov G. P., “Universalnye variatsionnye funktsionaly dlya postroeniya dvumernykh setok”, Preprinty IPM, 2001, 001, 36 pp.

[12] Godunov S. K., Prokopov G. P., “O raschetakh konformnykh otobrazhenii i postroenii raznostnykh setok”, Zh. vychisl. matem. i matem. fiz., 7:5 (1967), 1031–1059 | Zbl

[13] Prokopov G. P., “Variatsionnye metody rascheta dvumernykh setok pri reshenii nestatsionarnykh zadach”, Preprinty IPM, 2003, 004, 32 pp.

[14] Prokopov G. P., “Realizatsiya variatsionnogo podkhoda k raschetu dvumernykh setok v nestatsionarnykh zadachakh”, Preprinty IPM, 2005, 116, 36 pp.

[15] Prokopov G. P., “Vybor parametrov pri variatsionnom podkhode k raschetu regulyarnykh setok”, Preprinty IPM, 2006, 014, 32 pp.

[16] Azarenok B. N., “O postroenii strukturirovannykh setok v dvumernykh nevypuklykh oblastyakh s pomoschyu otobrazhenii”, Zh. vychisl. matem. i matem. fiz., 49:5 (2009), 826–839 | Zbl

[17] Knupp P., Luczak R., “Truncation error in grid generation: a case study”, Numerical Methods for Partial Differential Equations, 11 (1995), 561–571 | DOI | Zbl

[18] Knupp P., Steinberg S., Fundamentals of Grid Generation, CRC Press, Boca Raton, 1993, 286 pp.

[19] Ivanenko S. A., Adaptivno-garmonicheskie setki, VTs RAN, M., 1997, 181 pp.

[20] Lavrentev M. A., Shabat B. V., Metody teorii funktsii kompleksnogo peremennogo, Nauka, M., 1987, 716 pp. | Zbl

[21] Koppenfels V., Shtalman F., Praktika konformnykh otobrazhenii, IL, M., 1963, 607 pp. | Zbl

[22] Yanke E., Emde F., Lesh F., Spetsialnye funktsii, Nauka, M., 1977, 342 pp.

[23] Beitman G., Erdeii A., Vysshie transtsendentnye funktsii. Gipergeometricheskaya funktsiya. Funktsiya Lezhandra, Nauka, M., 1965, 295 pp.

[24] Vecheslavov V. V., Kokoulin V. I., “Opredelenie parametrov konformnogo preobrazovaniya odnosvyaznykh mnogougolnykh oblastei”, ZhVMiMF, 13:4 (1973), 865–872 | Zbl

[25] Prokopov G. P., “Ob organizatsii sravneniya algoritmov i programm postroeniya regulyarnykh dvumernykh raznostnykh setok”, Voprosy atomnoi nauki i tekhniki. Ser. Matem. modelir. fiz. protsessov, 1989, no. 3, 98–108

[26] Bezrodnykh S. I., Vlasov V. I., “Ob odnoi probleme konstruktivnoi teorii garmonicheskikh otobrazhenii”, Sovremennaya matematika. Fundamentalnye napravleniya, 46 (2012), 5–30 | Zbl

[27] Ivanenko S. A., “Upravlenie formoi yacheek v protsesse postroeniya setok”, Zh. vychisl. matem. i matem. fiz., 40:1 (2000), 1662–1684 | Zbl

[28] Ivanenko S. A., “Postroenie nevyrozhdennykh setok”, ZhVMiMF, 28:10 (1988), 1498–1506

[29] Garanzha V. A., Kaporin I. E., “Regulyarizatsiya barernogo variatsionnogo metoda postroeniya raschetnykh setok”, Zh. vychisl. matem. i matem. fiz., 39:9 (1999), 1489–1503 | Zbl

[30] Charakhchyan A. A., “Ob odnom podkhode k umensheniyu vychislitelnykh zatrat pri postroenii krivolineinykh setok”, Zh. vychisl. matem. i matem. fiz., 38:2 (1998), 344–350 | Zbl