LYNX (Lithosphere dYnamics Numerical toolboX) is another multiphysics modelling solution developed within the group of Basin Modelling at GFZ. Lynx is also based on the flexible, object-oriented numerical framework MOOSE (developed at the Idaho National Laboratories).

Lynx is a novel numerical simulator for modelling thermo-poromechanical coupled processes driving the deformation dynamics of the lithosphere. The formulation adopted in Lynx relies on an efficient implementation of a thermodynamically consistent visco-elasto-plastic rheology with anisotropic porous-visco-plastic damage feedback. The main target is to capture the multiphysics coupling responsible for semi-brittle and semi-ductile behaviour of porous rocks as also relevant to strain localization and faulting processes. More information on the governing equations, their derivation and their implementation together with a list of synthetic and real case applications can be found in two publications by Jacquey and Cacace (2019, a,b).

User group
International user community from geosciences


  • Object-orientation: flexible modular structure within easy to be extended modules by the user
  • Geometric agnosticism: 1D/2D/3D Finite Elements from the user required
  • Hybrid parallelism: multi-threading and MPI
  • Proved scalability on HPC architectures - JUWELS cluster Module at JSC

Realistic physics-based rheological description of lithosphere deformation dynamics based on:

  • Explicit incorporation of the lithosphere visco-elasto-plastic rheology including nonlinear feedback effects from the energetics of the system
  • its extension to account for time-dependent brittle behavior via an overstress (viscoplastic) formulation
  • Thermodynamically consistent formulation of semi-brittle semi-ductile deformation modes - brittle deformation via damage mechanics and for ductile deformation via a rate-dependent viscoplastic formulation
  • Poro-damage feedback via dynamic porosity (volumetric mechanical response)
  • Implicit and efficient numerical implementation within a limited amount of internal iterations
  • Use of Automatic Differentation techniques to compute the full Jacobian contribution of the system matrix

LYNX is open sourced under the GNU GENERAL PUBLIC LICENSE v3 and can be obtained upon request by contacting one of the authors.

back to top of main content