Equation of State reference
I download all these papers & resources and put them on PAS server.
EOS papers
MESA EOS
- https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/MESA%201.pdf
- https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/MESA%205.pdf
OPAL:
- https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/Rogers+96.pdf
- https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/Rogers+Nayfonov02.pdf
SCVH:
PTEH:
HELM:
Other useful resources
- EOS Presentation (very useful!) \\[https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/OPAL%20PPT.pdf]
- Statistical Mechanics testbook \\[https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/Callen,%20Herbert%20B%20-%20Thermodynamics%20and%20an%20Introduction%20to%20Thermostatistics%202nd%20Edition.pdf]
- Another statistical textbook (referred in EOS papers and used as undergrad textbook) \\[https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/Fundamentals-of-Statistical-And-Thermal-Physics.pdf]
- More OPAL EOS papers \\[https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/Brian%201995.pdf]
- More OPAL EOS papers \\[https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/Lin-D%c3%a4ppen2010_Article_EmulatingTheOPALEquationOfStat.pdf]
- Other MESA papers (they don't have information about EOS)
https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/MESA%202.pdf \\[https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/MESA%203.pdf] \\[https://www.pas.rochester.edu/~yishengtu/research_files/EOS_reference/MESA%204.pdf]
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
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
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:
- Ramp up Solar Wind; ramp down asteroid density and temperature (for absorbing boundary condition); ramp down resistivity.
- Introduce is a stronger Diffusion factor using fluxes which should have smoothed things out (it didn't).
- 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.
- 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
Current state of the Equations and example plots using Earth-like inputs
4x4 Grid of Temperature/Population vs time: Horizontal Plots vary by Orbital Distance, Vertical Plots vary by the Maximum Allowed Population (Carrying Capacity)
4x4 Grid of Population vs Temperature: Horizontal Plots vary by Orbital Distance, Vertical Plots vary by the Maximum Allowed Population (Carrying Capacity)