Context Navigation

Posts for the month of May 2022

CE

Computer time

Code is now running on Cori, in addition to Anvil
Enough time for 1 run on Anvil (need to finish quota by June 30), and maybe another run on Cori
Plan to apply for more time on Cori

Reduction of ambient pressure and density

Method 1 from last post: Starting from t=0 with low ambient pressure and density

Abandoned this method
- Eventually get Hypre error for tabular EOS run
- Need to set MinDensity = ambient density, which is very inefficient
- Get artifacts in ambient medium that slow code and force extra unwanted refinement

Method 3 from last post: Restarting from frame 86 of old runs and modifying ambient

Got this to work well with ideal gas EOS run
- Reduced ambient density and pressure gradually over ~1-2 frames (0.25-0.50 days)
- When done gradually, can retain momentum-conserving self-gravity
- Works well when done using either conservative (modify density, momentum, energy) or primitive (modify density and pressure) forms

Does not quite work yet for tabular EOS run
- If modify density, momentum and energy, then at some point pressure and temperature decrease to very low values, leading to artifacts and unwanted extra refinement
- If modify density and pressure, then code eventually crashes due to protections detecting NANs (energy)
- If modify density and momentum alone, code eventually crashes
- If modify density alone, code slows down catastrophically

Discussion

Solutions for reducing ambient in tabular EOS run
Strategy for paper

Posted: 3 years ago
Author: Luke
Categories: (none)
Comments (0)

Download in other formats:

RSS Feed