×

The stability and accuracy of EPIC algorithms. (English) Zbl 0667.76156

In the absence of particle discretisation effects (quadrature errors) EPIC algorithms give accurate and unconditionally stable schemes for integrating kinematic fluid and MHD equations. If particle discretisation is introduced as a consistently applied trapezium rule quadrature, the unconditional stability is retained. However, it is shown that not all quadratures lead to stable schemes. General results on linear stability are given. An analysis of dispersive and dissipative effects of finite particle number is presented, and the efficacy of the algorithm is illustrated by kinematics and dynamic test problems. The kinematic test with discrete particles showed unexpected stability and accuracy for large Courant numbers.

MSC:

76W05 Magnetohydrodynamics and electrohydrodynamics
Full Text: DOI

References:

[1] Eastwood, J. W., Particle Methods, (Winkler, K.-H.; Norman, M. L., Astrophysical Radiation Hydrodynamics (1986), Reidel: Reidel Dordrecht) · Zbl 0890.00021
[2] Eastwood, J. W.; Arter, W., (Morton, K. W.; Baines, M. J., Numerical Methods for Fluid Dynamics (1986), Oxford Univ. Press: Oxford Univ. Press Oxford), 581
[3] Eastwood, J. W., Comput. Phys. Commun., 43, 89 (1986)
[4] Hockney, R. W.; Eastwood, J. W., Computer Simulation using Particles (1981), McGraw-Hill: McGraw-Hill New York · Zbl 0662.76002
[5] Roache, P. J., Computational Fluid Dynamics, ((1976), Hermosa: Hermosa Albuquerque), 242
[6] Eastwood, J. W.; Arter, W., Phys. Rev. Lett., 57, 2528 (1986)
[7] Arter, W.; Eastwood, J. W., The Effect of Aliasing Error upon Numerical Solutions of the Hydrodynamic and Magnetohydrodynamic Equations, Transp. Th. Stat. Phys. (1987), to appear · Zbl 0622.76058
[8] Eastwood, J. W.; Arter, W., Culham Preprint CLM-P-782, IMA J. Numer. Anal (1986), to appear
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.