[FLASH-USERS] Magnetic resistivity SpitzerHighZ

[Yingchao Lu]

Hi FLASH users,

It seems the magnetic resistivity for Spitzer HighZ in FLASH is three
orders of magnitude smaller than those found in the literature.

The code in
source/physics/materialProperties/MagneticResistivity/MagneticResistivityMain/SpitzerHighZ/MagneticResistivity_fullState.F90
calculates the electron-ion equilibrium time
res_ieEquilTime(zbar, abar, tele, tion, nion, eqtime)

And the code in
source/physics/materialProperties/MagneticResistivity/MagneticResistivityMain/SpitzerHighZ/res_ieEquilTime.F90
  eqtime = res_ieTimeCoef * &
       3.0 * res_boltz**1.5 / (8.0 * sqrt(2*PI) * res_qele**4) * &
       (mion * tele + res_mele * tion)**1.5 / &
       ( sqrt(mion*res_mele) * zbar**2 * nion * ll)

When Ti=Te, we can make some approximation, shown in the attached file.

There are several references for Spitzer magnetic resistivity and they are
consistent with each other.
http://people.hao.ucar.edu/judge/homepage/PHSX515/fall2012/Braginskii1965.pdf
https://www.cfa.harvard.edu/~namurphy/Lectures/Ay253_04_BeyondIdealMHD.pdf
Physics Reports 417 (2005) 1–209

I think the (7) to (10) in the image should be correct. Could any one
double check?

Thanks,
Yingchao
[image: Screen Shot 2018-11-20 at 5.26.52 PM.png]
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20181120/4982d164/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2018-11-20 at 5.26.52 PM.png
Type: image/png
Size: 103621 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20181120/4982d164/attachment.png>