Posts for the month of January 2020

Equation of State reference

I download all these papers & resources and put them on PAS server.

EOS papers

MESA EOS

OPAL:

SCVH:

PTEH:

HELM:

Other useful resources

Jet updates 2020-01/27

Jet runs and status:

Model(SS) # Stampede run # Restart Mdot/M_Edd M2/Msun a_i/Rsun Ambient ptcl. buffer status
0 2 No 1000 1 49 High no completed (173 frames)
3 6 No 1 1 49 High no frame 0-72
4 7 Yes 1 1 a(frame 55) High no problem with restarting
5 8 Yes 1 1 a(frame 130) High no problem with restarting
6 9 No 1 ½ 49 High no frame 0-72
7 - Yes 1000 1 49 High no submit after debugging
8 - Yes 1000 1 49 High no submit after debugging
1 - No 1 1 73.5 Low yes maybe?
2 (no Jet) - No 0 1 73.5 Low yes maybe?

Update 1/27

Radiation pressure paper

Revision submitted. Hopefully it'll be accepted soon and I can move the data off BlueHive.

MHD

Ramp vs No-Ramp

http://www.pas.rochester.edu/~adebrech/MHD/Run2_ramp_comp.png

Slight differences. But the non-ramping version is working, so that's probably the better way to go. Working on some resolution tests now.

Charge Exchange

http://www.pas.rochester.edu/~adebrech/ChargeExchange/HD209458b_charge_exchange0000.png

I have 10 more frames of this, but no room on BlueHive currently. The low-pressure wind is working (have 5 frames or so), but I believe I've pretty much used up my allotted portion of the allocation on Stampede2.

I will start the postprocessing on BlueHive, so that we can start looking at synthetic observations.

AMR line transfer

This may become useful for the higher field cases. In order to decrease the discontinuity between the planet and the ambient, we need higher resolution around the planet - but we can force the rest of the simulation to stay at e.g. level 4 maximum while refining the planet to a much higher level. However, this would be practically impossible to run in the "projected to max level" version of the line transfer routine. I'm working on updating and merging this into all of my branches.

Currently waiting on a timing test for the buffered nonblocking send with ray stack version, to see if that's an improvement over the last looping version (which may or may not have used IBsend during the test - if this isn't an improvement, probably have to check again).

At some point I would like to make a graph of the dependencies and see if there's a best way to prioritize the rays.

HPC summer school

Statements are due midnight tonight, AoE.

Exoplanets III

Talk abstract for MHD submitted. Deadline is this Friday, I assume we'll hear back Feb-Mar. Working on flight and hotel.

Misc

I've updated the development branch to use the MPI_F08 module, which allows type and bounds checking, makes iErr optional, and in general is both better and the only way to use MPI that's completely compatible with the Fortran standard. Sent a merge request on Gitlab.

I now have my own professional domain: alexcomputes.com (was tired of using my much less professional free domain on my CV).

Convection project update (EOS)

Physics of MESA EOS

pdf: https://www.pas.rochester.edu/~yishengtu/research_files/CEE_gp_meeting_01272020/EOS.pdf
pptx: https://www.pas.rochester.edu/~yishengtu/research_files/CEE_gp_meeting_01272020/EOS.pptx

Integrating EOS into AstroBEAR

  • Finished getting EOS from MESA; Inverted MESA EOS tables for astrobear
  • Starting to add modules into EOS.f90

Next Step

  • continue putting EOS into Astrobear & run test modules
  • calculate whether new EOS will naturally bring convection

Issue With Conductive Shell Around Asteroid

The problem with enveloping the asteroid with a low resisitivity shell to simulate a mock atmosphere persists.

Things we have tried so far:

  1. Ramp up Solar Wind; ramp down asteroid density and temperature (for absorbing boundary condition); ramp down resistivity.
  2. Introduce is a stronger Diffusion factor using fluxes which should have smoothed things out (it didn't).
  3. Make a function for applying boundary conditions instead of copy-pasting it 12 places. Solved an asymmetry error in the code (still unsure how), but doesn't help otherwise.
  4. Trying different resisitivites. All with conductive shells have reached Min Timestep (I was too ambitious with the values perhaps).

Notes:

  • Simulations stop either by density protections at the night side or zones of high B in the day side.
  • The Oran wind mode (700 km/s, 35 particles/cc) seems to run a little longer than our wind (500 km/s, 300 particles/cc).

Planet-Civ Evolution Project Update 1/27/20

Planet-Civ Evolution Project Update 1/19