• Nuclear Engineering and Technology
• /
• v.41 no.5
• /
• pp.677-690
• /
• 2009

A previous study demonstrated that the two-row fuel assembly has much more favorable neutron-physical and thermal-hydraulic behavior than the conventional one-row fuel assemblies. Based on the newly developed two-row fuel assembly, an SCWR core is proposed and analyzed. The performance of the proposed core is investigated with 3-D coupled neutron-physical and thermal-hydraulic calculations. During the coupling procedure, the thermal-hydraulic behavior is analyzed using a sub-channel analysis code and the neutron-physical performance is computed with a 3-D diffusion code. This paper presents the main results achieved thus far related to the distribution of some neutronic and thermal-hydraulic parameters. It shows that with adjustment of the coolant and moderator mass flow in different assemblies, promising neutron-physical and thermal-hydraulic behavior of the SCWR core is achieved. A sensitivity study of the heat transfer correlation is also performed. Since the pin power in fuel assemblies can be non-uniform, a sub-channel analysis is necessary in order to investigate the detailed distribution of thermal-hydraulic parameters in the hottest fuel assembly. The sub-channel analysis is performed based on the bundle averaged parameters obtained with the core analysis. With the sub-channel analysis approach, more precise evaluation of the hot channel factor and maximum cladding surface temperature can be achieved. The difference in the results obtained with both the sub-channel analysis and the fuel assembly homogenized method confirms the importance of the sub-channel analysis.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2004.10a
• /
• pp.63-64
• /
• 2004

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

• Journal of Magnetics
• /
• v.13 no.3
• /
• pp.97-101
• /
• 2008

We investigated the exchange bias fields at the NiFe/FeMn and FeMn/CoFe interfaces in 18.9-nm NiFe/15.0-nm FeMn/17.6-nm CoFe trilayer thin films as the annealing temperature was varied from room temperature to $250^{\circ}C$ in a vacuum for 1 hour in a magnetic field of 150 Oe. Interestingly, magnetic hysteresis (M-H) measurements showed that NiFe/FeMn/CoFe trilayer thin films exhibited a completely contrasting variation of the exchange bias fields at both the NiFe/FeMn and FeMn/CoFe interfaces with annealing temperatures. High-angle X-ray diffraction (XRD) measurements indicated the absence of any discernible effect of thermal treatment on the NiFe(111) and FeMn(111) peaks. The compositional depth profile obtained from X-ray photoelectron spectroscopy (XPS) results presented the asymmetric compositional depth profiles of the Mn and Fe atoms throughout the FeMn layer. We contend that this asymmetric compositional depth profile and the preferential Mn diffusion into the NiFe layer, compared to that into the CoFe layer, are conclusive experimental evidence of the contrasting variation of the exchange bias fields at two interfaces having a common polycrystalline FeMn(111) layer.

• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.

• 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.13 no.4
• /
• pp.264-276
• /
• 1981

A simulation procedure which can represent time-dependent nuclear characteristics of TRIGA Mark-III reactor is developed. CITATION, a multi-group diffusion-depletion program, has been utilized as calculational tool. The group structure employed in this study consists of 7 groups: -3-fast and 4-thermal-which is conventionally utilized in TRIGA type reactor analysis. Three-dimensional nuclear characteristics are synthesized by combining results from two-dimensional plane calculation and two-dimensional cylinder calculation, since direct three-dimensional approach is not yet possible. An effort ia made to develope a method which can extract effective zone and group dependent bucklings by neutron diffusion theory rather than conventional zone and/or group independent Ducklings by neutron transport theory, since neutron leakage is quite high for small core such as research reactors. It is turned out that the method developed in this study gives satisfactory results. The calculation is performed under assumptions that all control rods are fully withdrawn, that no samples are inserted in the irradiation holes and that the core is located in the center of the reactor pool. Burnup-dependent variation of core excess reactivity, time dependent change of Xe-135 poisoning and reactivity worth of rotary specimen rack are calculated and compared with operation records. Neutron flux and power distribution as well as neutron spectrum in each irradiation .facility are presented.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1133-1138
• /
• 2003

The pressure tubes, which contain high temperature heavy water and fuel, are within the core of a CANDU nuclear reactor, and are thus subjected to high stresses, temperature gradient, and neutron flux. Further, it is well known that pressure tubes of cold-worked Zr-2.5Nb materials result in hydrogen diffusion, which create fully-hydrided regions (frequently called Blister). Thus a proper investigation of hydrogen diffusion within zirconium-alloy nuclear components, such as CANDU pressure tube and fuel channels is essential to predict the structural integrity of these components. In this respect, this paper presents numerical investigation of hydrogen diffusion to quantify the hydrogen concentration for blister growth of CANDU pressure tube. For this purpose, coupled temperature-hydrogen diffusion analyses are performed by means of two-dimensional finite element analysis. Comparison of predicted temperature field and blister with published test data shows good agreement.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.2
• /
• pp.189-195
• /
• 2004

The pressure tubes, which contain high temperature heavy water and fuel, are within the core of a CANDU nuclear reactor, and are thus subjected to high stresses, temperature gradient, and neutron flux. Further, it is well known that pressure tubes of cold-worked Zr-2.5Nb materials result in hydrogen diffusion, which create fully-hydrided regions (frequently called Blister). Thus a proper investigation of hydrogen diffusion within zirconium-alloy nuclear components, such as CANDU pressure tube and fuel channels is essential to predict the structural integrity of these components. In this respect, this paper presents numerical investigation of hydrogen diffusion to quantify the hydrogen concentration fur blister growth of CANDU pressure tube. For this purpose, coupled temperature-hydrogen diffusion analyses are performed by means of two-dimensional finite element analysis. Comparison of predicted temperature field and blister with published test data shows good agreement.

• Nuclear Engineering and Technology
• /
• v.46 no.4
• /
• pp.513-520
• /
• 2014

In this paper, three types of operational and industrial absorbers used at research reactors, including Ag-In-Cd alloy, $B_4C$, and Hf are selected for sensitivity analyses. Their integral effects on the main neutronic core parameters important to safety issues are investigated. These parameters are core excess reactivity, shutdown margin, total reactivity worth of control rods, thermal neutron flux, power density distribution, and Power Peaking Factor (PPF). The IAEA 10 MW benchmark core is selected as the case study to verify calculations. A two-dimensional, three-group diffusion model is selected for core calculations. The well-known WIMS-D4 and CITATION reactor codes are used to carry out these calculations. It is found that the largest shutdown margin is gained using the $B_4C$; also the lowest PPF is gained using the Ag-In-Cd alloy. The maximum point power densities belong to the inside fuel regions surrounding the central flux trap (irradiation position), surrounded by control fuel elements, and the peripheral fuel elements beside the graphite reflectors. The greatest and least fluctuation of the point power densities are gained by using $B_4C$ and Ag-In-Cd alloy, respectively.