• Nuclear Engineering and Technology
• /
• v.56 no.4
• /
• pp.1296-1309
• /
• 2024

The aim of this work is to explore the effect of the double subdiffusion on the stability in BWRs. A BWR novel reduced order model with double subdiffusion effects: reduced order fractional model (DS-F-ROM) to describe the neutron and heat transfer processes was proposed for this study. The double subdiffusion was developed with a fractional-order two-equation model, and with different fractional-orders and relaxation times. The stability analysis was carried out using the root-locus method and change from the s to the W domain and were confirmed using the time-domain evolution of neutron flux for a unit step change in reactivity. The results obtained using the reduced fractional-order model are presented for different anomalous diffusion coefficient values. Results are compared with normal diffusion and P1 equations, which are obtained straightforwardly with DS-ROM when relaxation time tends to zero, and when the anomalous diffusion coefficient tends to one, respectively.

• Nuclear Engineering and Technology
• /
• v.53 no.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
• /
• v.21 no.3
• /
• pp.193-204
• /
• 1989

This paper develops a new method for reconstructing neutron flux distribution, that is based on the maximum entropy Principle in information theory. The Probability distribution that maximizes the entropy Provides the most unbiased objective Probability distribution within the known partial information. The partial information are the assembly volume-averaged neutron flux, the surface-averaged neutron fluxes and the surface-averaged neutron currents, that are the results of the nodal calculation. The flux distribution on the boundary of a fuel assembly, which is the boundary condition for the neutron diffusion equation, is transformed into the probability distribution in the entropy expression. The most objective boundary flux distribution is deduced using the results of the nodal calculation by the maximum entropy method. This boundary flux distribution is then used as the boundary condition in a procedure of the imbedded heterogeneous assembly calculation to provide detailed flux distribution. The results of the new method applied to several PWR benchmark problem assemblies show that the reconstruction errors are comparable with those of the form function methods in inner region of the assembly while they are relatively large near the boundary of the assembly. The incorporation of the surface-averaged neutron currents in the constraint information (that is not done in the present study) should provide better results.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.10a
• /
• pp.103-108
• /
• 1995

By applying the P$_1$ and P$_2$ equations to the operator form of a synthetic acceleration, we derive the p,-acceleration (diffusion synthetic acceleration: DSA) and P$_2$-acceleration (p$_2$SA) schemes in one dimensional slab geometry. We Fourier-analyze the derived acceleration schemes with the discrete-ordinates transport equation and showed that the DSA outperforms the P$_2$SA. These results confirm that one cannot simply assume that replacement of the DSA with a higher order approximation will lead to a better acceleration performance.

• Proceedings of the Korea Society for Energy Engineering kosee Conference
• /
• 2000.11a
• /
• pp.115-120
• /
• 2000

The analysis of flux distribution under stead-state in large power reactors with assymetry reactivity insertions requires the use of three-dimensional diffusion calculations. For the purpose, consistently formulated modern nodal methods based on higher order interface techniques have become popular tools for flux distributions in large commercial nuclear reactors. Among the earlier developments, the nodal Green's function method obtains its nodal interface equation from the transverse-integrated integral diffusion equation using a finite-medium Green's function. In this method, the outgoing current from a node surface is formulated as a response of the incoming currents and the spatially integrated neutron source within the same node. The well-known nodal expansion method is also based on an interface partial current formulation. Nodal methods high-level interface variables, i.e., interface net current and flux, may be more computationally efficient than the nodal Green's function method because they have one fewer unknown per interface. The Analytic Nodal Method(ANM), which can be classified as an interface net current technique and, was faster in solving some standard benchmark problems than the other two methods.(omitted)

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.223-229
• /
• 2020

The present work proposes a solution to the static Boltzmann transport equation approximated by the simplified P₃ (SP₃) on angular, and the analytic coarse mesh finite difference (ACMFD) for spatial variables. Multi-group SP₃-ACMFD equations in 3D rectangular geometry are solved using the GMRES solution technique. As the core time dependent analysis necessitates the solution of an eigenvalue problem for an initial condition, this work is hence devoted to development and verification of the proposed static SP₃-ACMFD solver. A 3D multi-group static diffusion solver is also developed as a byproduct of this work to assess the improvement achieved using the SP₃ technique. Static results are then compared against transport benchmarks to assess the proximity of SP₃-ACMFD solutions to their full transport peers. Results prove that the approach can be considered as an acceptable interim approximation with outputs superior to the diffusion method, close to the transport results, and with the computational costs less than the full transport approach. The work would be further generalized to time dependent solutions in Part II.

• Nuclear Engineering and Technology
• /
• v.3 no.3
• /
• pp.121-128
• /
• 1971

A unified modal-nodal expansion of tile angular distribution of neutron flux in one spatial dimension is considered, following the proposal of Harms. Several standard nodal and/or modal methods of analysis are shown to be specializations of this technique. The modal-nodal moment from of the mono-energetic transport equation with isotropic sources and scattering is derived and the infinite medium eigenvalue problem solved. The technique is shown to yield results which approximate the exact value of the inverse diffusion length in non-multiplying media more accurately than standard methods of equal or somewhat greater computational complexity.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.05a
• /
• pp.27-32
• /
• 1997

In this paper, a new multigrid method is developed to solve the reactor eigenvalue problems. The new algorithm can be used in any matrix equation concerned with the eigenvalue problem. The finite difference neutron diffusion problem is considered demonstration of the performance of the new multigrid algorithm. The numerical results show that the new multigrid algorithm works well and requires much shorter (7~10 times) computing time compaired to the production code VENTURE.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1998.05a
• /
• 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.

• 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.