SE2-ANL Subchannel Analysis with New Heating Calculation Based on VARIANT Transport Code

Journal Article · Fri Jul 01 04:00:00 EDT 2016 · Transactions of the American Nuclear Society

OSTI ID:23042838

Yang, G.; Yang, W. S.; Deng, P. ^[1]

---

1. Purdue University, School of Nuclear Engineering, West Lafayette, IN 47906-2017 (United States)

The multi-assembly, steady-state sub-channel analysis code SE2-ANL has been widely used in thermal-hydraulic analysis of sodium cooled fast reactors (SFR). It is used to compute the core-wide temperature profiles for a given flow distribution and to determine the assembly flow rates yielding a desired (assembly-wise) peak temperature distribution. It is a modified version of SUPERENERGY-2 coupled with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In the current version of SE2-ANL, the neutron and gamma fluxes are calculated using the finite-difference diffusion theory option of the DIF3D code with six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. In order to determine the pin power distributions more accurately by eliminating these limitations, a new coupled heating calculation procedure has recently developed based on the neutron and gamma transport calculations with the VARIANT code. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at the center of each fuel pin and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The SE2-ANL code has been modified to interface with the new, coupled neutron and gamma heating calculation scheme. In addition, several improvements have been made on the computational methods and models of SE2-ANL. To eliminate the limit on the axial mesh size due to the numerical instability problem of the current explicit scheme, the numerical scheme for discretizing the energy equation in the axial direction has been changed to the {theta}-method of the SLTHEN code. An automated flow allocation scheme has also been implemented in SE2-ANL to determine assembly flow rates such that the peak cladding mid-wall temperatures of individual orifice zones are equalized over the burn cycle. The impacts of the new heating calculation scheme based on VARIANT transport solutions on the cladding and fuel temperatures have been tested against the existing method of SE2-ANL. Using the existing and new heating calculation schemes, full-core subchannel analyses were performed for the metal core design of the 1000 MWt Advanced Burner Reactor (ABR-1000), and the resulting coolant, cladding and fuel temperatures were compared. (authors)

• OSTI ID: 23042838
• Journal Information: Transactions of the American Nuclear Society, Journal Name: Transactions of the American Nuclear Society Vol. 115; ISSN 0003-018X
• Country of Publication: United States
• Language: English

Related Subjects

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS
97 MATHEMATICS AND COMPUTING
CLADDING
FLOW RATE
FUEL PINS
GAMMA SOURCES
HEAT
METALS
NEUTRONS
POWER DISTRIBUTION
RADIATION HEATING
SODIUM COOLED REACTORS
STEADY-STATE CONDITIONS
TEMPERATURE DISTRIBUTION
THERMAL HYDRAULICS