• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05a
• /
• pp.34-39
• /
• 1996

A consistent general order nodal method for solving the 3-D neutron diffusion equation in (x-y-z) geometry has ben derived by using a weighted integral technique and expanding the spatial variables 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 the analytic solutions of the transverse-integrated quasi -one dimensional equations and a consistent expansion for the spatial variables so that it renders the use of an approximation for the transverse leakages no necessary. Thus, we can expect extremely accurate solutions and the solution would converge exactly when the mesh width is decreased or the approximation order is increased since the equation set is consistent mathematically.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.16 no.2
• /
• pp.203-209
• /
• 2018

To correctly predict the neutron behavior based on diffusion calculations, it is necessary to adopt well-specified boundary conditions using suitable diffusion approximations to transport boundary conditions. Boundary conditions such as the zero net-current, the Marshak, the Mark, the zero scalar flux, and the Albedo condition have been used extensively in diffusion theory to approximate the reflective and vacuum conditions in transport theory. In this paper, we derive and analyze these conditions to prove their mathematical validity and to understand their physical implications, as well as their relationships with one another. To show the validity of these diffusion boundary conditions, we solve a sample problem. The results show that solutions of the diffusion equation with these well-formulated boundary conditions are very close to the solution of the transport equation with transport boundary conditions.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1996.05a
• /
• pp.131-136
• /
• 1996

An analytic nodal expansion method has been derived for the multigroup neutron diffusion equation in 2-D cylindrical(R-Z) coordinate. In this method we used the second order Legendre polynomials for source, and transverse leakage, and then the diffusion eqaution was solved analytically. This formalism has been applied to 2-D LWR model. $textsc{k}$$_{eff}$, power distribution, and computing time have been compared with those of ADEP code(finite difference method). The benchmark showed that the analytic nodal expansion method in R-Z coordinate has good accuracy and quite faster than the finite difference method. This is another merit of using R-Z coordinate in that the transverse integration over surfaces is better than the linear integration over length. This makes the discontinuity factor useless.s.

• Proceedings of the Korea Society for Energy Engineering kosee Conference
• /
• 1993.11a
• /
• 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
• /
• v.52 no.3
• /
• pp.485-498
• /
• 2020

A diffusion synthetic acceleration (DSA) technique for the S_N transport equation discretized with the linear discontinuous expansion method with subcell balance (LDEM-SCB) on unstructured tetrahedral meshes is presented. The LDEM-SCB scheme solves the transport equation with the discrete ordinates method by using the subcell balances and linear discontinuous expansion of the flux. Discretized DSA equations are derived by consistently discretizing the continuous diffusion equation with the LDEM-SCB method, however, the discretized diffusion equations are not fully consistent with the discretized transport equations. In addition, a fine mesh rebalance (FMR) method is devised to accelerate the discretized diffusion equation coupled with the preconditioned conjugate gradient (CG) method. The DSA method is applied to various test problems to show its effectiveness in speeding up the iterative convergence of the transport equation. The results show that the DSA method gives small spectral radii for the tetrahedral meshes having various minimum aspect ratios even in highly scattering dominant mediums for the homogeneous test problems. The numerical tests for the homogeneous and heterogeneous problems show that DSA with FMR (with preconditioned CG) gives significantly higher speedups and robustness than the one with the Gauss-Seidel-like iteration.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1259-1268
• /
• 2017

The alpha eigenvalue problem in multigroup neutron diffusion is studied with particular attention to the theoretical analysis of the model. Contrary to previous literature results, the existence of eigenvalue and eigenflux clustering is investigated here without the simplification of a unique fissile isotope or a single emission spectrum. A discussion about the negative decay constants of the neutron precursors concentrations as potential eigenvalues is provided. An in-hour equation is derived by a perturbation approach recurring to the steady state adjoint and direct eigenvalue problems of the effective multiplication factor and is used to suggest proper detection criteria of flux clustering. In spite of the prior work, the in-hour equation results give a necessary and sufficient condition for the existence of the eigenvalue-eigenvector pair. A simplified asymptotic analysis is used to predict bands of accumulation of eigenvalues close to the negative decay constants of the precursors concentrations. The resolution of the problem in one-dimensional heterogeneous problems shows numerical evidence of the predicted clustering occurrences and also confirms previous theoretical analysis and numerical results.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1479-1486
• /
• 2019

The purpose of the present study is to develop the 3D static and noise simulator based on Galerkin Finite Element Method (GFEM) using the unstructured hexahedral elements. The 3D, 2G neutron diffusion and noise equations are discretized using the unstructured hexahedral by considering the linear approximation of the shape function in each element. The validation of the static calculation is performed via comparison between calculated results and reported data for the VVER-1000 benchmark problem. A sensitivity analysis of the calculation to the element type (unstructured hexahedral or tetrahedron elements) is done. Finally, the neutron noise calculation is performed for the neutron noise source of type of variable strength using the Green function technique. It is shown that the error reduction in the static calculation is considerable when the unstructured tetrahedron elements are replaced with the hexahedral ones. Since the neutron flux distribution and neutron multiplication factor are appeared in the neutron noise equation, the more accurate calculation of these parameters leads to obtaining the neutron noise distribution with high accuracy. The investigation of the changes of the neutron noise distribution in axial direction of the reactor core shows that the 3D neutron noise analysis is required instead of 2D.

• 전기의세계
• /
• v.29 no.4
• /
• pp.256-259
• /
• 1980

A singular pertubation theory is applied to obtain an approximate solution for suboptimal control of nuclear reactors with spatially distributed parameters. The inverse of the neutron velocity is regarded as a small perturbing parameter, and the model, adopted for simplicity, is a cylindrically symmetrical reactor whose dynamics are described by the one group diffusion equation with one delayed neutron group. The Helmholtz mode expansion is used for the application of the optimal theory for lumped parameter systems to the spatially distributed parameter systems. An asymptotic expansion of the feedback gain matrix is obtained with construction of the boundary layer correction up to the first order.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2004.10a
• /
• pp.151-152
• /
• 2004

An acceleration technique combined with the discrete ordinates method which has been widely used in the solution of neutron transport phenomena is applied to the solution of radiative transfer equation. The self-adjoint form of the second order radiation intensity equation is used to enhance the stability of the solution, and a new linearization method is developed to avoid the nonlinearity of the material temperature equation. This new acceleration method is applied to the well known Marshak wave problem, and the numerical result is compared with that of a non-accelerated calculation

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