Numerical study of the generalised Klein-Gordon equations

24 pages, 10 figures, 56 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/International audienceIn this study, we discuss an approximate set of equations describing water wave propagating in deep water. These generalized Klein-Gordon (gKG) equations possess a variational formulation, as well as a canonical Hamiltonian and multi-symplectic structures. Periodic travelling wave solutions are constructed numerically to high accuracy and compared to a seventh-order Stokes expansion of the full Euler equations. Then, we propose an efficient pseudo-spectral discretisation, which allows to assess the stability of travelling waves and localised wave packets

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

Classical integrable field theories with defects and near-integrable boundaries

In the first part of this thesis algebro-geometric solutions for the sine-Gordon and KdV equations in the presence of a type I integrable defect are found, generalising the previously known soliton solutions. Elliptic (genus one) solutions where the defect induces only a phase shift are obtained via ansätze for the fields on each side of the defect. Algebro-geometric solutions for arbitrary genus and involving soliton emission by the defect are constructed using a Darboux transformation, exploiting the fact that the defect equations have the form of a Bäcklund transformation at a point. All the soliton and phase-shifted elliptic solutions to the defect equations are recovered as limits of the algebro-geometric solutions constructed in this way. Certain energy and momentum conserving defects for the Kadomtsev-Petviashvili equation are then presented as a first step towards the construction of integrable defects in higher dimensions. Algebro-geometric solutions to the sine-Gordon equation on the half-line with an integrable two parameter boundary condition are obtained by imposing a corresponding restriction on the Lax pair eigenfunction or, alternatively, as a Darboux transformation of the known algebro-geometric solution for the Dirichlet boundary. Finally, the collision of sine-Gordon solitons with a Robin type boundary is examined. This boundary is typically non-integrable but becomes an integrable Neumann or Dirichlet boundary for certain values of a boundary parameter. Depending on the boundary parameter and initial velocity an antikink may be reflected into various combinations of kinks, antikinks and breathers. The soliton content of the field after the collision is numerically determined by computing the discrete scattering data associated with the inverse scattering method. A highlight of this investigation is the discovery of an intricate structure of resonance windows caused by the production of a breather which can collide multiple times with the boundary before escaping as a lighter breather or antikink

エネルギー関数を持つ発展方程式に対する幾何学的数値計算法

学位の種別: 課程博士審査委員会委員 : (主査)東京大学教授 松尾 宇泰, 東京大学教授 中島 研吾, 東京大学准教授 鈴木 秀幸, 東京大学准教授 長尾 大道, 東京大学准教授 齋藤 宣一University of Tokyo(東京大学