• Contents V46
• Online Suppl
• CTAC 2004
• Part 1
• Part 2
• Part 3
• Part 4
• Author index
• Journal home
• RSS News
• Information
• for Readers
• for Authors
• for Referees
• Other e-journals
• Aust Math Soc
• 
  

ANZIAM J. 46(E) ppC458--C471, 2005.

A modified finite volume method incorporating radial basis functions for simulating diffusion

T. J. Moroney

I. W. Turner

(Received 2 December 2004, revised 15 April 2005)

Abstract

The finite volume method is the favoured numerical technique for solving (possibly coupled, nonlinear, anisotropic) diffusion equations. The method transforms differential equations into a system of nonlinear, algebraic equations through the process of discretisation. The accuracy of this discretisation determines to a large extent the accuracy of the final solution. We present a new method of discretisation which is designed to achieve high accuracy without imposing excessive computational requirements. The method employs radial basis functions as a means of local gradient interpolation. When combined with high order Gaussian quadrature integration methods, the interpolation based on radial basis functions produces an efficient and accurate discretisation. The resulting nonlinear, algebraic system is solved efficiently using a Jacobian-free Newton--Krylov method. Information obtained from the Newton--Krylov iterations is used to construct an effective preconditioner in order to reduce the number of nonlinear iterations required to achieve an accurate solution. Results to date have been promising, with the method giving accuracy several orders of magnitude better than simpler methods based on shape functions.

Download to your computer

• Click here for the PDF article (269 kbytes) We suggest printing 2up.
• Click here for its BiBTeX record

Authors

T. J. Moroney

School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia. mailto:t.moroney@qut.edu.au

I. W. Turner

School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia. mailto:i.turner@qut.edu.au

Published June 9, 2005. ISSN 1446-8735

References

1. K. Burrage and J. Erhel. On the performance of various adaptive preconditioned gmres strategies. Numerical Linear Algebra with Applications, 5(2):101--121, 1998. http://citeseer.ist.psu.edu/362843.html
2. B. Carpentieri, I. S. Duff, and L. Giraud. Sparse pattern selection strategies for robust frobenius-norm minimization preconditioners in electromagnetism. Numerical linear algebra with applications, 7:(7--8), 2000. http://citeseer.csail.mit.edu/carpentieri00sparse.html
3. R. D. Cook. Concepts and applications of finite element analysis. Wiley, 2001.
4. J. F. Epperson. An Introduction to Numerical Methods and Analysis. Wiley, 2002.
5. B. Fornberg, T. A. Driscoll, G. Wright, and R. Charles. Observations on the behavior of radial basis function approximations near boundaries. Computers & Mathematics with Applications, 43(3-5):473--490, 2002. URL
6. P. Jayantha and I. W. Turner. A new higher order control-volume finite-element least-squares strategy for simulating transport in highly anisotropic media. Submitted to Journal of Computational Physics, 2002.
7. M. N. Ozi\c sik. Heat Conduction. Wiley, 1980.
8. R. Morgan. gmres with deflated restarting. SIAM J. Sci. Comput., 24(1):20--37, 2002. http://epubs.siam.org/sam-bin/dbq/article/36465
9. M. Powell. The theory of radial basis function approximations in 1990, in Advances in Numerical Analysis, Vol. II: Wavelets, Subdicision Algorithms and Radial Basis Functions, pages 105--210. Oxford University Press, 1990.