Data driven non-equilibrium moist phase exchanges for atmospheric convection within a discontinuous Galerkin model of the compressible Euler equations
Abstract
A discrete formulation of the multiphase moist compressible Euler equations is presented for which the energy and tracer variance conserving dynamics is coupled to the non-equilibrium moist phase exchanges in a thermodynamically consistent manner. A neural network is then trained to learn the mass exchanges between vapour, liquid and ice phases in atmospheric convection. The network is trained on convection resolving output from a regional configuration of the LFRic model with a multi-moment microphysics parameterisation (CASIM). The loss function for learning these phase exchanges is formulated under the assumptions of thermal and mechanical equilibrium (same temperature and pressure for all phases), and chemical dis-equilibrium (different Gibbs free energies for all phases). The network outputs determine the exchanges of vapour, liquid and ice so as to conserve mass, and the resulting change in entropy is determined from the network outputs so as to conserve energy. The neural network is implemented in a thermodynamically consistent manner within the discontinuous Galerkin moist compressible Euler model in order to simulate the formation of three-phase clouds for convection at sub-km resolution. The results are compared to those from a physics based representation of three-phase moist processes at thermodynamic equilibrium, as well as to an idealised configuration of LFRic model.
Related
- A Fully Discrete Local Discontinuous Galerkin Method for Quasilinear Stochastic Convection-Diffusion-Type Equations
- Entropy-Stable and Well-Balanced Discontinuous Galerkin Methods for the Compressible Euler Equations in Vector-Invariant Form
- A Massively Parallel Hybridizable Discontinuous Galerkin Solver for Direct Numerical Simulation of Compressible Flows on GPUs
- Well-balanced second-order approximation of the compressible atmospheric Euler equations
- GPU-Accelerated Nested Schur-Complement Solvers for Atmospheric Flow Simulations Using Horizontally-Explicit, Vertically-Implicit Schemes and Discontinuous Galerkin Methods
BibTeX
@misc{lee2026data,
title = {{Data driven non-equilibrium moist phase exchanges for atmospheric convection within a discontinuous Galerkin model of the compressible Euler equations}},
author = {David Lee and Kieran Ricardo and Junwei Lyu},
howpublished = {arXiv},
year = {2026},
month = aug,
doi = {10.48550/arxiv.2607.13360},
eprint = {2607.13360},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.13360},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Data driven non-equilibrium moist phase exchanges for atmospheric convection within a discontinuous Galerkin model of the compressible Euler equations
AU - David Lee
AU - Kieran Ricardo
AU - Junwei Lyu
PY - 2026
DA - 2026/08/28
PB - arXiv
DO - 10.48550/arxiv.2607.13360
UR - https://arxiv.org/abs/2607.13360
AB - A discrete formulation of the multiphase moist compressible Euler equations is presented for which the energy and tracer variance conserving dynamics is coupled to the non-equilibrium moist phase exchanges in a thermodynamically consistent manner. A neural network is then trained to learn the mass exchanges between vapour, liquid and ice phases in atmospheric convection. The network is trained on convection resolving output from a regional configuration of the LFRic model with a multi-moment microphysics parameterisation (CASIM). The loss function for learning these phase exchanges is formulated under the assumptions of thermal and mechanical equilibrium (same temperature and pressure for all phases), and chemical dis-equilibrium (different Gibbs free energies for all phases). The network outputs determine the exchanges of vapour, liquid and ice so as to conserve mass, and the resulting change in entropy is determined from the network outputs so as to conserve energy. The neural network is implemented in a thermodynamically consistent manner within the discontinuous Galerkin moist compressible Euler model in order to simulate the formation of three-phase clouds for convection at sub-km resolution. The results are compared to those from a physics based representation of three-phase moist processes at thermodynamic equilibrium, as well as to an idealised configuration of LFRic model.
ER -