5,090 research outputs found
Collocation methods for a class of second order initial value problems with oscillatory solutions
We derive and analyse two families of multistep collocation methods for periodic initial-value problems of the form y" = f(x, y); y((^x)o) = yo, y(^1)(xo) = zo involving ordinary differential equations of second order in which the first derivative does not appear explicitly. A survey of recent results and proposed numerical methods is given in chapter 2. Chapter 3 is devoted to the analysis of a family of implicit Chebyshev methods proposed by Panovsky k Richardson. We show that for each non-negative integer r, there are two methods of order 2r from this family which possess non-vanishing intervals of periodicity. The equivalence of these methods with one-step collocation methods is also established, and these methods are shown to be neither P-stable nor symplectic. In chapters 4 and 5, two families of multistep collocation methods are derived, and their order and stability properties are investigated. A detailed analysis of the two-step symmetric methods from each class is also given. The multistep Runge-Kutta-Nystrom methods of chapter 4 are found to be difficult to analyse, and the specific examples considered are found to perform poorly in the areas of both accuracy and stability. By contrast, the two-step symmetric hybrid methods of chapter 5 are shown to have excellent stability properties, in particular we show that all two-step 27V-point methods of this type possess non-vanishing intervals of periodicity, and we give conditions under which these methods are almost P-stable. P-stable and efficient methods from this family are obtained and demonstrated in numerical experiments. A simple, cheap and effective error estimator for these methods is also given
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
A multigrid perspective on the parallel full approximation scheme in space and time
For the numerical solution of time-dependent partial differential equations,
time-parallel methods have recently shown to provide a promising way to extend
prevailing strong-scaling limits of numerical codes. One of the most complex
methods in this field is the "Parallel Full Approximation Scheme in Space and
Time" (PFASST). PFASST already shows promising results for many use cases and
many more is work in progress. However, a solid and reliable mathematical
foundation is still missing. We show that under certain assumptions the PFASST
algorithm can be conveniently and rigorously described as a multigrid-in-time
method. Following this equivalence, first steps towards a comprehensive
analysis of PFASST using block-wise local Fourier analysis are taken. The
theoretical results are applied to examples of diffusive and advective type
A fractional B-spline collocation method for the numerical solution of fractional predator-prey models
We present a collocation method based on fractional B-splines for the solution of fractional differential problems. The key-idea is to use the space generated by the fractional B-splines, i.e., piecewise polynomials of noninteger degree, as approximating space. Then, in the collocation step the fractional derivative of the approximating function is approximated accurately and efficiently by an exact differentiation rule that involves the generalized finite difference operator. To show the effectiveness of the method for the solution of nonlinear dynamical systems of fractional order, we solved the fractional Lotka-Volterra model and a fractional predator-pray model with variable coefficients. The numerical tests show that the method we proposed is accurate while keeping a low computational cost
Numerical Solution of ODEs and the Columbus' Egg: Three Simple Ideas for Three Difficult Problems
On computers, discrete problems are solved instead of continuous ones. One
must be sure that the solutions of the former problems, obtained in real time
(i.e., when the stepsize h is not infinitesimal) are good approximations of the
solutions of the latter ones. However, since the discrete world is much richer
than the continuous one (the latter being a limit case of the former), the
classical definitions and techniques, devised to analyze the behaviors of
continuous problems, are often insufficient to handle the discrete case, and
new specific tools are needed. Often, the insistence in following a path
already traced in the continuous setting, has caused waste of time and efforts,
whereas new specific tools have solved the problems both more easily and
elegantly. In this paper we survey three of the main difficulties encountered
in the numerical solutions of ODEs, along with the novel solutions proposed.Comment: 25 pages, 4 figures (typos fixed
- …