Coastal flow models and boundary conditions

24 oct. 2022, 10:00 → 27 oct. 2022, 14:00 Europe/Paris

Amphithéâtre Laurent Schwartz, building 1R3 (Institut de Mathématiques de Toulouse)

Aims and scope

 

This workshop will address the interplay between boundary conditions and transport or dispersive phenomena. A central topic will be the modelling of physical processes in which wave propagation is mainly governed by dispersion, be it encoded either in the underlying governing equations or in the chosen numerical discretization, and wave breaking that comes from the nonlinear features of the physical process. The model equations that will be considered range from toy linear models such as the transport, Schrödinger or Airy equations, to complex, nonlinear systems such as the Euler-Korteweg, Boussinesq, Green-Naghdi or the full water wave equations, as arising in coastal oceanography. We wish to bring together experts in modelling issues, analysis of partial differential equations, numerical analysis and scientific computing in order to foster exchanges between various view points.

 

The workshop is part of the thematic semester Modelling, numerical analysis and scientific computing in PDEs organized by the Labex CIMI and will also serve as the closing conference of the ANR project NABUCO.

 

List of speakers

 

• Mini-courses :
• Jean-François Coulombel & Grégory Faye (CNRS, Toulouse)
• Vincent Duchêne (CNRS, Rennes)
• Anne Mangeney (Université Paris Cité et Institut Universitaire de France)

 

• Confirmed talks
• Corentin Audiard (Sorbonne Université, Paris)
• Geoffrey Beck (INRIA, Rennes)
• Eric Blayo (Université Grenoble Alpes)
• David Chiron (Université Côte d'Azur, Nice)
• Ali Haidar (Université de Montpellier)
• Maria Kazakova (Université Savoie Mont-Blanc, Chambéry)
• Pierre Le Barbenchon (Université Rennes 1)
• Martin Parisot (INRIA, Bordeaux)
• Aric Wheeler (Indiana University, Bloomington)

 

Program, titles and abstracts

 

All details can be found in the document attached at the bottom of this page (Programme_workshop.pdf).

 

Practical information

 

The workshop will take place within the Institut de Mathématiques de Toulouse, in the Amphitheater Laurent Schwartz. A detailed map of campus is available here : Access to campus. The amphitheater is located in building no 23 on the campus map (rue Sébastienne Guyot).

From downtown, take the Metro B line up to "Université Paul Sabatier" station. The 1R3 mathematics building is two hundred meters away from the Metro station.

Programme_Workshop.pdf

lundi 24 octobre

12:00 Lunch

1 Large time stability issues for finite difference schemes (session 1)

We prove generalized Gaussian bounds for the large time behavior of finite difference approximations of the transport equation. The goal of the talk is to review the connection with the central limit theorem in probability theory, some history of the problem and to present the general methodology leading to sharp estimates for the Green's function. The talk is based on a joint work with Grégory Faye.

Orateur: Jean-François Coulombel (Institut de Mathématiques de Toulouse)

Illustrations.pdf

15:30 Coffee break

2 Derivation and justification of shallow water models for surface gravity waves and stratified flows (session 1)

We shall discuss how standard models in oceanography can be rigorously justified as asymptotic models, here in the shallow water regime, focusing on the general strategy and mathematical tools which are involved. Most of the discussion will be devoted to the justification of the Saint-Venant system for homogeneous potential flows with a free surface. Yet going through this somewhat standard procedure will allow us to unveil, by contrast, how much is yet to be understood when we venture into non-hydrostatic models or heterogeneous flows.

Orateur: Vincent Duchêne (CNRS & Université de Rennes 1)

mardi 25 octobre

3 Waves - floating structure interactions in Boussinesq regime

In the context of nearshore wave energy facilities, we have tackled, with David Lannes and Lisl Weynans, the interaction of waves with a floating structure immersed in a 2D fluid. Some difficulties come from the presence of several surfaces: the surface of the sea and the contact surface between the structure and the fluid. The horizontal plane is decomposed into two regions: the exterior region where the surface of the fluid is in contact with the air, and the interior region where it is in contact with the bottom of the object. In the exterior region, we have the standard wave equations, where the surface is free but the pressure is constrained (equal to the atmospheric pressure). In the interior region, the opposite happens: the pressure is free but the surface is constrained, which changes the structure of the equations. Finally, coupling conditions between both regions are needed. We show how to implement this program in the case where the waves are described by the nonlinear dispersive Boussinesq equations. The difficulties related to the computation of the contact surface are overcome by considering an augmented formulation: in addition to the usual equations, we find two hidden ODEs of the water column at the contact points between the waves and the structure. Finally, we propose a numerical method that exploits the added-mass effect, the dispersive boundary layer and this two hidden ODEs.

Orateur: Dr Geoffrey Beck (INRIA Rennes)

Slides-Nabuco.pdf

4 Numerical simulation of shallow-water and floating object nonlinear interactions

In this work, a novel numerical algorithm is introduced for the study of nonlinear interactions between free-surface shallow-water flows and a partly immersed floating object. The object's motion may be either prescribed, or computed as a response to the hydrodynamic forcing. Away from the object, the nonlinear hyperbolic shallow-water equations are used, while the description of the flow beneath the floating object reduces to a nonlinear algebraic equation for the free-surface, together with a nonlinear ordinary differential equation for the computation of the discharge. In both domains, the model accounts for the possible topography variations. When a free-motion is allowed (with heaving, surging and pitching in the horizontal one-dimensional case), this set of equations has to be supplemented with the Newton's second law for the object's motion, involving the force and torque applied over the object by the surrounding fluid, and a part of this external forcing is regarded as an added-mass effect, in order to benefit from its stabilizing influence. At the discrete level, we introduce a discontinuous Galerkin (DG) approximation relying on some arbitrary-order polynomials. This DG method is stabilized with a recent \textit{a posteriori} Local Subcell Correction method through Finite-Volume, in the vicinity of the solution's singularities. The motion of the water-object contact-points is described with the help of an Arbitrary-Lagrangian-Eulerian description. We show that the discontinuous nature of the chosen DG approximation leads to a very natural coupling between the exterior and interior flows, resulting in a global method which ensures the preservation of the water-height positivity at the sub-cell level, preserves the class of motionless steady-states (well-balancing) and preserves the Geometric Conservation Law. Several numerical computations, involving well-balancing, prescribed motions and their impact on the surrounding fluid, or the nonlinear interactions between the object and surface waves, are investigated and highlight that our numerical model effectively allows to model the wave-floating body interactions, with a robust computation of the air-water-body contact-points evolution, as well as of the possible strong flow singularities, and retains the highly accurate sub-cell resolution of DG schemes.

Orateur: M. Ali Haidar (Université de Montpellier)

Haidar.pdf

10:40 Coffee break

5 Quelques aspects mathématiques et numériques dans la représentation des interactions océan-atmosphère

Les interactions océan-atmosphère (OA) jouent un rôle important dans de nombreux phénomènes, tels que cyclones tropicaux ou dynamique du climat. La représentation de ces interactions au sein d’un système de modélisation OA consiste principalement à évaluer les flux échangés entre les deux milieux, et à les imposer à l’interface air-mer.
Cet exposé abordera quelques questions mathématiques et numériques liées à ces deux aspects et présentera une synthèse des travaux récents menés sur ces thématiques dans notre équipe au Laboratoire Jean Kuntzmann : modélisation des couches limites de surface de part et d’autre de l’interface air-mer, schémas de discrétisation au voisinage de cette interface, algorithmique de couplage.

Orateur: Prof. Eric Blayo (Université Grenoble Alpes)

talk-CL-Toulouse22.pdf

12:15 Lunch

6 Seismic waves: a unique source of information on glaciers and landslides (session 1)

I will show here how the coupling between numerical modelling of glaciers and landslides associated with the analysis of the generated seismic waves makes it possible to recover unique information on the source processes. These waves provide a unique tool to constrain the models by providing detailed information on the dynamics of landslides and iceberg calving. This approach allows quantifying in particular the location and volume of events and relate landslide or iceberg calving activity with external forcing such as seismic, volcanic or climatic activity, with strong implication on hazard assessment.

Orateur: Prof. Anne Mangeney (Université Paris Cité)

anr-nabucco_24-10-2022_mangeney-cours1.pptx

15:30 Coffee break

7 Seismic waves: a unique source of information on glaciers and landslides (session 2)

I will show here how the coupling between numerical modelling of glaciers and landslides associated with the analysis of the generated seismic waves makes it possible to recover unique information on the source processes. These waves provide a unique tool to constrain the models by providing detailed information on the dynamics of landslides and iceberg calving. This approach allows quantifying in particular the location and volume of events and relate landslide or iceberg calving activity with external forcing such as seismic, volcanic or climatic activity, with strong implication on hazard assessment.

Orateur: Prof. Anne Mangeney (Université Paris Cité)

anr-nabuco-24-10-2022_mangeney-cours2.pptx

19:30 Social dinner

mercredi 26 octobre

8 Sharp regularity results for solutions of boundary value problems

We study the well-posedness of initial boundary value problems for the linear Schrödinger equations on a half space. The boundary data lie in a (allegedly optimal) Bourgain type Sobolev space, which allows to include Neuman and transparent boundary conditions in the analysis. Strichartz estimates (in $L^2$) are obtained thanks to an explicit solution formula. In the case of Dirichlet boundary data, the regularity of solutions is obtained provided natural compatibility conditions are satisfied. The regularity results concern fractional regularity, and include the more delicate case where the initial data are in $H^{1/2}$. The proof of regularity uses an interpolation argument that can be applied to other boundary value problems.

Orateur: Corentin Audiard (Sorbonne Université)

toulouse_2022.pdf

9 Nonlinear damping vs. formation of singularities

We explore the boundaries of damping estimates by comparing and contrasting two closely related models of combustion, the Majda and ZND models. We show that singularities form in the unweighted Lipschitz norm in finite time on both sides of the shock for both models, extending classical results of John and Liu to suitable variable coefficient systems. On the other hand, we show some energy estimates on the Majda model in exponentially weighted norms, which allows us to obtain an orbital asymptotic stability result, and also that the ZND model does not admit such estimates. This work is joint with Paul Blochas.

Orateur: M. Aric Wheeler (Indiana University)

Wheeler.pdf

10:40 Coffee break

10 Smooth branches of travelling waves for the 2D Nonlinear Schrödinger equation

We shall present some results on the existence of smooth branches of travelling waves for the 2D nonlinear Schrödinger equation parametrized by the speed. In the limit of small speed (joint works with E. Pacherie), the travelling wave has two well separated vortices and we prove that these are the only minimizers of the energy for fixed momentum. In the limit where the speed is close to the speed of sound, we obtain rarefaction pulses described by rational lump solutions to the KP-I equation.

Orateur: Dr David Chiron (Université Côte d'Azur)

Chiron.pdf

12:15 Lunch

11 Large time stability issues for finite difference schemes (session 2)

Following session 1, we shall present some sharp estimates for the Green's function of finite difference schemes that approximate the incoming transport equation on a half-line. The stability analysis relies on a precise description of the possible eigenvalues arising from the interplay between the finite difference scheme and the numerical boundary conditions. The talk is based on a joint work with Jean-François Coulombel.

Orateur: Gregory FAYE

LectureNotesGFaye.pdf

slides_sharp.pdf

15:30 Coffee break

12 Derivation and justification of shallow water models for surface gravity waves and stratified flows (session 2)

We shall discuss how standard models in oceanography can be rigorously justified as asymptotic models, here in the shallow water regime, focusing on the general strategy and mathematical tools which are involved. Most of the discussion will be devoted to the justification of the Saint-Venant system for homogeneous potential flows with a free surface. Yet going through this somewhat standard procedure will allow us to unveil, by contrast, how much is yet to be understood when we venture into non-hydrostatic models or heterogeneous flows.

Orateur: Vincent Duchêne (CNRS & Université de Rennes 1)

jeudi 27 octobre

13 Numerical proof of stability for finite difference schemes in a finite space domain

The goal of this talk is to study the stability of finite differences schemes for scalar hyperbolic initial boundary value problem. It is based on the GKS theory (introduced by Gustafsson, Kreiss and Sundström) and it deals with the Kreiss-Lopantinskii determinant. We will use complex analysis and geometric consideration to find zeros of this determinant and conclude on the stability of the scheme.

Orateur: M. Pierre Le Barbenchon (Université Rennes 1)

LeBarbenchon.pdf

14 PML methods for mixed hyperbolic-dispersive equations

The classical PML approach is first applied to the linearised Korteweg-de Vries equation. These equations are not always stable, the main obstruction being the classical condition found in the literature on PMLs that we recover in our analysis. We introduce two alternative strategies to design absorbing boundary conditions. We start from studying hyperbolic relaxation of the Korteweg-de Vries equation. In this case, the complete PML equations are not, again, completely stable. However, a version of the PML equations for this system derived without the source term is found to be stable and can absorb outgoing wave. Finally, we consider BBM-Boussinesq system that model bi-directional waves at the surface of an inviscid fluid layer. We show that the PML equations are always stable in this case. We illustrate numerically stability properties of diferents PML models. This talk is based on recent joint work with Christophe Besse, Sergey Gavrilyuk and Pascal Noble.

Orateur: Dr Maria Kazakova (Université Savoie Mont Blanc)

Kazakova_Waves.pdf

10:40 Coffee break

15 On the "projection" structure of the Green-Naghdi equations

Many reduced models of the water wave equations, in particular the Green-Naghdi equations, can be presented as a "projection" of the shallow water equations onto a set of admissible functions. We will see why we use quotation marks for this property, and how to use it to its advantage for many numerical and modeling problems (entropic stability, boundary conditions, balanced, HPC, coupling, dispersion relation).

Orateur: Martin Parisot (Inria)

Parisot.pdf

12:15 Lunch