<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
Hello FLASH users,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
I have two questions for you.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
#1</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
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.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
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.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1">
Using the formalism from this paper <a href="https://iopscience.iop.org/article/10.1088/1741-4326/ac25c1" id="LPlnk552811">
https://iopscience.iop.org/article/10.1088/1741-4326/ac25c1</a>: </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1">
<span> </span>gamma_perp >= kappa_perp</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1">
In addition, in the paper in Figure 2, there is clearly visible, that the ratio of these transport coefficients, <span style="display: inline !important; background-color: rgb(255, 255, 255);" class="ContentPasted2">kappa_perp /
<span style="display: inline !important; background-color: rgb(255, 255, 255);" class="ContentPasted3">
gamma_perp</span>, </span>should be always greater than 3.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1">
Therefore, I suspect an implementation error in the computation of the kappa_pepr coefficient.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1">
When reading through the two files that compute kappa_perp and gamma_perp coefficients:</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted6 ContentPasted7">
/FLASH4.7/source/physics/materialProperties/Conductivity/ConductivityMain/JiHeld/Conductivity_anisoFullState.F90<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5">
/FLASH4.7/source/physics/materialProperties/Thermoelectric/ThermoelectricMain/DaviesWen/Thermoelectric_fullState.F90<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5">
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.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5 ContentPasted8">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5 ContentPasted8">
Finally, I discovered that both codes compute different values of magnetization, with a difference of sqrt(4.0*PI).</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5 ContentPasted8 ContentPasted9">
I believe that in <span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">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...<br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5 ContentPasted8 ContentPasted9">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5 ContentPasted8 ContentPasted9">
I propose a quick fix:</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5 ContentPasted8 ContentPasted9 ContentPasted10">
add line 'xmag = xmag / sqrt(4.0*PI)' after line number 133, in file <span style="display: inline !important; background-color: rgb(255, 255, 255);" class="ContentPasted11">Conductivity_anisoFullState.F90.</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted1 ContentPasted4 ContentPasted5 ContentPasted8 ContentPasted9">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; color: black;">I am happy to provide further details and the setup of the simulation where the incorrect fraction appears.</span><br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0">2#</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0">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></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0 ContentPasted12">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></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0 ContentPasted12">Do you have any references or experience with this?</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted13">
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="display: inline !important; background-color: rgb(255, 255, 255);">diff_eleFlCoef = </span>0.05-0.1...</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0"><br>
</span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<div style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0">
Regards,</div>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0">Jiří Löffelmann</span><br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<span style="font-size: 12pt; margin: 0px; color: black; background-color: white;" class="ContentPasted0"><br>
</span></div>
</body>
</html>