<div dir="ltr">Hello Benjamin,<div><br></div><div>Yes, we are aware of this bug and will release a patch very soon to correct this issue. Thank you !</div><div><br></div><div>Best,</div><div><br></div><div>Ben.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le ven. 6 sept. 2019 à 11:13, Benjamin Spiers <<a href="mailto:benjamin.spiers@physics.ox.ac.uk">benjamin.spiers@physics.ox.ac.uk</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">





<div lang="EN-GB">
<div class="gmail-m_-3835434501878239766WordSection1">
<p class="MsoNormal">Hi All,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I have recently been trying to use the source term implementation of the Biermann battery in LaserSlab-like simulations in order to study magnetic field generation in long-pulse laser interactions.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I found that with laser powers of 450GW I ran into issues with the timestep in my simulations becoming impractically small (<1fs), which I believe to be because the ions in these simulations reached temperatures of 10-100MeV, and so had
 enormous associated sound speeds. I also found that this issue still arose when I added a multiplicative factor to all updates to U(*_VAR,…) quantities in hy_uhd_biermannSource.F90, even when said factor was set to zero (disabling the source terms altogether).
 However, setting hy_useBiermann = .false. stopped the issue.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I believe that the issue is caused by the calls to `Eos_wrapped(MODE_DENS_EI, blkLimits, blockID)` at the start and end of the main loop in this file. This EOS mode doesn’t properly perform a 3T update, and somehow flip-flopping between
 EOS modes like this seems to cause significant artificial heating of plasma ions. Switching to MODE_DENS_EI_GATHER seems to have fixed the issue for me: magnetic fields are generated as expected and the plasma remains at a sensible temperature with sensibly-sized
 time-steps throughout. <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I believe what is required to solve this issue is to change the two cases of MODE_DENS_EI to the corresponding _GATHER mode, and also to move the<i>
</i>#ifdef FLASH_UHD_3T from the next line to encompass the whole loop, so as to protect 1T simulations from the inverse issue (It already disables field generation in the 1T case so this shouldn’t be a regression).<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Best wishes,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Ben<u></u><u></u></p>
</div>
</div>

</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><span style="background-color:rgb(255,255,255)"><font color="#3d85c6">Benjamin Khiar<br></font></span></div><div><span style="background-color:rgb(255,255,255)"><font color="#3d85c6">Postdoctoral researcher</font></span></div><div><span style="background-color:rgb(255,255,255)"><font color="#3d85c6">FLASH Center for Computational Science</font></span></div><div><span style="background-color:rgb(255,255,255)"><font color="#3d85c6">University of Chicago</font></span></div></div></div></div></div>