julesg/uw3project repository overview

⁠underworld3

Building

Refer to the Dockerfile for uw3 build instructions.

For development, building inplace will prob be preferable. Remove any existing installations, then run.

python setup.py build_ext --inplace

For in place usage, you will need to set an appropriate PYTHONPATH.

Milestones

• Spherical stokes
• Buoyancy drive stokes
• Advection diffusion
• High Ra, constant viscosity
• Highly temp-dep viscosity

Objectives

Ingredients in achieving the above

[T⁠] Topology & Meshing

• spherical annulus - https://github.com/julesghub/cubie⁠
• Cartesian
• Different element types (at least Linear / Quadratic & Hex, Tet)

[D⁠] Disc

• Cont Galerkin
• Disc Galerkin
• Semi-lagrangian
• Free-slip BC on surface

[P⁠] Physics

• Stokes-Boussinesq
• Temp-dep rheology
• Buoyancy driven convection
• Non-linear viscosity (Jacobian ?) and yielding in particular
• Viscoelasticity
• Energy equation, resolve bdry layers
• kermit the 🐸

[S⁠] Solvers

• Block Stokes solvers
• Semi-lagrangian
• TS (address this later)

PIC for composition

• Viscosity, buoyancy, ...
• Nearest neighbour (k-d tree ? 🌳 )
• 2D - L2 projection into FEM space (Petsc shall provide)
• 3D - L2 projection into FEM space (Petsc shall provide but not in 3D)

[O1⁠ O2⁠] Output

• HDF5 -> XDMF -> Paraview
• LavaVu

[V⁠] Exact solutions

• MMS
• Analytical
  • https://www.solid-earth-discuss.net/se-2017-71/se-2017-71.pdf⁠ -https://www.researchgate.net/publication/304784132_Benchmark_solutions_for_Stokes_flows_in_cylindrical_and_spherical_geometry⁠

Tasks

• Investigate cython for functions

• Petsc4py compatibility

• alias .data and .array

• rejig a few models top down ...

• audit the petsc / petsc4py features corresponding to uw ones

• Can we get early into TS (TS-lite)

• (Build a gLucifer equivalent for the DM - as generic as possible)

• Evaluate Julian's C-based prototype