Volume 5, Issue 2-4
Spectral Elliptic Solvers in a Finite Cylinder

F. Auteri & L. Quartapelle

DOI:

Commun. Comput. Phys., 5 (2009), pp. 426-441.

Published online: 2009-02

Preview Full PDF 176 1738
Export citation
  • Abstract

New direct spectral solvers for the 3D Helmholtz equation in a finite cylindrical region are presented. A purely variational (no collocation) formulation of the problem is adopted, based on Fourier series expansion of the angular dependence and Legendre polynomials for the axial dependence. A new Jacobi basis is proposed for the radial direction overcoming the main disadvantages of previously developed bases for the Dirichlet problem. Nonhomogeneous Dirichlet boundary conditions are enforced by a discrete lifting and the vector problem is solved by means of a classical uncoupling technique. In the considered formulation, boundary conditions on the axis of the cylindrical domain are never mentioned, by construction. The solution algorithms for the scalar equations are based on double diagonalization along the radial and axial directions. The spectral accuracy of the proposed algorithms is verified by numerical tests.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-5-426, author = {}, title = {Spectral Elliptic Solvers in a Finite Cylinder}, journal = {Communications in Computational Physics}, year = {2009}, volume = {5}, number = {2-4}, pages = {426--441}, abstract = {

New direct spectral solvers for the 3D Helmholtz equation in a finite cylindrical region are presented. A purely variational (no collocation) formulation of the problem is adopted, based on Fourier series expansion of the angular dependence and Legendre polynomials for the axial dependence. A new Jacobi basis is proposed for the radial direction overcoming the main disadvantages of previously developed bases for the Dirichlet problem. Nonhomogeneous Dirichlet boundary conditions are enforced by a discrete lifting and the vector problem is solved by means of a classical uncoupling technique. In the considered formulation, boundary conditions on the axis of the cylindrical domain are never mentioned, by construction. The solution algorithms for the scalar equations are based on double diagonalization along the radial and axial directions. The spectral accuracy of the proposed algorithms is verified by numerical tests.

}, issn = {1991-7120}, doi = {https://doi.org/}, url = {http://global-sci.org/intro/article_detail/cicp/7741.html} }
TY - JOUR T1 - Spectral Elliptic Solvers in a Finite Cylinder JO - Communications in Computational Physics VL - 2-4 SP - 426 EP - 441 PY - 2009 DA - 2009/02 SN - 5 DO - http://doi.org/ UR - https://global-sci.org/intro/article_detail/cicp/7741.html KW - AB -

New direct spectral solvers for the 3D Helmholtz equation in a finite cylindrical region are presented. A purely variational (no collocation) formulation of the problem is adopted, based on Fourier series expansion of the angular dependence and Legendre polynomials for the axial dependence. A new Jacobi basis is proposed for the radial direction overcoming the main disadvantages of previously developed bases for the Dirichlet problem. Nonhomogeneous Dirichlet boundary conditions are enforced by a discrete lifting and the vector problem is solved by means of a classical uncoupling technique. In the considered formulation, boundary conditions on the axis of the cylindrical domain are never mentioned, by construction. The solution algorithms for the scalar equations are based on double diagonalization along the radial and axial directions. The spectral accuracy of the proposed algorithms is verified by numerical tests.

F. Auteri & L. Quartapelle. (2020). Spectral Elliptic Solvers in a Finite Cylinder. Communications in Computational Physics. 5 (2-4). 426-441. doi:
Copy to clipboard
The citation has been copied to your clipboard