geodynamics / aspect

A parallel, extensible finite element code to simulate convection in both 2D and 3D models.
https://aspect.geodynamics.org/
Other
228 stars 237 forks source link

A talk for subduction model #6071

Open 3Faker3 opened 1 month ago

3Faker3 commented 1 month ago

Hi everyone,

I am trying to run a subduction model but have been facing some difficulties. Despite changing various parameters like the solver settings, discretization, mesh refinement, compositional fields, initial composition and temperature model, boundary velocity (both free subduction and fixed velocity), visco-plastic parameters, and everything else I could think of, the model still fails and stops at step 1.

For a successful subduction model run, what are the critical aspects I should focus on? What could I be missing? I also noticed that someone mentioned this earlier.https://github.com/geodynamics/aspect/issues/1828

I would appreciate any guidance you can provide. Thank you in advance!

(version:2.6.0 pre deal.ii: 9.5.2) The uploaded zip file contains all the required files.

subdutcion.zip

gassmoeller commented 1 month ago

Hi @3Faker3, Your model stops at a curious point, since it looks like it is stopping before even starting the timestep. I also do not see anything obviously wrong in your parameter file. Do you have any error message on screen? Also please rerun the model in debug mode, it is generally advisable to start a new model in debug mode first, because it will often give you much more detailed error messages (see here.

Generally in cases like yours we recommend to start from a simple model that is already in ASPECT, make sure that runs properly, and then step-by-step modify the model in the direction you want to go. So after every modification (change of material model, introduction of compositions, setting initial temperature ,...) you can check if the model still works. You could start from this cookbook for example and then modify it step-by-step towards your intended setup. Or if you want to start from a very simple case, use this one.

Hope that helps.

3Faker3 commented 1 month ago

Hi @3Faker3, Your model stops at a curious point, since it looks like it is stopping before even starting the timestep. I also do not see anything obviously wrong in your parameter file. Do you have any error message on screen? Also please rerun the model in debug mode, it is generally advisable to start a new model in debug mode first, because it will often give you much more detailed error messages (see here.

Generally in cases like yours we recommend to start from a simple model that is already in ASPECT, make sure that runs properly, and then step-by-step modify the model in the direction you want to go. So after every modification (change of material model, introduction of compositions, setting initial temperature ,...) you can check if the model still works. You could start from this cookbook for example and then modify it step-by-step towards your intended setup. Or if you want to start from a very simple case, use this one.

Hope that helps. hi, Thank you first. The log.txt contains all the screen output, and the two cookbooks you mentioned have successfully run on my supercomputer. I am so confused that I still can't figure out what is wrong.