• Nuclear Engineering and Technology
• /
• 제53권10호
• /
• pp.3133-3150
• /
• 2021

This work deals with the validation of a high-fidelity multiphysics system coupling the Serpent 2 Monte Carlo neutron transport code with SUBCHANFLOW, a subchannel thermalhydraulics code, and TRANSURANUS, a fuel-performance analysis code. The results for a full-core pin-by-pin burnup calculation for the ninth operating cycle of the Temelín II VVER-1000 plant, which starts from a fresh core, are presented and assessed using experimental data. A good agreement is found comparing the critical boron concentration and a set of pin-level neutron flux profiles against measurements. In addition, the calculated axial and radial power distributions match closely the values reported by the core monitoring system. To demonstrate the modeling capabilities of the three-code coupling, pin-level neutronic, thermalhydraulic and thermomechanic results are shown as well. These studies are encompassed in the final phase of the EU Horizon 2020 McSAFE project, during which the Serpent-SUBCHANFLOW-TRANSURANUS system was developed.

• Nuclear Engineering and Technology
• /
• 제53권5호
• /
• pp.1676-1685
• /
• 2021

We develop a new high-fidelity multiphysics model to simulate boron chemistry in the porous Chalk River Unidentified Deposit (CRUD) deposits. Heat transfer, capillary flow, solute transport, and chemical reactions are fully coupled. The evaporation of coolant in the deposits is included in governing equations modified by the volume-averaged assumption of wick boiling. The axial offset anomaly (AOA) of the Seabrook nuclear power plant is simulated. The new model reasonably predicts the distributions of temperature, pressure, velocity, volumetric boiling heat density, and chemical concentrations. In the thicker CRUD regions, 60% of the total heat is removed by evaporative heat transfer, causing boron species accumulation. The new model successfully shows the quantitative effect of coolant evaporation on the local distributions of boron. The total amount of boron in the CRUD layer increases by a factor of 1.21 when an evaporation-driven increase of soluble and precipitated boron concentrations is reflected. In addition, the concentrations of B(OH)₃ and LiBO₂ are estimated according to various conditions such as different CRUD thickness and porosity. At the end of the cycle in the AOA case, the total mass of boron incorporated in CRUD deposits of a reference single fuel rod is estimated to be about 0.5 mg.

• Nuclear Engineering and Technology
• /
• 제51권6호
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.