Versions of the collocation and least residuals method for solving problems of mathematical physics in the convex quadrangular domains

V. A. Belyaev ; V. P. Shapeev

Modelirovanie i analiz informacionnyh sistem, Tome 24 (2017) no. 5, pp. 629-648.

Résumé

The new versions of the collocations and least residuals (CLR) method of high-order accuracy are proposed and implemented for the numerical solution of the boundary value problems for PDE in the convex quadrangular domains. Their implementation and numerical experiments are performed by the examples of solving the biharmonic and Poisson equations. The solution of the biharmonic equation is used for simulation of the stress-strain state of an isotropic plate under the action of the transverse load. Differential problems are projected into the space of fourth-degree polynomials by the CLR method. The boundary conditions for the approximate solution are put down exactly on the boundary of the computational domain. The versions of the CLR method are implemented on the grids, which are constructed by two different ways. In the first version, a “quasiregular” grid is constructed in the domain, the extreme lines of this grid coincide with the boundaries of the domain. In the second version, the domain is initially covered by a regular grid with rectangular cells. Herewith, the collocation and matching points that are situated outside the domain are used for approximation of the differential equations in the boundary cells that had been crossed by the boundary. In addition the “small” irregular triangular cells that had been cut off by the domain boundary from rectangular cells of the initial regular grid are joined to adjacent quadrangular cells. This technique allowed to essentially reduce the conditionality of the system of linear algebraic equations of the approximate problem in comparison with the case when small irregular cells together with other cells were used as independent ones for constructing an approximate solution of the problem. It is shown that the approximate solution of problems converges with high order and matches with high accuracy with the analytical solution of the test problems in the case of the known solution in numerical experiments on the convergence of the solution of various problems on a sequence of grids.

Détail

Keywords: collocations and least residuals method, boundary value problem, non-canonical domain, irregular grid, high order approximation, biharmonic equation.
Mots-clés : Poisson’s equation

Bibliographie

[1] Sleptsov A. G., “Kollokatsionno-setochnoe resheniya ellepticheskih zadach”, Modelirovanie v mekhanike, 5(22):2 (1991), 101–126 (in Russian) | MR | Zbl

[2] Sleptsov A. G., Shokin Yu. I., “An adaptive grid-projection method for elliptic problems”, Comput. Math. Math. Phys., 37 (1997), 558–571 | MR | Zbl

[3] Belyaev V. V., Shapeev V. P., “Metod kollokaziy i naimenshikh kvadratov na adaptivnykh setkakh v oblasti s krivolineinoi granitsei”, Vychislitelnye tekhnologii, 5:4 (2000), 12–21 (in Russian) | MR

[4] Shapeev V. P., Belyaev V. A., “Varianty metoda kollokatsiy i naimenshikh nevyazok povushennoy tochnosti v oblasti s krivolineinoi granitsei”, Vychislitelnye tekhnologii, 21:5 (2016), 95–110 (in Russian)

[5] Golushko S. K., Idimeshev S. V., Shapeev V. P., “Metod kollokatsiy i naimenshikh nevyazok v prilozhenii k zadacham mekhaniki izotropnykh plastin”, Vychislitelnye tekhnologii, 18:6 (2013), 31–43 (in Russian)

[6] Semin L. G., Sleptsov A. G., Shapeev V. P., “Metod kollokatsiy-naimenshikh kvadratov dlya uravneniy Stoksa”, Vychislitelnye tekhnologii, 1:2 (1996), 90–98 (in Russian) | MR | Zbl

[7] Isaev V. I., Shapeev V. P., “Development of the collocations and least squares method”, Proc. Inst. Math. Mech., 261 (2008), 87–106 | MR

[8] Isaev V. I., Shapeev V. P., “High-accuracy versions of the collocations and least squares method for the numerical solution of the Navier–Stokes equations”, Comput. Math. and Math. Phys., 50 (2010), 1670–1681 | DOI | MR | Zbl

[9] Shapeev V. P., Vorozhtsov E. V., Isaev V. I., Idimeshev S. V., “Metod kollokatsiy i naimenshikh nevyazok dlya trekhmernykh uravneniy Navie–Stoksa”, Vychislitelnye metody i programmirovanie, 14 (2013), 306–322 (in Russian)

[10] Shapeev V., “Collocation and least residuals method and its applications”, EPJ Web of Conferences, 108 (2016), 01009 | DOI

[11] Vorozhtsov E. V., Shapeev V. P., “O kombinirovanii sposobov uskoreniya iteracionnyh processov pri chislennom reshenii uravneniy Navie–Stoksa”, Vychislitelnye metody i programmirovanie, 18 (2017), 80–102 (in Russian)

[12] Sleptsov A. G., “Ob uskorenii skhodimosti lineinykh iteratsiy II”, Modelirovanie v mekhanike, 3(20):5 (1989), 118–125 (in Russian) | MR | Zbl

[13] Saad Y., Numerical Methods for Large Eigenvalue Problems, Manchester University Press, Manchester, 1991 | MR

[14] Isaev V. I., Shapeev V. P., Eremin S. A., “Issledovanie svoistv metoda kollokatsii i naimenshikh kvadratov resheniya kraevykh zadach dlya uravneniya Puassona i uravneniy Navie–Stoksa”, Vychislitelnye tekhnologii, 12:3 (2007), 1–19 (in Russian)

[15] Timoshenko S. P. Woinowsky-Krieger S., Theory of Plates and Shells, 2dn edn., McGraw-Hill Book Company, 1959