A Spectral Method in Time for Initial-Value Problems
- 1 Division of Fusion Plasma Physics, Association EURATOM-VR, Alfven Laboratory,School of Electrical Engineering, KTH Royal Institute of Technology, Stockholm, Sweden
Abstract
A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this generalized weighted residual method (GWRM). The approximate solutions obtained are thus analytical, finite order multivariate polynomials. The method avoids time step limitations. To determine the spectral coefficients, a system of algebraic equations is solved iteratively. A root solver, with excellent global convergence properties, has been developed. Accuracy and efficiency are controlled by the number of included Chebyshev modes and by use of temporal and spatial subdomains. As examples of advanced application, stability problems within ideal and resistive magnetohydrodynamics (MHD) are solved. To introduce the method, solutions to a stiff ordinary differential equation are demonstrated and discussed. Subsequently, the GWRM is applied to the Burger and forced wave equations. Comparisons with the explicit Lax-Wendroff and implicit Crank-Nicolson finite difference methods show that the method is accurate and efficient. Thus the method shows potential for advanced initial value problems in fluid mechanics and MHD.
- R. D. Richtmyer and K. W. Morton, ”Difference Methods for Initial-value Problems”, Krieger Publishing, 1994.
- D. S. Harned and W. Kerner, ”Semi-implicit method for three-dimensional compressible magnetohydrodynamic simulation”, J. Comp. Phys. Vol. 60, No. 1, 1985, pp. 62-75. doi:10.1016/0021-9991(85)90017-8
- D. S. Harned and D. D. Schnack, ”Semi-implicit method for long time scale magnetohydrodynamic computations in three dimensions”, J. Comp. Phys. Vol. 65, No. 1, 1986, pp. 57-70. doi:10.1016/0021-9991(86)90004-5
- J. Scheffel, ”Semi-analytical solution of initial-value problems”, TRITA-ALF-2004-03, KTH Royal Institute of Technology, Stockholm, Sweden, 2004.
- B. A. Finlayson and L. E. Scriven, ”The Method of Weighted Residuals – a Review”, Appl. Mech. Rev. Vol. 19, No. 9, 1966, pp. 735-748.
- C. A. J. Fletcher, ”Computational Techniques for Fluid Dynamics”, Vols 1-2, Springer, 2000.
- R. Peyret and T. D. Taylor, ”Computational Methods for Fluid Flow”, Springer, New York 1983. doi:10.1007/978-3-642-85952-6
- H. Tal-Ezer, ”Spectral Methods in Time for Hyperbolic Equations”, SIAM J. Numer. Anal. Vol. 23, No. 1, 1986, pp. 11-26. doi:10.1137/0723002
- H. Tal-Ezer, ”Spectral Methods in Time for Parabolic Problems”, SIAM J. Numer. Anal. Vol. 26, No. 1, 1989, pp. 1-11. doi:10.1137/0726001
- G. Ierley, B. Spencer and R. Worthing, ”Spectral Methods in Time for a Class of Parabolic Partial Differential Equations”, J. Comput. Phys. Vol. 102, No. 1, 1992, pp. 88-97. doi:10.1016/S0021-9991(05)80008-7
- J.-G. Tang and H.-P. Ma, ”Single and Multi-interval Legendre ?-methods in Time for Parabolic Equations”, Adv. in Comput. Math. Vol. 17, No. 4, 2002, pp. 349-367. doi:10.1023/A:1016273820035
- M. Dehghan and A. Taleei, ”Numerical Solution of Nonlinear Schr?dinger Equation by Using Time-Space Pseudo-Spectral Method”, Num. Meth. for Part. Diff. Eq. Vol. 26, No. 4, 2010, pp. 979-992.
- D. Gottlieb and S. A. Orszag, ”Numerical Analysis of Spectral Methods: Theory and Applications”, SIAM, Philadelphia, 1987.
- J. C. Mason and D. C. Handscomb, ”Chebyshev Polynomials”, Chapman & Hall, 2003.
- W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery, ”Numerical Recipes”, Cambridge University Press, 1992.