Difference between revisions of "ET Workshop Fall 2011 EOS"

From Einstein Toolkit Documentation
Jump to: navigation, search
(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
 +
 +
Christian:
 +
* To fix: single NS case for hot EOS
 +
* handle multiple composition evolutions, requires EOS tweaks
 +
* more important working EOS + MHD
 +
 
== MHD Con2Prim ==
 
== 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.
 
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.

Revision as of 16:06, 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

Christian:

  • To fix: single NS case for hot EOS
  • handle multiple composition evolutions, requires EOS tweaks
  • more important working EOS + MHD

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

Code: File:GRHydro Con2PrimM pt EOSOmni.c