Meeting Update 1118

  • ART problem, I think Baowei is getting a mismatch between the bottom temperatures of his profile and the given bottom heat flux (the bottom temperature should give a flux that is approx close to q0 given by Riccardo).

  • Wrote a routine that reads from Jacques' Argon and SiO2 ionization fraction Tables, and wired it to the resistivity routine. Current features:
  1. read the Table.txt, from cell density and temperature, infer the ionization fraction Z by averaging the lower and higher bounds (should be better with Jonathan's extrapolation routine, I haven't yet received it btw).
  2. look at the tracer field and determine which table to use (which material has higher fraction). It may be worthwhile to do another average.
  3. from Z, calculate desired physics quantity on that cell using Spitzer equation (resistivity and heat conduction).

I also modified heat conduction routine to do the same thing, so that only ionized electrons can conduct heat, and only along B field directions (the anisotropicity is determined by gyroradius). So I think we are more sophisticatd than before with the explicit routines, except they don't do subcycles.

Preliminary run (beta = 10):

new nluf run

  • Got contacted by UIUC comp astro group, read some of their papers (they do GRMHD, I found the magnetized disk simulations particularly interesting).

  • TSF paper (short version) should be out to Adam and Eric this week. The isothermal runs got a memory error, see below.

​http://www.pas.rochester.edu/~shuleli/TSF.e3864249

  • Started high-res magnetized TSF (Krumholz) runs on bluestreak. Until the end of this year, I'll be more dedicated to writing (short letter for the TSF, resistive clumps paper, and also potentially the longer TSF paper) while waiting for these magnetized runs.

Attachments (1)

Download all attachments as: .zip

Comments

No comments.