<div dir="ltr">Hello Petr,<div><br></div><div>are you using the </div><div><br></div><div>GRHydro::reconstruct_Wv = yes<br></div><div><br></div><div>parameter? </div><div>I need to use it in my accretion disk runs when using the Ideal_Fluid EOS...Using the polytrope this parameter wasn't needed if I recall correctly..</div>
<div><br></div><div>hope this helps..</div><div><br></div><div>best wishes,</div><div><br></div><div>Vassili</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Nov 26, 2013 at 6:00 PM, Petr Tsatsin <span dir="ltr"><<a href="mailto:ptsatsin@my.fau.edu" target="_blank">ptsatsin@my.fau.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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" target="_blank">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.<span class="HOEnZb"><font color="#888888"><br>
</font></span></div><span class="HOEnZb"><font color="#888888"><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><a href="tel:%28561%29-297-3386" value="+15612973386" target="_blank">(561)-297-3386</a></div>
<div>
cell: <a href="tel:1-561-299-8048" value="+15612998048" target="_blank">1-561-299-8048</a></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></font></span></div></div>
<br>_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br></blockquote></div><br></div>