Anna Kalogirou

Variational Water Wave Modelling: from Continuum to Experiment.

Variational methods are investigated asymptotically and numerically to model water waves in tanks with wave generators. As a validation, our modelling results using (dis)continuous Galerkin finite element methods will be compared to a soliton splash event resulting after a sluice gate is removed during a finite time in a long water channel with a contraction at its end.

An in-depth numerical study of the two-dimensional Kuramoto–Sivashinsky equation

The Kuramoto–Sivashinsky equation in one spatial dimension (1D KSE) is one of the most well-known and well-studied partial differential equations. It exhibits spatio-temporal chaos that emerges through various bifurcations as the domain length increases. There have been several notable analytical studies aimed at understanding how this property extends to the case of two spatial dimensions. In this study, we perform an extensive numerical study of the Kuramoto–Sivashinsky equation (2D KSE) to complement this analytical work. We explore in detail the statistics of chaotic solutions and classify the solutions that arise for domain sizes where the trivial solution is unstable and the long-time dynamics are completely two-dimensional. While we find that many of the features of the 1D KSE, including how the energy scales with system size, carry over to the 2D case, we also note several differences including the various paths to chaos that are not through period doubling.

Nonlinear dynamics of a dispersive anisotropic Kuramoto–Sivashinsky equation in two space dimensions

A Kuramoto–Sivashinsky equation in two space dimensions arising in thin film flows is considered on doubly periodic domains. In the absence of dispersive effects, this anisotropic equation admits chaotic solutions for sufficiently large length scales with fully two-dimensional profiles; the one-dimensional dynamics observed for thin domains are structurally unstable as the transverse length increases. We find that, independent of the domain size, the characteristic length scale of the profiles in the streamwise direction is about 10 space units, with that in the transverse direction being approximately three times larger. Numerical computations in the chaotic regime provide an estimate for the radius of the absorbing ball in L2 in terms of the length scales, from which we conclude that the system possesses a finite energy density. We show the property of equipartition of energy among the low Fourier modes, and report the disappearance of the inertial range when solution profiles are two-dimensional. Consideration of the high-frequency modes allows us to compute an estimate for the analytic extensibility of solutions in C2. We also examine the addition of a physically derived third-order dispersion to the problem; this has a destabilizing effect, in the sense of reducing analyticity and increasing amplitude of solutions. However, sufficiently large dispersion may regularize the spatio-temporal chaos to travelling waves. We focus on dispersion where chaotic dynamics persist, and study its effect on the interfacial structures, absorbing ball and properties of the power spectrum.

Instability of two-layer film flows due to the interacting effects of surfactants, inertia and gravity

We consider a two-fluid shear flow where the interface between the two fluids is coated with an insoluble surfactant. An asymptotic model is derived in the thin-layer approximation, consisting of a set of nonlinear partial differential equations describing the evolution of the film and surfactant disturbances at the interface. The model includes important physical effects such as Marangoni forces (caused by the presence of surfactant), inertial forces arising in the thick fluid layer, as well as gravitational forces. The aim of this study is to investigate the effect of density stratification or gravity—represented through the Bond number Bo—on the flow stability and the interplay between the different (de)stabilisation mechanisms. It is found that gravity can either stabilise or destabilise the interface (depending on fluid properties) but not always as intuitively anticipated. Different traveling-wave branches are presented for varying Bo, and the destabilising mechanism associated with the Marangoni fo...

Mathematical and Numerical Modelling of Wave Impact on Wave-Energy Buoys

We report on the mathematical and numerical modelling of amplified rogue waves driving a wave-energy device in a contraction. This wave-energy device consists of a floating buoy attached to an AC-induction motor and constrained to move upward only in a contraction, for which we have realised a working scale-model. A coupled Hamiltonian system is derived for the dynamics of water waves and moving wave-energy buoys. This nonlinear model consists of the classical water wave equations for the free surface deviation and velocity potential, coupled to a set of equations describing the dynamics of a wave-energy buoy. As a stepping stone, the model is solved numerically for the case of linear shallow water waves causing the motion of a simple buoy structure with V-shaped cross-sections, using a variational (dis)continuous Galerkin finite element method.

Modelling of Nonlinear Wave-Buoy Dynamics Using Constrained Variational Methods

We consider a comprehensive mathematical and numerical strategy to couple water-wave motion with rigid ship dynamics using variational principles. We present a methodology that applies to three-dimensional potential flow water waves and ship dynamics. For simplicity, in this paper we demonstrate the method for shallow-water waves coupled to buoy motion in two dimensions, the latter being the symmetric motion of a crosssection of a ship. The novelty in the presented model is that it employs a Lagrange multiplier to impose a physical restriction on the water height under the buoy in the form of an inequality constraint. A system of evolution equations can be obtained from the model and consists of the classical shallow-water equations for shallow, incompressible and irrotational waves, and relevant equations for the dynamics of the wave-energy buoy. One of the advantages of the variational approach followed is that, when combined with symplectic integrators, it eliminates any numerical damping and preserves the discrete energy; this is confirmed in our numerical results.