Essential spectrum of systems of singular differential equations

In this paper we develop a new method to determine the essential spectrum of coupled systems of singular differential equations. Applications to problems from magnetohydrodynamics and astrophysics are given.

Fully Adaptive Newton--Galerkin Methods for Semilinear Elliptic Partial Differential Equations

In this paper we develop an adaptive procedure for the numerical solution of general, semilinear elliptic problems with possible singular perturbations. Our approach combines both prediction-type adaptive Newton methods and a linear adaptive finite element discretization (based on a robust a posteriori error analysis), thereby leading to a fully adaptive Newton–Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for various examples

Computational Comparison of Continuous and Discontinuous Galerkin Time-Stepping Methods for Nonlinear Initial Value Problems

This article centers on the computational performance of the continuous and discontinuous Galerkin time stepping schemes for general first-order initial value problems in R n , with continuous nonlinearities. We briefly review a recent existence result for discrete solutions from [6], and provide a numerical comparison of the two time discretization methods.

Terahertz macrospin dynamics in insulating ferrimagnets

We investigate numerically the excitation of nonlinear magnetic interactions in a ferrite material by an energetic pump pulse of terahertz (THz) radiation. The calculations are performed by solving the coupled Maxwell and Landau-Lifshitz-Gilbert differential equations. In a time-resolved THz pump/THz probe scheme, it is demonstrated that Faraday rotation of a delayed THz probe pulse can be used to map these interactions. Our study is motivated by the ability of soft x-ray free electron lasers to perform time-resolved imaging of the magnetization process at the submicrometer and subpicosecond length and time scales.

An adaptive Newton-method based on a dynamical systems approach

The traditional Newton method for solving nonlinear operator equations in Banach spaces is discussed within the context of the continuous Newton method. This setting makes it possible to interpret the Newton method as a discrete dynamical system and thereby to cast it in the framework of an adaptive step size control procedure. In so doing, our goal is to reduce the chaotic behavior of the original method without losing its quadratic convergence property close to the roots. The performance of the modified scheme is illustrated with various examples from algebraic and differential equations.