|
Cococubed.com
|
| Pressure from the LS and an NSE EOS |
Home
Commercial:
Artwork
Software
Teaching materials
Bicycle sag support
Free:
Family Album
Pretty astronomy pictures
Some astronomy codes
... Stellar equation of states
... EOS with ionization
... EOS for supernovae
... Chemical potentials
... Stellar atmospheres
... Voigt Function
... Polytropic stars
... Cold white dwarfs
... Hotter white dwarfs
... Cold neutron stars
... Stellar opacities
... Neutrino energy loss rates
... Ephemeris routines
... Fermi-Dirac functions
... Galactic chemical evolution
... Nuclear reaction networks
... Nuclear statistical equilibrium
... Laminar deflagrations
... CJ detonations
... ZND detonations
... Fitting to conic sections
... Unusual linear algebra
... Derivatives on uneven grids
... Pentadiagonal solver
... Quadratics, Cubics, Quartics
... Supernova light curves
... Exact Riemann solutions
... 1D PPM Hydrodynamics
... Verification problems
... EZ stellar evolution
... FLASH code
... Mesa code
Some astronomy talks
Some research
Bicycle adventures
Contact us:
J.D. Maldonado
F.X.Timmes, my vitae
The electron-positron portion of the two equations of state are identical. Differences between the two EOS routines originate from the model used for nucleons (interacting nucleons in a liquid dropish model versus non-interacting Boltzmann nucleons), and the model used for the composition.
LS Baryon Pressure:
Pressure from baryons Ye=0.5 |
Pressure from baryons Ye=0.4 |
Pressure from baryons Ye=0.3 |
Pressure from baryons Ye=0.2 |
In certain regions of the rho-T plane, contributions from the coulomb lattice terms cause the baryon pressure to become negative.
LS Total Pressure:
Total pressure Ye=0.5 |
Total pressure Ye=0.4 |
Total pressure Ye=0.3 |
Total pressure Ye=0.2 |
NSE Baryon Pressure:
Pressure from baryons Ye=0.5 |
Pressure from baryons Ye=0.4 |
Pressure from baryons Ye=0.3 |
Pressure from baryons Ye=0.2 |
Even though the nuclei in the NSE-based model are a perfect gas, the surface isn't planar because of the changing composition. This causes the ripples in the surface.
NSE Total Pressure:
Total pressure Ye=0.5 |
Total pressure Ye=0.4 |
Total pressure Ye=0.3 |
Total pressure Ye=0.2 |
LS and NSE Pressures Compared:
Pressure Ye=0.5 |
Pressure Ye=0.4 |
Pressure Ye=0.3 |
Pressure Ye=0.2 |
Overall, i think the agreement between the two is acceptable given the behavior of the coulomb lattice terms discussed above.
Pressure differences Ye=0.5 |
Pressure differences Ye=0.4 |
Pressure differences Ye=0.3 |
Pressure differences Ye=0.2 |