Convection project update 09/30/2019

New Work

Interpolated the OPAL EOS with the provided interpolation function. The result agrees with MESA in the region applicable to OPAL EOS.

Pressure
https://www.pas.rochester.edu/~yishengtu/research_files/CEE_gp_meeting_09302019/comparison2.png https://www.pas.rochester.edu/~yishengtu/research_files/CEE_gp_meeting_09302019/comparison1.png Internal energy
https://www.pas.rochester.edu/~yishengtu/research_files/CEE_gp_meeting_09302019/comparison3.png https://www.pas.rochester.edu/~yishengtu/research_files/CEE_gp_meeting_09302019/comparison4.png

Next step

  1. Now I am using a bash script that coordinates a Python script and two fortran scripts to do the interpolation, which is very non-efficient (computation time ~5 hours). Modifying the fortran script to get a one-step function may save a lot of time (It looks like MESA only use less then 1 second to do the calculations).
  2. Inverting the tables to get other quantities as a function of density and internal energy.

Comments

No comments.