Difference between revisions of "FINIX PB2 example input"

From Serpent Wiki
Jump to: navigation, search
(Main Serpent input)
Line 80: Line 80:
  
 
mesh 10 2 1000 3000 0 -7.62 7.62 -7.62 7.62 -250.0 250.0</nowiki>
 
mesh 10 2 1000 3000 0 -7.62 7.62 -7.62 7.62 -250.0 250.0</nowiki>
 +
 +
== ''materials'' file  ==
  
 
== References ==
 
== References ==
 +
 
<references/>
 
<references/>

Revision as of 14:18, 1 December 2015

The input model is based on the Peach Bottom 2 BWR assembly from the UAM benchmark [1] with the segmented absorber rod replaced by a full length one.

Main Serpent input

% --- Asymmetric BWR assembly with Gd-pins -------------------

set title "Peach Bottom 2 from UAM, FINIX calculation"

% --- Remember to set TMS limits for materials

include materials

include pins

include geometry

% --- Reflective boundary condition in XY, black in Z:

set bc 2 2 1

% --- Link FINIX input files

set finrodfile      "./finixfiles/rod_library.txt"
set finoptionsfile  "./finixfiles/options_library.txt"
set finscenariofile "./finixfiles/scenario_library.txt"

% --- Define sub-pin-level power binning
% Name bin1, 6 radial zones in pellet, 24 axial zones

finbin bin1 0.0 0.60579 6 -182.88 182.88 24

% --- Link finix solvers for each pin type

finrod 1 rod_pb2_nogd  options_pb2_all scenario_pb2_all bin1
finrod 2 rod_pb2_nogd  options_pb2_all scenario_pb2_all bin1
finrod 3 rod_pb2_nogd  options_pb2_all scenario_pb2_all bin1
finrod 4 rod_pb2_nogd  options_pb2_all scenario_pb2_all bin1
finrod 5 rod_pb2_yesgd options_pb2_all scenario_pb2_all bin1

% --- Neutron population and criticality cycles:
% --- This is simulated on each iteration

set pop 5000 10 100

% --- This is the number of iterations

set ccmaxiter 10

% --- Use a smaller number of inactive cycles after each FINIX solution
% --- (Fission source passing on, number of alternative cycles is 10)

set fsp 1 10

% --- Total power

set power 3.9429e+06

% --- Geometry and mesh plots:

plot 3 1000 1000 0.0 
plot 3 1000 1000 0.0 -22.86 7.62 -22.86 7.62
plot 2 500 1500
plot 1 500 1500

% --- xy-plot of fission power / thermal flux

mesh 3 1000 1000 0 -7.62 7.62 -7.62 7.62 -250.0 250.0

% --- xz-plot of fission power / thermal flux

mesh 2 1000 3000 0 -7.62 7.62 -7.62 7.62 -250.0 250.0

% --- xy-plot of temperature distribution / thermal flux

mesh 10 3 1000 1000 0 -7.62 7.62 -7.62 7.62 -250.0 250.0

% --- xz-plot of temperature distribution / thermal flux

mesh 10 2 1000 3000 0 -7.62 7.62 -7.62 7.62 -250.0 250.0

materials file

References

  1. ^ Ivanov et al. "BENCHMARKS FOR UNCERTAINTY ANALYSIS IN MODELLING (UAM) FOR THE DESIGN, OPERATION AND SAFETY ANALYSIS OF LWRs, Volume I: Specification and Support Data for Neutronics Cases (Phase I), Version 2.1", NEA/NSC/DOC(2013)7