<html 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=utf-8">
<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:0in;
font-size:10.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.m-6764176604076760465contentpasted2
{mso-style-name:m_-6764176604076760465contentpasted2;}
span.m-6764176604076760465contentpasted3
{mso-style-name:m_-6764176604076760465contentpasted3;}
span.m-6764176604076760465contentpasted11
{mso-style-name:m_-6764176604076760465contentpasted11;}
span.m-6764176604076760465contentpasted0
{mso-style-name:m_-6764176604076760465contentpasted0;}
span.EmailStyle22
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Thank you for pointing that out Mathieu. There may be a units issue for transport coefficients when using non-default UnitSystem (something other than ‘none’). The transport routines assume that the B-field
has 1/sqrt(4*pi) built-in already. We will look into this more carefully for all transport processes and include any fixes in a patch and/or in a subsequent minor release version.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Jiri,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I also noticed a mistake in the inequality in your email. The paper you referenced indeed shows that kappa_perp / gamma_perp > 3. So the general relation should be:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">kappa_perp > gamma_perp.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;color:black">--</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:#181A1B">Eddie Hansen</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:#181A1B">Research Scientist</span><span style="color:black"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:#181A1B">Flash Center for Computational Science</span><span style="color:black"><o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Mathieu Bailly-Grandvaux <mbaillygrandvaux@eng.ucsd.edu><br>
<b>Date: </b>Monday, May 15, 2023 at 3:51 PM<br>
<b>To: </b>Hansen, Eddie <ehansen@pas.rochester.edu><br>
<b>Cc: </b>Loffelmann, Jiri <loffejir@fjfi.cvut.cz>, flash-users@flash.rochester.edu <flash-users@flash.rochester.edu><br>
<b>Subject: </b>[EXT] Re: [FLASH-USERS] E-MHD Aniso conductivity BUG<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">Hi Eddie,<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">I chime in here about the B-field units, in case it is a relevant point of concern with the new routines for heat transport in the last released version of FLASH. What I used to do is to declare the flag "UnitSystem
= 'cgs'" to avoid the 4pi normalization factor such as B can be expressed directly in Gauss (compared to the default "UnitSystem='none'", for which the 4pi factor is needed). Is it still correct and compatible with the new routines?<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">Best,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">Mathieu<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">On Mon, May 15, 2023, 12:41 PM Hansen, Eddie <<a href="mailto:ehansen@pas.rochester.edu">ehansen@pas.rochester.edu</a>> wrote:<o:p></o:p></span></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">Hello,<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">The ‘chi’ that you are referring to is different between those two routines on purpose, but it’s not because of different B-field units…<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">The fits for the kappa’s from Ji & Held, Phys. Plasmas 20, 042144, 2013 use a collision time tau that is multiplied by an extra factor of Z. So with
chi = omega*tau, the chi‘s in JiHeld/Conductivity_anisoFullState vs DaviesWen/Thermoelectric_fullState are different by a factor of Z.
<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt">The flux limiters for the Thermoelectric effects (e.g., Nernst) are not well-tested yet. They were implemented because we knew they might be needed
in the future. This topic is addressed in Davies, et al., Phys. Plasmas 22, 112703 (2015).<o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;color:black">--</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:#181A1B">Eddie Hansen</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:#181A1B">Research Scientist</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:#181A1B">Flash Center for Computational Science</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt"> <o:p></o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">flash-users <<a href="mailto:flash-users-bounces@flash.rochester.edu" target="_blank">flash-users-bounces@flash.rochester.edu</a>> on behalf of Loffelmann, Jiri <<a href="mailto:loffejir@fjfi.cvut.cz" target="_blank">loffejir@fjfi.cvut.cz</a>><br>
<b>Date: </b>Monday, May 15, 2023 at 1:43 PM<br>
<b>To: </b><a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a> <<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a>><br>
<b>Subject: </b>[FLASH-USERS] E-MHD Aniso conductivity BUG</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">Hello FLASH users,</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">I have two questions for you.</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">#1</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">I noticed an error in the computation of anisotropic thermal conduction coefficients inside the new FLASH version 4.7 and I would like
to check that my belief is correct.</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">In my E-MHD simulations, using all the newly added MHD terms (Seebeck, Nernst, CrossNernst, aniso. thermal cond.), the magnetic advection
velocity corresponding to the Nernst term is in some situations similar or even greater than the advection of temperature due to grad(T) in the direction perpendicular to the magnetic field.</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">Using the formalism from this paper
<a href="https://urldefense.com/v3/__https:/iopscience.iop.org/article/10.1088/1741-4326/ac25c1__;!!Mih3wA!DuALZDn2Aiv1otRU3Bm64LhM3F-Oxfd7veokgy1ikyGQzwuH1akCr1srtyZRzD8LEimISKKkDrSKzSaG-q-LbgQm$" target="_blank">
https://iopscience.iop.org/article/10.1088/1741-4326/ac25c1</a>: </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> gamma_perp >= kappa_perp</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">In addition, in the paper in Figure 2, there is clearly visible, that the ratio of these transport coefficients, <span class="m-6764176604076760465contentpasted2"><span style="background:white">kappa_perp
/ </span></span><span class="m-6764176604076760465contentpasted3"><span style="background:white">gamma_perp</span></span><span class="m-6764176604076760465contentpasted2"><span style="background:white">, </span></span>should be always greater than 3.</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">Therefore, I suspect an implementation error in the computation of the kappa_pepr coefficient.</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">When reading through the two files that compute kappa_perp and gamma_perp coefficients:</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">/FLASH4.7/source/physics/materialProperties/Conductivity/ConductivityMain/JiHeld/Conductivity_anisoFullState.F90</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">/FLASH4.7/source/physics/materialProperties/Thermoelectric/ThermoelectricMain/DaviesWen/Thermoelectric_fullState.F90</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">I noticed that the computation is very similar and truly, for the same values of magnetization ('chi'), it gives the coefficients with
the correct ratio.</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">Finally, I discovered that both codes compute different values of magnetization, with a difference of sqrt(4.0*PI).</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">I believe that in Conductivity_anisoFullState, the magnetic field is in Gaussian units and in Thermoelectric_fullState it is in different
physical units. And since the computation of the magnetization is precisely the same, it differs in the end. I am not really sure why the units are different, as I am not as familiar with all of the code, but either way, I believe the magnetization should
be the same...</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">I propose a quick fix:</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">add line 'xmag = xmag / sqrt(4.0*PI)' after line number 133, in file <span class="m-6764176604076760465contentpasted11"><span style="background:white">Conductivity_anisoFullState.F90.</span></span></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">I am happy to provide further details and the setup of the simulation where the incorrect fraction appears.</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span class="m-6764176604076760465contentpasted0"><span style="font-size:12.0pt;color:black;background:white">2#</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span class="m-6764176604076760465contentpasted0"><span style="font-size:12.0pt;color:black;background:white">Secondly, I would like to ask about the usage of flux limiters for
the Nernst and CrossNernst terms in the E-MHD model inside FLASH.</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span class="m-6764176604076760465contentpasted0"><span style="font-size:12.0pt;color:black;background:white">I would like to use E-MHD for the laser-target problems and I am wondering
what values to use for 'nernstflcoef' and 'crossfieldflcoef'?</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span class="m-6764176604076760465contentpasted0"><span style="font-size:12.0pt;color:black;background:white">Do you have any references or experience with this?</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black">Should I use the same value as for thermal diffusion (diff_eleFlCoef)? When for laser-target problems, it is usually reasonable to set <span style="background:white">diff_eleFlCoef
= </span>0.05-0.1...</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;background:white">
<span style="font-size:12.0pt;color:black">Regards,</span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span class="m-6764176604076760465contentpasted0"><span style="font-size:12.0pt;color:black;background:white">Jiří Löffelmann</span></span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:12.0pt;color:black"> </span><span style="font-size:11.0pt"><o:p></o:p></span></p>
</div>
</div>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt">_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://urldefense.com/v3/__https:/flash.rochester.edu/mailman/listinfo/flash-users__;!!Mih3wA!DuALZDn2Aiv1otRU3Bm64LhM3F-Oxfd7veokgy1ikyGQzwuH1akCr1srtyZRzD8LEimISKKkDrSKzSaG-o0CDY7M$" target="_blank">https://urldefense.com/v3/__https://flash.rochester.edu/mailman/listinfo/flash-users__;!!Mih3wA!DuALZDn2Aiv1otRU3Bm64LhM3F-Oxfd7veokgy1ikyGQzwuH1akCr1srtyZRzD8LEimISKKkDrSKzSaG-o0CDY7M$</a>
<o:p></o:p></span></p>
</blockquote>
</div>
</div>
</body>
</html>