A UNIFORMLY AND OPTIMALLY ACCURATE METHOD FOR THE KLEIN-GORDON-ZAKHAROV SYSTEM IN SIMULTANEOUS HIGH-PLASMA-FREQUENCY AND SUBSONIC LIMIT REGIME *

We present a uniformly and optimally accurate numerical method for solving the Klein-Gordon-Zakharov (KGZ) system with two dimensionless parameters 0 < ε ≤ 1 and 0 < γ ≤ 1, which are inversely proportional to the plasma frequency and the acoustic speed, respectively. In the simultaneous high-plasma-frequency and subsonic limit regime, i.e. ε < γ → 0 + , the KGZ system collapses to a cubic Schrödinger equation, and the solution propagates waves with O(ε 2)-wavelength in time and meanwhile contains rapid outgoing initial layers with speed O(1/γ) in space due to the incompatibility of the initial data. By presenting a multiscale decomposition of the KGZ system, we propose a multiscale time integrator Fourier pseduospectral method which is explicit, efficient and uniformly accurate for solving the KGZ system for all 0 < ε < γ ≤ 1. Numerical results are reported to show the efficiency and accuracy of scheme. Finally, the method is applied to investigate the convergence rates of the KGZ system to its limiting models when ε < γ → 0 +

Multiscale Methods and Analysis for Highly Oscillatory Differential Equations

Ph.DDOCTOR OF PHILOSOPH

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

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

Master'sMASTER OF SCIENC

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described