• Proceedings of the Korean Nuclear Society Conference
• /
• 2003.10a
• /
• pp.127.1-127.1
• /
• 2003

• Proceedings of the Korean Nuclear Society Conference
• /
• 2005.10a
• /
• pp.514-515
• /
• 2005

• Nuclear Engineering and Technology
• /
• v.19 no.4
• /
• pp.317-324
• /
• 1987

The purpose of this work is to present a spatial neutron kinetics computational scheme for the analysis of space-dependent transients like rod ejection accident of pressurized water reactors. In this work modified Borresen's 1.5 group coarse mesh scheme was formulated for the neutronic computation of the space-dependent transients and applied to the analysis of hypothetical rod ejection accident of KNU no. 1 PWR core at BOC, HZP. The computational accuracy of the modified Borresen's scheme is examined by comparing calculations for core power and control rod worths with startup core physics test results. Effects of such parameters as ejected rod worths and number of delayed neution group ell transient results as well as computational efficiency are also examined. OB this basis it is suggested that the modified Borresen's method is a useful scheme for the analysis of spatial neutron transients of PWR's.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2004.10a
• /
• pp.57-58
• /
• 2004

• Proceedings of the Korean Nuclear Society Conference
• /
• 1999.10a
• /
• pp.65.2-65
• /
• 1999

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3301-3312
• /
• 2023

This study incorporates a high-fidelity transient analysis solver based on multigroup CMFD in the MOC code STREAM. Transport modeling with heterogeneous geometries of the reactor core increases computational cost in terms of memory and time, whereas the multigroup CMFD reduces the computational cost. The reactor condition does not change at every time step, which is a vital point for the utilization of CMFD. CMFD correction factors are updated from the transport solution whenever the reactor core condition changes, and the simulation continues until the end. The transport solution is adjusted once CMFD achieves the solution. The flux-weighted method is used for rod decusping to update the partially inserted control rod cell material, which maintains the solution's stability. A smaller time-step size is needed to obtain an accurate solution, which increases the computational cost. The adaptive step-size control algorithm is robust for controlling the time step size. This algorithm is based on local errors and has the potential capability to accept or reject the solution. Several numerical problems are selected to analyze the performance and numerical accuracy of parallel computing, rod decusping, and adaptive time step control. Lastly, a typical pressurized LWR was chosen to study the rod-ejection accident.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05c
• /
• pp.311-316
• /
• 1996

제어봉이탈사고시의 핵연료봉 거동을 연구로에서의 반응도사고 모사실험 결과와 기존의 핵연료 손상기준을 비교하여 분석하였다. 반응도사고시 고연소도 핵연료의 손상은 주로 PCMI 기구로 발생하는데, 고연소도에서의 피복관의 부식 및 수소화 그리고 방사선조사에 의한 연성감소와 산화층 박리로 인한 수소화합물의 국부적인 집중화로 인한 피복관의 현저한 연성감소가 주요 원인이었다. 기존의 핵연료 손상 기준에서 DNB가 일어날때 핵연료 손상이 발생한다는 가정은 낮은 핵연료엔탈피에서 핵연료 손상이 일어나는 것과 동일함을 확인하였으며, 현재까지 발표된 실험자료와 핵연료손상기구의 분석을 통해 연소도에 따른 반응도사고시의 핵연료손상기준을 예비적으로 유도하였다. 핵연료손상은 낮은 연소도에는 DNB로 발생하고 고연소도에서는 PCMI로 발생할수 있기 때문에, 과도상태에서의 고연소도 핵연료의 건전성 유지를 위해서는 피복관 산화층의 박리로 인한 수소화합물의 집중화로 피복관의 연성이 감소되는 것을 방지할 필요가 있다.

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3494-3515
• /
• 2022

This paper presents the development of a nodal diffusion code, RAST-V, and its verification and validation for VVER (vodo-vodyanoi energetichesky reactor) analysis. A VVER analytic solver has been implemented in an in-house nodal diffusion code, RAST-K. The new RAST-K version, RAST-V, uses the triangle-based polynomial expansion nodal method. The RAST-K code provides stand-alone and two-step computation modes for steady-state and transient calculations. An in-house lattice code (STREAM) with updated features for VVER analysis is also utilized in the two-step method for cross-section generation. To assess the calculation capability of the formulated analysis module, various verification and validation studies have been performed with Rostov-II, and X2 multicycles, Novovoronezh-4, and the Atomic Energy Research benchmarks. In comparing the multicycle operation, rod worth, and integrated temperature coefficients, RAST-V is found to agree with measurements with high accuracy which RMS differences of each cycle are within ±47 ppm in multicycle operations, and ±81 pcm of the rod worth of the X2 reactor. Transient calculations were also performed considering two different rod ejection scenarios. The accuracy of RAST-V was observed to be comparable to that of conventional nodal diffusion codes (DYN3D, BIPR8, and PARCS).

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.258-270
• /
• 2020

A new task of using the Jacobian-Free-Newton-Krylov (JFNK) method for the PWR core transient simulations involving neutronics, thermal hydraulics and mechanics is conducted. For the transient scenario of PWR, normally the Picard iteration of the coupled coarse-mesh nodal equations and parallel channel TH equations is performed to get the transient solution. In order to solve the coupled equations faster and more stable, the Newton Krylov (NK) method based on the explicit matrix was studied. However, the NK method is hard to be extended to the cases with more physics phenomenon coupled, thus the JFNK based iteration scheme is developed for the nodal method and parallel-channel TH method. The local gap conductance is sensitive to the gap width and will influence the temperature distribution in the fuel rod significantly. To further consider the local gap conductance during the transient scenario, a 1D mechanics model is coupled into the JFNK scheme to account for the fuel thermal expansion effect. To improve the efficiency, the physics-based precondition and scaling technique are developed for the JFNK iteration. Numerical tests show good convergence behavior of the iterations and demonstrate the influence of the fuel thermal expansion effect during the rod ejection problems.

• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.