A splitting of a third-order partial differential equation into a first-order and a second-order one is proposed as the basis for a mixed finite element method to approximate its solution. A time-continuous numerical method is described and error estimates for its solution are demonstrated. Finally,
A family of fourth-order parallel splitting methods for parabolic partial differential equations
โ Scribed by M. S. A. Taj; E. H. Twizell
- Publisher
- John Wiley and Sons
- Year
- 1997
- Tongue
- English
- Weight
- 149 KB
- Volume
- 13
- Category
- Article
- ISSN
- 0749-159X
No coin nor oath required. For personal study only.
โฆ Synopsis
A family of numerical methods which are L-stable, fourth-order accurate in space and time, and do not require the use of complex arithmetic is developed for solving second-order linear parabolic partial differential equations. In these methods, second-order spatial derivatives are approximated by fourth-order finitedifference approximations, and the matrix exponential function is approximated by a rational approximation consisting of three parameters. Parallel algorithms are developed and tested on the one-dimensional heat equation, with constant coefficients, subject to homogeneous and time-dependent boundary conditions. These methods are also extended to two-and three-dimensional heat equations, with constant coefficients, subject to homogeneous boundary conditions.
๐ SIMILAR VOLUMES
Approximations where the derivatives are corrected so as to satisfy linear completeness on the derivatives are investigated. These techniques are used in particle methods and other mesh-free methods. The basic approximation is a Shepard interpolant which possesses only constant completeness. The der
A generalized and more complete methodology for treating boundary conditions in the Differential Quadrature Method (DQM) is presented. This improved approach eliminates the deficiencies of the -type grid arrangement, which represents an approximation, by applying the boundary conditions exactly. Two
Elliptic PDEs with variable coefficients in a domain with complex geometry occur in many ocean models. The parallelization of the elliptic solver by the Shur complement method is presented for the ice-ocean model BRIOS. The Schur complement method is usually employed as an iterative solver, but for