<html>#2272: GRHydro C2P inconsistency with and without magnetic field evolution
<table style='border-spacing: 1ex 0pt; '>
<tr><td style='text-align:right'> Reporter:</td><td>Lorenzo Sala</td></tr>
<tr><td style='text-align:right'> Status:</td><td>new</td></tr>
<tr><td style='text-align:right'>Milestone:</td><td></td></tr>
<tr><td style='text-align:right'> Version:</td><td>ET_2018_09</td></tr>
<tr><td style='text-align:right'> Type:</td><td>bug</td></tr>
<tr><td style='text-align:right'> Priority:</td><td>major</td></tr>
<tr><td style='text-align:right'>Component:</td><td>EinsteinToolkit thorn</td></tr>
</table>
<p>I have found an inconsistency in simulating static TOVs with tabulated EoS with GRHydro. These TOV are not magnetized, thus evolving with or without magnetic fields evolution should not change the results. Unfortunately, evolving magnetic fields leads to crash. <br />
I’m importing initial data from Lorene/Nrotstar using my thorn ID_Nrotstar, setting up T and Y_e in beta-equilibrium, with uniform initial specific entropy at 1 kb/baryon; this creates a gradient in temperature, from a central ~17MeV to an atmosphere value of 0.01MeV. <br />
Attached you can find plots of maximum rho vs time: if magnetic field is initially set up to 0 and magnetic field evolution is activated, as you can see in the attached .par, <strong>rho gets stuck and the simulation crashes after a while</strong> due to NaN in metric. If I don’t evolve magnetic field (thus, using a different C2P), I don’t have this issue. Attached you can find the .par files, as you can see they only differ in two lines concerning Bvec evolution.<br />
Thank you for your attention,<br />
Lorenzo Sala</p>
<p><table><tr><td>attachment:</td><td><a href="Comparison.png">https://api.bitbucket.org/2.0/repositories/einsteintoolkit/tickets/issues/2272/attachments/Comparison.png</a></td></tr>
<tr><td>attachment:</td><td><a href="tov_LS220_entropyslice_EvolBvec.par">https://api.bitbucket.org/2.0/repositories/einsteintoolkit/tickets/issues/2272/attachments/tov_LS220_entropyslice_EvolBvec.par</a></td></tr>
<tr><td>attachment:</td><td><a href="tov_LS220_entropyslice_NOEvolBvec.par">https://api.bitbucket.org/2.0/repositories/einsteintoolkit/tickets/issues/2272/attachments/tov_LS220_entropyslice_NOEvolBvec.par</a></td></tr>
</table></p>
<p>--<br/>
Ticket URL: <a href='https://bitbucket.org/einsteintoolkit/tickets/issues/2272/grhydro-c2p-inconsistency-with-and-without'>https://bitbucket.org/einsteintoolkit/tickets/issues/2272/grhydro-c2p-inconsistency-with-and-without</a></p>
</html>