Almost all currently available methods are based on iterative techniques using multigrid [7,23], domain decom-We present a direct, adaptive solver for the Poisson equation which can achieve any prescribed order of accuracy. It is based on position [11], or some other preconditioning strategy. Una do
A Fast 3D Poisson Solver of Arbitrary Order Accuracy
β Scribed by E. Braverman; M. Israeli; A. Averbuch; L. Vozovoi
- Publisher
- Elsevier Science
- Year
- 1998
- Tongue
- English
- Weight
- 482 KB
- Volume
- 144
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
β¦ Synopsis
We present a direct solver for the Poisson and Laplace equations in a 3D rectangular box. The method is based on the application of the discrete Fourier transform accompanied by a subtraction technique which allows reducing the errors associated with the Gibbs phenomenon and achieving any prescribed rate of convergence. The algorithm requires O(N 3 log N ) operations, where N is the number of grid points in each direction. We show that our approach allows accurate treatment of singular cases which arise when the boundary function is discontinuous or incompatible with the differential equation.
π SIMILAR VOLUMES
We present a fourth order numerical solution method for the singular Neumann boundary problem of Poisson equations. Such problems arise in the solution process of incompressible Navier-Stokes equations and in the time-harmonic wave propagation in the frequence space with the zero wavenumber. The equ
In many realistic fluid-dynamical simulations the specification of the boundary conditions, the error sources, and the number of time steps to reach a steady state are important practical considerations. In this paper we study these issues in the case of the lattice-BGK model. The objective is to pr
An implicit split-operator FFT algorithm for the numerical solution of the time-dependent Schrodinger equation is implemented for the electronic structure of αΈ¦q and and H . The covalent versus separated-atoms behavior is described by two 2 2 distinct steady states to which the imaginary-time Schrodi