Hostname: page-component-78c5997874-xbtfd Total loading time: 0 Render date: 2024-11-10T19:28:46.087Z Has data issue: false hasContentIssue false
  • English
  • Français

Second Order Multigrid Methods for Elliptic Problems with Discontinuous Coefficients on an Arbitrary Interface, I: One Dimensional Problems

Published online by Cambridge University Press:  28 May 2015

Armando Coco  and
Giovanni Russo
Armando Coco
Affiliation:
Dipartimento di Matematica e Informatica, Università di Catania, Viale Andrea Doria, 6, 95125, Catania, Italy
Giovanni Russo*
Affiliation:
Dipartimento di Matematica e Informatica, Università di Catania, Viale Andrea Doria, 6, 95125, Catania, Italy
*
*Corresponding author.Email address:russo@dmi.unict.it
Get access

Abstract

In this paper we present a one dimensional second order accurate method to solve Elliptic equations with discontinuous coefficients on an arbitrary interface. Second order accuracy for the first derivative is obtained as well. The method is based on the Ghost Fluid Method, making use of ghost points on which the value is defined by suitable interface conditions. The multi-domain formulation is adopted, where the problem is split in two sub-problems and interface conditions will be enforced to close the problem. Interface conditions are relaxed together with the internal equations (following the approach proposed in [10] in the case of smooth coefficients), leading to an iterative method on all the set of grid values (inside points and ghost points). A multigrid approach with a suitable definition of the restriction operator is provided. The restriction of the defect is performed separately for both sub-problems, providing a convergence factor close to the one measured in the case of smooth coefficient and independent on the magnitude of the jump in the coefficient. Numerical tests will confirm the second order accuracy. Although the method is proposed in one dimension, the extension in higher dimension is currently underway [12] and it will be carried out by combining the discretization of [10] with the multigrid approach of [11] for Elliptic problems with non-eliminated boundary conditions in arbitrary domain.

Keywords

35J2565N0665N55Elliptic equationdiscontinuous coefficientsecond order accuracymultigridarbitrary interfacejump conditions
Type
Research Article
Information
Numerical Mathematics: Theory, Methods and Applications , Volume 5 , Issue 1 , February 2012 , pp. 19 - 42
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1] Alcouffe, R. E., Brandt, A., Dendy, J. J. E., and Painter, J. W. The multigrid method for the diffusion equation with strongly discontinuous coefficients. Journal on Scientific and Statistical Computing, 2: 430–454, 1981.CrossRefGoogle Scholar
[2] Angot, P., Bruneau, C.-H., and Fabrie, P. A penalization method to take into account obstacles in incompressible viscous flows. Numer. Math., 81, 1999.CrossRefGoogle Scholar
[3] Babuška, I. The finite element method for elliptic equations with discontinuous coefficients. Computing, 5:207–213, 1970.CrossRefGoogle Scholar
[4] Bao, G., Wei, G., and Zhao, S. Numerical solution of the Helmholtz equation with high wave numbers. J. Numer. Methods Engng., 59:389–408, 2004.CrossRefGoogle Scholar
[5] Bramble, J. and King, J. A finite element method for interface problems in domains with smooth boundaries and interfaces. Adv. Comput. Math., 6:109–138, 1996.CrossRefGoogle Scholar
[6] Bramble, J. H. and Hubbard, B. E. Approximation of solutions of mixed boundary value problems for Poisson’s equation by finite differences. J. Assoc. Comput. Mach., 12:114–123, 1965.CrossRefGoogle Scholar
[7] Briggs, W. L., V Henson, E., and McCormick, S. F. A Multigrid Tutorial. SIAM, 2000.CrossRefGoogle Scholar
[8] Chantalat, F., Bruneau, C.-H., Galusinski, C., and Iollo, A. Level-set, penalization and cartesian meshes: A paradigm for inverse problems and optimal design. Journal of Computational Physics, 228:6291–6315, 2009.CrossRefGoogle Scholar
[9] Chen, H., Min, C., and Gibou, F. A supra-convergent finite difference scheme for the Poisson and heat equations on irregular domains and non-graded adaptive Cartesian grids. Journal of Scientific Computing, 31:19–60, 2007.CrossRefGoogle Scholar
[10] Coco, A. and Russo, G. A fictitious time method for the solution of Poisson equation in an arbitrary domain embedded in a square grid. Journal of Computation Physics. Under revision.Google Scholar
[11] Coco, A. and Russo, G. Multigrid approach for Poisson’s equation with mixed boundary condition in an arbitrary domain.Google Scholar
[12] Coco, A. and Russo, G. Second order multigrid methods for elliptic problems with discontinuous coefficients on an arbitrary interface, II: higher dimensional problems.Google Scholar
[13] Donea, J. An arbitrary Lagrangian-Eulerian finite element method for transient fluid-structure interactions. Computer Methods in Applied Mechanics and Engineering, 33:689–723, 1982.CrossRefGoogle Scholar
[14] Fedkiw, R., Aslam, T, Merriman, B., and Osher, S. A Non-Oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (The Ghost Fluid Method). Journal of Computational Physics, 152:457–492, 1999.CrossRefGoogle Scholar
[15] Formaggia, L. and Nobile, F. Stability analysis of second-order time accurate schemes for ALE-FEM. Computer Methods in Applied Mechanics and Engineering, 193:4097–4116, 2004.CrossRefGoogle Scholar
[16] Gibou, F. and Fedkiw, R. A second-order-accurate symmetric discratization of the poisson equation on irregular domains. Journal of Computational Physics, 176:205–227, 2002.CrossRefGoogle Scholar
[17] Gibou, F. and Fedkiw, R. A fourth order accurate discretization for the laplace and heat equations on arbitary domains, with applications to the stefan problem. Journal of Computational Physics, 202:577–601, 2005.CrossRefGoogle Scholar
[18] Glowinski, R., Pan, T. W., Hesla, T. I., and Joseph, D. D. A distributed Lagrange multiplier/fictitious domain method for particulate flows. International Journal of Multiphase Flow, 25:755–794, 1999.CrossRefGoogle Scholar
[19] Hackbusch, W. Multi-grid methods and applications. Springer, 1985.CrossRefGoogle Scholar
[20] Hackbusch, W. Elliptic Differential Equations: Theory and Numerical Treatment. Springer, 2003.Google Scholar
[21] Dendy, J. J. E. Black Box Multigrid. Journal of Computational Physics, 48:366–386, 1982.CrossRefGoogle Scholar
[22] Johansen, H. and Colella, P. A Cartesian Grid Embedded Boundary Method for Poisson Equation on Irregular Domains. Journal of Computational Physics, 147:60–85, 1998.CrossRefGoogle Scholar
[23] LeVeque, R. and Li, Z. The immersed interface method for elliptic equations with discontinuous coefficients and singular sources. SIAM J. Numer. Anal., 31:1019–1044, 1994.CrossRefGoogle Scholar
[24] Mayo, A. The fast solution of Poisson’s and the biharmonic equations on irregular regions. SIAM J. Numer. Anal., 21:285–299, 1984.CrossRefGoogle Scholar
[25] Papac, J., Gibou, F., and Ratsch, C. Efficient Symmetric Discretization for the Poisson, Heat and Stefan-Type Problems with Robin Boundary Conditions. Journal of Computational Physics, 229:875–889, 2010.CrossRefGoogle Scholar
[26] Peskin, C. S. Numerical analysis of blood flow in the heart. Journal of Computational Physics, 25:220–252, 1977.CrossRefGoogle Scholar
[27] Quarteroni, A. and Valli, A. Domain Decomposition Methods for Partial Differential Equations. Numerical Mathematics and Scientific Computation, 1999.CrossRefGoogle Scholar
[28] Ruge, J. W. and Stüben. Algebraic multigrid. Multigrid methods. [29] Sarthou, A., Vincent, S., Caltagirone, J., and Angot, P. Eulerian-Lagrangian grid coupling and penalty methods for the simulation of multiphase flows interacting with complex objects. International Journal for Numerical Methods in Fluids, 00:1–6, 2007.Google Scholar
[30] Shortley, G. H. and Weller, R. The numerical solution of laplace’s equation. J. Appl. Phys., 9:334–348, 1938.CrossRefGoogle Scholar
[31] U.Trottemberg, Oosterlee, C., and Schuller, A. Multigrid. Academic Press, 2000.Google Scholar
[32] Yu, S., Zhou, Y., and Wei, G. Matched Interface and Boundary (MIB) method for elliptic problems with sharp-edged interfaces. Journal of Computational Physics, 224:729–756, 2007.CrossRefGoogle Scholar