Summary for the current status of the Ablative RT project

  • Jonathan suggested to test the flux and with hydro-off. Ideally we can extend the test to 3-cases

  1. currently we still have problem with the hydro-off test result will be posted in the third part.
  • When doing the test we find a mismatch of for the bottom flux and the equation:
    1. Putting everything in computational units, the bottom flux calculated from Kappa1 and Temperature is 2.32E-4 and the bottom flux converted from Betti's data is 3.48E-4 (5.876e18 W/m2 → 5.876e21 erg/s/cm2 and divided by fluxScale=pScale*velScale). The difference is a factor of 1.5 which is , or gamma7 in the code, or as in Jonathan's blog post.

2.One way of understanding this is that Kappa1~K0/Cv where from Betti, but in AstroBEAR we have

and

  1. One easy fix for this is to multiply both size of the equation by gamma7: Since in the code there was already a gamma7 in the left side of the equation (which is a bug), we just include gamma7 in the Kappa1 as Jonathan mentioned in his blog. So now the Kappa1 becomes

which is 1.5 time larger than before.

  1. Since we still have problem with the hydro-off test, I can only show the effect of this 1.5 factor with bug-buried results
without gamma7
with gamma7
  • Results for flux test: the way I ran the temperature limiting case tests is multiplying a small factor to the Temperature) with hydro on and off.
    1. limiting case with hydro off, multiply by 1e-17

  1. limiting case with hydro on, multiply by 1e-13 — limited by the BCs in BeforeStep, cannot go too small
  1. normal hydro-off test

Attachments (9)

Comments

No comments.