[Commits] [svn:einsteintoolkit] EOS_Omni/trunk/doc/ (Rev. 106)

rhaas at tapir.caltech.edu rhaas at tapir.caltech.edu
Wed Mar 12 22:02:12 CDT 2014


User: rhaas
Date: 2014/03/12 10:02 PM

Modified:
 /trunk/doc/
  documentation.tex

Log:
 EOS_Omni: update documentation
 
 From: Roland Haas <rhaas at tapir.caltech.edu>

File Changes:

Directory: /trunk/doc/
======================

File [modified]: documentation.tex
Delta lines: +16 -3
===================================================================
--- trunk/doc/documentation.tex	2014-03-13 03:02:09 UTC (rev 105)
+++ trunk/doc/documentation.tex	2014-03-13 03:02:11 UTC (rev 106)
@@ -113,9 +113,10 @@
   interface and implements multiple analytic EOS, and also provides table
   reader and interpolation routines for finite-temperature
   microphysical EOS available from
-    \url{http://www.stellarcollapse.org}. Currently, the implemented
-  analytic EOS are the polytropic EOS, the gamma-law EOS, and a hybrid EOS
-  consisting of a 2-piece piecewise-polytrope with an a thermal,
+    \url{http://www.stellarcollapse.org}. In addition cold and barotropic
+  tabulated EOS are provided for use in initial data thorns. Currently, the
+  implemented analytic EOS are the polytropic EOS, the gamma-law EOS, and a
+  hybrid EOS consisting of a 2-piece piecewise-polytrope with a thermal,
   gamma-law component.
 \end{abstract}
 
@@ -173,6 +174,7 @@
       for stellar core collapse simulations.
     \item \texttt{eoskey = 4}: Finite-temperature microphysical EOS
     \item \texttt{eoskey = 5}: Cold tabulated EOS with Gamma-Law
+    \item \texttt{eoskey = 6}: Tabulated barotropic EOS
   \end{itemize}
 \item \texttt{havetemp} determines whether the EOS is to be called as
   a function of $(\rho,\epsilon,Y_e)$ (\texttt{havetemp = 0}), or as a
@@ -211,9 +213,20 @@
 
 \begin{verbatim}
 eos_omni::poly_gamma 
+eos_omni::poly_gamma_initial
 eos_omni::poly_k 
 \end{verbatim}
 
+The only non-obvious parameter here is \texttt{poly\_gamma\_initial}.
+In most simulations it should not be set.
+In simulations that are run with a different adiabatic index than
+what was used to set up the initial data, \texttt{poly\_gamma} should
+be the evolution value, and \texttt{poly\_gamma\_initial} should be the
+initial data value. \texttt{EOS\_Omni} then rescales \texttt{poly\_k} such
+that the \emph{cgs} value of \texttt{poly\_k} is the same for initial data and
+evolution. Since the units of \texttt{poly\_k}
+($[K] = \frac{\rm N}{\rm m^2} \left(\frac{\rm m^3}{\rm kg}\right)^{\Gamma}$)
+depend on \texttt{poly\_gamma} this is not fully trivial.
 
 Check \texttt{param.ccl} for parameters for the other EOS\@.
 



More information about the Commits mailing list