clfryer
commented
3 years ago My setup codes typically are reworked constantly to set up different conditions. The setup code I gave you starts an explosion at the edge of a red giant star so that I can model the shock progression down a steep density gradient. I model the very outer layer of the shock starting at a radius of rb (and a mass of 0.001Msun).
uint is the internal energy. In an ideal gas, uint is propotional to rho and Temperatue. To keep temperature constant, I set it proportional to 1/rho.
ycc is the abundance fraction. In some calculations, this is followed. Not so in these but I still read them in. They should probably be set to zero in the code, but my compiler does that for undefined values and I was doing the setup quickly so didn't worry about it. It is read in 1dsimp, but not used.
abar is the average mass... Here I set it to the value for pure hydrogen.
chris
From: Sam Slocum @.***> Sent: Wednesday, August 11, 2021 8:05:31 AM To: SNeModel-Dev/SNeRadHydro Cc: Fryer, Christopher Lee; Assign Subject: [EXTERNAL] [SNeModel-Dev/SNeRadHydro] Basic Questions (#1)
What is rb and what is its relation to the radius?
Is x[] the radius at that grid point? It looks like it is changing depending on how much mass is remaining to be positioned in the grid. What is the motivation for this?
Looking at v[] it seems that there is velocity being placed toward the center, but not further out. Is that to represent the flow of material into a stationary surrounding material?
uint[]? What is it? It seems like a weighted value of some magic numbers I couldn't easily identify.
abar[]?
ycc[][]? yccave[][]?
— You are receiving this because you were assigned. Reply to this email directly, view it on GitHubhttps://github.com/SNeModel-Dev/SNeRadHydro/issues/1, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AC367UX23QLWY3XYZVHNLYDT4J7SXANCNFSM5B6S4KUA. Triage notifications on the go with GitHub Mobile for iOShttps://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Androidhttps://play.google.com/store/apps/details?id=com.github.android&utm_campaign=notification-email.
What is rb and what is its relation to the radius?
Is x[] the radius at that grid point?
It looks like it is changing depending on how much mass is remaining to be positioned in the grid. What is the motivation for this?
Looking at v[] it seems that there is velocity being placed toward the center, but not further out. Is that to represent the flow of material into a stationary surrounding material?
uint[]? What is it? It seems like a weighted value of some magic numbers I couldn't easily identify.
abar[]?
ycc[][]? yccave[][]?