Discretization of boundary integral equations leads, in general, to fully populated complex valued non-Hermitian systems of equations. In this paper we consider the e cient solution of these boundary element systems by preconditioned iterative methods of Krylov subspace type. We devise preconditione
Iterative solution of systems of equations in the dual reciprocity boundary element method for the diffusion equation
✍ Scribed by V. Bulgakov; B. Šarler; G. Kuhn
- Publisher
- John Wiley and Sons
- Year
- 1998
- Tongue
- English
- Weight
- 233 KB
- Volume
- 43
- Category
- Article
- ISSN
- 0029-5981
No coin nor oath required. For personal study only.
✦ Synopsis
In this paper the diffusion equation is solved in two-dimensional geometry by the dual reciprocity boundary element method (DRBEM). It is structured by fully implicit discretization over time and by weighting with the fundamental solution of the Laplace equation. The resulting domain integral of the diffusive term is transformed into two boundary integrals by using Green's second identity, and the domain integral of the transience term is converted into a finite series of boundary integrals by using dual reciprocity interpolation based on scaled augmented thin plate spline global approximation functions. Straight line geometry and constant field shape functions for boundary discretization are employed. The described procedure results in systems of equations with fully populated unsymmetric matrices. In the case of solving large problems, the solution of these systems by direct methods may be very time consuming. The present study investigates the possibility of using iterative methods for solving these systems of equations. It was demonstrated that Krylov-type methods like CGS and GMRES with simple Jacobi preconditioning appeared to be efficient and robust with respect to the problem size and time step magnitude. This paper can be considered as a logical starting point for research of iterative solutions to DRBEM systems of equations.
📜 SIMILAR VOLUMES
In this paper we analyze convergence of basic iterative Jacobi and Gauss-Seidel type methods for solving linear systems which result from finite element or finite volume discretization of convection-diffusion equations on unstructured meshes. In general the resulting stiffness matrices are neither M
This paper describes the solution of a steady state natural convection problem in porous media by the dual reciprocity boundary element method (DRBEM). The boundary element method (BEM) for the coupled set of mass, momentum, and energy equations in two dimensions is structured by the fundamental sol
The Laplace transform is applied to remove the time-dependent variable in the di usion equation. For nonharmonic initial conditions this gives rise to a non-homogeneous modiÿed Helmholtz equation which we solve by the method of fundamental solutions. To do this a particular solution must be obtained