A Lagrangian method for solving the spherical shallow water equations using power diagrams

2025-08-11

AI Summary

Lagrangian methods for global atmospheric simulation suffer from excessive solution smoothing, limiting their competitiveness against mainstream Eulerian approaches. Method: This paper proposes a novel Lagrangian framework based on spherical power diagrams for solving the spherical shallow-water equations. Fluid is discretized using particles; an adaptive spherical power diagram mesh is dynamically constructed to ensure exact mass conservation. For the first time in a Lagrangian setting, numerical stability is achieved without artificial viscosity. Momentum and energy conservation are enforced via semi-discrete optimal transport and a semi-implicit time integrator. Results: Experiments demonstrate that billion-scale meshes can be partitioned within two minutes, with conservation accuracy matching the state-of-the-art Lagrangian methods. The approach significantly improves both physical fidelity and computational efficiency.

Technology Category

• Planning, Routing, and Scheduling: Optimization of Spatio-temporal Systems
• Search and Optimization: Sampling/Simulation-based Search
• Constraint Satisfaction and Optimization: Distributed CSP/Optimization

Application Category

• Graph Algorithms and Modeling for the Web: Efficient manipulation of static and dynamic Web-related graphs
• Systems and Infrastructure for Web, Mobile and WoT: Experiences and lessons learnt from Web-based algorithms and system deployments
• Search and Retrieval-Augmented AI: Efficiency and scalability of Web search engines

Abstract

Numerical simulations of the air in the atmosphere and water in the oceans are essential for numerical weather prediction. The state-of-the-art for performing these fluid simulations relies on an Eulerian viewpoint, in which the fluid domain is discretized into a mesh, and the governing equations describe the fluid motion as it passes through each cell of the mesh. However, it is unclear whether a Lagrangian viewpoint, in which the fluid is discretized by a collection of particles, can outperform Eulerian simulations in global atmospheric simulations. To date, Lagrangian approaches have shown promise, but tend to produce smoother solutions. In this work, a new Lagrangian method is developed to simulate the atmosphere in which particles are represented with spherical power cells. We introduce an efficient algorithm for computing these cells which are then used to discretize the spherical shallow water equations. Mass conservation is enforced by solving a semi-discrete optimal transport problem and a semi-implicit time stepping procedure is used to advance the solution in time. We note that, in contrast to previous work, artificial viscosity is not needed to stabilize the simulation. The performance of the spherical Voronoi diagram calculation is first assessed, which shows that spherical Voronoi diagrams of 100 million sites can be computed in under 2 minutes on a single machine. The new simulation method is then evaluated on standard benchmark test cases, which shows that momentum and energy conservation of this new method is comparable to the latest Lagrangian approach for simulating the spherical shallow water equations.

Problem

• Develops Lagrangian method for spherical shallow water equations
• Compares Lagrangian and Eulerian atmospheric simulation performance
• Ensures mass conservation without artificial viscosity

Innovation

• Lagrangian method with spherical power cells
• Efficient algorithm for spherical Voronoi diagrams
• Semi-discrete optimal transport for mass conservation