[Mesa-users] Nuclear binding energy released

[Bill Paxton]

for energy from burning with nuclear reaction neutrinos subtracted, you could use cumulative_nuclear_heating

where
           s% total_nuclear_heating = dt*dot_product(s% dm(1:nz), s% eps_nuc(1:nz))
and
           s% cumulative_nuclear_heating = &
               s% cumulative_nuclear_heating_old + s% total_nuclear_heating

you can add cumulative_nuclear_heating to your history columns list.

 -b


> On Jan 24, 2019, at 3:44 PM, Josiah Schwab via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Hi Shiv,
> 
> I want to compute the nuclear binding energy released by burning over the whole simulation.
> 
> In that case, you also might consider computing the binding energy at the first step, storing that value, calculating it at subsequent steps, and then taking the difference.
> 
> Note that the mass fractions are in s% xa(i,k) and there is also a binding energy value on chem_isos, i.e. chem_isos% binding_energy(ci).
> 
> The bookkeeping might be complicated if your models have mass loss though.
> 
> Hope that helps,
> Josiah 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190124/dbd947e9/attachment.html>