search

underworldcode / underworld2

underworld2: A parallel, particle-in-cell, finite element code for Geodynamics.
http://www.underworldcode.org/
Other
168 stars 58 forks source link

Get no result when using both Stokes and AdvectionDiffusion systems in update() #364

Closed jerryjohn42 closed 5 years ago

jerryjohn42 commented 5 years ago

Hi all, I am a new user of underworld code, and now I am trying to produce a model with a subducting plate and mantle convection. The model should be affected by both density and temperature difference, so I use Stokes and AdvectionDiffusion systems to get the results. It works while just running advector.integrate(dt) in update(). But while I run advector.integrate(dt) and advDiff.integrate(dt) together in update(), the code has been running for long but I get no result, not even figure of step equal to zero. Can anyone help me to check the code and point it out whether it is my code error or because the limit of my computer calculation ability. Thanks. Jerry subduction slab test.zip

jmansour commented 5 years ago

Hi Jerry

I was able to run your code fine using UW2.7. Note that I had to lower the resolution and reduce the maxSteps in order to obtain an output in reasonable time.

I recommend updating to 2.7 if you aren't already using it. Note that I don't believe this is what made a difference. Also try doing as I've done, and lowering your resolution significantly (I ran 32x8) and setting maxSteps=1 to ensure you are getting the correct execution. Once it all appears to be operating as you would expect, then start increasing the resolution again.

Let us know how you progress. John

jerryjohn42 commented 5 years ago

Hi Jerry

I was able to run your code fine using UW2.7. Note that I had to lower the resolution and reduce the maxSteps in order to obtain an output in reasonable time.

I recommend updating to 2.7 if you aren't already using it. Note that I don't believe this is what made a difference. Also try doing as I've done, and lowering your resolution significantly (I ran 32x8) and setting maxSteps=1 to ensure you are getting the correct execution. Once it all appears to be operating as you would expect, then start increasing the resolution again.

Let us know how you progress. John

Hi John, Thanks for your reply, I lowered the resolution as your suggestion and get the expected results. Next time I will use a lower resolution in my PC to test first and try higher resolution in server. Thanks for your help. : ) Jerry