Fix evolve base state with new time-stepping scheme

AMReX-Astro / MAESTROeX

A C++ low Mach number stellar hydrodynamics code

https://amrex-astro.github.io/MAESTROeX/

BSD 3-Clause "New" or "Revised" License

40 stars 22 forks source link

Fix evolve base state with new time-stepping scheme #27

Open doreenfan opened 6 years ago

doreenfan commented 6 years ago

Issue can be reproduced using inputs_3d_C.128.avgbase.dr for uneven base state spacing, and inputs_3d_C.128.avgbase.dx for even base state spacing. Both tests fail after approx. 20 time steps, near the base_cutoff_radius. Backtrace files show that error is due to either a failure to converge during MAC projection, or eos failure to compute enthalpy outside base_cutoff_radius in UpdateScal().

ajnonaka commented 5 years ago

Mostly fixed; still using the original make_w0 for the velocity splitting. If this approach would be simplified it would be nice.