• Nuclear Engineering and Technology
• /
• 제49권6호
• /
• pp.1240-1249
• /
• 2017

An essential component of the neutron transport solver is the resonance self-shielding calculation used to determine equivalence cross sections. The neutron transport code, MPACT, is currently using the subgroup self-shielding method, in which the method of characteristics (MOC) is used to solve purely absorbing fixed-source problems. Recent efforts incorporating multigroup kernels to the MOC solvers in MPACT have reduced runtime by roughly $2{\times}$. Applying the same concepts for self-shielding and developing a novel lumped parameter approach to MOC, substantial improvements have also been made to the self-shielding computational efficiency without sacrificing any accuracy. These new multigroup and lumped parameter capabilities have been demonstrated on two test cases: (1) a single lattice with quarter symmetry known as VERA (Virtual Environment for Reactor Applications) Progression Problem 2a and (2) a two-dimensional quarter-core slice known as Problem 5a-2D. From these cases, self-shielding computational time was reduced by roughly $3-4{\times}$, with a corresponding 15-20% increase in overall memory burden. An azimuthal angle sensitivity study also shows that only half as many angles are needed, yielding an additional speedup of $2{\times}$. In total, the improvements yield roughly a $7-8{\times}$ speedup. Given these performance benefits, these approaches have been adopted as the default in MPACT.

• Nuclear Engineering and Technology
• /
• 제53권10호
• /
• pp.3264-3274
• /
• 2021

In a severe accident of light water reactor (LWR), molten core material (corium) can be released into the wet cavity, and a fuel-coolant interaction (FCI) can occur. The molten jet with high speed is broken and fragmented into small debris, which may cause a steam explosion or a molten core concrete interaction (MCCI). Since the premixing stage where the jet breakup occurs has a large impact on the severe accident progression, the understanding and evaluation of the jet breakup phenomenon are highly important. Therefore, in this study, the jet breakup simulations were performed using the Smoothed Particle Hydrodynamics (SPH) method which is a particle-based Lagrangian numerical method. For the multi-fluid system, the normalized density approach and improved surface tension model (CSF) were applied to the in-house SPH code (single GPU-based SOPHIA code) to improve the calculation accuracy at the interface of fluids. The jet breakup simulations were conducted in two cases: (1) jet breakup without structures, and (2) jet breakup with structures (control rod guide tubes). The penetration depth of the jet and jet breakup length were compared with those of the reference experiments, and these SPH simulation results are qualitatively and quantitatively consistent with the experiments.

• Environmental Engineering Research
• /
• 제24권3호
• /
• pp.513-520
• /
• 2019

In this investigation, nano ZnO was sonochemically synthesized by a novel method using a methionine precursor. A narrow size distribution (41-50 nm) of nano ZnO was achieved that was immobilized on perlite and applied as a catalyst in catalytic ozonation. The catalyst was characterized by fourier transform infrared spectroscopy, BET surface area, and field emission scanning electron microscope. The ozonation of recalcitrant Remazol black 5 (RB5) di-azo dye solution by means of the synthesized catalyst was investigated in a bubble column slurry reactor. The influence of pH values (7, 9, 11), catalyst dosage (8, 12, 15, $20g\;L^{-1}$) and reaction time (10, 20, 30, 60 min) was investigated. Although the dye color was completely removed by single ozonation at a higher reaction time, the applied nanocatalyst improved the dye declorination kinetics. Also, the degradation of the hazardous aromatic fraction of the dye was enhanced five-times by catalytic ozonation at a low reaction time (10 min) and a neutral pH. The second-order kinetics was best fitted in terms of both RB5 color and its aromatic fraction removal. The total organic carbon analysis indicated a significant improvement in the mineralization of RB5 by catalytic ozonation using the nano-ZnO/perlite catalyst.

• Nuclear Engineering and Technology
• /
• 제51권6호
• /
• pp.1525-1531
• /
• 2019

The indigenous nuclear power program of India is based mainly on a series of Pressurised Heavy Water Reactors (PHWRs). A burst correlation for Indian PHWR fuel claddings has been developed and empirical burst parameters are determined. The burst correlation is developed from data available in literature for single-rod transient burst tests performed on Indian PHWR claddings in inert environment. The heating rate and internal overpressure were in the range of 7 K/s-73 K/s and 3 bar-80 bar, respectively, during the burst tests. A burst criterion for inert environment, which assumes that deformation is controlled by steady state creep, has been developed using the empirical burst parameters. The burst criterion has been validated with experimental data reported in literature and the prediction of burst parameters is in a fairly good agreement with the experimental data. The burst criterion model reveals that increasing the heating rate increases the burst temperature. However, at higher heating rates, burst strain is decreased considerably and an early rupture of the claddings without undergoing considerable ballooning is observed. It is also found that the degree of anisotropy has significant influence on the burst temperature and burst strain. With increasing degree of anisotropy, the burst temperature for claddings increases but there is a decrease in the burst strain. The effect of anisotropy in the ${\alpha}$-phase is carried over to ${\alpha}+{\beta}$-phase and its effect on the burst strain in the ${\alpha}+{\beta}$-phase too can be observed.

• Nuclear Engineering and Technology
• /
• 제51권1호
• /
• pp.293-302
• /
• 2019

Many researches have been done to develop and improve the performance of personal (individual) dosimeter response to cover a wide of neutron energy range (from thermal to fast). Depending on the individual category of the dosimeter, the semiconductor sensor has been used to simplify and lightweight. In this plan, it's very important to have a fairly accurate counting of doses rate in different energies. With a general design and single-sensor simulations, all optimal thicknesses have been extracted. The performance of the simulation scheme has been compared with the commercial and laboratory samples in the world. Due to the deviation of all dosimeters with a flat energy response, in this paper, has been used an idea of one semi-conductor sensor to have the flat energy-response in the entire neutron energy range. Finally, by analyzing of the sensors data as arrays for the first time, we have reached a nearly flat and acceptable energy-response. Also a comparison has been made between Lucite-PMMA ($H_5C_5O_2$) and polyethylene-PE ($CH_2$) as a radiator and $B_4C$ has been studied as absorbent. Moreover, in this paper, the effect of gamma dose in the dosimeter has been investigated and shown around the standard has not been exceeded.

• Nuclear Engineering and Technology
• /
• 제54권10호
• /
• pp.3803-3810
• /
• 2022

OpenMC is a community-driven open-source Monte Carlo neutron and photon transport simulation code. The Weight Window Mesh (WWM) function and an automatic Global Variance Reduction (GVR) method was recently developed and implemented in a developmental branch of OpenMC. This WWM function and GVR method broaden OpenMC's usage in general purposes deep penetration shielding calculations. However, the Local Variance Reduction (LVR) method, which suits the source-detector problem, is still missing in OpenMC. In this work, the Weight Window Generator (WWG) function has been developed and benchmarked for the same branch. This WWG function allows OpenMC to generate the WWM for the source-detector problem on its own. Single-material cases with varying shielding and sources were used to benchmark the WWG function and investigate how to set up the particle histories utilized in WWG-run and WWM-run. Results show that there is a maximum improvement of WWM generated by WWG. Based on the above results, instructions on determining the particle histories utilized in WWG-run and WWM-run for optimal computation efficiency are given and tested with a few multi-material cases. These benchmarks demonstrate the ability of the OpenMC WWG function and the above instructions for the source-detector problem. This developmental branch will be released and merged into the main distribution in the future.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1996년도 춘계학술발표회논문집(2)
• /
• pp.429-434
• /
• 1996

K-REM［1］, which is under development as a realistic evaluation model of large break LOCA, is applied to the analysis of cold leg guillotine break of Ulchin 3&4. Fuel parameters on which statistical analysis of their effects on the peak cladding temperature (PCT) are made and system parameters on which the concept of limiting value approach (LVA) are applied, are determined from the single parameter sensitivity study. 3 parameters of fuel gap conductance, fuel thermal conductivity and power peaking factor are selected as fuel related ones and 4 parameters of axial power shape, reactor power, decay heat and the gas pressure of safety injection tank (SIT) are selected as plant system related ones. Response surface of PCT is generated from the plant calculation results and on which Monte Carlo sampling is made to get plant application uncertainty which is statistically combined with code uncertainty to produce the 95th percentile PCT. From the break spectrum analysis, blowdown PCT of 1350.23 K and reflood PCT of 1195.56 K are obtained for break discharge coefficients of 0.8 and 0.5, respectively.

• Nuclear Engineering and Technology
• /
• 제55권6호
• /
• pp.1999-2010
• /
• 2023

Undesirable flapping motion of discs can cause the failure of swing check valves in nuclear passive safety systems. Time-resolved particle image velocimetry (PIV) was employed to investigate the flow characteristics around a free-to-rotate plate and the motion response, with the Reynolds numbers, based on the hydraulic diameter of the channel, from 1.32 × 10⁴ to 3.95 × 10⁴. Appreciable flapping motion (±3.52°) appeared at the Reynolds number of 2.6 × 10⁴ with the frequency of 5.08 Hz. In the low-Reynolds-number case, the plate showed negligible flapping. In the high-Reynolds-number case, the deflection angle increased with reduced flapping amplitude. The torque from the fluid determined the flapping amplitude. In the low-Reynolds-number case, Karman vortices were absent. With increasing Reynolds numbers, Karman vortices developed behind the plate with larger deflection angles. Strong interaction between the wake flow from the leading and trailing edge of the plate was observed. Based on power spectrum density (PSD) analysis, the vortex shedding frequency coincided with the flapping frequency, and the amplitude was positively correlated to the strength of the vortices. Proper orthogonal decomposition (POD) modes evince that, in the case of appreciable motion, coherent structures exhibited a larger spatial scale, enhancing the magnitude of the external torque on the plate.

• Nuclear Engineering and Technology
• /
• 제55권7호
• /
• pp.2454-2465
• /
• 2023

Film boiling may lead to burnout of the heating element. Even though burnout does not occur, the heating element is subject to deformation because it is not sufficiently strong to withstand external loads. In particular, the ballooning and rupture of a tube under film boiling are important phenomena in the field of nuclear reactor safety. If the tube-type cladding of nuclear fuel ruptures owing to high internal pressure and thermal load, radioactive materials inside the cladding are released to the coolant. Therefore, predicting the ballooning and rupture is important. This study presents numerical simulations to predict the ballooning behavior and rupture time of a horizontal tube at high internal pressure under saturated film boiling. To do so, a multi-step coupled simulation of conjugated film boiling heat transfer and ballooning using creep model is adopted. The numerical methods and models are validated against experimental values. Two different nonuniform heat flux distributions and four different internal pressures are considered. The three-step simulation is enough to obtain a convergent result. However, the single-step simulation also successfully predicts the rupture time. This is because the film boiling heat transfer characteristics are slightly affected by the tube geometry related to creep ballooning.

• 방사성폐기물학회지
• /
• 제21권3호
• /
• pp.397-410
• /
• 2023

This study examined the efficacy of various chlorinating agents in partitioning light water reactor spent fuel, with the aim of optimizing the chlorination process. Through thermodynamic equilibrium calculations, we assessed the outcomes of employing MgCl₂, NH₄Cl, and Cl₂ as chlorinating agents. A comparison was drawn between using a single agent and a sequential approach involving all three agents (MgCl₂, NH₄Cl, and Cl₂). Following heat treatment, the utilization of MgCl₂ as the sole chlorinating agent resulted in a moderate separation. Specifically, this method yielded a solid separation with 96.9% mass retention, 31.7% radioactivity, and 44.2% decay heat, relative to the initial spent fuel. In contrast, the sequential application of the chlorinating agents following heat treatment led to a final solid separation characterized by 93.1% mass retention, 5.1% radioactivity, and 15.4% decay heat, relative to the original spent fuel. The findings underscore the potential effectiveness of a sequential chlorination strategy for partitioning spent fuel. This approach holds promise as a standalone technique or as a complementary process alongside other partitioning processes such as pyroprocessing. Overall, our findings contribute to the advancement of spent fuel management strategies.