Dear all,
I found Heister’s model and change it with ASPECT-version2.4.0.
Here is heister’s paper High accuracy mantle convection simulation through modern numerical methods – II: realistic models and problems | Geophysical Journal International | Oxford Academic, and it’s input file paper-aspect-methods-2-data/application/global_convection_3d.prm at master · tjhei/paper-aspect-methods-2-data · GitHub
I found it’s hard to slove after 30+Myrs, the time step become about 3000 yrs after 30+Myrs, it seems that this model wouldn’t end after 37600 CPU hours which is described in his paper(even won’t be end after 3700 time steps, because I run it about 1800 time steps only 39Myrs). I even used the “project q1 only viscosity” compared to his model, it should be more faster. it seems hard to solve after the plume were generated.
The addition files are my input file and logfile.
heister3d.prm (6.3 KB)
log.txt (2.8 MB)