Source author record

V. A. Dougalis

V. A. Dougalis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis math.AP

Catalog footprint

What is connected

7works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Notes on the Boussinesq-Full dispersion systems for internal waves: Numerical solution and solitary waves

In this paper we study some theoretical and numerical issues of the Boussinesq/Full dispersion system. This is a a three-parameter system of pde's that models the propagation of internal waves along the interface of two-fluid layers with rigid lid condition for the upper layer, and under a Boussinesq regime for the upper layer and a full dispersion regime for the lower layer. We first discretize in space the periodic initial-value problem with a Fourier-Galerkin spectral method and prove error estimates for several ranges of values of the parameters. Solitary waves of the model systems are then studied numerically in several ways. The numerical generation is analyzed by approximating the ode system with periodic boundary conditions for the solitary-wave profiles with a Fourier spectral scheme, implemented in a collocation form, and solving iteratively the corresponding algebraic system in Fourier space with the Petviashvili method accelerated with the minimal polynomial extrapolation technique. Motivated by the numerical results, a new result of existence of solitary waves is proved. In the last part of the paper, the dynamics of these solitary waves is studied computationally, To this end, the semidiscrete systems obtained from the Fourier-Galerkin discretization in space are integrated numerically in time by a Runge-Kutta Composition method of order four. The fully discrete scheme is used to explore numerically the stability of solitary waves, their collisions, and the resolution of other initial conditions into solitary waves.

preprint2020arXiv

Galerkin finite element methods for the numerical solution of two classical-Boussinesq type systems over variable bottom topography

We consider two `Classical' Boussinesq type systems modelling two-way propagation of long surface waves in a finite channel with variable bottom topography. Both systems are derived from the 1-d Serre-Green-Naghdi (SGN) system; one of them is valid for stronger bottom variations, and coincides with Peregrine's system, and the other is valid for smaller bottom variations. We discretize in the spatial variable simple initial-boundary-value problems (ibvp's) for both systems using standard Galerkin-finite element methods and prove $L^2$ error estimates for the ensuing semidiscrete approximations. We couple the schemes with the 4th order-accurate, explicit, classical Runge-Kutta time-stepping procedure and use the resulting fully discrete methods in numerical simulations of dispersive wave propagation over variable bottoms with several kinds of boundary conditions, including absorbing ones. We describe in detail the changes that solitary waves undergo when evolving under each system over a variety of variable-bottom environments. We assess the efficacy of both systems in approximating these flows by comparing the results of their simulations with each other, with simulations of the SGN-system, and with available experimental data from the literature.

preprint2016arXiv

Notes on Galerkin-finite element methods for the Shallow Water equations with characteristic boundary conditions

We consider the Shallow Water equations in the supercritical and subcritical cases in one space variable,posed in a finite spatial interval with characteristic boundary conditions at the endpoints, which, as is well known, are transparent,i.e. allow outgoing waves to exit without generating spurious reflected waves. Assuming that the resulting initial-boundary-value problems have smooth solutions,we approximate them in space using standard Galerkin-finite element methods and prove L^2 error estimates for the semidiscrete problems on quasiuniform meshes.We discretize the problems in the temporal variable using an explicit,fourth-order accurate Runge-Kutta scheme and check, by means of numerical experiment, that the resulting fully discrete schemes have excellent absorption properties.

preprint2014arXiv

Notes on error estimates for the standard Galerkin-finite element method for the Shallow Water equations

We consider a simple initial-boundary-value problem for the shallow water equations in one space dimension, and also the analogous problem for a symmetric variant of the system. Assuming smoothness of solutions, we discretize these problems in space using standard Galerkin-finite element methods and prove $L^{2}$-error estimates for the semidiscrete problems for quasiuniform and uniform meshes. In particular we show that in the case of spatial discretizations with piecewise linear continuous functions on a uniform mesh, suitable compatibility conditions at the boundary and superaccuracy properties of the $L^{2}$ projection on the finite element subspaces lead to an optimal-order $O(h^{2})$ $L^{2}$-error estimate. We also examine temporal discretizations of the semidiscrete problems by three explicit Runge-Kutta methods (the Euler, improved Euler, and the Shu-Osher scheme) and prove $L^{2}$-error estimates, which are of optimal order in the temporal variable, under appropriate stability conditions. In a final section of remarks we prove optimal-order $L^{2}$-error estimates for smooth spline spatial discretizations of the periodic initial-value problem for the systems. We also prove that small-amplitude, appropriately transformed solutions of the symmetric system are close to the corresponding solutions of the usual system while they are both smooth, thus providing a justification of the symmetric system.

preprint2014arXiv

Notes on the numerical solution of the Benjamin equation

In this paper we consider the Benjamin equation, a partial differential equation that models one-way propagation of long internal waves of small amplitude along the interface of two fluid layers under the effects of gravity and surface tension. We solve the periodic initial-value problem for the Benjamin equation numerically by a new fully discrete hybrid finite-element / spectral scheme, which we first validate by pinning down its accuracy and stability properties. After testing the evolution properties of the scheme in a study of propagation of single - and multi-pulse solitary waves of the Benjamin equation, we use it in an exploratory mode to illuminate phenomena such as overtaking collisions of solitary waves, and the stability of single-, multi-pulse and 'depression' solitary waves.

preprint2011arXiv

A Finite Difference method for the Wide-Angle `Parabolic' equation in a waveguide with downsloping bottom

We consider the third-order wide-angle `parabolic' equation of underwater acoustics in a cylindrically symmetric fluid medium over a bottom of range-dependent bathymetry. It is known that the initial-boundary-value problem for this equation may not be well posed in the case of (smooth) bottom profiles of arbitrary shape if it is just posed e.g. with a homogeneous Dirichlet bottom boundary condition. In this paper we concentrate on downsloping bottom profiles and propose an additional boundary condition that yields a well posed problem, in fact making it $L^2$-conservative in the case of appropriate real parameters. We solve the problem numerically by a Crank-Nicolson-type finite difference scheme, which is proved to be unconditionally stable and second-order accurate, and simulates accurately realistic underwater acoustic problems.

preprint2010arXiv

Notes on error estimates for Galerkin approximations of the 'classical' Boussinesq system and related hyperbolic problems

We consider the `classical' Boussinesq system in one space dimension and its symmetric analog. These systems model two-way propagation of nonlinear, dispersive long waves of small amplitude on the surface of an ideal fluid in a uniform horizontal channel. We discretize an initial-boundary-value problem for these systems in space using Galerkin-finite element methods and prove error estimates for the resulting semidiscrete problems and also for their fully discrete analogs effected by explicit Runge-Kutta time-stepping procedures. The theoretical orders of convergence obtained are consistent with the results of numerical experiments that are also presented.