[Users] Binary Neutron Star evolution with Gamma-Law EOS

[Petr Tsatsin]

Hello,
I've recently tried to perform simulation of the binary neutron stars
started from Meudon initial data. I've encountered a following problem. I
can run the simulations without any problems when I use polytropic EOS.
However, when I switch to Gamma-Law simulations crushed  soon after start
somewhere in the conversion of primitive quantities to conserved
quantities. The typical log is look like this one:

WARNING[L1,P28] (GRHydro): Con2Prim failed; stopping the code.
WARNING[L1,P28] (GRHydro): Even with mesh refinement, this point is not
restricted from a finer level, so this is really an error
WARNING[L1,P28] (GRHydro):             on carpet reflevel:  2
WARNING[L1,P28] (GRHydro):                   xyz location:
12.00000        0.000000       -2.000000
WARNING[L1,P28] (GRHydro):                         radius:     12.16553
WARNING[L1,P28] (GRHydro): hydro_excision_mask, C2Pfailed:   0   1.0000000
WARNING[L1,P28] (GRHydro):        rho eps press w_lorentz:
0.1047685E-08   0.5280317E-01   0.5532109E-10    1.013257
WARNING[L1,P28] (GRHydro):                     velocities:
-0.3371566E-02   0.1176264      -0.3492353E-01
WARNING[L1,P28] (GRHydro):                       dens tau:
0.2037191E-08   0.1388378E-09
WARNING[L1,P28] (GRHydro):                           scon:
-0.2835004E-10   0.4723716E-09  -0.1051106E-09
WARNING[L1,P28] (GRHydro):                       3-metric:
1.420906      -0.7104738E-01  -0.2077218E-01    1.770598
0.4372303E-01    1.468061
WARNING[L1,P28] (GRHydro):                   lapse, shift:
0.8118283       0.1691328E-01  -0.3824454E-01  -0.5210642E-03
WARNING level 0 in thorn GRHydro processor 28 host
c559-402.stampede.tacc.utexas.edu
  (line 2404 of GRHydro_Con2Prim.F90):
  -> Aborting.
cactus_SBNS:
/work/02566/tsatsin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:314:
int Carpet::Abort(const _cGH *, int): Assertion `0' failed.

The parameter files for "polytropic" runs and "Gamma-Law" are identical
except, obviously, the parameters controlling the EOS.
I was wondering if somebody experienced similar problems before? I would be
very grateful for any help and suggestions.
I also include the some relevant parameters I used for "Gamma-law" EOS
below. May be I'm missing some important parameters or values are not
properly set?

 # Hydro paramters
ActiveThorns = "EOS_Omni"
EOS_Omni::poly_k = 123.6
EOS_Omni::poly_gamma = 2.0
EOS_Omni::gl_k = 123.6
EOS_Omni::gl_gamma = 2.0

ActiveThorns = "GRHydro"

GRHydro::initial_Gamma = 2.0

HydroBase::evolution_method      = "GRHydro"

#GRHydro::riemann_solver            = "Marquina"
GRHydro::riemann_solver     = "HLLE"
GRHydro::GRHydro_eos_type   = "General"
GRHydro::GRHydro_eos_table  = "Ideal_Fluid"
GRHydro::recon_method       = "ppm"
GRHydro::GRHydro_stencil    = 3
GRHydro::bound              = "flat"
GRHydro::rho_abs_min        = 1.e-10
#GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau
#GRHydro::GRHydro    = 10    # gij(6), alpha, beta(3)

GRHydro::c2p_reset_pressure        = "yes"
grhydro::grhydro_c2p_warnlevel         =  0
grhydro::grhydro_c2p_warn_from_reflevel      =  2


ActiveThorns = "GenericFD NewRad"
ActiveThorns = "ML_BSSN ML_BSSN_Helper"
ActiveThorns = "ML_ADMConstraints"


Thank you and Happy Thanksgiving.
Petr.
-- 
Petr Tsatsin
Graduate Student - Department of Physics
Charles E. Schmidt College of Science
Florida Atlantic University
(561)-297-3386
cell: 1-561-299-8048
ptsatin at fau.edu, ptsatsin at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20131126/85332fa0/attachment.html