2,644 research outputs found
Numerical approximations for population growth model by Rational Chebyshev and Hermite Functions collocation approach: A comparison
This paper aims to compare rational Chebyshev (RC) and Hermite functions (HF)
collocation approach to solve the Volterra's model for population growth of a
species within a closed system. This model is a nonlinear integro-differential
equation where the integral term represents the effect of toxin. This approach
is based on orthogonal functions which will be defined. The collocation method
reduces the solution of this problem to the solution of a system of algebraic
equations. We also compare these methods with some other numerical results and
show that the present approach is applicable for solving nonlinear
integro-differential equations.Comment: 18 pages, 5 figures; Published online in the journal of "Mathematical
Methods in the Applied Sciences
A fast and well-conditioned spectral method for singular integral equations
We develop a spectral method for solving univariate singular integral
equations over unions of intervals by utilizing Chebyshev and ultraspherical
polynomials to reformulate the equations as almost-banded infinite-dimensional
systems. This is accomplished by utilizing low rank approximations for sparse
representations of the bivariate kernels. The resulting system can be solved in
operations using an adaptive QR factorization, where is
the bandwidth and is the optimal number of unknowns needed to resolve the
true solution. The complexity is reduced to operations by
pre-caching the QR factorization when the same operator is used for multiple
right-hand sides. Stability is proved by showing that the resulting linear
operator can be diagonally preconditioned to be a compact perturbation of the
identity. Applications considered include the Faraday cage, and acoustic
scattering for the Helmholtz and gravity Helmholtz equations, including
spectrally accurate numerical evaluation of the far- and near-field solution.
The Julia software package SingularIntegralEquations.jl implements our method
with a convenient, user-friendly interface
An approximation algorithm for the solution of the nonlinear Lane-Emden type equations arising in astrophysics using Hermite functions collocation method
In this paper we propose a collocation method for solving some well-known
classes of Lane-Emden type equations which are nonlinear ordinary differential
equations on the semi-infinite domain. They are categorized as singular initial
value problems. The proposed approach is based on a Hermite function
collocation (HFC) method. To illustrate the reliability of the method, some
special cases of the equations are solved as test examples. The new method
reduces the solution of a problem to the solution of a system of algebraic
equations. Hermite functions have prefect properties that make them useful to
achieve this goal. We compare the present work with some well-known results and
show that the new method is efficient and applicable.Comment: 34 pages, 13 figures, Published in "Computer Physics Communications
Spectral method for matching exterior and interior elliptic problems
A spectral method is described for solving coupled elliptic problems on an
interior and an exterior domain. The method is formulated and tested on the
two-dimensional interior Poisson and exterior Laplace problems, whose solutions
and their normal derivatives are required to be continuous across the
interface. A complete basis of homogeneous solutions for the interior and
exterior regions, corresponding to all possible Dirichlet boundary values at
the interface, are calculated in a preprocessing step. This basis is used to
construct the influence matrix which serves to transform the coupled boundary
conditions into conditions on the interior problem. Chebyshev approximations
are used to represent both the interior solutions and the boundary values. A
standard Chebyshev spectral method is used to calculate the interior solutions.
The exterior harmonic solutions are calculated as the convolution of the
free-space Green's function with a surface density; this surface density is
itself the solution to an integral equation which has an analytic solution when
the boundary values are given as a Chebyshev expansion. Properties of Chebyshev
approximations insure that the basis of exterior harmonic functions represents
the external near-boundary solutions uniformly. The method is tested by
calculating the electrostatic potential resulting from charge distributions in
a rectangle. The resulting influence matrix is well-conditioned and solutions
converge exponentially as the resolution is increased. The generalization of
this approach to three-dimensional problems is discussed, in particular the
magnetohydrodynamic equations in a finite cylindrical domain surrounded by a
vacuum
Spectral methods for modeling supersonic chemically reacting flow fields
A numerical algorithm was developed for solving the equations describing chemically reacting supersonic flows. The algorithm employs a two-stage Runge-Kutta method for integrating the equations in time and a Chebyshev spectral method for integrating the equations in space. The accuracy and efficiency of the technique were assessed by comparison with an existing implicit finite-difference procedure for modeling chemically reacting flows. The comparison showed that the procedure presented yields equivalent accuracy on much coarser grids as compared to the finite-difference procedure with resultant significant gains in computational efficiency
High-order, Dispersionless "Fast-Hybrid" Wave Equation Solver. Part I: Sampling Cost via Incident-Field Windowing and Recentering
This paper proposes a frequency/time hybrid integral-equation method for the
time dependent wave equation in two and three-dimensional spatial domains.
Relying on Fourier Transformation in time, the method utilizes a fixed
(time-independent) number of frequency-domain integral-equation solutions to
evaluate, with superalgebraically-small errors, time domain solutions for
arbitrarily long times. The approach relies on two main elements, namely, 1) A
smooth time-windowing methodology that enables accurate band-limited
representations for arbitrarily-long time signals, and 2) A novel Fourier
transform approach which, in a time-parallel manner and without causing
spurious periodicity effects, delivers numerically dispersionless
spectrally-accurate solutions. A similar hybrid technique can be obtained on
the basis of Laplace transforms instead of Fourier transforms, but we do not
consider the Laplace-based method in the present contribution. The algorithm
can handle dispersive media, it can tackle complex physical structures, it
enables parallelization in time in a straightforward manner, and it allows for
time leaping---that is, solution sampling at any given time at
-bounded sampling cost, for arbitrarily large values of ,
and without requirement of evaluation of the solution at intermediate times.
The proposed frequency-time hybridization strategy, which generalizes to any
linear partial differential equation in the time domain for which
frequency-domain solutions can be obtained (including e.g. the time-domain
Maxwell equations), and which is applicable in a wide range of scientific and
engineering contexts, provides significant advantages over other available
alternatives such as volumetric discretization, time-domain integral equations,
and convolution-quadrature approaches.Comment: 33 pages, 8 figures, revised and extended manuscript (and now
including direct comparisons to existing CQ and TDIE solver implementations)
(Part I of II
- …