Auteur/autrice : davidlannes

Mathocean du 21/01/19

(Bordeaux, salle 2 de l’IMB)

  • 14h00-15h00 Nicolas Barral (Imperial College London) : Framework pour des simulations côtières avec adaptation de maillage anisotrope
  • 15h00-16h00 Kevin Martins (EPOC): Wave transformation in the surf zone: new insights from 2D LiDAR scanners

Nicolas Barral: Framework pour des simulations côtières avec adaptation de maillage  anisotrope. La plupart des simulations numériques pour l’ingénierie marine requièrent la modélisation de larges zones géographiques, alors que la quantité d’intérêt est très localisée (par exemple les systèmes de génération d’énergie marémotrice, ou l’étude de l’impact de tsunamis sur des zones habitées ou à risque). La taille des domaines considérés, combinée à la haute précision requise conduisent généralement à des simulations très coûteuses en temps de calcul.
L’adaptation de maillage a démontré depuis une vingtaine d’années son efficacité pour améliorer la précision numérique tout en réduisant le temps de calcul en permettant un contrôle fin de la répartition des degrés de liberté. En particulier, l’adaptation de maillage anisotrope tire parti de la flexibilité des maillages non structurés, et permet d’optimiser à la fois la taille et l’orientation des éléments du maillage. La génération de ces maillages adaptés est basée sur des champs de métriques Riemanniennes, eux-même dérivés d’un estimateur d’erreur, a priori ou a posteriori. La construction de cet estimateur d’erreur, propre au modèle étudié, est essentiel au processus d’adaptation. Dans cette présentation, on considère des métriques « multi-échelles » basée sur la Hessienne des champs de solution, qui sont capables de capturer les différentes échelles des dynamiques océaniques et des structures d’ingénierie étudiées.
L’efficacité de ces estimateurs d’erreur peut être considérablement améliorée en considérant des approches dites « goal oriented », qui sont en cours de développement dans notre framework: en considérant l’adjoint de la solution dans le modèle d’erreur, seuls les phénomènes physiques ayant une influence sur une certaine quantité cible sont maillées et résolues avec précision.
Dans cet exposé, je présenterai l’implémentation d’un framework permettant de réaliser de telles simulations adaptatives dans notre code de simulations côtières Thetis, s’appuyant sur son intégration sous-jacente avec PETSc. On considéra des applications pour la modélisation de la propagation de tsunamis, ainsi que la modélisation de structures de générations d’énergie marémotrice (barrages ou turbines).

Kevin Martins: Wave transformation in the surf zone: new insights from 2D LiDAR scanners: Over the last few decades, remote-sensing technology became an essential tool for improving our understanding of nearshore processes. Very recently, the application of lidar scanners (Light Detection And Ranging) to study swash zone hydrodynamics and morphodynamics was made possible thanks to the initiatives of few researchers. In this talk, I will present my PhD work, in which we explore the use of 2D lidar scanners to monitor the time-varying surface elevation of breaking waves in the surf zone. The surf zone constitutes one of the most challenging environments where to deploy instruments, due to the energetic wave conditions often found there. Hence, obtaining direct measurements of wave profiles at a high sampling rate represents a unique opportunity to better understand wave transformation in the surf zone.
In the first part, I will describe the working principles of 2D lidar scanners, and its potential applications to study nearshore hydro- and morphodynamics. In a second part, we use an innovative laboratory dataset to validate a CFD model of breaking wave in a wave flume at prototype scale (BARDEXII experiments). This numerical model is used to study and quantify the influence of wave reflection, generated at the inner surf/swash zone boundary, on surf zone hydrodynamics at various temporal scales. Then, the main field experiments performed during my PhD is presented. An array of three lidar scanners was deployed from a nearshore pier at Saltburn-by-the-Sea (UK) to monitor waves from the shoaling area and up to the runup limit. New insights on individual wave properties at breaking and on energy dissipation rates in the inner surf zone were obtained and will be presented. I will conclude by briefly describing my current work and how lidar data can be used to improve the parameterizations of physical processes in the surf zone.

MathOcean du 10/12/18

(Toulouse, 14h à l’IMT)

Eliane Becache (INRIA et ENSTA Paris Tech) : minicours sur les PML

Eliane Becache : La méthode des couches PMLs est très populaire pour borner le domaine de calcul dans les problèmes d’ondes et a été  appliquée  à de nombreux modèles  (acoustique, électromagnétisme, élasticité…)  transitoires et fréquentiels.
Cette présentation  commencera par rappeler la définition des  PMLs et leurs propriétés.  Nous nous intéresserons ensuite à leur analyse mathématique, en particulier à des questions de stabilité. Nous présenterons d’abord l’analyse pour des modèles d’ondes transitoires hyperboliques linéaires non-dispersifs  pour lesquels on peut exhiber une condition nécessaire de stabilité. Ces résultats permettent de comprendre pourquoi les PMLs peuvent devenir instables dans certains milieux anisotropes et donnent des pistes  pour construire des PMLs stabilisées. Nous aborderons ensuite des questions de convergence des PMLs pour des problèmes fréquentiels dans des guides d’ondes, qui mettent en évidence des différences étonnantes par rapport au cas transitoire. Nous finirons par des résultats  plus récents concernant des modèles linéaires dispersifs ainsi que par les problèmes qui restent ouverts.

MathOcean du 26/11/2018

(Bordeaux, salle 2 de l’IMB)

  • 14h00-15h00 Anne Mangeney (Institut de Physique du Globe de Paris): Recovering ice mass loss in Greenland from seismic data and mechanical modelling.
  • 15h00-16h00 Lisl Weynans et David Lannes (IMB): Generating boundary conditions and floating objects in Boussinesq systems.

Anne Mangeney (Institut de Physique du Globe de Paris) Recovering ice mass loss in Greenland from seismic data and mechanical modelling (joint work with Amandine Sergeant, Pauline Bonnet, Vladislav Yastrebov, Olivier Castelnau, Eléonore Stutzmann, Jean-Paul Montagner, and Patrick Quetey). Quantification of ice mass loss is a key issue to better understand the impact of climate change and in particular to constrain regional and global circulation models. The Greenland ice sheet is losing mass as a result of both increased surface melting and runoff and increased ice discharge from marine-terminating outlet glaciers. Ice discharge at these glaciers includes submarine ice melting and iceberg production. Calving of icebergs generate glacial earthquakes of magnitudes up to 5 that can be recorded at 100’s of km from the source. These waves result from the forces applied by km-scale icebergs against the terminus face, as they slowly capsize. Long-period inversion of the recorded seismic signal makes it possible to calculate these forces. Former studies failed to recover iceberg volume from the earthquake magnitude or from the maximum force due to the complex physical processes involved. We recently showed that numerical modelling of the iceberg capsize taking into account the interplay between these processes help constraining seismic inversion, making it possible to recover icebergs volume. This coupled seismic and modelling approach allows to monitor ice-mass loss from iceberg capsize in Greenland over a 20-year period of recorded glacial seismicity and to investigate its relation to climate change, ocean temperature, or glacier dynamics.

Lisl Weynans et David Lannes: Generating boundary conditions and floating objects in Boussinesq systems. Boussinesq systems are dispersive perturbations of the (hyperbolic) nonlinear shallow water equations, and are widely used for applications in coastal oceanography and more recently for the description of wave-structure interactions (eg for marine renewable energies). We shall consider two different phenomenons in this talk, in which we exhibit a dispersive boundary layer that plays an important role:
– Generating boundary conditions: for numerical simulations, one quantity (eg the surface elevation) is known at the entrance of the domain, and one is let to solve an initial boundary value problem. The theory is known for hyperbolic systems but there is no general theory in the presence of dispersion. Several techniques have been proposed, which are often difficult to implement and/or time consuming. We present here a simple method based on the analysis of the dispersive boundary layer, and that can be implemented without extra computational cost.
– Nonlinear dispersive waves with a floating object (joint work with D. Bresch and G. Métivier). We consider here a the interaction of wave described by a Boussinesq system with a floating object. We show that the problem can be reduced to a simple transmission problem (and even to an ODE) and show how to solve it by analysing the dispersive boundary layer.

A nonlinear weakly dispersive method for recovering the elevation of irrotational surface waves from pressure measurements

P. Bonneton, D. Lannes, K. Martins, H. Michallet, Coastal Engineering 138 (2018), 1–8. Download

We present the derivation of a nonlinear weakly dispersive formula to reconstruct, from pressure measurements, the surface elevation of nonlinear waves propagating in shallow water. The formula is simple and easy to use as it is local in time and only involves first and second order time derivatives of the measured pressure. This novel approach is evaluated on laboratory and field data of shoaling waves near the breaking point. Unlike linear methods, the nonlinear formula is able to reproduce at the individual wave scale the peaked and skewed shape of nonlinear waves close to the breaking point. Improvements in the frequency domain are also observed as the new method is able to accurately predict surface wave elevation spectra over four harmonics. The nonlinear weakly dispersive formula derived in this paper represents an economic and easy to use alternative to direct wave elevation measurement methods (e.g. acoustic surface tracking and LiDAR scanning).

Reconstruction of water depth time series of a group of waves observed
in the field.