PhD defence announcement
The PhD defence of Adrien Garinet (LEGOS) will take place on Tuesday, December 10th at 9:30 am (UTC+1)* in Coriolis room, Observatoire Midi-Pyrénées, where he will be presenting his work : “Evidence and reduction of spurious numerical mixing in a numerical model including internal tides : a case study on the South-East Asian seas”.
Figure extraite de Garinet, A., Herrmann, M., Marsaleix, P., and Pénicaud, J.: Spurious numerical mixing under strong tidal forcing: a case study in the south-east Asian seas using the Symphonie model (v3.1.2), Geosci. Model Dev., 17, 6967–6986, https://doi.org/10.5194/gmd-17-6967-2024, 2024.
The defense will be held in english, and the jury will be composed of :
Guillaume ROULLET, rapporteur, UBO
Knut KLINGBEIL, rapporteur, IOW
Florian LEMARIÉ, examinateur, INRIA
Julien LE SOMMER, examinateur, IGE
Aurore Voldoire, examinatrice, CNRM
Marine Herrmann, directrice de thèse, IRD
Patrick Marsaleix, co-directeur de thèse, CNRS
Abstract:
Because they attempt to represent the complexity of continuous processes on discrete grids, computer simulations used to study the ocean have a fundamental tendency to generate numerical inaccuracies that must be controlled in order for the output to be realistic. Among these, spurious numerical mixing between layers of different densities is particularly troublesome. Indeed, physical mixing in this direction, though of major importance, is quite subtle, still poorly understood, and can easily be overwhelmed by numerical errors. In numerical models that use fixed levels to represent the vertical physics of the ocean, transport across these levels – and errors associated with their discretization – is the main source of such spurious cross-isopycnal mixing.
This work focuses first on showing that since these errors scale with velocities, explicitly resolving the strong vertical displacements generated by internal tides propagating in the interior of the ocean can lead to intense levels of such spurious mixing if care is not taken in how vertical advection is implemented. Using the Symphonie ocean model, and taking the simulation of the South-East Asian Seas as a realistic case study, we show that while internal tides are essential for a correct representation of the water masses in the region, numerical errors generated by using an overly diffusive vertical advection scheme can severely bias the output.
Investigating the choice of a better scheme for handling vertical transport, we then recognize that the reduction of numerical mixing cannot be properly achieved without taking into account the properties of the time-stepping algorithm. We thus propose a new way to design advection schemes aiming at counteracting its spurious effect. Overall, this work offers a deep look into the way vertical advection is handled in the fixed coordinate model and suggests ways to improve it.