Numerical study of the sharp stratification limit towards bilayer models
Fuente:
arXiv
Saved in:
| Main Author: | |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866915867029143552 |
|---|---|
| author | Fradin, Théo |
| author_facet | Fradin, Théo |
| contents | In the study of oceanic flows at the geophysical scale, the phenomenon of density stratification plays a central role in the dynamics of the system. Two categories of mathematical models are commonly used to describe the role played by the density stratification: on the one hand, continuously stratified models - such as the stratified Euler equations in a strip, considered in the present article - offer an accurate description of vertical effects, but come with a high level of complexity, both at the theoretical and numerical levels. On the other hand, bilayer models approximate the stratification by a piecewise constant profile. In the latter case, the main point is to study the evolution of the free interface between both layers, which leads to a substantially simplified model. In the present article, we compare both approaches in the framework of the linearized stratified Euler equations around density profiles that are close to piecewise constant profiles, and prove the convergence towards the bilayer Euler equations. However, in the presence of a shear flow, bilayer models have a range of validity limited by the presence of Kelvin-Helmholtz instabilities. In this case, we use a suitable normal modes decomposition to compute numerically the dispersion relation of this linearized model, and provide numerical evidence that the Kelvin-Helmholtz instabilities limit the applicability of two widely used bilayer models, namely the bilayer Euler equations and the bilayer shallow-water equations. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2603_15287 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Numerical study of the sharp stratification limit towards bilayer models Fradin, Théo Fluid Dynamics Analysis of PDEs 35Q35 (Primary), 76B70, 76M22, 86A05 (Secondary) In the study of oceanic flows at the geophysical scale, the phenomenon of density stratification plays a central role in the dynamics of the system. Two categories of mathematical models are commonly used to describe the role played by the density stratification: on the one hand, continuously stratified models - such as the stratified Euler equations in a strip, considered in the present article - offer an accurate description of vertical effects, but come with a high level of complexity, both at the theoretical and numerical levels. On the other hand, bilayer models approximate the stratification by a piecewise constant profile. In the latter case, the main point is to study the evolution of the free interface between both layers, which leads to a substantially simplified model. In the present article, we compare both approaches in the framework of the linearized stratified Euler equations around density profiles that are close to piecewise constant profiles, and prove the convergence towards the bilayer Euler equations. However, in the presence of a shear flow, bilayer models have a range of validity limited by the presence of Kelvin-Helmholtz instabilities. In this case, we use a suitable normal modes decomposition to compute numerically the dispersion relation of this linearized model, and provide numerical evidence that the Kelvin-Helmholtz instabilities limit the applicability of two widely used bilayer models, namely the bilayer Euler equations and the bilayer shallow-water equations. |
| title | Numerical study of the sharp stratification limit towards bilayer models |
| topic | Fluid Dynamics Analysis of PDEs 35Q35 (Primary), 76B70, 76M22, 86A05 (Secondary) |
| url | https://arxiv.org/abs/2603.15287 |