<span style="font-family:SimSun;">
<p>
<span style="font-family:SimSun;">Hi all!</span>
</p>
<p>
<span style="font-family:SimSun;"> </span>
</p>
<p>
<span style="font-family:SimSun;">I met some problems when trying to add the resistivity η to the simulation
of <strong>magnetoHD/NohCylindrical(user manual 32.2.4)</strong>,</span>
</p>
<p>
<span style="font-family:SimSun;"> </span>
</p>
<p>
<span style="font-family:SimSun;">First, the units of resistivity troubles me. In the CGS UnitSystem,
whether setting </span><i><span style="font-family:SimSun;"> <strong>resistivity = 1 </strong></span></i><span style="font-family:SimSun;">in </span><i><span style="font-family:SimSun;">“flash.par” </span></i><span style="font-family:SimSun;">means set magnetic resistivity η=1
cm</span><sup><span style="font-family:SimSun;">2</span></sup><span style="font-family:SimSun;">/s=10</span><sup><span style="font-family:SimSun;">-11 </span></sup><span style="font-family:SimSun;">Ω·m and magnetic viscosity
υ</span><sub><span style="font-family:SimSun;">m</span></sub><span style="font-family:SimSun;">=η·c·c/(4·pi), where c=3·10</span><sup><span style="font-family:SimSun;">10 </span></sup><span style="font-family:SimSun;">cm/s?</span>
</p>
<p>
<span style="font-family:SimSun;"> </span>
</p>
<p>
<span style="font-family:SimSun;">Second, the calculation about resistivity troubles me. Looking through
the kernels and helpers of Diffuse Unit, I can’t find any calculation about resistivity.
Is it in the kernels and helpers of <strong>Hydro Unit</strong>? I can’t find the calculation about <strong>magnetic diffusion</strong> in the Hydro Unit. Is it needed to </span><em><span style="font-family:Times New Roman;"><strong>REQUIRE Diffuse</strong></span></em><span style="font-family:SimSun;">?</span>
</p>
<p>
<span style="font-family:SimSun;"> </span>
</p>
<p>
<span style="font-family:SimSun;">Finally, When I set these parameter in </span><i><span style="font-family:SimSun;">“flash.par”</span></i><span style="font-family:SimSun;"> as below and setup with </span><i><span style="font-family:Times New Roman;"><strong>+usm</strong></span></i>
</p>
<blockquote style="padding:0px;margin-top:0px;margin-right:0px;margin-bottom:0px;" dir="ltr">
<p>
<i><span style="font-family:Times New Roman;"><strong># Units
(none, SI or CGS)</strong></span></i>
</p>
<p>
<i><span style="font-family:Times New Roman;"><strong>UnitSystem = "CGS"</strong></span></i>
</p>
<p>
<i><span style="font-family:Times New Roman;"><strong>##
Resistive MHD parameters ##</strong></span></i>
</p>
<p>
<i><span style="font-family:Times New Roman;"><strong>useDiffuse =
.true.</strong></span></i>
</p>
<p>
<i><span style="font-family:Times New Roman;"><strong> dt_diff_factor
= 0.8</strong></span></i>
</p>
<p>
<i><strong><span style="font-family:Times New Roman;">useMagneticResistivity = .true.</span></strong></i>
</p>
<p>
<i><strong><span style="font-family:Times New Roman;"><span class="ke-content-forecolor" style="color:#E53333;"> resistivity = 0.0</span> # magnetic Resistivity</span></strong></i>
</p>
</blockquote>
<p>
mpirun stopped with the same message, even if setting a large range of <i>resistivity</i>, from 10<sup>20</sup> to 10<sup>-40</sup>
</p>
<p>
<b><i> DRIVER_ABORT:
[Hydro]: Computed dt is not positive! Aborting!</i></b>
</p>
<p>
<b><i> dtCheck= 0.0000000000000000</i></b><b></b>
</p>
<p>
<br>
</p>
<p>
How to solve these problem? Really need someone to help!
</p>
<p>
</p>
<p>
Zhang
</p>
</span><br>