×

Multigrid solution of high order discretisation for three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind. (English) Zbl 1102.65125

Summary: We propose two compact finite difference approximations for three-dimensional biharmonic equation with Dirichlet boundary conditions of the second kind. In these methods there is no need to define special formulas near the boundaries and boundary conditions are incorporated with these techniques. The unknown solution and its second derivatives are carried as unknowns at grid points. We derive second-order and fourth-order approximations on a 27 point compact stencil. Classical iteration methods such as Gauss-Seidel and successive overrelaxation for solving the linear system arising from the second-order and fourth-order discretisation suffer from slow convergence. In order to overcome this problem we use the multigrid method which exhibit grid-independent convergence and solve the linear system of equations in small amount of computer time. The fourth-order finite difference approximations are used to solve several test problems and produce high accurate numerical solutions.

MSC:

65N55 Multigrid methods; domain decomposition for boundary value problems involving PDEs
65N06 Finite difference methods for boundary value problems involving PDEs
35J40 Boundary value problems for higher-order elliptic equations
65N12 Stability and convergence of numerical methods for boundary value problems involving PDEs
65F10 Iterative numerical methods for linear systems

Software:

Wesseling
Full Text: DOI

References:

[1] Altas, I.; Dym, J.; Gupta, M. M.; Manohar, R., Multigrid solution of automatically generated high order discretisation for the biharmonic equation, SIAM J. Sci. Comput., 19, 1575-1585 (1998) · Zbl 0918.65079
[2] Altas, I.; Erhel, J.; Gupta, M. M., High accuracy solution of three dimensoinal biharmonic equations, Numer. Alg., 29, 1-19 (2002) · Zbl 0992.65115
[3] Bakhvalov, N. S., On the convergence of a relaxation method with natural constraints on the elliptic operator, USSR Comp. Math. Math. Phys., 6, 101-135 (1966) · Zbl 0154.41002
[4] Bauer, L.; Riess, E. L., Block five diagonal matrices and the fast numerical solution of the biharmonic equation, Math. Comput., 26, 311-326 (1972) · Zbl 0257.65034
[5] Brandt, A., Multi-level adaptive technique (MLAT) for fast numerical solution of boundary value problems, (Proceedings of the Third International Conference on Numerical Methods in Fluid Mechanics, Paris, 1972. Proceedings of the Third International Conference on Numerical Methods in Fluid Mechanics, Paris, 1972, Lecture Notes in Physics, vol. 18 (1973), Springer-Verlag), 82-89 · Zbl 0259.76013
[6] A. Brandt, Multigrid techniques: 1984 guide with applications to fluid dynamics, in: GMD-Studien Nr.85, Gesselschaft für Mathematik und Datenverarbeitung MBH, St. Augustin, 1984.; A. Brandt, Multigrid techniques: 1984 guide with applications to fluid dynamics, in: GMD-Studien Nr.85, Gesselschaft für Mathematik und Datenverarbeitung MBH, St. Augustin, 1984. · Zbl 0581.76033
[7] Brandt, A., Multi-level adaptive solutions to boundary value problems, Math. Comput., 31, 333-390 (1977) · Zbl 0373.65054
[8] Briggs, W., A Multigrid Tutorial (1984), SIAM: SIAM Philadelphia
[9] Conte, S. D.; Dames, R. T., On an alternating direction method for solving the plate problem with mixed boundary conditions, J. Assoc. Comput. Mach., 7, 264-273 (1960) · Zbl 0098.09605
[10] de Vries, R. W.; Zandbergen, P. J., The numerical solution of the biharmonic equation using a multigrid method, (Bellhaus, W. F.; Hussaini, M. Y., Advances in Fluid Dynamics (1989), Springer-Verlag: Springer-Verlag New York), 25-35
[11] Fedorenko, R. P., A relaxation method for solving elliptic difference equations, USSR Comput. Math. Math. Phys., 1, 5, 1092-1096 (1962) · Zbl 0163.39303
[12] Fedorenko, R. P., The speed of convergence of one iterative process, USSR Comp. Math. Math. Phys., 4, 3, 227-235 (1964) · Zbl 0148.39501
[13] Gupta, M. M.; Manohar, R., Direct solution of biharmonic equation using noncoupled approach, J. Comput. Phys., 33, 236-248 (1979) · Zbl 0426.65058
[14] Gupta, M. M., Discretization error estimates for certain splitting procedures for solving first biharmonic boundary value problem, SIAM J. Numer. Anal., 12, 364-377 (1975) · Zbl 0272.65094
[15] Gupta, M. M.; Ehrlich, L. W., Some difference schemes for the biharmonic equation, SIAM J. Numer. Anal., 12, 773-790 (1975) · Zbl 0331.65061
[16] Kelmanson, M. A., Modified integral equation solution of viscous flows near sharp corners, Comput. Fluids, 11, 307-314 (1983) · Zbl 0559.76032
[17] Meurant, G., Computer Solution of Large Linear Systems (1999), North-Holland · Zbl 0934.65032
[18] Stephenson, J. W., Single cell discretisations of order two and four for biharmonic problems, J. Comput. Phys., 55, 65-80 (1984) · Zbl 0542.65051
[19] Trottenberg, U.; Oosterlee, C. W.; Schüller, A., Multigrid (2001), Academic Press · Zbl 0976.65106
[20] Wesseling, P., An Introduction to Multigrid Methods (1992), John Wiley & Sons: John Wiley & Sons Chichester · Zbl 0760.65092
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.