8,579 research outputs found
A multigrid continuation method for elliptic problems with folds
We introduce a new multigrid continuation method for computing solutions of nonlinear elliptic eigenvalue problems which contain limit points (also called turning points or folds). Our method combines the frozen tau technique of Brandt with pseudo-arc length continuation and correction of the parameter on the coarsest grid. This produces considerable storage savings over direct continuation methods,as well as better initial coarse grid approximations, and avoids complicated algorithms for determining the parameter on finer grids. We provide numerical results for second, fourth and sixth order approximations to the two-parameter, two-dimensional stationary reaction-diffusion problem: Δu+λ exp(u/(1+au)) = 0.
For the higher order interpolations we use bicubic and biquintic splines. The convergence rate is observed to be independent of the occurrence of limit points
Stabilization of Unstable Procedures: The Recursive Projection Method
Fixed-point iterative procedures for solving nonlinear parameter dependent problems can converge for some interval of parameter values and diverge as the parameter changes. The Recursive Projection Method (RPM), which stabilizes such procedures by computing a projection onto the unstable subspace is presented. On this subspace a Newton or special Newton iteration is performed, and the fixed-point iteration is used on the complement. As continuation in the parameter proceeds, the projection is efficiently updated, possibly increasing or decreasing the dimension of the unstable subspace. The method is extremely effective when the dimension of the unstable subspace is small compared to the dimension of the system. Convergence proofs are given and pseudo-arclength continuation on the unstable subspace is introduced to allow continuation past folds. Examples are presented for an important application of the RPM in which a “black-box” time integration scheme is stabilized, enabling it to compute unstable steady states. The RPM can also be used to accelerate iterative procedures when slow convergence is due to a few slowly decaying modes
Evolution of the L1 halo family in the radial solar sail CRTBP
We present a detailed investigation of the dramatic changes that occur in the
halo family when radiation pressure is introduced into the
Sun-Earth circular restricted three-body problem (CRTBP). This
photo-gravitational CRTBP can be used to model the motion of a solar sail
orientated perpendicular to the Sun-line. The problem is then parameterized by
the sail lightness number, the ratio of solar radiation pressure acceleration
to solar gravitational acceleration. Using boundary-value problem numerical
continuation methods and the AUTO software package (Doedel et al. 1991) the
families can be fully mapped out as the parameter is increased.
Interestingly, the emergence of a branch point in the retrograde satellite
family around the Earth at acts to split the halo family
into two new families. As radiation pressure is further increased one of these
new families subsequently merges with another non-planar family at
, resulting in a third new family. The linear stability of
the families changes rapidly at low values of , with several small
regions of neutral stability appearing and disappearing. By using existing
methods within AUTO to continue branch points and period-doubling bifurcations,
and deriving a new boundary-value problem formulation to continue the folds and
Krein collisions, we track bifurcations and changes in the linear stability of
the families in the parameter and provide a comprehensive overview of
the halo family in the presence of radiation pressure. The results demonstrate
that even at small values of there is significant difference to the
classical CRTBP, providing opportunity for novel solar sail trajectories.
Further, we also find that the branch points between families in the solar sail
CRTBP provide a simple means of generating certain families in the classical
case.Comment: 31 pages, 17 figures, accepted by Celestial Mechanics and Dynamical
Astronom
Systematic experimental exploration of bifurcations with non-invasive control
We present a general method for systematically investigating the dynamics and
bifurcations of a physical nonlinear experiment. In particular, we show how the
odd-number limitation inherent in popular non-invasive control schemes, such as
(Pyragas) time-delayed or washout-filtered feedback control, can be overcome
for tracking equilibria or forced periodic orbits in experiments. To
demonstrate the use of our non-invasive control, we trace out experimentally
the resonance surface of a periodically forced mechanical nonlinear oscillator
near the onset of instability, around two saddle-node bifurcations (folds) and
a cusp bifurcation.Comment: revised and extended version (8 pages, 7 figures
Continuation of connecting orbits in 3D-ODEs: (I) Point-to-cycle connections
We propose new methods for the numerical continuation of point-to-cycle
connecting orbits in 3-dimensional autonomous ODE's using projection boundary
conditions. In our approach, the projection boundary conditions near the cycle
are formulated using an eigenfunction of the associated adjoint variational
equation, avoiding costly and numerically unstable computations of the
monodromy matrix. The equations for the eigenfunction are included in the
defining boundary-value problem, allowing a straightforward implementation in
AUTO, in which only the standard features of the software are employed.
Homotopy methods to find connecting orbits are discussed in general and
illustrated with several examples, including the Lorenz equations. Complete
AUTO demos, which can be easily adapted to any autonomous 3-dimensional ODE
system, are freely available.Comment: 18 pages, 10 figure
- …