Fast Poisson, Fast Helmholtz and fast linear elastostatic solvers on rectangular parallelepipeds
FFT-based fast Poisson and fast Helmholtz solvers on rectangular parallelepipeds for periodic boundary conditions in one-, two and three space dimensions can also be used to solve Dirichlet and Neumann boundary value problems. For non-zero boundary conditions, this is the special, grid-aligned case of jump corrections used in the Explicit Jump Immersed Interface method. Fast elastostatic solvers for periodic boundary conditions in two and three dimensions can also be based on the FFT. From the periodic solvers we derive fast solvers for the new 'normal' boundary conditions and essential boundary conditions on rectangular parallelepipeds. The periodic case allows a simple proof of existence and uniqueness of the solutions to the discretization of normal boundary conditions. Numerical examples demonstrate the efficiency of the fast elastostatic solvers for non-periodic boundary conditions. More importantly, the fast solvers on rectangular parallelepipeds can be used together with the Immersed Interface Method to solve problems on non-rectangular domains with general boundary conditions. Details of this are reported in the preprint The Explicit Jump Immersed Interface Method for 2D Linear Elastostatics by the author.
- Research Organization:
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
- Sponsoring Organization:
- Computing Sciences Directorate
- DOE Contract Number:
- DE-AC02-05CH11231
- OSTI ID:
- 982430
- Report Number(s):
- LBNL-43565; TRN: US201014%%33
- Country of Publication:
- United States
- Language:
- English
Similar Records
A fast parallel Poisson solver on irregular domains applied to beam dynamics simulations
A fast parallel 3D Poisson solver with longitudinal periodic and transverse open boundary conditions for space-charge simulations