Thieulot - Computational Infrastructure for Geodynamics

ASPECT: FROM
BENCHMARKING TO SUBDUCTION
C. Thieulot, A. Glerum, M. Fraters, W. Spakman
ASPECT @ THE UU
•
ASPECT on Cartesius supercomputer
•
ASPECT: visco-plastic thermal convection
•
ASPECT: Aegean subduction (A. Glerum)
•
ASPECT: Caribbean region (M. Fraters)
ASPECT ON CARTESIUS
•
it is the Dutch supercomputer
https://surfsara.nl/nl/systems/cartesius/description
ASPECT ON CARTESIUS
•
Pilot study to •
carry out scaling tests
•
identify (& fix) bottlenecks
•
With Thomas Geenen @
•
if good scaling results, one can ask for more cpu-hours
https://surfsara.nl/nl/systems/cartesius/description
ASPECT ON CARTESIUS
•
realistic test case: the 3D slab detachment paper
sticky air
lithosphere (power-law, n=4)
mantle (power-law, n=3)
0
5
x
0
100
m
k
0
0x66
ASPECT ON CARTESIUS
PRELIMINARY
~3 Mdofs
256
ASPECT ON CARTESIUS
•
Status: delay •
MPI versions compatibility issues (impi, xmpi)
•
(restart, visualisation, …)
•
thanks to Thomas for help
VISCO-PLASTIC BENCHMARK
submitted March, 3rd.
VISCO-PLASTIC BENCHMARK
T=0
free slip on sides
T=1
VISCO-PLASTIC BENCHMARK
Measurements:
Avg. Temp.
Nusselt Nb.
Vrms
Work
Visc. Dissip.
Many codes (FVM/FEM):
- YACC - Plaatjes
- CHIC
- GAIA
- StreamV
- StagYY
- FEniCS
- Fluidity
- MC3D
- ELEFANT
- ASPECT
VISCO-PLASTIC BENCHMARK
Cases:
Case 1
VISCO-PLASTIC BENCHMARK
Case 1 stagnant lid
VISCO-PLASTIC BENCHMARK
Case 5a Bifurcation test
periodic
mobile periodic stagnant
AEGEAN SUBDUCTION
AEGEAN SUBDUCTION
AEGEAN SUBDUCTION
AEGEAN SUBDUCTION
Simplified test case of instantaneous modelling
subduction
wedge
No
r th
Sub
duc
ti
on z
o
ne
-
L
350
0x1
700
x70
kinematic b.c. + open b.c.
Ana
toli
an
faul
t
nt
a
ev
0km
lt
u
fa
-
-
AEGEAN SUBDUCTION
Casper Pranger, MSc Thesis @ UU
“Numerical modelling of the instantaneous subduction dynamics of the Banda Arc region”
oceanic crust
➜ meshing tool in Mathematica
wedge
slab
AEGEAN SUBDUCTION
Workflow
Struct. Geology
Mantle Tomography
+
Meshing tool
Triangulated interfaces between domains
(Point coordinates + Connectivity)
Ascii file read by ASPECT
Anne’s plugin
At startup, ci computed for each node
+
Use geometrical info to set temperature
(in progress)
AEGEAN SUBDUCTION
3D slab
Mesh points of triangulated surfaces
Proof of concept
N
AEGEAN SUBDUCTION
Composition isocontours
crust
lithosphere
slab