Coupling flow and CPO

The coupled evolution of flow and CPO can be solved by joining different components of specfab together with a momentum balance solver (e.g., Stokes solver):

That is, specfab can model:

• CPO evolution given bulk stress, strain-rate and temperature fields (green piece)
• Viscous anisotropy induced by the CPO (orange piece)
• Bulk anisotropic rheologies given the local viscous anisotropy (blue piece)

Elmer/Ice

See Lilien et al. (2023) for ice flow modelling.

FEniCS

See Rathmann et al. (2024) for olivine modelling.