A Galerkin finite element formulation of diffusion processes based on a diagonal capacity matrix is analysed from the standpoint of local stability and convergence. The theoretical analysis assumes that the conductance matrix is locally diagonally dominant, and it is shown that one can always construct a finite element network of linear triangles satisfying this condition. Time derivatives are replaced by finite differences, leading to a mixed explicit‐implicit system of algebraic equations which can be efficiently solved by a point iterative technique. In this work the accelerated point iterative method is adopted and is shown to converge when the conductance matrix is locally diagonally dominant. Several examples are included in Part II of this paper to demonstrate the efficiency of the new approach.
|Original language||English (US)|
|Number of pages||15|
|Journal||International Journal for Numerical Methods in Engineering|
|State||Published - 1977|
ASJC Scopus subject areas
- Numerical Analysis
- Applied Mathematics