Theodoros Katsaounis

Denys Dutykh, Theodoros Katsaounis, Dimitrios Mitsotakis

Finite volume schemes are commonly used to construct approximate solutions to conservation laws. In this study we extend the framework of the finite volume methods to dispersive water wave models, in particular to Boussinesq type systems. We focus mainly on the application of the method to bidirectional nonlinear, dispersive wave propagation in one space dimension. Special emphasis is given to important nonlinear phenomena such as solitary waves interactions, dispersive shock wave formation and the runup of breaking and non-breaking long waves.

Denys Dutykh, Theodoros Katsaounis, Dimitrios Mitsotakis

We extend the framework of the finite volume method to dispersive unidirectional water wave propagation in one space dimension. In particular we consider a KdV-BBM type equation. Explicit and IMEX Runge-Kutta type methods are used for time discretizations. The fully discrete schemes are validated by direct comparisons to analytic solutions. Invariants conservation properties are also studied. Main applications include important nonlinear phenomena such as dispersive shock wave formation, solitary waves and their various interactions.

Denys Dutykh, Theodoros Katsaounis, Dimitrios Mitsotakis

Historically the finite volume methods have been developed for the numerical integration of conservation laws. In this study we present some recent results on the application of such schemes to dispersive PDEs. Namely, we solve numerically a representative of Boussinesq type equations in view of important applications to the coastal hydrodynamics. Numerical results of the runup of a moderate wave onto a non-uniform beach are presented along with great lines of the employed numerical method (see D. Dutykh et al. (2011) for more details).

Theodoros Katsaounis, Irene Kyza

We provide a posteriori error estimates in the $L^\infty(L^2)-$norm for relaxation time discrete and fully discrete schemes for a class of evolution nonlinear Schrödinger equations up to the critical exponent. In particular for the discretization in time we use the relaxation Crank-Nicolson-type scheme introduced by Besse in \cite{Besse}. For the discretization in space we use finite element spaces that are allowed to change between time steps. The estimates are obtained using the reconstruction technique. Through this technique the problem is converted to a perturbation of the original partial differential equation and this makes it possible to use nonlinear stability arguments as in the continuous problem. In particular, main ingredients we use in our analysis are the Gagliardo-Nirenberg inequality and the two conservation laws (mass and energy conservation) of the continuous problem. Numerical results illustrate that the estimates are indeed of optimal order of convergence.

Theodoros Katsaounis, Irene Kyza

We derive optimal order a posteriori error estimates for fully discrete approximations of linear Schrödinger-type equations, in the $L^\infty(L^2)-$norm. For the discretization in time we use the Crank-Nicolson method, while for the space discretization we use finite element spaces that are allowed to change in time. The derivation of the estimators is based on a novel elliptic reconstruction that leads to estimates which reflect the physical properties of Schrödinger equations. The final estimates are obtained using energy techniques and residual-type estimators. Various numerical experiments for the one-dimensional linear Schrödinger equation in the semiclassical regime, verify and complement our theoretical results. The numerical implementations are performed with both uniform partitions and adaptivity in time and space. For adaptivity, we further develop and analyze an existing time-space adaptive algorithm to the cases of Schrödinger equations. The adaptive algorithm reduces the computational cost substantially and provides efficient error control for the solution and the observables of the problem, especially for small values of the Planck constant.

Theodoros Katsaounis, Dimitrios Mitsotakis

In this paper we perform a numerical study on the interesting phenomenon of soliton reflection of solid walls. We consider the 2D nonlinear Schrodinger equation as the underlying mathematical model and we use an implicit-explicit type Crank-Nicolson finite element scheme for its numerical solution. After verifying the perfect reflection of the solitons on a vertical wall, we present the imperfect reflection of a dark soliton on a diagonal wall.