This package has some tools for finite volume methods in JAX. We use the staggered grids (i.e. Arakawa Grids) to define variables an the nodes, faces, and centers. This package includes all of the key operations necessary for interacting with these variables including: 1) differences 2) interpolations, and 3) reconstructions. In addition, we provide a way to mask boundaries which are consistent with Arakawa Grids.
Operators.
This package has simple operators that are useful for calculating differences.
These include the standard difference and laplacian operators.
It also has some geostrophic specific operators like the geostrophic_gradient, the divergence, and vorticity.
Masks. This package includes some useful masks that are necessary for interacting with the Arakawa C-/D-Grids.
Interpolation.
This package includes some simple interpolation schemes for moving between grid spaces.
It has linear interpolation schemes based on the mean which include the arithmetic, geometric, harmonic, and quadratic.
We also have some general purpose interpolation schemes for the Cartesian grid including the linear interpolation scheme and splines.
Reconstructions. This package includes some reconstruction methods to calculate the fluxes typically found within the advection terms. For example, the flux found in the vector invariant formulation for the Shallow Water equations (example). Another example, the flux found in the advection term for the QG equations (example).
We can install it directly through pip
pip install git+https://github.com/jejjohnson/finitevolXWe also use poetry for the development environment.
git clone https://github.com/jejjohnson/finitevolX.git
cd finitevolX
conda create -n finitevolx python=3.11 poetry
poetry installWe can also use conda for the development environment.
git clone https://github.com/jejjohnson/finitevolX.git
cd finitevolX
conda env create -f environment.yaml
conda activate finitevolxAll of these examples go for correctness and readability.
Linear Shallow Water Model.
In the scripts located in scripts/swm_linear.py, we have an example using a linear shallow water model.
This script was taken from dionhaefner/shallow-water/shallow_water_simple.py.
It has been rewritten with the finitevolx helper functions.
NonLinear Shallow Water Model.
In the scripts located in scripts/swm.py, we have an example using a nonlinear shallow water model.
This script was taken from dionhaefner/shallow-water/shallow_water_nonlinear.py.
This uses the vector invariant formulation which involves the potential vorticity and kinetic energy.
It has been rewritten with the finitevolx and fieldx helper functions.
Note: there is no lateral viscosity term. It only uses an implicit diffusion reconstruction method, i.e. improved weno method, which prevents the simulation from blowing up. However, a 3 pt scheme is sufficient. The 5pt schemes results in weird oscillations along the boundaries.
To run the scripts, run poetry and ensure that it installs with
poetry install --with expor alternatively we can simply use the conda environment from above.
Software
- kernex - differentiable stencils
- FiniteDiffX - finite difference tools in JAX
- PyFVTool - Finite Volume Tool in Python
- jaxinterp2d - CartesianGrid interpolator for JAX
- ndimsplinejax - SplineGrid interpolator for JAX
Algorithms
- Thiry et al, 2023 | MQGeometry 1.0 - the WENO reconstructions applied to the multilayer Quasi-Geostrophic equations, Arakawa Grid masks
- Roullet & Gaillard, 2021 | pyRSW - the WENO reconstructions applied to the shallow water equations