Spectral methods for hyperbolic initial boundary value problems on parallel computers

• Published in: Journal of computational and applied mathematics Vol. 134; no. 1; pp. 165 - 190
• Main Authors: Dutt, Pravir K., Bedekar, Smita
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.09.2001; Elsevier
• Subjects: Domain decomposition methods; Energy estimates on multidomains; Error estimates for the analytic case; Exact sciences and technology; Exponential accuracy; Galerkin-Collocation method; Mathematics; Numerical analysis; Numerical analysis. Scientific computation; p and h– p versions of FEM; Parallel pre-conditioners; Partial differential equations, boundary value problems; Partial differential equations, initial value problems and time-dependant initial-boundary value problems; Sciences and techniques of general use; Symmetric formulation; Well posedness of continuity boundary conditions; 65N35; Primary 65M06; 65MI5; Domain decomposition methods; Exponential accuracy; Well posedness of continuity boundary conditions; Parallel pre-conditioners; Energy estimates on multidomains; Error estimates for the analytic case; Symmetric formulation; p and h– p versions of FEM; 65M12; Galerkin-Collocation method; Energy analysis; Spectral method; Error estimation; Continuous; Initial value problem; Well posed problem; Parallel computer; Boundary condition; Finite element method; Accuracy; Boundary value problem; Hyperbolic equation; Galerkin method; Collocation method; Domain decomposition; Preconditioning; Symmetry
• ISSN: 0377-0427; 1879-1778
• DOI: 10.1016/S0377-0427(00)00535-5

In this paper a parallel spectral algorithm is developed for hyperbolic initial boundary value problems in one space dimension. The Galerkin-Collocation method, which is spectrally accurate in both space and time, is parallelized by using domain decomposition. This procedure leads to a minimization problem in which there is coupling at inter-domain boundaries. We construct a decoupled preconditioner which can be used to iteratively solve the minimization problem. Symmetric formulation of the problem, which is needed to compute the residual for the normal equations, is discussed. The methodology outlined for computing the normal equations applies equally well to computation of the residual for the p and h– p versions of the finite element method. There is, therefore, no need to compute the mass and stiffness matrices to obtain the residual, as is normally done. This leads to a great saving in time and memory particularly for solving nonlinear problems using the p and h– p versions of the finite element method. The method we discuss in this paper generalizes to hyperbolic initial boundary value problems in multidimensions too, provided the computational boundaries we have introduced are noncharacteristic and the system is symmetrizable. Finally, we show that for the case of analytic coefficients and data, satisfying all the required compatibility conditions so that the solution is analytic, the numerical solution is exponentially accurate in N, where N is proportional to the number of subdomains and the number of degrees of freedom in each element.