1,943 research outputs found
Spectral Methods for Numerical Relativity. The Initial Data Problem
Numerical relativity has traditionally been pursued via finite differencing.
Here we explore pseudospectral collocation (PSC) as an alternative to finite
differencing, focusing particularly on the solution of the Hamiltonian
constraint (an elliptic partial differential equation) for a black hole
spacetime with angular momentum and for a black hole spacetime superposed with
gravitational radiation. In PSC, an approximate solution, generally expressed
as a sum over a set of orthogonal basis functions (e.g., Chebyshev
polynomials), is substituted into the exact system of equations and the
residual minimized. For systems with analytic solutions the approximate
solutions converge upon the exact solution exponentially as the number of basis
functions is increased. Consequently, PSC has a high computational efficiency:
for solutions of even modest accuracy we find that PSC is substantially more
efficient, as measured by either execution time or memory required, than finite
differencing; furthermore, these savings increase rapidly with increasing
accuracy. The solution provided by PSC is an analytic function given
everywhere; consequently, no interpolation operators need to be defined to
determine the function values at intermediate points and no special
arrangements need to be made to evaluate the solution or its derivatives on the
boundaries. Since the practice of numerical relativity by finite differencing
has been, and continues to be, hampered by both high computational resource
demands and the difficulty of formulating acceptable finite difference
alternatives to the analytic boundary conditions, PSC should be further pursued
as an alternative way of formulating the computational problem of finding
numerical solutions to the field equations of general relativity.Comment: 15 pages, 5 figures, revtex, submitted to PR
The automatic solution of partial differential equations using a global spectral method
A spectral method for solving linear partial differential equations (PDEs)
with variable coefficients and general boundary conditions defined on
rectangular domains is described, based on separable representations of partial
differential operators and the one-dimensional ultraspherical spectral method.
If a partial differential operator is of splitting rank , such as the
operator associated with Poisson or Helmholtz, the corresponding PDE is solved
via a generalized Sylvester matrix equation, and a bivariate polynomial
approximation of the solution of degree is computed in
operations. Partial differential operators of
splitting rank are solved via a linear system involving a block-banded
matrix in operations. Numerical
examples demonstrate the applicability of our 2D spectral method to a broad
class of PDEs, which includes elliptic and dispersive time-evolution equations.
The resulting PDE solver is written in MATLAB and is publicly available as part
of CHEBFUN. It can resolve solutions requiring over a million degrees of
freedom in under seconds. An experimental implementation in the Julia
language can currently perform the same solve in seconds.Comment: 22 page
Spectral methods for partial differential equations
Origins of spectral methods, especially their relation to the Method of Weighted Residuals, are surveyed. Basic Fourier, Chebyshev, and Legendre spectral concepts are reviewed, and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic, and mixed type. Fluid dynamical applications are emphasized
Spectral methods for CFD
One of the objectives of these notes is to provide a basic introduction to spectral methods with a particular emphasis on applications to computational fluid dynamics. Another objective is to summarize some of the most important developments in spectral methods in the last two years. The fundamentals of spectral methods for simple problems will be covered in depth, and the essential elements of several fluid dynamical applications will be sketched
Black hole evolution by spectral methods
Current methods of evolving a spacetime containing one or more black holes
are plagued by instabilities that prohibit long-term evolution. Some of these
instabilities may be due to the numerical method used, traditionally finite
differencing. In this paper, we explore the use of a pseudospectral collocation
(PSC) method for the evolution of a spherically symmetric black hole spacetime
in one dimension using a hyperbolic formulation of Einstein's equations. We
demonstrate that our PSC method is able to evolve a spherically symmetric black
hole spacetime forever without enforcing constraints, even if we add dynamics
via a Klein-Gordon scalar field. We find that, in contrast to
finite-differencing methods, black hole excision is a trivial operation using
PSC applied to a hyperbolic formulation of Einstein's equations. We discuss the
extension of this method to three spatial dimensions.Comment: 20 pages, 17 figures, submitted to PR
A second-order continuity domain-decomposition technique based on integrated Chebyshev polynomials for two-dimensional elliptic problems
This paper presents a second-order continuity
non-overlapping domain decomposition (DD) technique for
numerically solving second-order elliptic problems in
two-dimensional space. The proposed DD technique uses integrated Chebyshev polynomials to represent the solution in subdomains. The constants of integration are utilized to impose continuity of the second-order normal derivative of the solution at the interior points of subdomain interfaces. To also achieve a C2 (C squared) function
at the intersection of interfaces, two additional unknowns are introduced at each intersection point. Numerical results show that the present DD method yields a higher level of accuracy than conventional DD techniques based on differentiated Chebyshev polynomials
Spectral methods in general relativistic astrophysics
We present spectral methods developed in our group to solve three-dimensional
partial differential equations. The emphasis is put on equations arising from
astrophysical problems in the framework of general relativity.Comment: 51 pages, elsart (Elsevier Preprint), 19 PostScript figures,
submitted to Journal of Computational & Applied Mathematic
- …