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
© Copyright 2024