Difference between revisions of "ET Workshop Fall 2011 EOS"
(→EOS) |
(→EOS) |
||
| Line 15: | Line 15: | ||
* Josh: implementing should not be a big problem | * Josh: implementing should not be a big problem | ||
* Josh: Some functions are not be fully implemented for the non-polytropic EOS_Omni code | * Josh: Some functions are not be fully implemented for the non-polytropic EOS_Omni code | ||
| + | == MHD Con2Prim == | ||
| + | To perform Conservative to Primitive variable conversion in MHD for a generic EOS, we need to use code that is not dependent on the adiabatic index, since it will not be constant. Some notes describing a possible way to implement a more generic Newton-Raphson routine have been uploaded, along with some untested code to perform this task. | ||
| + | |||
| + | Notes: [[Image:Con2prim_eosomni.pdf]] | ||
| + | |||
| + | Code: [[Image:GRHydro_Con2PrimM_pt_EOSOmni.c]] | ||
Revision as of 14:39, 31 October 2011
EOS
Discussion about which EOS interfaces Cactus has/had:
- EOS_Base
- EOSG_Base
- EOS_Omni
Who uses what:
- EOS_Omni: ET
- EOSG_Base: "nobody" -> to be archived
- EOS_Base: GaTech code, Whisky_Exp -> to remain, but will be dropped from the ET thornlist
- some version of Whisky has "own" interface to tabulated EOS
MHD
- MHD-Con2Prim assumes polytrope EOS
- Code in Wiki: is general, works for polytrope EOS, but is not tested otherwise
- Josh: implementing should not be a big problem
- Josh: Some functions are not be fully implemented for the non-polytropic EOS_Omni code
MHD Con2Prim
To perform Conservative to Primitive variable conversion in MHD for a generic EOS, we need to use code that is not dependent on the adiabatic index, since it will not be constant. Some notes describing a possible way to implement a more generic Newton-Raphson routine have been uploaded, along with some untested code to perform this task.
Notes: File:Con2prim eosomni.pdf
