<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi All,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></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).<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Best wishes,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Ben<o:p></o:p></p>
</div>
</body>
</html>