Burn_nseAtDens.F90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	Burn_nseAtDens (qbar_nse, sumyi_nse, approxtemp, edot, Yedot, Ye, dens, emq)

Function/Subroutine Documentation

Burn_nseAtDens()

subroutine Burn_nseAtDens	(	real	qbar_nse,
		real	sumyi_nse,
		real	approxtemp,
		real	edot,
		real	Yedot,
		real	Ye,
		real	dens,
		real	emq
	)

This function gives the NSE final state for a constant density burn via table lookup. qbar is returned in MeV/nucleon. Only the compositional information (qbar, 1/Abar, edot, Yedot) should be used for hydrodynamics purposes. Accurate thermodynamic prperties of the NSE final state (rho, T) should be obtained by solving eint - qbar_nse = emq . This is important due to the limited accuracy of the table interpolation (i.e., the interpolated values cannot satisfy any constraints on their own.)

qbar_nse – Average binding energy per nucleon (averaged over nse) sumyi_nse – Sum over the abundances. (Y = X/A = X/(#of nucleons)) approxtemp – ?? Ask Dean edot – Neutrino and anti-neutrino loss rates [ergs/g] Yedot – Change in electron fraction wrt time (neutronization rate) Ye – Electron fraction dens – Density in/out emq – E-Q = Internal energy minus binding energy (per volume?)

NOTES Takes place after burning is complete. NSE expands after burning and the equilibrium state is a function of the density and the temperature. Been burned, flame is done, but then the hot stuff expands. This changes the composition and releases energy.

Definition at line 46 of file Burn_nseAtDens.F90.