1,720,998 research outputs found
Spectral Matrix Conditioning in Cylindrical and Spherical Elliptic Equations
No full textIn the spectral solution of 3-D Poisson equations in cylindrical and spherical coordinates including the axis or the center, it is convenient to employ radial basis functions that depend on the Fourier wavenumber or on the latitudinal mode. This idea has been adopted by Matsushima and Marcus and by Verkley for planar problems and pursued by the present authors for spherical ones. For the Dirichlet boundary value problem in both geometries, original bases have been introduced built upon Jacobi polynomials which lead to a purely diagonal representation of the radial second-order differential operator of all spectral modes. This note details the origin of such a diagonalization which extends to cylindrical and spherical regions the properties of the Legendre basis introduced by Jie Shen for Cartesian domains. Closed form expressions are derived for the diagonal elements of the stiffness matrices as well as for the elements of the tridiagonal mass matrices occurring in evolutionary problems. Furthermore, the bound on the condition number of the spectral matrices associated with the Helmholtz equation are determined, proving in a rigorous way one of the main advantages of the proposed radial bases
On the Approximation of the Unsteady Navier-Stokes Equations by Finite Element Projection Methods
No full textNavier-Stokes Spectral Solver in a Sphere
No full textThe paper presents the first implementation of a primitive variable spectral method for calculating viscous flows inside a sphere. A variational formulation of the Navier–Stokes equations is adopted using a fractional-step time discretization with the classical second-order backward difference scheme combined with explicit extrapolation of the nonlinear term. The resulting scalar and vector elliptic equations are solved by means of the direct spectral solvers developed recently by the authors. The spectral matrices for radial operators are characterized by a minimal sparsity – diagonal stiffness and tridiagonal mass matrix. Closed-form expressions of their nonzero elements are provided here for the first time, showing that the condition number of the relevant matrices grows as the second power of the truncation order. A new spectral elliptic solver for the velocity unknown in spherical coordinates is also described that includes implicitly the Coriolis force in a rotating frame, but requires a minimal coupling between the modal velocity components in the Fourier space. The numerical tests confirm that the proposed method achieves spectral accuracy and ensures infinite differentiability to all orders of the numerical solution, by construction. These results indicate that the new primitive variable spectral solver is an effective alternative to the spectral method recently proposed by Kida and Nakayama, where the velocity field is represented in terms of poloidal and toroidal functions
Galerkin-Laguerre Spectral Solution of Self-Similar Boundary Layer Problems
Get PDFIn this work the Laguerre basis for the biharmonic equation introduced by Jie Shen is employed in the spectral solution of self-similar problems of the boundary layer theory. An original Petrov-Galerkin formulation of the Falkner-Skan equation is presented which is based on a judiciously chosen special basis function to capture the asymptotic behaviour of the unknown. A spectral method of remarkable simplicity is obtained for computing Falkner-Skan-Cooke boundary layer flows. The accuracy and efficiency of the Laguerre spectral approximation is illustrated by determining the linear stability of nonseparated and separated flows according to the Orr-Sommerfeld equation. The pentadiagonal matrices representing the derivative operators are explicitly provided in an Appendix to aid an immediate implementation of the spectral solution algorithms
Spectral Elliptic Solvers in a Finite Cylinder
No full textNew 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 newJacobi basis is proposed for the radial direction overcoming themain 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
On the spectral solution of the three-dimensional Navier-Stokes equations in spherical and cylindrical regions
No full text- …
