<div dir="ltr"><div>Hello,<br></div>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: <br>
<br clear="all"><div><div>WARNING[L1,P28] (GRHydro): Con2Prim failed; stopping the code.<br>WARNING[L1,P28] (GRHydro): Even with mesh refinement, this point is not restricted from a finer level, so this is really an error<br>
WARNING[L1,P28] (GRHydro): on carpet reflevel: 2<br>WARNING[L1,P28] (GRHydro): xyz location: 12.00000 0.000000 -2.000000<br>WARNING[L1,P28] (GRHydro): radius: 12.16553<br>
WARNING[L1,P28] (GRHydro): hydro_excision_mask, C2Pfailed: 0 1.0000000<br>WARNING[L1,P28] (GRHydro): rho eps press w_lorentz: 0.1047685E-08 0.5280317E-01 0.5532109E-10 1.013257<br>WARNING[L1,P28] (GRHydro): velocities: -0.3371566E-02 0.1176264 -0.3492353E-01<br>
WARNING[L1,P28] (GRHydro): dens tau: 0.2037191E-08 0.1388378E-09<br>WARNING[L1,P28] (GRHydro): scon: -0.2835004E-10 0.4723716E-09 -0.1051106E-09<br>WARNING[L1,P28] (GRHydro): 3-metric: 1.420906 -0.7104738E-01 -0.2077218E-01 1.770598 0.4372303E-01 1.468061<br>
WARNING[L1,P28] (GRHydro): lapse, shift: 0.8118283 0.1691328E-01 -0.3824454E-01 -0.5210642E-03<br>WARNING level 0 in thorn GRHydro processor 28 host <a href="http://c559-402.stampede.tacc.utexas.edu">c559-402.stampede.tacc.utexas.edu</a><br>
(line 2404 of GRHydro_Con2Prim.F90):<br> -> Aborting.<br>cactus_SBNS: /work/02566/tsatsin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:314: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.<br><br>
</div><div>The parameter files for "polytropic" runs and "Gamma-Law" are identical except, obviously, the parameters controlling the EOS.<br></div><div>I was wondering if somebody experienced similar problems before? I would be very grateful for any help and suggestions.<br>
</div><div>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?<br><br> # Hydro paramters<br>ActiveThorns = "EOS_Omni"<br>
EOS_Omni::poly_k = 123.6<br>EOS_Omni::poly_gamma = 2.0<br>EOS_Omni::gl_k = 123.6<br>EOS_Omni::gl_gamma = 2.0<br><br>ActiveThorns = "GRHydro"<br><br>GRHydro::initial_Gamma = 2.0<br><br>HydroBase::evolution_method = "GRHydro"<br>
<br>#GRHydro::riemann_solver = "Marquina"<br>GRHydro::riemann_solver = "HLLE"<br>GRHydro::GRHydro_eos_type = "General"<br>GRHydro::GRHydro_eos_table = "Ideal_Fluid"<br>
GRHydro::recon_method = "ppm"<br>GRHydro::GRHydro_stencil = 3<br>GRHydro::bound = "flat"<br>GRHydro::rho_abs_min = 1.e-10<br>#GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau<br>
#GRHydro::GRHydro = 10 # gij(6), alpha, beta(3)<br><br>GRHydro::c2p_reset_pressure = "yes"<br>grhydro::grhydro_c2p_warnlevel = 0<br>grhydro::grhydro_c2p_warn_from_reflevel = 2<br><br>
<br>ActiveThorns = "GenericFD NewRad"<br>ActiveThorns = "ML_BSSN ML_BSSN_Helper"<br>ActiveThorns = "ML_ADMConstraints"<br><br><br></div><div>Thank you and Happy Thanksgiving.<br></div><div>Petr.<br>
</div><div>-- <br><div dir="ltr"><div>Petr Tsatsin</div><div>Graduate Student - Department of Physics</div><div>Charles E. Schmidt College of Science</div><div>Florida Atlantic University</div><div>(561)-297-3386</div><div>
cell: 1-561-299-8048</div><div><a href="mailto:ptsatin@fau.edu" target="_blank">ptsatin@fau.edu</a>, <a href="mailto:ptsatsin@gmail.com" target="_blank">ptsatsin@gmail.com</a></div><div></div><div></div><div></div><div></div>
<div></div></div>
</div></div></div>