A Semi-Discrete Optimal Transport Scheme for the 3D Incompressible Semi-Geostrophic Equations

Théo Lavier
2024-11-01
Abstract:We describe a mesh-free three-dimensional (3D) numerical scheme for solving the incompressible semi-geostrophic equations, based on semi-discrete optimal transport techniques. These results generalise previous two-dimensional (2D) implementations. The optimal transport methods we adopt are known for their structural preservation and energy conservation qualities and achieve an excellent level of efficiency and numerical energy-conservation. We use this scheme to generate numerical simulations of an important benchmark problem. To our knowledge, this is the first fully 3D simulation of this particular cyclone, evidencing the model's applicability to atmospheric and oceanic phenomena and offering a novel, robust tool for meteorological and oceanographic modelling.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a grid - free three - dimensional numerical scheme for solving the incompressible semi - geostrophic equations. Specifically, the research aims to: 1. **Extension to three - dimensional simulation**: Previous studies have mainly focused on semi - geostrophic flow simulation in two - dimensional cases, while this paper realizes a fully three - dimensional numerical simulation for the first time. This provides a new and powerful tool for modeling atmospheric and oceanic phenomena. 2. **Utilization of the optimal transport method**: The semi - discrete optimal transport technique is adopted to construct the numerical scheme. This method is known for its structure - preserving and energy - conserving properties, and can ensure numerical energy conservation while computing efficiently. 3. **Simulation of isolated semi - geostrophic cyclones**: Use the newly developed numerical scheme to simulate the formation and development process of isolated semi - geostrophic cyclones. This is the first complete three - dimensional simulation of this specific cyclone, demonstrating the application potential of the model in meteorological and oceanic modeling. 4. **Overcoming the limitations of traditional methods**: Through the optimal transport method, the problem that the transformation from geostrophic space to physical space fails in a finite time is solved, ensuring that the transformation is always effective throughout the simulation process. ### Mathematical background The semi - geostrophic system (SG system) mentioned in the paper is a set of equations derived under the assumptions of hydrostatic equilibrium and geostrophic equilibrium, which is suitable for describing the motion of shallow, rotation - dominated fluids. The mathematical form of this system is as follows: \[ \partial_t\alpha_t+\text{div}(\alpha_t v[\alpha_t]) = 0, \] where, \[ v[\alpha_t]=J(\text{id}-T^{-1}),\quad J=\begin{pmatrix} 0&-1&0\\ 1&0&0\\ 0&0&0 \end{pmatrix}, \] $\alpha_t$ is a probability measure - valued mapping, $T:X\rightarrow Y$ is the optimal transport mapping from $X$ to $Y$, and satisfies: \[ T = \arg\min_{T:X\rightarrow Y,T^\#1_X=\alpha_t}\int_X\|x - T(x)\|^2dx. \] ### Initial conditions The initial conditions include a symmetric baroclinic jet, combined with a uniform baroclinic shear component, which is controlled by the parameter $A\in\mathbb{R}$, and the non - dimensionalized pressure $\Phi(x)$ is defined as: \[ \Phi(x)=\frac{1}{2}\left(\arctan\left(\frac{x_2}{1 + x_3}\right)-\arctan\left(\frac{x_2}{1 - x_3}\right)\right)-0.12x_2x_3-\frac{1}{2}A(x_2^2 - x_3^2). \] In addition, the temperature is perturbed by the perturbation function $h(x_1,x_2)$ to generate initial instability. ### Simulation results Through numerical simulation, the author shows the formation and development process of cyclones and verifies the effectiveness and accuracy of the method. In particular, in long - term simulations, this method can keep the total energy stable, which reflects its advantages in modeling complex atmospheric and oceanic phenomena. ### Conclusion This research has not only made important progress in numerical methods, but also provided a new perspective for understanding atmospheric and oceanic phenomena. Through three - dimensional simulation, the behavior of cyclones and other large - scale weather systems can be more comprehensively understood, providing strong support for weather forecasting and climate research.