• Proceedings of the Korea Society for Energy Engineering kosee Conference
• /
• 1994.11a
• /
• pp.141-145
• /
• 1994

A computer code for solving the 3-D time-dependent multigroup neutron diffusion equation by a coupled reactor kinetics method recently developed has been developed and for evaluating its applicability in CANDU transient analysis applied to a 3-D kinetics benchmark problem which reveals non-uniform loss of coolant accident followed by an asymmetric insertion of shutdown devices. The performance of the method and code has been compared with the CANDU design code, CERBERUS, employing a finite difference improved quasistatic method.

• Nuclear Engineering and Technology
• /
• v.47 no.5
• /
• pp.608-616
• /
• 2015

The migration of silver (Ag) in silicon carbide (SiC) and $^{110m}Ag$ through SiC of irradiated tristructural isotropic (TRISO) fuel has been studied for the past three to four decades. However, there is no satisfactory explanation for the transport mechanism of Ag in SiC. In this work, the diffusion coefficients of Ag measured and/or estimated in previous studies were reviewed, and then pre-exponential factors and activation energies from the previous experiments were evaluated using Arrhenius equation. The activation energy is $247.4kJ{\cdot}mol^{-1}$ from Ag paste experiments between two SiC layers produced using fluidized-bed chemical vapor deposition (FBCVD), $125.3kJ{\cdot}mol^{-1}$ from integral release experiments (annealing of irradiated TRISO fuel), $121.8kJ{\cdot}mol^{-1}$ from fractional Ag release during irradiation of TRISO fuel in high flux reactor (HFR), and $274.8kJ{\cdot}mol^{-1}$ from Ag ion implantation experiments, respectively. The activation energy from ion implantation experiments is greater than that from Ag paste, fractional release and integral release, and the activation energy from Ag paste experiments is approximately two times greater than that from integral release experiments and fractional Ag release during the irradiation of TRISO fuel in HFR. The pre-exponential factors are also very different depending on the experimental methods and estimation. From a comparison of the pre-exponential factors and activation energies, it can be analogized that the diffusion mechanism of Ag using ion implantation experiment is different from other experiments, such as a Ag paste experiment, integral release experiments, and heating experiments after irradiating TRISO fuel in HFR. However, the results of this work do not support the long held assumption that Ag release from FBCVD-SiC, used for the coating layer in TRISO fuel, is dominated by grain boundary diffusion. In order to understand in detail the transport mechanism of Ag through the coating layer, FBCVD-SiC in TRISO fuel, a microstructural change caused by neutron irradiation during operation has to be fully considered.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1157-1164
• /
• 2017

The Monte Carlo criticality simulation of decoupled systems, as for instance in large reactor cores, has been a challenging issue for a long time. In particular, due to limited computer time resources, the number of neutrons simulated per generation is still many order of magnitudes below realistic statistics, even during the start-up phases of reactors. This limited number of neutrons triggers a strong clustering effect of the neutron population that affects Monte Carlo tallies. Below a certain threshold, not only is the variance affected but also the estimation of the eigenvectors. In this paper we will build a time-dependent diffusion equation that takes into account both spatial correlations and population control (fixed number of neutrons along generations). We will show that its solution obeys a traveling wave dynamic, and we will discuss the mechanism that explains this biasing of local tallies whenever leakage boundary conditions are applied to the system.

• Nuclear Engineering and Technology
• /
• v.53 no.11
• /
• pp.3653-3664
• /
• 2021

In PWR three-dimensional transient coupling calculation software CORCA-K, the nodal Green's function method and diagonal implicit Runge Kutta method are used to solve the spatiotemporal neutron dynamic diffusion equation, and the single-phase closed channel model and one-dimensional cylindrical heat conduction transient model are used to calculate the coolant temperature and fuel temperature. The LMW, NEACRP and PWR MOX/UO₂ benchmarks and FangJiaShan (FJS) nuclear power plant (NPP) transient control rod move cases are used to verify the CORCA-K. The effects of burnup, fuel effective temperature and ejection rate on the control rod ejection process of PWR are analyzed. The conclusions are as follows: (1) core relative power and fuel Doppler temperature are in good agreement with the results of benchmark and ADPRES, and the deviation between with the reference results is within 3.0% in LMW and NEACRP benchmarks; 2) the variation trend of FJS NPP core transient parameters is consistent with the results of SMART and ADPRES. And the core relative power is in better agreement with the SMART when weighting coefficient is 0.7. Compared with SMART, the maximum deviation is -5.08% in the rod ejection condition and while -5.09% in the control rod complex movement condition.

• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1213-1224
• /
• 2024

The liquid lead-bismuth cooled fast reactor has been in a single-phase, low-pressure, and high-temperature state for a long time during operation. Considering the requirement of calculation efficiency for long-term transient accident calculation, based on a homogeneous hydrodynamic model, one-dimensional heat conduction model, coolant flow and heat transfer model, neutron kinetics model, coolant and material properties model, this study used the fully implicit difference scheme algorithm of the convection-diffusion term to solve the basic conservation equation, to develop the transient analysis program NUSOL-LMR 2.0 for the lead-bismuth fast reactor system. The steady-state and typical design basis accidents (including reactivity introduction, loss of flow caused by main pump idling, excessive cooling, and plant power outage accidents) for the ABR have been analyzed. The results are compared with the international system analysis software ATHENA. The results indicate that the developed program can stably, accurately, and efficiently predict the transient accident response and safety characteristics of the lead-bismuth fast reactor system.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.230-237
• /
• 2020

In this part, an implicit time dependent solution is presented for the Boltzmann transport equation discretized by the analytic coarse mesh finite difference method (ACMFD) over the spatial domain as well as the simplified P₃ (SP₃) for the angular variable. In the first part of this work we proposed a SP₃-ACMFD approach to solve the static eigenvalue equations which provide the initial conditions for temp dependent equations. Having solved the 3D multi-group SP₃-ACMFD static equations, an implicit approach is resorted to ensure stability of time steps. An exponential behavior is assumed in transverse integrated equations to establish a relationship between flux moments and currents. Also, analytic integration is benefited for the time-dependent solution of precursor concentration equations. Finally, a multi-channel one-phase thermal hydraulic model is coupled to the proposed methodology. Transient equations are then solved at each step using the GMRES technique. To show the sufficiency of proposed transient SP₃-ACMFD approximation for a full core analysis, a comparison is made using transport peers as the reference. To further demonstrate superiority, results are compared with a 3D multi-group transient diffusion solver developed as a byproduct of this work. Outcomes confirm that the idea can be considered as an economic interim approach which is superior to the diffusion approximation, and comparable with transport in results.