From 76aab72a062c8a918102412f4ebdc5bc616973e4 Mon Sep 17 00:00:00 2001 From: Austin Harris Date: Tue, 2 Jul 2024 09:41:26 -0400 Subject: [PATCH] Add test for log(ft) to test driver --- test/Data_SN231_logft/ab_co | 59 ++++++++++++++++++++++++++++ test/Test_Problems/setup_logft_sn231 | 20 ++++++++++ test/test_settings_logft | 32 +++++++++++++++ test/test_xnet.sh | 15 +++++++ 4 files changed, 126 insertions(+) create mode 100644 test/Data_SN231_logft/ab_co create mode 100644 test/Test_Problems/setup_logft_sn231 create mode 100644 test/test_settings_logft diff --git a/test/Data_SN231_logft/ab_co b/test/Data_SN231_logft/ab_co new file mode 100644 index 00000000..324f8189 --- /dev/null +++ b/test/Data_SN231_logft/ab_co @@ -0,0 +1,59 @@ +50% C12, 49% O16, 1% Ne22 (Ye=0.4995) + n 0.0000000E+00 p 0.0000000E+00 d 0.0000000E+00 t 0.0000000E+00 + he3 0.0000000E+00 he4 0.0000000E+00 li6 0.0000000E+00 li7 0.0000000E+00 + be7 0.0000000E+00 be9 0.0000000E+00 b8 0.0000000E+00 b10 0.0000000E+00 + b11 0.0000000E+00 c11 0.0000000E+00 c12 4.1250000E-02 c13 0.0000000E+00 + c14 0.0000000E+00 n12 0.0000000E+00 n13 0.0000000E+00 n14 0.0000000E+00 + n15 0.0000000E+00 o14 0.0000000E+00 o15 0.0000000E+00 o16 3.0937500E-02 + o17 0.0000000E+00 o18 0.0000000E+00 o19 0.0000000E+00 f17 0.0000000E+00 + f18 0.0000000E+00 f19 0.0000000E+00 f20 0.0000000E+00 f21 0.0000000E+00 + ne18 0.0000000E+00 ne19 0.0000000E+00 ne20 0.0000000E+00 ne21 0.0000000E+00 + ne22 4.5454545E-04 ne23 0.0000000E+00 ne24 0.0000000E+00 na20 0.0000000E+00 + na21 0.0000000E+00 na22 0.0000000E+00 na23 0.0000000E+00 na24 0.0000000E+00 + na25 0.0000000E+00 na26 0.0000000E+00 mg22 0.0000000E+00 mg23 0.0000000E+00 + mg24 0.0000000E+00 mg25 0.0000000E+00 mg26 0.0000000E+00 mg27 0.0000000E+00 + mg28 0.0000000E+00 mg29 0.0000000E+00 al24 0.0000000E+00 al25 0.0000000E+00 + al26 0.0000000E+00 al27 0.0000000E+00 al28 0.0000000E+00 al29 0.0000000E+00 + al30 0.0000000E+00 si26 0.0000000E+00 si27 0.0000000E+00 si28 0.0000000E+00 + si29 0.0000000E+00 si30 0.0000000E+00 si31 0.0000000E+00 si32 0.0000000E+00 + si33 0.0000000E+00 p28 0.0000000E+00 p29 0.0000000E+00 p30 0.0000000E+00 + p31 0.0000000E+00 p32 0.0000000E+00 p33 0.0000000E+00 p34 0.0000000E+00 + p35 0.0000000E+00 s30 0.0000000E+00 s31 0.0000000E+00 s32 0.0000000E+00 + s33 0.0000000E+00 s34 0.0000000E+00 s35 0.0000000E+00 s36 0.0000000E+00 + s37 0.0000000E+00 cl32 0.0000000E+00 cl33 0.0000000E+00 cl34 0.0000000E+00 + cl35 0.0000000E+00 cl36 0.0000000E+00 cl37 0.0000000E+00 cl38 0.0000000E+00 + cl39 0.0000000E+00 ar34 0.0000000E+00 ar35 0.0000000E+00 ar36 0.0000000E+00 + ar37 0.0000000E+00 ar38 0.0000000E+00 ar39 0.0000000E+00 ar40 0.0000000E+00 + ar41 0.0000000E+00 ar42 0.0000000E+00 k36 0.0000000E+00 k37 0.0000000E+00 + k38 0.0000000E+00 k39 0.0000000E+00 k40 0.0000000E+00 k41 0.0000000E+00 + k42 0.0000000E+00 k43 0.0000000E+00 k44 0.0000000E+00 ca38 0.0000000E+00 + ca39 0.0000000E+00 ca40 0.0000000E+00 ca41 0.0000000E+00 ca42 0.0000000E+00 + ca43 0.0000000E+00 ca44 0.0000000E+00 ca45 0.0000000E+00 ca46 0.0000000E+00 + ca47 0.0000000E+00 ca48 0.0000000E+00 sc41 0.0000000E+00 sc42 0.0000000E+00 + sc43 0.0000000E+00 sc44 0.0000000E+00 sc45 0.0000000E+00 sc46 0.0000000E+00 + sc47 0.0000000E+00 sc48 0.0000000E+00 sc49 0.0000000E+00 ti42 0.0000000E+00 + ti43 0.0000000E+00 ti44 0.0000000E+00 ti45 0.0000000E+00 ti46 0.0000000E+00 + ti47 0.0000000E+00 ti48 0.0000000E+00 ti49 0.0000000E+00 ti50 0.0000000E+00 + ti51 0.0000000E+00 v45 0.0000000E+00 v46 0.0000000E+00 v47 0.0000000E+00 + v48 0.0000000E+00 v49 0.0000000E+00 v50 0.0000000E+00 v51 0.0000000E+00 + v52 0.0000000E+00 v53 0.0000000E+00 cr47 0.0000000E+00 cr48 0.0000000E+00 + cr49 0.0000000E+00 cr50 0.0000000E+00 cr51 0.0000000E+00 cr52 0.0000000E+00 + cr53 0.0000000E+00 cr54 0.0000000E+00 cr55 0.0000000E+00 mn49 0.0000000E+00 + mn50 0.0000000E+00 mn51 0.0000000E+00 mn52 0.0000000E+00 mn53 0.0000000E+00 + mn54 0.0000000E+00 mn55 0.0000000E+00 mn56 0.0000000E+00 mn57 0.0000000E+00 + fe51 0.0000000E+00 fe52 0.0000000E+00 fe53 0.0000000E+00 fe54 0.0000000E+00 + fe55 0.0000000E+00 fe56 0.0000000E+00 fe57 0.0000000E+00 fe58 0.0000000E+00 + fe59 0.0000000E+00 fe60 0.0000000E+00 co52 0.0000000E+00 co53 0.0000000E+00 + co54 0.0000000E+00 co55 0.0000000E+00 co56 0.0000000E+00 co57 0.0000000E+00 + co58 0.0000000E+00 co59 0.0000000E+00 co60 0.0000000E+00 co61 0.0000000E+00 + ni55 0.0000000E+00 ni56 0.0000000E+00 ni57 0.0000000E+00 ni58 0.0000000E+00 + ni59 0.0000000E+00 ni60 0.0000000E+00 ni61 0.0000000E+00 ni62 0.0000000E+00 + ni63 0.0000000E+00 ni64 0.0000000E+00 cu57 0.0000000E+00 cu58 0.0000000E+00 + cu59 0.0000000E+00 cu60 0.0000000E+00 cu61 0.0000000E+00 cu62 0.0000000E+00 + cu63 0.0000000E+00 cu64 0.0000000E+00 cu65 0.0000000E+00 zn59 0.0000000E+00 + zn60 0.0000000E+00 zn61 0.0000000E+00 zn62 0.0000000E+00 zn63 0.0000000E+00 + zn64 0.0000000E+00 zn65 0.0000000E+00 zn66 0.0000000E+00 zn67 0.0000000E+00 + zn68 0.0000000E+00 ga62 0.0000000E+00 ga63 0.0000000E+00 ga64 0.0000000E+00 + ga65 0.0000000E+00 ga66 0.0000000E+00 ga67 0.0000000E+00 ga68 0.0000000E+00 + ga69 0.0000000E+00 ge65 0.0000000E+00 ge66 0.0000000E+00 ge67 0.0000000E+00 + ge68 0.0000000E+00 ge69 0.0000000E+00 ge70 0.0000000E+00 diff --git a/test/Test_Problems/setup_logft_sn231 b/test/Test_Problems/setup_logft_sn231 new file mode 100644 index 00000000..018c858a --- /dev/null +++ b/test/Test_Problems/setup_logft_sn231 @@ -0,0 +1,20 @@ + +## Output Controls +3 Diagnostic Output Level +2 Per Timestep Output Level +# ASCII output filename root, network will append zone number +Test_Results/ev_logft_sn160_ +# Binary output filename root, network will append zone number +Test_Results/ts_logft_sn160_ +# Species to output in ASCII output (format 14a5): 14 + p he4 c12 o16 ne20 mg24 si28 ca40 ti44 cr48 fe52 fe54 ni56 zn60 +## Input Controls +# Nuclear Data Directory +Data_SN231_logft +# Initial Abundance and Thermodynamic Trajectory Files +Data_SN231_logft/ab_co +Test_Problems/th_co_burn_1 +Data_SN231_logft/ab_co +Test_Problems/th_co_burn_2 +Data_SN231_logft/ab_co +Test_Problems/th_co_burn_3 diff --git a/test/test_settings_logft b/test/test_settings_logft new file mode 100644 index 00000000..7e6c6ea2 --- /dev/null +++ b/test/test_settings_logft @@ -0,0 +1,32 @@ +## Problem Description +2000 timesteps +diagnostic level 1 +New parameters +## Job Controls +1 Initial Zone +3 # of Zones +2 Include Weak Reactions (yes=1,no=0,only=-1) +1 Include Screening (yes=1) +1 Process Nuclear Data at Run Time (yes=1,no=0) +## Neutrinos +0 Include Neutrino Reactions (yes=1, no=0) +## NSE Initial Conditions +11.0 Temperature in GK to use NSE initial conditions instead of file +## Integration Controls +1 Choice of integration Scheme (1=Backward Euler, 2= Bader-Deufelhard) +9999 Max. number of timesteps before quit +5 Max. iterations per step +1 Rebuild the jacobian every ijac iterations after the first +0 Convergence Condition (Mass Cons.=0, (dY/Y small)=1) +1.00E-01 Max. Abundance Change per timestep +1.00E-07 Smallest Abundance used in timestep calculation +1.00E-06 Mass Conservation Limit +1.00E-04 Convergence Criterion +1.00E-30 Lower Abundance limit, smaller abundances = 0 +2.00E+00 Max. Factor to change dt in a timestep +## Self-heating Controls +1 Include self-heating (yes=1,no=0) +1.00E-02 Max. Temperature Change per timestep +1.00E-04 Temperature Convergence Criterion +## Zone Batching Controls +2 Blocking size for zone loop diff --git a/test/test_xnet.sh b/test/test_xnet.sh index f33cc70f..e8ad307d 100755 --- a/test/test_xnet.sh +++ b/test/test_xnet.sh @@ -45,6 +45,7 @@ ## torch47 (47) 52 ## SN160 (160) 53 ## SN160 (160) (BDF integrator) 54 +## SN231 (231) (log-ft rates) 55 ## ===================================================================== ## All parallel test problems 30 ## 4 different SN alpha (16) 31 @@ -109,6 +110,13 @@ function do_test_bdf { test_diff Test_Results/net_diag_$2 Test_Problems/Results/net_diag_$2 $2 } +function do_test_logft { + cat test_settings_logft Test_Problems/setup_$2 >| control + $1 + mv -f net_diag01 Test_Results/net_diag_$2 + test_diff Test_Results/net_diag_$2 Test_Problems/Results/net_diag_$2 $2 +} + function do_test_batch { cat test_settings_batch Test_Problems/setup_$2 >| control $1 @@ -297,6 +305,13 @@ for xnet in ${xnet_list[@]}; do do_test_bdf $xnet $test_name fi + # Self-heating test using SN160 (explosive burning of degenerate C/O) (with BDF integrator) + if [ $itest -eq 50 -o $itest -eq 55 ]; then + echo "Test: Self-heating from explosive burning of degenerate C/O with SN231 network and log(ft) rates" + test_th="logft"; test_net="sn231"; test_name=${test_th}_${test_net} + do_test_logft $xnet $test_name + fi + # Zone batching test using alpha (explosive burning of degenerate C/O) if [ $itest -eq 60 -o $itest -eq 61 ]; then echo "Test: Zone batching w/ self-heating from explosive burning of degenerate C/O with alpha network"