• Nuclear Engineering and Technology
• /
• 제49권5호
• /
• pp.920-927
• /
• 2017

With ever-advancing computer technology, the Monte Carlo (MC) neutron transport calculation is expanding its application area to nuclear reactor transient analysis. Dynamic MC (DMC) neutron tracking for transient analysis requires efficient algorithms for delayed neutron generation, neutron population control, and initial condition modeling. In this paper, a new MC steady-state simulation method based on time-dependent MC neutron tracking is proposed for steady-state initial condition modeling; during this process, prompt neutron sources and delayed neutron precursors for the DMC transient simulation can easily be sampled. The DMC method, including the proposed time-dependent DMC steady-state simulation method, has been implemented in McCARD and applied for two-dimensional core kinetics problems in the time-dependent neutron transport benchmark C5G7-TD. The McCARD DMC calculation results show good agreement with results of a deterministic transport analysis code, nTRACER.

• Nuclear Engineering and Technology
• /
• 제55권7호
• /
• pp.2712-2722
• /
• 2023

For an accurate and efficient time-dependent Monte Carlo (TDMC) neutron transport analysis, several advanced methods are newly developed and implemented in the Seoul National University Monte Carlo code, McCARD. For an efficient control of the neutron population, a dynamic weight window method is devised to adjust the weight bounds of the implicit capture in the time bin-by-bin TDMC simulations. A moving geometry module is developed to model a continuous insertion or withdrawal of a control rod. Especially, the history-based batch method for the TDMC calculations is developed to predict the unbiased variance of a bin-wise mean estimate. The developed methods are verified for three-dimensional problems in the C5G7-TD benchmark, showing good agreements with results from a deterministic neutron transport analysis code, nTRACER, within the statistical uncertainty bounds. In addition, the TDMC analysis capability implemented in McCARD is demonstrated to search the optimum detector positions for the pulsed-neutron-source experiments in the Kyoto University Critical Assembly and AGN201K.

• Nuclear Engineering and Technology
• /
• 제53권12호
• /
• pp.3861-3878
• /
• 2021

The steady-state simplified spherical harmonics equations (SP_N equations) are a higher order approximation to the neutron transport equations than the neutron diffusion equation that also have reasonable computational demands. This work extends these results for the analysis of transients by comparing of two formulations of time-dependent SP_N equations considering different treatments for the time derivatives of the field moments. The first is the full system of equations and the second is a diffusive approximation of these equations that neglects the time derivatives of the odd moments. The spatial discretization of these methodologies is made by using a high order finite element method. For the time discretization, a semi-implicit Euler method is used. Numerical results show that the diffusive formulation for the time-dependent simplified spherical harmonics equations does not present a relevant loss of accuracy while being more computationally efficient than the full system.

• Nuclear Engineering and Technology
• /
• 제56권3호
• /
• pp.916-921
• /
• 2024

The control rod worth sensitivity to the neutron detector position in the rod drop experiment is studied by the time-dependent Monte Carlo (TDMC) neutron transport calculations for AGN-201K educational reactor and the Kyoto University Critical Assembly. The TDMC simulations of the rod drop experiments are conducted by the Seoul National University Monte Carlo (MC) code, McCARD, yielding time-dependent neutron densities at detector positions. The detector-position-dependent results of the total control rod worth calculated by the extrapolation, the integral counting, and the inverse methods are compared with the numerical reference using the MC eigenvalue calculations and the experimental results. From these comparisons, it is observed that the total control rod worth can be estimated with a considerable difference depending on the detector position through the rod drop experiment. The proposed TDMC simulation of the rod drop experiment can be applied for searching a better detector position or quantifying a bias for the control rod worth measurement.

• Nuclear Engineering and Technology
• /
• 제43권3호
• /
• pp.243-256
• /
• 2011

This paper considers the concept of analog computing based on a cellular neural network (CNN) paradigm to simulate nuclear reactor dynamics using a time-dependent second order form of the neutron transport equation. Instead of solving nuclear reactor dynamic equations numerically, which is time-consuming and suffers from such weaknesses as vulnerability to transient phenomena, accumulation of round-off errors and floating-point overflows, use is made of a new method based on a cellular neural network. The state-of-the-art shows the CNN as being an alternative solution to the conventional numerical computation method. Indeed CNN is an analog computing paradigm that performs ultra-fast calculations and provides accurate results. In this study use is made of the CNN model to simulate the space-time response of scalar flux distribution in steady state and transient conditions. The CNN model also is used to simulate step perturbation in the core. The accuracy and capability of the CNN model are examined in 2D Cartesian geometry for two fixed source problems, a mini-BWR assembly, and a TWIGL Seed/Blanket problem. We also use the CNN model concurrently for a typical small PWR assembly to simulate the effect of temperature feedback, poisons, and control rods on the scalar flux distribution.

• Nuclear Engineering and Technology
• /
• 제4권2호
• /
• pp.102-108
• /
• 1972

등방성이고 단면적이 상수인 경우의 지발 중성자를 가진 시간 종속 중성자 수송 방정식이 해석적으로 풀어지고 있다. 두 개로 구분된 시간 영역에 있어서의 방정식이 점근적 방법에 의하여 원래의 수송 방정식으로부터 얻어 지고 있다. 각 시간 영역에 있어서의 근사해는 중성자 속도의 역수 정도로 시간에 있어서 균일하게 유용하다는 것이 보여 지고 있다.

• Nuclear Engineering and Technology
• /
• 제2권4호
• /
• pp.263-268
• /
• 1970

시간과 에너지에 종속된 중성자 전도 방정식에 나타나는 연속 에너지 전도 연산자의 스펙트럼(Spectrum)을 분석했다. 스펙트럼에 관한 4가지 정리를 증명하고 일반화된 Mellin 에너지변화의 Convolution 정리를 얻었다. 또한 최종해에 필요한 완전성정리를 증명하고 점근적으로 가장 우세한 시간붕괴상수 1 - c를 발견하였다.

• 에너지공학
• /
• 제15권4호
• /
• pp.250-255
• /
• 2006

원자력 분야에서 중성자 수송계산을 위해 개발되어 널리 사용되는 방향차분법을 시간 종속 복사 열전달식의 해를 구하는데 적용하였다. 광자의 방향별 밀도를 자체수반형 2계 편미분방정식으로 나타내어 해의 안정성을 높였고 매질의 온도방정식의 비선형성은 다단계 선형화법을 사용하여 근사하였다. 본 연구에서 개발된 해법을 전형적인 Marshak wave 문제에 적용하였고 계산 결과를 기존의 Monte Carlo의 계산결과와 비교하여 그 우월성을 보였다.

• Nuclear Engineering and Technology
• /
• 제2권4호
• /
• pp.249-254
• /
• 1970

고 에너지 영역에서 비탄성핵을 가진 충격중성자 전도문제에서 시간과 에너지의 과도적인 특성은 연속 에너지 전도연산자로서 면밀히 해결했다. 점근적으로 가장 우세한 하나의 고정된 시간고유치를 발견했다. 완전한 해는 다음 세 부분으로 되어 있다. 1. 점근적으로 가장 우세한 시간과 에너지의 분리 가능형. 2. 과도적으로 감속하는 순 에너지에 의해서 결정된 분리 불가능형. 3. 점근적으로 무시할 수 있는 시간과 에너지의 혼합에 의하여 결정된 분리 불가능형.

• Nuclear Engineering and Technology
• /
• 제55권5호
• /
• pp.1593-1603
• /
• 2023

In deterministic and Monte Carlo transport codes, b-delayed emission is included using a group structure where all of the precursors are grouped together in 6 groups or families, but given the increase in computational power, nowadays there is no reason to keep this structure. Furthermore, there have been recent efforts to compile and evaluate all the available b-delayed neutron emission data and to measure new and improved data on individual precursors. In order to be able to perform a transient Monte Carlo simulation, data from individual precursors needs to be implemented in a transport code. This work is the first step towards the development of a tool to explore the effect of individual precursors in a fissile system. In concrete, individual precursor data is included by expanding the capabilities of the open source Monte Carlo code OpenMC. In the modified code - named Time Dependent OpenMC or OpenMC(TD)- time dependency related to β-delayed neutron emission was handled by using forced decay of precursors and combing of the particle population. The data for continuous energy neutron cross-sections was taken from JEFF-3.1.1 library. Regarding the data needed to include the individual precursors, cumulative yields were taken from JEFF-3.1.1 and delayed neutron emission probabilities and delayed neutron spectra were taken from ENDF-B/VIII.0. OpenMC(TD) was tested in a monoenergetic system, an energy dependent unmoderated system where the precursors were taken individually or in a group structure, and in a light-water moderated energy dependent system, using 6-groups, 50 and 40 individual precursors. Neutron flux as a function of time was obtained for each of the systems studied. These results show the potential of OpenMC(TD) as a tool to study the impact of individual precursor data on fissile systems, thus motivating further research to simulate more complex fissile systems.