• Nuclear Engineering and Technology
• /
• v.56 no.3
• /
• pp.772-784
• /
• 2024

The hexagonal nodal code RENUS has been enhanced to handle irregularly deformed hexagonal assemblies. The underlying RENUS methods involving triangle-based polynomial expansion nodal (T-PEN) and corner point balance (CPB) were extended in a way to use line and surface integrals of polynomials in a deformed hexagonal geometry. The nodal calculation is accelerated by the coarse mesh finite difference (CMFD) formulation extended to unstructured geometry. The accuracy of the unstructured nodal solution was evaluated for a group of 2D SFR core problems in which the assembly corner points are arbitrarily displaced. The RENUS results for the change in nuclear characteristics resulting from fuel deformation were compared with those of the reference McCARD Monte Carlo code. It turned out that the two solutions agree within 18 pcm in reactivity change and 0.46% in assembly power distribution change. These results demonstrate that the proposed unstructured nodal method can accurately model heterogeneous thermal expansion in hexagonal fueled cores.

• Journal of Environmental Science International
• /
• v.17 no.6
• /
• pp.623-631
• /
• 2008

Many recent studies have concentrated upon the radiative effects of atmospheric aerosols. Though their scattering and absorption of radiation, aerosols can also induce some other important environment effects. In this study, new radiation code and aerosol data within Atmosphere General Circulation Model (AGCM) is used to assess the aerosol radiative forcing and to analyze relative climate effects. The new Kangnung National University AGCM Stratospheric-15 (KNU AGCM ST15) was integrated by using two sets of radiative effect of aerosols: CTRL as not a radiative effect of aerosols and AERO as a radiative effect of aerosols. Two cases show the difference of net shortwave radiation budget at top-of-atmosphere (TOA) is found to be about $-3.4Wm^{-2}$, at the surface (SFC) is about $-5.6Wm^{-2}$. Consequently the mean atmospheric absorption due to aerosol layer in global is about $2.2Wm^{-2}$. This result confirms the existence of a negative forcing due to the direct effect of aerosols at the surface and TOA in global annual mean. In addition, it is found that cooling over at the surface air temperature due to radiative effect of aerosols is about $0.17^{\circ}C$. It is estimated that radiative forcing of the net upward longwave radiation taken as the indirect effect of aerosol is much smaller than that of the direct effect as there is about $0.2Wm^{-2}$ of positive forcing both at TOA and at SFC. From this study, It made an accurate estimation of considering effect of aerosols that is negative effect. This may slow the rate of projected global warming during the $21^{st}$ century.

• Tunnel and Underground Space
• /
• v.13 no.1
• /
• pp.64-75
• /
• 2003

This study presents a newly suggested test method of Mode II fracture toughness measurement called "Punch Through Shear Test" which was originally proposed by Backers and Stephansson in 2001. The purpose of this study is to check the validity of the suggested testing method by performing Mode II fracture toughness tests for Daejeon Granite. In addition, the optimal specimen geometry for the testing and the relation between Mode II fracture toughness and confining pressure were also investigated. Fractured surface was observed to be very smooth with lots of rock debris which came off fracture surface which obviously implies that the surface was sheared off. This confirms that Mode II fracturing actually occurred. In addition, numerical analyses including continuum analysis, particle flow code analysis and crack propagation simulations were performed. Results of these numerical analyses indicated that the cracks occurred in the specimen were predominantly in Mode II and these cracks led to failure of the test specimen. From this investigation, it can be concluded that the newly suggested "Punch Through Shear Test" method provides a reliable means of determining the Mode II fracture toughness. fracture toughness.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.10
• /
• pp.779-787
• /
• 2013

The effects of various nozzle configurations on infrared signature are investigated for the purpose of analysing the infrared signature level of aircraft propulsion system. A virtual subsonic aircraft is selected and then a circular convergent nozzle, which meets the mission requirements, is designed. Convergent nozzles of different configurations are designed with different geometric profiles. Using a compressible Navier-Stokes-Fourier CFD code, an analysis of thermal flow field and nozzle surface temperature distribution is conducted. From the information of plume flow field and nozzle surface temperature distribution, IR signature of plume and nozzle surface is calculated through the narrow-band model and the RadThermIR code. Finally, qualitative information for IR signature reduction is obtained through the analysis of the effects of various nozzle configurations on IR signature.

• Journal of the Korean Wood Science and Technology
• /
• v.41 no.6
• /
• pp.497-506
• /
• 2013

The optimal condition for the torrefaction of eucalyptus (Eucalyptus globulus) was investigated by response surface methodology. The carbon content in the torrefied biomass increased with the severity factor (SF), while hydrogen and oxygen contents decreased. The calorific value of torrefied biomass ranged from 20.23 to 21.29 MJ/kg, depending on the torrefaction conditions. This implied that the energy contained in the torrefied biomass increased by 1.6 to 6.9%, when compared with that of the untreated biomass. The weight loss of biomass increased as the SF increased. The Code level of reaction temperature had the highest impact on the energy yield of torrefied biomass, while the effect of Code level of reaction time was considerably low. The highest energy yield was obtained at low SF.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• Journal of Radiation Protection and Research
• /
• v.48 no.3
• /
• pp.131-143
• /
• 2023

Background: This study examined the detection limit of thyroid screening monitoring conducted at the time of the Fukushima Daiichi Nuclear Power Plant (FDNPP) accident in 2011 using a Monte Carlo simulation. Materials and Methods: We calculated the detection limit of a NaI(Tl) survey meter to measure ¹³¹I accumulation in the thyroid gland of children. Mathematical phantoms of 1- and 5-year-old children were developed in the simulation of the Particle and Heavy Ion Transport code System code. Contamination of the body surface with eight radionuclides found after the FDNPP accident was assumed to have been deposited on the neck and shoulder area. Results and Discussion: The detection limit was calculated as a function of ambient dose rate. In the case of 40 Bq/cm² contamination on the body surface of the neck, the present simulations showed that residual thyroid radioactivity corresponding to thyroid dose of 100 mSv can be detected within 21 days after intake at the ambient dose rate of 0.2 µSv/hr and within 11 days in the case of 2.0 µSv/hr. When a time constant of 10 seconds was used at the dose rate of 0.2 µSv/hr, the estimated survey meter output error was 5%. Evaluation of the effect of individual differences in the location of the thyroid gland confirmed that the measured value would decrease by approximately 6% for a height difference of ±1 cm and increase by approximately 65% for a depth of 1 cm. Conclusion: In the event of a nuclear disaster, simple measurements carried out using a NaI(Tl) scintillation survey meter remain effective for assessing ¹³¹I intake. However, it should be noted that the presence of short-half-life radioactive materials on the body surface affects the detection limit.

• Journal of Navigation and Port Research
• /
• v.32 no.1
• /
• pp.9-14
• /
• 2008

VOF method is known as one of the most effective numerical techniques handling two-fluid domains of different density simultaneously. Present study deals with the numerical analysis of flow field around bodies steadily moving near free-surface using FLUENT-VOF method. Validations were made by applying to three typical examples ; 2-D submerged hydrofoil, 3-D surface piercing body and container ship. It was found that the commercial software, FLUENT, is useful in practical use, and VOF method is capable of handling free-surface around moving bodies although discussions are limited to the analysis in qualitative sense.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.683-688
• /
• 2007

Surface heat treatment of SM45C round bar by diode laser was simulated to find it's condition by using commercial finite element code MARC. Due to axisymmetric geometry, a quarter of model for SM45C round bar was considered and user subroutines were applied to boundary condition for the heat transfer. Material properties such as conductivity, specific heat and mass density were given as a function of temperature. Rotation speed of round bar and feed rate of beam were considered to design heat source model. Shape parameter values of heat source were determined by beam profile. As results, Three dimensional heat source model for diode laser beam conditions of surface hardening has been designed by the comparison between the finite element analysis results and experimental data on SM45C round bar. Diode laser surface hardening for SM45C round bar was successfully simulated and it should be useful to determine optimal heat treatment condition.