• Nuclear Engineering and Technology
• /
• v.28 no.3
• /
• pp.320-330
• /
• 1996

The streaming rays(SR) method is improved and extended to multigroup, anisotropic scattering, and three-dimensional angular space(x-y-z(infinite))problems. This method is applied to the shielding problems in which the ray effect occurs seriously. For verification, the results of MORSE-CG code are used as reference solution and the results of TWODANT code are compared. The results show that solutions of the SR method are much better than those of the TWODANT code and are in good agreement with those of the MORSE-CG code. Also, to reduce computing time, two acceleration algorithms are implemented in the SR method : the standard coarse-mesh rebalance and a new angular two-grid acceleration.

• Nuclear Engineering and Technology
• /
• v.55 no.1
• /
• pp.324-329
• /
• 2023

Different types of deterministic solution methods were used to solve neutron transport equations corresponding to half-space and slab albedo problems. In these types of solution methods, in addition to the error of the numerical solutions, the obtained results contain truncation and discretization errors. In the present work, a non-analog Monte Carlo method is provided to simulate the half-space and slab albedo problems with Inönü and Anlı-Güngör strongly anisotropic scattering functions. For each scattering function, the sampling method of the direction of the scattered neutrons is presented. The effects of different beams with different angular dependencies and the effects of different scattering parameters on the reflection probability are investigated using the developed Monte Carlo method. The validity of the Monte Carlo method is also confirmed through the comparison with the published data.

• Nuclear Engineering and Technology
• /
• v.44 no.2
• /
• pp.161-176
• /
• 2012

McCARD is a Monte Carlo (MC) neutron-photon transport simulation code. It has been developed exclusively for the neutronics design of nuclear reactors and fuel systems. It is capable of performing the whole-core neutronics calculations, the reactor fuel burnup analysis, the few group diffusion theory constant generation, sensitivity and uncertainty (S/U) analysis, and uncertainty propagation analysis. It has some special features such as the anterior convergence diagnostics, real variance estimation, neutronics analysis with temperature feedback, $B_1$ theory-augmented few group constants generation, kinetics parameter generation and MC S/U analysis based on the use of adjoint flux. This paper describes the theoretical basis of these features and validation calculations for both neutronics benchmark problems and commercial PWR reactors in operation.

• Nuclear Engineering and Technology
• /
• v.42 no.5
• /
• pp.520-545
• /
• 2010

This paper presents the neutronics modeling capabilities of the fast reactor simulation system SHARP, which ANL is developing as part of the U.S. DOE's NEAMS program. We discuss the three transport solvers (PN2ND, SN2ND, and MOCFE) implemented in the UNIC code along with the multigroup cross section generation code $MC^2$-3. We describe the solution methods and modeling capabilities, and discuss the improvement needs for each solver, focusing on massively parallel computation. We present the performance test results against various benchmark problems and ZPR-6 and ZPPR critical experiments. We also discuss weak and strong scalability results for the SN2ND solver on the ZPR-6 critical assembly benchmarks.

• Journal of Radiation Protection and Research
• /
• v.15 no.2
• /
• pp.113-122
• /
• 1990

After four years of planning, equipment acquisition, facility construction and beam testing, the KCCH cyclotron facility was put into operation in November1986. Now the KCCH cyclotron(MC-50) has been used for four years in neutron therapy and radioisotope production. Up to December 1989, 179(1852 sessions) patient have undergone neutron therapy. Radioisotope production for nuclear medicine use was started from March 1989 after extensive work to overcome target transport, target melting, beam diagnostic and chemical processing problems. This status report introduces the cyclotron facility, and the experiences of neutron therapy and isotope production with the MC-50 cyclotron. Besides, the operation results and the general troubles of the MC-50 during 1989 are summarized. Total operation time was 1252.5 hours. Four hundred hours were used for neutron therapy of 599 treatment sessions and 832.5 hours for radioisotope production. Total amount of produced raioisotope was 1695 mCi(Ga-67 : 1478mCi, Tl-201 : 107 mCi, I-123 : 25mCi, In-111 : 85mCi). Twenty hours were used for scheduled beam testing. In 1989, 882% of the planned operation were performed on schedule and this rats is improved remarkably compared to 71.0% in 1988.

• Nuclear Engineering and Technology
• /
• v.53 no.4
• /
• pp.1079-1087
• /
• 2021

NTP-ERSN is a package developed for solving the multigroup form of the discrete ordinates, characteristics and collision probability of the Boltzmann transport equation in one-dimensional cartesian geometry, by combining pin cells. In this work, C5G7 MOX benchmark is used to verify the accuracy and efficiency of NTP-ERSN package, by treating reactor core problems without spatial homogenization. This benchmark requires solutions in the form of normalized pin powers as well as the vectors and the eigenvalue. All NTP-ERSN simulations are carried out with appropriate spatial and angular approximations. A good agreement between NTP-ERSN results with those obtained with OpenMC calculation code for seven energy groups. In addition, our studies about angular and mesh refinements are carried out to produce better quality solution. Moreover, NTP-ERSN GUI has also been updated and adapted to python 3 programming language.

• Nuclear Engineering and Technology
• /
• v.41 no.3
• /
• pp.279-286
• /
• 2009

Improving over a previous study [1], this paper provides a Monte Carlo method for the heat conduction analysis of problems with complicated geometry (such as a pebble with dispersed fuel particles). The method is based on the theoretical results of asymptotic analysis of neutron transport equation. The improved method uses an appropriate boundary layer correction (with extrapolation thickness) and a scaling factor, rendering the problem more diffusive and thus obtaining a heat conduction solution. Monte Carlo results are obtained for the randomly distributed fuel particles of a pebble, providing realistic temperature distributions (showing the kernel and graphite-matrix temperatures distinctly). The volumetric analytic solution commonly used in the literature is shown to predict lower temperatures than those of the Monte Carlo results provided in this paper.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1326-1338
• /
• 2017

The continuous energy Monte Carlo neutron transport code, MC21, was coupled to the CTF subchannel thermal-hydraulics code using a combination of Consortium for Advanced Simulation of Light Water Reactors (CASL) tools and in-house Python scripts. An MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 6 demonstrated good agreement with MC21/COBRA-IE and VERA solutions. The MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 7, Watts Bar Unit 1 at beginning of cycle hot full power equilibrium xenon conditions, is the first published coupled Monte Carlo neutronics/subchannel T-H solution for this problem. MC21/CTF predicted a critical boron concentration of 854.5 ppm, yielding a critical eigenvalue of $0.99994{\pm}6.8E-6$ (95% confidence interval). Excellent agreement with a VERA solution of Problem 7 was also demonstrated for integral and local power and temperature parameters.

• Nuclear Engineering and Technology
• /
• v.46 no.6
• /
• pp.783-796
• /
• 2014

The coarse mesh finite difference (CMFD) method is applied to the discontinuous finite element method based discrete ordinate calculation for source convergence acceleration. The three-dimensional (3-D) DFEM-Sn code FEDONA is developed for general geometry applications as a framework for the CMFD implementation. Detailed methods for applying the CMFD acceleration are established, such as the method to acquire the coarse mesh flux and current by combining unstructured tetrahedron elements to rectangular coarse mesh geometry, and the alternating calculation method to exchange the updated flux information between the CMFD and DFEM-Sn. The partial current based CMFD (p-CMFD) is also implemented for comparison of the acceleration performance. The modified p-CMFD method is proposed to correct the weakness of the original p-CMFD formulation. The performance of CMFD acceleration is examined first for simple two-dimensional multigroup problems to investigate the effect of the problem and coarse mesh sizes. It is shown that smaller coarse meshes are more effective in the CMFD acceleration and the modified p-CMFD has similar effectiveness as the standard CMFD. The effectiveness of CMFD acceleration is then assessed for three-dimensional benchmark problems such as the IAEA (International Atomic Energy Agency) and C5G7MOX problems. It is demonstrated that a sufficiently converged solution is obtained within 7 outer iterations which would require 175 iterations with the normal DFEM-Sn calculations for the IAEA problem. It is claimed that the CMFD accelerated DFEM-Sn method can be effectively used in the practical eigenvalue calculations involving general geometries.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.732-740
• /
• 2021

Scattering resonance of medium mass nuclides leads complex spectrum in the fast reactor, which requires thousands of energy groups in the spectrum calculation. When the broad-group cross sections are collapsed, reaction rate cannot be completely conserved. To eliminate the error from energy collapsing, the Super-homogenization method in energy collapsing (ESPH) was employed in the fast reactor code SARAX. An ESPH factor was derived based on the ESPH-corrected SN transport equation. By applying the factor in problems with reflective boundary condition, both the effective multiplication factor and reaction rate were conserved. The fixed-source iteration was used to ensure the stability of ESPH iteration. However, in the energy collapsing process of SARAX, the vacuum boundary condition was adopted, which was necessary for fast reactors with strong heterogeneity. To further reduce the error caused by leakage, an additional conservation factor was proposed to correct the neutron current in energy collapsing. To evaluate the performance of ESPH with conservation factor, numerical benchmarks of fast reactors were calculated. The results of broad-group calculation agreed well with the direct full-core Monte-Carlo calculation, including the effective multiplication factor, radial power distribution, total control rod worth and sodium void worth.