<div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Zhang! <div><br></div><div><span style="font-family:SimSun">>>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 style="font-family:SimSun"><span style="font-family:"Times New Roman""><strong>REQUIRE Diffuse</strong></span></em><span style="font-family:SimSun">? </span><br></div><div><font face="arial, helvetica, sans-serif">The calculation of resistivity is made in source/physics/materialProperties/MagneticResistivity/MagneticResistivityMain/Constant (for constant resistivity). </font></div><div><font face="arial, helvetica, sans-serif">Also, in your config file, have you included a require statament like this?:  </font></div><div><font face="arial, helvetica, sans-serif">   REQUIRES physics/materialProperties/MagneticResistivity/MagneticResistivityMain </font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Best!</font></div><div><span style="font-family:SimSun"><br></span></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El vie., 29 mar. 2019 a las 11:52, zdy0903 (<<a href="mailto:zdy0903@stu.xjtu.edu.cn" target="_blank">zdy0903@stu.xjtu.edu.cn</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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>^{<span style="font-family:SimSun">2</span>}<span style="font-family:SimSun">/s=10</span>^{<span style="font-family:SimSun">-11 </span>}<span style="font-family:SimSun">Ω·m and magnetic viscosity
υ</span>_{<span style="font-family:SimSun">m</span>}<span style="font-family:SimSun">=η·c·c/(4·pi), where c=3·10</span>^{<span style="font-family:SimSun">10 </span>}<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="gmail-m_-2944167736361421641gmail-m_4544061255065082111ke-content-forecolor" style="color:rgb(229,51,51)"> 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²⁰ to 10^-40
</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></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_-2944167736361421641gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font face="monospace, monospace" style="background-color:rgb(234,209,220)"><i><font color="#000000">Dr. Ernesto Zurbriggen</font></i></font><div><font face="monospace, monospace"><i style="background-color:rgb(234,209,220)"><br></i></font><div><span style="background-color:rgb(234,209,220)"><font face="monospace, monospace"><font color="#444444"><b>*</b><i> Consejo Nacional de Investigaciones Científicas y Técnicas </i></font></font><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13"><b>(CONICET)</b></font></i><i style="font-family:monospace,monospace;font-size:12.8px"><font color="#444444">, Instituto de Astronomía Teórica y Experimental </font><font color="#274e13"><b>(IATE)</b>, Argentina</font><font color="#444444">.</font></i></span></div><div><span style="background-color:rgb(234,209,220)"><span style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace"><b>*</b><i> Universidad Nacional de Córdoba </i></span><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13"><b>(UNC)</b></font></i><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">, </i><font face="monospace, monospace" style="font-size:12.8px"><i><font color="#444444">Observatorio Astronómico de Córdoba </font></i></font><i style="font-size:12.8px;font-family:monospace,monospace"><font color="#274e13"><b>(OAC)</b>, Argentina</font></i><i style="font-size:12.8px;color:rgb(68,68,68);font-family:monospace,monospace">.</i><br></span></div></div><div><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px;background-color:rgb(234,209,220)"><br></i></div><div><i style="color:rgb(68,68,68);font-family:monospace,monospace;font-size:12.8px;background-color:rgb(234,209,220)">Teléfono: +54 0351 4331064-5, interno 222.</i></div><div><span style="font-size:12.8px"><i><span style="background-color:rgb(234,209,220)"><font color="#0000ff"><u>Laprida 854, Bº Observatorio (CP 5000), Córdoba, Argentina</u></font><font color="#444444">.</font></span><br></i></span><font color="#444444"><br><br><br></font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>