A conservative, discontinuous Galerkin, tracer transport scheme using compatible finite elements
This addresses mass conservation issues in geophysical fluid modeling for climate and weather simulations, but it is incremental as it builds on existing compatible finite element methods.
The paper tackled the problem of conserving mass in scalar tracer transport within geophysical fluid models by solving a conservative transport equation for tracer density instead of the advective equation for mixing ratio, ensuring mass conservation in both continuous and discrete equations with a discontinuous Galerkin scheme. Results demonstrated accurate mass conservation in tests like terminator toy chemistry and a moist rising bubble.
This paper outlines a conservative transport scheme for scalar tracers within a compatible finite element model for geophysical fluid equations. Instead of using the advective transport equation for a mixing ratio, a conservative transport equation is solved for the tracer density of the mixing ratio multiplied by the dry density. This ensures mass conservation in the continuous equations, which can be preserved in the discrete equations with a discontinuous Galerkin transport scheme. Our method is designed to work for two placements of the mixing ratio in a Charney-Phillips vertical staggering: either co-located with the dry density or vertically staggered from it. The new scheme is designed to conserve the tracer density and ensure consistency by maintaining a constant mixing ratio. Additionally, a mass-conserving limiter is developed to ensure non-negativity in the co-located configuration. Tests with terminator toy chemistry and a moist rising bubble show the use of the new transport scheme with physics terms and its ability to accurately model mass conservation of moisture species in a dynamical core setup.