Skip to content

Commit a36a0d1

Browse files
authored
add NSE network compatibility docs (#1852)
1 parent 4db027a commit a36a0d1

File tree

1 file changed

+25
-0
lines changed

1 file changed

+25
-0
lines changed

Docs/source/nse_net.rst

Lines changed: 25 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -14,6 +14,31 @@ under the assumption of NSE, following the procedure outlined in :cite:`Calder_2
1414
Self-consistent NSE does not support the templated C++ networks
1515
(like ``aprox13``). You should use a pynucastro-generated network.
1616

17+
Additionally, the pynucastro network must meet the following requirements:
18+
19+
* Every rate should have an inverse rate, recomputed with detailed balancing
20+
21+
* The "Q" value should be recomputed using the masses when doing detailed
22+
balance.
23+
24+
* Approximate rates like $(\alpha, p)(p, \gamma)$ are okay, but modified
25+
rates that change the stoichiometry are not supported.
26+
27+
.. tip::
28+
29+
The test in ``Microphysics/nse_solver/nse_compatibility/`` can be used to
30+
check if a network is supported. It will integrate at fixed $(rho, T)$
31+
and compare the resulting mass fractions to the NSE solution.
32+
33+
If a network has weak rates, then there will be a $Y_e$ evolution which
34+
will prevent this test from giving perfect agreement, but if the $Y_e$
35+
evolution is slow, this test can still give useful results.
36+
37+
.. tip::
38+
39+
Presently it is recommended to use the ``ase`` network.
40+
41+
1742
There are two methods for solving the chemical potentials:
1843

1944
1. Hybrid Powell's method. This is the default method.

0 commit comments

Comments
 (0)