Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

create an NSE table with finer spacing in rho and Ye #1427

Merged
merged 3 commits into from
Dec 29, 2023
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
224,106 changes: 177,905 additions & 46,201 deletions networks/aprox19/nse_aprox19.tbl

Large diffs are not rendered by default.

8 changes: 4 additions & 4 deletions networks/aprox19/nse_table_size.H
Original file line number Diff line number Diff line change
Expand Up @@ -12,20 +12,20 @@ namespace nse_table_size {
const std::string table_name{"nse_aprox19.tbl"};

constexpr int ntemp = 101;
constexpr int nden = 31;
constexpr int nye = 15;
constexpr int nden = 61;
constexpr int nye = 29;

constexpr Real logT_min = 9.4;
constexpr Real logT_max = 10.4;
constexpr Real dlogT = 0.01;

constexpr Real logrho_min = 7.0;
constexpr Real logrho_max = 10.0;
constexpr Real dlogrho = 0.1;
constexpr Real dlogrho = 0.05;

constexpr Real ye_min = 0.43;
constexpr Real ye_max = 0.5;
constexpr Real dye = 0.005;
constexpr Real dye = 0.0025;

}
#endif
4 changes: 2 additions & 2 deletions nse_tabular/make_nse_table.py
Original file line number Diff line number Diff line change
Expand Up @@ -128,8 +128,8 @@ def generate_table():
nse_net = make_nse_network()

Ts = np.logspace(9.4, 10.4, 101)
rhos = np.logspace(7, 10, 31)
yes = np.linspace(0.43, 0.5, 15)
rhos = np.logspace(7, 10, 61)
yes = np.linspace(0.43, 0.5, 29)

#Ts = np.logspace(9.4, 10.4, 26)
#rhos = np.logspace(7, 10, 16)
Expand Down
86 changes: 43 additions & 43 deletions unit_test/burn_cell/ci-benchmarks/aprox19_nse_unit_test.out
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
Initializing AMReX (23.11-5-gd36463103dae)...
AMReX (23.11-5-gd36463103dae) initialized
Initializing AMReX (23.12-21-gef38229189e3)...
AMReX (23.12-21-gef38229189e3) initialized
starting the single zone burn...
reading the NSE table (C++) ...
Maximum Time (s): 0.01
Expand Down Expand Up @@ -48,53 +48,53 @@ burn entered NSE during integration (after 280 steps), zone = (0, 0, 0)
recovering burn failure in NSE, zone = (0, 0, 0)
------------------------------------
successful? 1
- Hnuc = 2.576395456e+21
- added e = 2.576395456e+19
- final T = 1.063285386e+10
- Hnuc = 2.578034812e+21
- added e = 2.578034812e+19
- final T = 1.063446163e+10
------------------------------------
e initial = 1.516256085e+18
e final = 2.728021064e+19
e final = 2.72966042e+19
------------------------------------
new mass fractions:
H1 3.698300896e-05
He3 7.832720283e-11
He4 0.06033012727
C12 9.999991743e-31
N14 9.875775959e-14
O16 2.423296446e-15
Ne20 9.999991743e-31
Mg24 3.915447813e-20
Si28 2.173376783e-22
S32 9.999991743e-31
Ar36 9.999991743e-31
Ca40 9.999991743e-31
Ti44 4.768494533e-30
Cr48 9.999991743e-31
Fe52 2.410423559e-30
Fe54 9.999991743e-31
Ni56 9.999991743e-31
n 0.4681726624
p 0.4714602272
H1 5.812961393e-07
He3 9.821888481e-12
He4 0.06008117992
C12 9.999992174e-31
N14 1.534220194e-15
O16 2.310011368e-15
Ne20 1.352206382e-19
Mg24 5.363834622e-22
Si28 4.61640554e-23
S32 9.999992174e-31
Ar36 9.999992174e-31
Ca40 4.564517527e-29
Ti44 9.999992174e-31
Cr48 5.163875678e-30
Fe52 9.999992174e-31
Fe54 3.381529446e-28
Ni56 9.999992174e-31
n 0.4683168262
p 0.4716014125
------------------------------------
species creation rates:
omegadot(H1): -9.996301699
omegadot(He3): -2.499999992
omegadot(He4): -73.96698727
omegadot(H1): -9.99994187
omegadot(He3): -2.499999999
omegadot(He4): -73.99188201
omegadot(C12): -2.5
omegadot(N14): -2.5
omegadot(O16): -2.5
omegadot(Ne20): -8.256539967e-35
omegadot(Mg24): 3.915447813e-18
omegadot(Si28): 2.173376773e-20
omegadot(S32): -8.256539967e-35
omegadot(Ar36): -8.256539967e-35
omegadot(Ca40): -8.256539967e-35
omegadot(Ti44): 3.768494533e-28
omegadot(Cr48): -8.256539967e-35
omegadot(Fe52): 1.410423559e-28
omegadot(Fe54): -8.256539967e-35
omegadot(Ni56): -8.256539967e-35
omegadot(n): 46.81726624
omegadot(p): 47.14602272
number of steps taken: 38629
AMReX (23.11-5-gd36463103dae) finalized
omegadot(Ne20): 1.352206382e-17
omegadot(Mg24): 5.363834612e-20
omegadot(Si28): 4.61640544e-21
omegadot(S32): -7.826200767e-35
omegadot(Ar36): -7.826200767e-35
omegadot(Ca40): 4.464517527e-27
omegadot(Ti44): -7.826200767e-35
omegadot(Cr48): 4.163875678e-28
omegadot(Fe52): -7.826200767e-35
omegadot(Fe54): 3.371529446e-26
omegadot(Ni56): -7.826200767e-35
omegadot(n): 46.83168262
omegadot(p): 47.16014125
number of steps taken: 27367
AMReX (23.12-21-gef38229189e3) finalized
7 changes: 4 additions & 3 deletions unit_test/test_nse/ci-benchmarks/aprox21_ci.out
Original file line number Diff line number Diff line change
@@ -1,6 +1,7 @@
AMReX (22.10-25-g735c3513153f) initialized
Initializing AMReX (23.12-21-gef38229189e3)...
AMReX (23.12-21-gef38229189e3) initialized
starting the single zone burn...
0x7ffd73ad6d20
0x7ffcb5676e40
State Density (g/cm^3): 277355338.4
State Temperature (K): 5197769252
Mass Fraction (H1): 0.4
Expand Down Expand Up @@ -50,4 +51,4 @@ Fe56 : 1.686482626e-05
Ni56 : 0.8982137413
n : 9.442515318e-10
p : 0.001666851761
AMReX (22.10-25-g735c3513153f) finalized
AMReX (23.12-21-gef38229189e3) finalized
112 changes: 56 additions & 56 deletions unit_test/test_nse_interp/ci-benchmarks/aprox19.out
Original file line number Diff line number Diff line change
@@ -1,78 +1,78 @@
Initializing AMReX (23.11-5-gd36463103dae)...
AMReX (23.11-5-gd36463103dae) initialized
Initializing AMReX (23.12-21-gef38229189e3)...
AMReX (23.12-21-gef38229189e3) initialized
starting the single zone burn...
reading the NSE table (C++) ...
rho, T, Ye = 1230000000 5180000000 0.472
density value brackets: 9 < 9.089905111 < 9.1
density value brackets: 9.05 < 9.089905111 < 9.1
temperature value brackets: 9.71 < 9.71432976 < 9.72
Ye value brackets: 0.475 < 0.472 < 0.47
Ye value brackets: 0.4725 < 0.472 < 0.47

4 rho values:
8.9 55.16599379
9 55.28915823
9.1 55.39796653
9.2 55.49591413
9 55.5223665
9.05 55.57080116
9.1 55.61707675
9.15 55.66151541

cubic interpolated value: 55.38753351
cubic interpolated value: 55.6078922

4 T values:
9.7 55.16599379
9.71 54.97411678
9.72 54.71072288
9.73 54.35084601
9.7 55.5223665
9.71 55.39489494
9.72 55.22387568
9.73 54.99343581

cubic interpolated value: 54.87031639
cubic interpolated value: 55.32712406

4 Ye values:
0.48 55.16599379
0.475 55.41756427
0.47 55.68822983
0.465 55.93040523
0.475 55.5223665
0.4725 55.65227675
0.47 55.78142906
0.4675 55.90244371

cubic interpolated value: 55.58071765
cubic interpolated value: 55.678404

tricubic interpolated values:
abar = 55.60621897
bea = 8.755258739
dyedt = -0.004462626048
dabardt = 2.472895091e-16
dbeadt = 0.0001077443794
e_nu = 1.081686275e+16
abar = 55.60652462
bea = 8.755119523
dyedt = -0.004527880871
dabardt = 1.310223024e-16
dbeadt = 0.000111945508
e_nu = 1.101641642e+16
X(H1) = 0
X(He3) = 2.646187332e-14
X(He4) = 0.0003164651992
X(C12) = 7.949931897e-10
X(N14) = 3.087926554e-14
X(O16) = 1.98556354e-09
X(Ne20) = 2.420691526e-11
X(Mg24) = 1.01690623e-08
X(Si28) = 7.957002233e-06
X(S32) = 1.148035178e-05
X(Ar36) = 1.188979307e-05
X(Ca40) = 1.11725697e-05
X(Ti44) = 0.0001435646298
X(Cr48) = 0.01141230313
X(Fe52) = 0.07731922235
X(Fe54) = 0.9104458693
X(Ni56) = 0.0003104223608
X(n) = 2.004738472e-08
X(p) = 9.620335872e-06
X(He3) = 2.723457525e-14
X(He4) = 0.0003164137294
X(C12) = 7.975878903e-10
X(N14) = 3.102372259e-14
X(O16) = 1.998161464e-09
X(Ne20) = 2.444604777e-11
X(Mg24) = 1.029036116e-08
X(Si28) = 8.173012779e-06
X(S32) = 1.194076593e-05
X(Ar36) = 1.243414613e-05
X(Ca40) = 1.259531492e-05
X(Ti44) = 0.0001438001413
X(Cr48) = 0.01146726617
X(Fe52) = 0.07746625779
X(Fe54) = 0.9098424639
X(Ni56) = 0.0007087774945
X(n) = 2.078852742e-08
X(p) = 9.843630258e-06

testing temperature derivatives of cubic
first finite-difference derivatives
dAbar/dT = -1.070778042e-09
dbea/dT = -6.886297725e-12
dAbar/dT = -1.072562355e-09
dbea/dT = -6.867059575e-12
now using derivative of the interpolant
dAbar/dT = -1.070778122e-09
dbea/dT = -6.886272826e-12
dAbar/dT = -1.072562604e-09
dbea/dT = -6.867048589e-12

EOS e consistency check (old method): 1.395273834e+18 1.388449367e+18
updated T: 6394612134
change in abar: 55.60621897 50.27196857
EOS e consistency check (new method): 1.388449367e+18 1.388449367e+18
EOS e consistency check (old method): 1.395278886e+18 1.38844906e+18
updated T: 6394534499
change in abar: 55.60652462 50.26831386
EOS e consistency check (new method): 1.38844906e+18 1.38844906e+18

EOS p consistency check (old method): 6.622132093e+26 6.57785215e+26
updated T: 6466848772
change in abar: 55.60621897 49.50670173
EOS p consistency check (new method): 6.57785215e+26 6.57785215e+26
AMReX (23.11-5-gd36463103dae) finalized
EOS p consistency check (old method): 6.622159603e+26 6.577850616e+26
updated T: 6466757500
change in abar: 55.60652462 49.50320619
EOS p consistency check (new method): 6.577850616e+26 6.577850616e+26
AMReX (23.12-21-gef38229189e3) finalized
Loading