A Parameterized multi-step Newton method for solving systems of nonlinear equations

We construct a novel multi-step iterative method for solving systems of nonlinear equations by introducing a parameter. to generalize the multi-step Newton method while keeping its order of convergence and computational cost. By an appropriate selection of theta, the new method can both have faster convergence and have larger radius of convergence. The new iterative method only requires one Jacobian inversion per iteration, and therefore, can be efficiently implemented using Krylov subspace methods. The new method can be used to solve nonlinear systems of partial differential equations, such as complex generalized Zakharov systems of partial differential equations, by transforming them into systems of nonlinear equations by discretizing approaches in both spatial and temporal independent variables such as, for instance, the Chebyshev pseudo-spectral discretizing method. Quite extensive tests show that the new method can have significantly faster convergence and significantly larger radius of convergence than the multi-step Newton method.Peer ReviewedPostprint (author's final draft

Numerical solution of the Klein-Gordon equation in an unbounded domain.

Master of Science in Applied Mathematics. University of KwaZulu-Natal, Westville, 2018.Abstract available in PDF file

NUMERICAL STUDIES ON THE KLEIN-GORDON-SCHRODINGER EQUATIONS IN THE SINGULAR LIMIT REGIME

Master'sMASTER OF SCIENC

Uniformly Accurate Methods for Klein-Gordon type Equations

The main contribution of this thesis is the development of a novel class of uniformly accurate methods for Klein-Gordon type equations. Klein-Gordon type equations in the non-relativistic limit regime, i.e., $c\gg 1$, are numerically very challenging to treat, since the solutions are highly oscillatory in time. Standard Gautschi-type methods suffer from severe time step restrictions as they require a CFL-condition $c^2\tau<1$ with time step size $\tau$ to resolve the oscillations. Within this thesis we overcome this difficulty by introducing limit integrators, which allows us to reduce the highly oscillatory problem to the integration of a non-oscillatory limit system. This procedure allows error bounds of order $\mathcal{O}(c^{-2}+\tau^2)$ without any step size restrictions. Thus, these integrators are very efficient in the regime $c\gg 1$. However, limit integrators fail for small values of $c$. In order to derive numerical schemes that work well for small as well as for large $c$, we use the ansatz of "twisted variables", which allows us to develop uniformly accurate methods with respect to $c$. In particular, we introduce efficient and robust uniformly accurate exponential-type integrators which resolve the solution in the relativistic regime as well as in the highly oscillatory non-relativistic regime without any step size restriction. In contrast to previous works, we do not employ any asymptotic nor multiscale expansion of the solution. Compared to classical methods our new schemes allow us to reduce the regularity assumptions as they converge under the same regularity assumptions required for the integration of the corresponding limit system. In addition, the newly derived first- and second-order exponential-type integrators converge to the classical Lie and Strang splitting schemes for the limit system. Moreover, we present uniformly accurate schemes for the Klein-Gordon-Schrödinger and the Klein-Gordon-Zakharov system. For all uniformly accurate integrators we establish rigorous error estimates and underline their uniform convergence property numerically

Nonlinear Waves and Dispersive Equations

[no abstract available