Skip to content

Commit

Permalink
Merge pull request #163 from smpark7/inscatter-160
Browse files Browse the repository at this point in the history 
Fix InScatter Jacobian term for Serpent cross sections
  • Loading branch information
@gwenchee
gwenchee authored Sep 16, 2021
2 parents 58332b5 + 1cfc2d0 commit 9dc170e
Show file tree
Hide file tree
Showing 30 changed files with 428 additions and 6 deletions.
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_BETA_EFF.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 2.37389E-04 1.54596E-03 1.38220E-03 2.80828E-03 8.26765E-04 2.85805E-04
972 2.37210E-04 1.54495E-03 1.38131E-03 2.80634E-03 8.26142E-04 2.85674E-04
1022 2.37036E-04 1.54397E-03 1.38044E-03 2.80445E-03 8.25534E-04 2.85552E-04
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_CHIP.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 1.00000E+00 0.00000E+00
972 1.00000E+00 0.00000E+00
1022 1.00000E+00 0.00000E+00
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_CHIT.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 1.00000E+00 0.00000E+00
972 1.00000E+00 0.00000E+00
1022 1.00000E+00 0.00000E+00
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_DIFFCOEF.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 1.42796148502 1.26251148885
972 1.44313263688 1.2777364643
1022 1.45794060058 1.29262094408
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_FISS.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.00130417 0.0217635
972 0.00129635 0.0211145
1022 0.00128962 0.0205105
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_INVV.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 9.68949E-08 2.05657E-06
972 9.79604E-08 2.01960E-06
1022 9.91412E-08 1.98497E-06
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_KAPPA.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 194.31036598 194.286017874
972 194.309021483 194.28532525
1022 194.309660985 194.284939421
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_LAMBDA.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 1.24378E-02 3.06300E-02 1.11474E-01 3.02248E-01 1.17229E+00 3.07474E+00
972 1.24376E-02 3.06293E-02 1.11471E-01 3.02241E-01 1.17212E+00 3.07438E+00
1022 1.24374E-02 3.06286E-02 1.11468E-01 3.02234E-01 1.17194E+00 3.07400E+00
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_NSF.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.00318268 0.0530312
972 0.00316355 0.0514497
1022 0.0031471 0.0499778
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_NUBAR.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 2.44038737281 2.4367036552
972 2.44035175686 2.43669989817
1022 2.44033126037 2.43669340094
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_REMXS.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.00572091 0.0286412
972 0.00577444 0.02804167
1022 0.0058392 0.02751448
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_fuel_SP0.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.266445 0.00202371 0.0015943 0.249011
972 0.26352 0.00209278 0.00178607 0.246277
1022 0.260705 0.00217227 0.00199628 0.243624
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_BETA_EFF.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
972 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
1022 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_CHIP.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.00000E+00 0.00000E+00
972 0.00000E+00 0.00000E+00
1022 0.00000E+00 0.00000E+00
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_CHIT.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.00000E+00 0.00000E+00
972 0.00000E+00 0.00000E+00
1022 0.00000E+00 0.00000E+00
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_DIFFCOEF.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.985027580772 0.769666518891
972 0.985732507684 0.770022115035
1022 0.985875363665 0.770014999892
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_FISS.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.0 0.0
972 0.0 0.0
1022 0.0 0.0
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_INVV.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 9.98046E-08 2.07824E-06
972 1.00883E-07 2.03998E-06
1022 1.02080E-07 2.00420E-06
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_KAPPA.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.0 0.0
972 0.0 0.0
1022 0.0 0.0
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_LAMBDA.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
972 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
1022 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_NSF.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.0 0.0
972 0.0 0.0
1022 0.0 0.0
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_NUBAR.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.0 0.0
972 0.0 0.0
1022 0.0 0.0
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_REMXS.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.0042598939 0.003109062
972 0.0044414982 0.003517643
1022 0.004648825 0.003971745
3 changes: 3 additions & 0 deletions property_file_dir/newt_msre_converted_to_serpent/serpent_msre_mod_SP0.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
922 0.393826 0.00424578 0.00296472 0.45548
972 0.393324 0.00442734 0.0033761 0.454857
1022 0.393018 0.0046346 0.00383269 0.454403
8 changes: 6 additions & 2 deletions src/kernels/InScatter.C
View file
Original file line number Diff line number Diff line change
Expand Up @@ -91,8 +91,12 @@ InScatter::computeQpOffDiagJacobian(unsigned int jvar)
{
if (i == _group)
continue;
jac += -_test[_i][_qp] * _d_gtransfxs_d_temp[_qp][i + _group * _num_groups] * _phi[_j][_qp] *
computeConcentration((*_group_fluxes[i]), _qp);
if (_sss2_input)
jac += -_test[_i][_qp] * _d_gtransfxs_d_temp[_qp][i * _num_groups + _group] * _phi[_j][_qp] *
computeConcentration((*_group_fluxes[i]), _qp);
else
jac += -_test[_i][_qp] * _d_gtransfxs_d_temp[_qp][i + _group * _num_groups] * _phi[_j][_qp] *
computeConcentration((*_group_fluxes[i]), _qp);
}
}

Expand Down
6 changes: 4 additions & 2 deletions tests/twod_axi_coupled/auto_diff_rho.i
View file
Original file line number Diff line number Diff line change
Expand Up @@ -226,7 +226,7 @@ nt_scale=1e13

[Executioner]
type = Transient
end_time = 10000
num_steps = 80

nl_rel_tol = 1e-6
nl_abs_tol = 1e-6
Expand All @@ -241,6 +241,8 @@ nt_scale=1e13

nl_max_its = 30
l_max_its = 100
steady_state_detection = true
steady_state_tolerance = 1e-10

dtmin = 1e-5
# dtmax = 1
Expand All @@ -250,7 +252,7 @@ nt_scale=1e13
dt = 1e-3
cutback_factor = 0.4
growth_factor = 1.2
optimal_iterations = 20
optimal_iterations = 5
[../]
[]

Expand Down
Loading

0 comments on commit 9dc170e

Please sign in to comment.