The hp-version of the least-squares collocation method with integral collocation for solving a biharmonic equation

Cover Page


Cite item

Full Text

Abstract

A new algorithm for the numerical solution of the biharmonic equation is developed. It is based on the first implemented hp-version of the least-squares collocation method (hp-LSCM) with integral collocations for a fourth-order elliptic equation in combination with modern methods of accelerating iterative processes for solving systems of linear algebraic equations (SLAE). The hp-LSCM provides the possibilities to refine the grid (h-version) and increase the degree of polynomials to the arbitrary order (p-approach). The convergence of approximate solutions obtained by the implemented version of the method is analyzed using an example of a numerical simulation of the bending of a hinged isotropic plate. The high accuracy and the increased order of convergence using polynomials up to the tenth order in the hp-LSCM are shown.

The effectiveness of the combined application of algorithms for accelerating iterative processes to solve SLAE that are combined with LSCM is investigated. In this paper, we consider the application of the following algorithms: preconditioning of SLAE matrices; the iteration acceleration algorithm based on Krylov subspaces; the prolongation operation on a multigrid complex; parallelization using OpenMP; a modified algorithm for solving local SLAEs. The latter is implemented with iterations over subdomains (which are cells) and makes it possible to more effectively solve overdetermined SLAEs in the LSCM in the case of solving a linear differential equation. The form of the matrices does not change at each iteration. Only the elements of the vectors of their right parts corresponding to the matching conditions are modified. The calculation time on a personal computer is reduced by more than 350 times with the combined use of all acceleration techniques compared to the case when only preconditioning was used.

About the authors

Vasily P. Shapeev

Khristianovich Institute of Theoretical and Applied Mechanics, Siberian Branch of the Russian Academy of Sciences;
Novosibirsk State University

Email: shapeev.vasily@mail.ru
ORCID iD: 0000-0001-6761-7273
SPIN-code: 7128-5536
Scopus Author ID: 8948291400
ResearcherId: J-7816-2018
http://www.mathnet.ru/person48027

Dr. Phys. & Math. Sci., Professor; Chief Researcher; Lab. of Thermomechanics and Strength of New Materials1; Professor; Dept. of Mathematical Modeling2

4/1, Institutskaya st., Novosibirsk, 630090, Russian Federation; 1, Pirogova st., Novosibirsk, 630090, Russian Federation

Luka S. Bryndin

Khristianovich Institute of Theoretical and Applied Mechanics, Siberian Branch of the Russian Academy of Sciences;
Novosibirsk State University

Email: l.bryndin@g.nsu.ru
ORCID iD: 0000-0002-0211-5800
SPIN-code: 4155-0327
Scopus Author ID: 57209026180
ResearcherId: U-5845-2018
http://www.mathnet.ru/person151108

Junior Research Scientist; Lab. of Thermomechanics and Strength of New Materials1; Postgraduate Student; Dept. of Mathematical Modeling2

4/1, Institutskaya st., Novosibirsk, 630090, Russian Federation; 1, Pirogova st., Novosibirsk, 630090, Russian Federation

Vasily A. Belyaev

Khristianovich Institute of Theoretical and Applied Mechanics, Siberian Branch of the Russian Academy of Sciences

Author for correspondence.
Email: belyaevasily@mail.ru
ORCID iD: 0000-0001-5901-2880
SPIN-code: 3276-0460
Scopus Author ID: 57197726181
ResearcherId: L-2637-2016
http://www.mathnet.ru/person132321

Junior Research Scientist; Lab. of Thermomechanics and Strength of New Materials1

Russian Federation, 4/1, Institutskaya st., Novosibirsk, 630090, Russian Federation

References

  1. Timoshenko S., Woinowsky–Krieger S. Theory of Plates and Shells. New York, St. Louis, McGraw-Hill Book Comp., 1959, xiv+580 pp.
  2. Reddy J. N. Mechanics of Laminated Composite Plates and Shells. Theory and Analysis. Boca Raton, CRC Press, 2011, xxiv+832 pp. DOI: https://doi.org/10.1201/b12409.
  3. Chen G., Li Z., Lin P. A fast finite difference method for biharmonic equations on irregular domains and its application to an incompressible Stokes flow, Adv. Comput. Math., 2008, vol. 29, no. 2, pp. 113–133. DOI: https://doi.org/10.1007/s10444-007-9043-6.
  4. McLaurin G. W. A general coupled equation approach for solving the biharmonic boundary value problem, SIAM J. Numer. Anal., 1974, vol. 11, no. 1, pp. 14–33. DOI: https://doi.org/10.1137/0711003.
  5. Mai-Duy N., See H., Tran-Cong T. A spectral collocation technique based on integrated Chebyshev polynomials for biharmonic problems in irregular domains, Appl. Math. Model., 2009, vol. 33, no. 1, pp. 284–299. DOI: https://doi.org/10.1016/j.apm.2007.11.002.
  6. Shao W., Wu X., Chen S. Chebyshev tau meshless method based on the integration-differentiation for Biharmonic-type equations on irregular domain, Engin. Anal. Bound. Elem., 2012, vol. 36, no. 12, pp. 1787–1798. DOI: https://doi.org/10.1016/j.enganabound.2012.06.005.
  7. Shao W., Wu X. An effective Chebyshev tau meshless domain decomposition method based on the integration-differentiation for solving fourth order equations, Appl. Math. Model., 2015, vol. 39, no. 9, pp. 2554–2569. DOI: https://doi.org/10.1016/j.apm.2014.10.048.
  8. Guo H., Zhang Z., Zou Q. A C 0 linear finite element method for biharmonic problems, J. Sci. Comput., 2018, vol. 74, no. 3, pp. 1397–1422. DOI: https://doi.org/10.1007/s10915-017-0501-0.
  9. Golushko S. K., Idimeshev S. V., Shapeev V. P. Application of collocations and least residuals method to problems of the isotropic plates theory, Vychisl. Tekhnol., 2013, vol. 18, no. 6, pp. 31–43 (In Russian). EDN: RUDTZF.
  10. Belyaev V. A., Bryndin L. S., Shapeev V. P., Golushko S. K. The least squares collocation method with the integral form of collocation equations for bending analysis of isotropic and orthotropic thin plates, AIP Conf. Proc., 2020, vol. 2288, 030084. EDN: FWKXNV. DOI: https://doi.org/10.1063/5.0028298.
  11. Belyaev V. A., Bryndin L. S., Golushko S. K., Semisalov B. V., Shapeev V. P. h-, p-, and hp-versions of the least-squares collocation method for solving boundary value problems for biharmonic equation in irregular domains and their applications, Comput. Math. Math. Phys., 2022, vol. 62, no. 4, pp. 517–537. EDN: VXTBZN. DOI: https://doi.org/10.1134/S0965542522040029.
  12. Babuška I., Guo B. Q. The h, p and h-p version of the finite element method; basis theory and applications, Adv. Eng. Softw., 1992, vol. 15, no. 3–4, pp. 159–174. DOI: https://doi.org/10.1016/0965-9978(92)90097-Y.
  13. Ramšak M., Škerget L. A subdomain boundary element method for high-Reynolds laminar flow using stream function-vorticity formulation, Int. J. Numer. Meth. Fluids, 2004, vol. 46, no. 8, pp. 815–847. DOI: https://doi.org/10.1002/fld.776.
  14. Saad Y. Numerical Methods for Large Eigenvalue Problems. Philadelphia, SIAM, 2011, xvi+276 pp. DOI: https://doi.org/10.1137/1.9781611970739.
  15. Fedorenko R. P. Vvedenie v vychislitel’nuiu fiziku [Introduction to Computational Physics]. Moscow, Moscow Institute of Physics and Technology, 1994, 528 pp. (In Russian). EDN: QJUAEP.
  16. Vorozhtsov E. V., Shapeev V. P. A divergence-free method of collocations and least squares for the computation of incompressible fluid flows and its efficient implementation, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2020, vol. 24, no. 3, pp. 542–573 (In Russian). EDN: OSPEBL. DOI: https://doi.org/10.14498/vsgtu1758.
  17. Belyaev V. A. On the effective implementation and capabilities of the least-squares collocation method for solving second-order elliptic equations, Vychisl. Meth. Prog. [Num. Meth. Prog.], 2021, vol. 22, no. 3, pp. 211–229 (In Russian). EDN: VBQEPE. DOI: https://doi.org/10.26089/NumMet.v22r313.
  18. Ortega J. M. Introduction to Parallel and Vector Solution of Linear Systems. New York, Springer, 1988, xi+305 pp. DOI: https://doi.org/10.1007/978-1-4899-2112-3.
  19. Karamyshev V. B., Kovenya V. M., Sleptsov A. G., Cherny S. G. Variational method of accelerating linear iterations and its applications, Comput. Fluids., 1996, vol. 25, no. 5, pp. 467–484. EDN: LDTBIH. DOI: https://doi.org/10.1016/0045-7930(96)00011-4.
  20. Il’in V. P., Knysh D. V. Parallel decomposition methods in the trace spaces, Vychisl. Meth. Prog. [Num. Meth. Prog.], 2011, vol. 12, no. 1, pp. 110–119 (In Russian). EDN: OJAYQV.
  21. Tang J. M., Nabben R., Vuik C., Erlangga Y. A. Comparison of two-level preconditioners derived from deflation, domain decomposition and multigrid methods, J. Sci. Comput., 2009, vol. 39, no. 3, pp. 340–370. DOI: https://doi.org/10.1007/s10915-009-9272-6.
  22. Godunov S. K., Ryabenkii V. S Difference Schemes: An Introduction to the Underlying Theory, Studies in Mathematics and its Applications, vol. 19. Amsterdam, North-Holland, 1987, xvii+489 pp.
  23. Demmel J. W. Applied Numerical Linear Algebra. Philadelphia, SIAM, 1997, xi+416 pp. DOI: https://doi.org/10.1137/1.9781611971446.
  24. Sleptsov A. G. Grid-projection solution of an elliptic problem for an irregular grid, Russ. J. Numer. Anal. Math. Model., 1993, vol. 8, no. 6, pp. 519–543. DOI: https://doi.org/10.1515/rnam.1993.8.6.519.
  25. 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. Math. Phys., 2010, vol. 50, no. 10, pp. 1670–1681. EDN: OHNDSJ. DOI: https://doi.org/10.1134/S0965542510100040.
  26. Drozdov G. M., Shapeev V. P. CAS application to the construction of high-order difference schemes for solving Poisson equation, In: Computer Algebra in Scientific Computing, Lecture Notes in Computer Science, 8660; eds. V. P. Gerdt, W. Koepf, W. M. Seiler, E. V. Vorozhtsov. Cham, Springer, 2014, pp. 99–110. EDN: UEYEZH. DOI: https://doi.org/10.1007/978-3-319-10515-4_8.
  27. Golushko S. K., Nemirovskii Yu. V. Priamye i obratnye zadachi mekhaniki uprugikh kompozitnykh plastin i obolochek vrashcheniia [Direct and Inverse Problems in the Mechanics of Composite Plates and Shells]. Moscow, Fizmatlit, 2008, 432 pp. (In Russian). EDN: MUWRZD
  28. Idimeshev S. V. Modified method of collocations and least residuals and its application in the mechanics of multilayer composite beams and plates, Thesis of Dissertation (Cand. Phys. & Math. Sci.). Novosibirsk, 2016, 179 pp. (In Russian). EDN: VPHODY.

Supplementary files

Supplementary Files
Action
1. JATS XML

Copyright (c) 2022 Authors; Samara State Technical University (Compilation, Design, and Layout)

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies