• 한국에너지공학회:학술대회논문집
• /
• 한국에너지공학회 1993년도 추계학술발표회 초록집
• /
• pp.99-102
• /
• 1993

A consistent general order nodal method for solving the three-dimensional neutron diffusion equation in (x-y-z) geometry has been derived by using a weighted integral technique and expanding the spatial variable by the Legendre orthogonal series function. The equation set derived can be converted into any order nodal schemes. It forms a compact system for general order of nodal moments. The method utilizes fewer unknown variables in the schemes for iterative-convergence solution than other nodal methods listed in the literatures, and because the method utilizes the analytic solutions of the transverse-integrated one dimensional equations and a consistent approximation for a given spatial variable through all the solution procedures, which renders the use of an approximation for the transverse leakages no longer necessary, we can expect extremely accurate solutions and the solution would converge exactly when the mesh width is decreased or the approximation order is increased.

• Nuclear Engineering and Technology
• /
• 제32권6호
• /
• pp.605-617
• /
• 2000

Nodal transport methods are studied for the solution of two dimensional discrete-ordinates, simplified even-parity transport equation(SEP) which is known to be an approximation to the true transport equation. The polynomial expansion nodal method(PEN) and the analytic function expansion nodal method(AFEN）which have been developed for the diffusion theory are used for the solution of the discrete-ordinates form of SEP equation. Our study shows that while the PEN method in diffusion theory can directly be converted without complication, the AFEN method requires a theoretical modification due to the nonhomogeneous property of the transport equation. The numerical results show that the proposed two methods work well with the SEP transport equation with higher accuracies compared with the conventional finite difference method.

• Nuclear Engineering and Technology
• /
• 제54권5호
• /
• pp.1789-1803
• /
• 2022

RAST-F is a new full-core analysis code based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo code MCS and a multi-group nodal diffusion solver. To demonstrate the feasibility of using MCS/RAST-F for fast reactor analysis, this paper presents the coupled nodal code verification results for the MET-1000 and CAR-3600 benchmark cores. Three different multi-group cross-section calculation schemes are employed to improve the agreement between the nodal and reference solutions. The reference solution is obtained by the MCS code using continuous-energy nuclear data. Additionally, the MCS/RAST-F nodal solution is verified with results based on cross-section generated by collision probability code TULIP. A good agreement between MCS/RAST-F and reference solution is observed with less than 120 pcm discrepancy in k_eff and less than 1.2% root-mean-square error in power distribution. This study confirms the two-step approach MCS/RAST-F as a reliable tool for the three-dimensional simulation of reactor cores with fast spectrum.

• Nuclear Engineering and Technology
• /
• 제16권4호
• /
• pp.202-216
• /
• 1984

Nodal method가 소격격자 해석방법의 하나로 정립됨으로써, 계산격자가 비교적 크더라도 각 격자의 평균출력분포를 정확히 계산할 수 있게 하는 균질화변수틀 찾는 방법이 중요하게 되었다. 본 연구에서는 simplified equivalence theory와 approximate node equivalence theory의 두가지 근사방법을 가압경수형 원자로 문제에 적응하여 시험하여 보았다. 균질화계산과 노심분석계산 방법으로서 analytic nodal method에 기초를 둔 ANM 코드를 개발하였다. 여러 균질화 방법외 정확성을 KTDD 코드에 의한 reference solution과 비교하여 본 결과, 균질화 계산은 핵연료영역에서는 영역별 핵연료집합체 계산으로, baffle과 reflector의 공존 격자영역은 이들을 포함하는 color set 계산으로 수행할 수 있음을 알았다. Approximate node equivalence theory에 입각해서 approximate homogenized cross-section들과 approximate discontinuity factor들의 균질화 변수를 사용하면 출력분포와 임계도가 각각 0.8%, 0,1% 오차 범위내에서 예측되었다.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1998년도 춘계학술발표회논문집(1)
• /
• pp.79-86
• /
• 1998

The nonlinear analytic function expansion nodal (AFEN) method is applied to the solution of the time-dependent neutron diffusion equation. Since the AFEN method requires both the particular solution and the homogeneous solution to the transient fixed source problem, the derivation solution method is focused on finding the particular solution efficiently. To avoid complicated particular solutions, the source distribution is approximated by quadratic polynomials and the transient source is constructed such that the error due to the quadratic approximation is minimized. In addition, this paper presents a new two-node solution scheme that is derived by imposing the constraint of current continuity at the interface corner points. The method is verified through a series of applications to the NEACRP PWR rod ejection benchmark problem.

• Nuclear Engineering and Technology
• /
• 제53권9호
• /
• pp.2788-2802
• /
• 2021

Multigroup cross section (MG XS) generation by the UNIST in-house Monte Carlo (MC) code MCS for fast reactor analysis using nodal diffusion codes is reported. The feasibility of the approach is quantified for two sodium fast reactors (SFRs) specified in the OECD/NEA SFR benchmark: a 1000 MW_th metal-fueled SFR (MET-1000) and a 3600 MW_th oxide-fueled SFR (MOX-3600). The accuracy of a few-group XSs generated by MCS is verified using another MC code, Serpent 2. The neutronic steady-state whole-core problem is analyzed using MCS/RAST-K with a 24-group XS set. Various core parameters of interest (core k_eff, power profiles, and reactivity feedback coefficients) are obtained using both MCS/RAST-K and MCS. A code-to-code comparison indicates excellent agreement between the nodal diffusion solution and stochastic solution; the error in the core k_eff is less than 110 pcm, the root-mean-square error of the power profiles is within 1.0%, and the error of the reactivity feedback coefficients is within three standard deviations. Furthermore, using the super-homogenization-corrected XSs improves the prediction accuracy of the control rod worth and power profiles with all rods in. Therefore, the results demonstrate that employing the MCS MG XSs for the nodal diffusion code is feasible for high-fidelity analyses of fast reactors.

• Nuclear Engineering and Technology
• /
• 제53권5호
• /
• pp.1391-1402
• /
• 2021

The aim of this work is to prepare a neutron noise calculator based on the second order of average current nodal expansion method (ACNEM). Generally, nodal methods have the ability to fulfill the neutronic analysis with adequate precision using coarse meshes as large as a fuel assembly size. But, for the zeroth order of ACNEM, the accuracy of neutronic simulations may not be sufficient when coarse meshes are employed in the reactor core modeling. In this work, the capability of second order ACNEM is extended for solving the neutron diffusion equation in the frequency domain using coarse meshes. For this purpose, two problems are modeled and checked including a slab reactor and 2D BIBLIS PWR. For validating of results, a semi-analytical solution is utilized for 1D test case, and for 2D problem, the results of both forward and adjoint neutron noise calculations are exploited. Numerical results indicate that by increasing the order of method, the errors of frequency dependent coarse mesh solutions are considerably decreased in comparison to the reference. Accordingly, the accuracy of second order ACNEM can be acceptable for the neutron noise calculations by using coarse meshes in the nuclear reactor core.

• Structural Engineering and Mechanics
• /
• 제7권2호
• /
• pp.203-216
• /
• 1999

A previously proposed finite element formulation method is refined and modified to generate a new type of elements. The method is based on selecting a set of general solution modes for element formulation. The constant strain modes and higher order modes are selected and the formulation method is designed to ensure that the element will pass the basic single element test, which in turn ensures the passage of the basic patch test. If the element is to pass the higher order patch test also, the element stiffness matrix is in general asymmetric. The element stiffness matrix depends only on a nodal displacement matrix and a nodal force matrix. A symmetric stiffness matrix can be obtained by either modifying the nodal displacement matrix or the nodal force matrix. It is shown that both modifications lead to the same new element, which is demonstrated through numerical examples to be more robust than an assumed stress hybrid element in plane stress application. The method of formulation can also be used to arrive at the conforming displacement and hybrid stress formulations. The convergence of the latter two is explained from the point of view of the proposed method.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 1996년도 한국자동제어학술회의논문집(국내학술편); 포항공과대학교, 포항; 24-26 Oct. 1996
• /
• pp.1-4
• /
• 1996

It is difficult that a non-translational motion in a block is estimated by the block matching algorithm (BMA). In this paper, a nodal-displacement-based deformation model is used for this reason. This model assumes that a selected number of control nodes move freely in a block and that displacement of any interior point can be interpolated from nodal displacements. As a special case with a single node this model is equivalent to a translational model. And this model can represent more complex deformation using more nodes. We used an iterative gradient based search algorithm to estimate nodal displacement. Each iteration involves the solution of a simple linear equation. This method is called the deformable block matching algorithm (DBMA).