• Journal of Environmental Science International
• /
• v.25 no.1
• /
• pp.67-83
• /
• 2016

In order to improve the prediction of the regional air quality modeling in the Seoul metropolitan area, a sensitivity analysis using two PBL and microphysics (MP) options of the WRF model was performed during four seasons. The results from four sets of the simulation experiments (EXPs) showed that meteorological variables (especially wind field) were highly sensitive to the choice of PBL options (YSU or MYJ) and no significant differences were found depending on MP options (WDM6 or Morrison) regardless of specific time periods, i.e. day and night, during four seasons. Consequently, the EXPs being composed of YSU PBL option were identified to produce better results for meteorological elements (especially wind field) regardless of seasons. On the other hand, the accuracy of all simulations for summer and winter was somewhat lower than those for spring and autumn and the effect according to physics options was highly volatile by geographical characteristics of the observation site.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.185.2-185.2
• /
• 2014

We have studied the atomic and electronic structure of graphene nanoribbons (GNRs) on a hexagonal boron nitride (h-BN) sheet with intercalated atoms using first-principles calculations. The h-BN sheet is an insulator with the band gap about 6 eV and then it may a good candidate as a supporting dielectric substrate for graphene-based nanodevices. Especially, the h-BN sheet has the similar bond structure as graphene with a slightly longer lattice constant. For the computation, we use the Vienna ab initio simulation package (VASP). The generalized gradient approximation (GGA) in the form of the PBE-type parameterization is employed. The ions are described via the projector augmented wave potentials, and the cutoff energy for the plane-wave basis is set to 400 eV. To include weak van der Waals (vdW) interactions, we adopt the Grimme's DFT-D2 vdW correction based on a semi-empirical GGA-type theory. Our calculations reveal that the localized states appear at the zigzag edge of the GNR on the h-BN sheet due to the flat band of the zigzag edge at the Fermi level and the localized states rapidly decay into the bulk. The open-edged graphene with a large corrugation allows some space between graphene and h-BN sheet. Therefore, atoms or molecules can be intercalated between them. We have considered various types of atoms for intercalation. The atoms are initially placed at the edge of the GNR or inserted in between GNR and h-BN sheet to find the effect of intercalated atoms on the atomic and electronic structure of graphene. We find that the impurity atoms at the edge of GNR are more stable than in between GNR and h-BN sheet for all cases considered. The nickel atom has the lowest energy difference of ~0.2 eV, which means that it is relatively easy to intercalate the Ni atom in this structure. Finally, the magnetic properties of intercalated atoms between GNR and h-BN sheet are investigated.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.56 no.2
• /
• pp.83-93
• /
• 2020

In fisheries, the importance of designing efficient fish cages is being emphasized as aquaculture has become more production than capture fishing. Particularly, the gravity cage system is one of the popular fish cage system in Korea. Currently, gravity cages of various shapes and sizes are being widely designed and installed in offshore and inland seas. The cage is subject to external forces, such as currents and waves, and the shape of the structure and tension on the ropes changes according to these external forces. Thus, it is important to accurately calculate these dynamic behavior, including the external forces and tension on the structure during the design stage. In this study, three types of cage systems with an equal internal volume of 8000 ㎥ were analyzed using mass-spring models and their behavior was interpreted through simulations. These simulations were used to analyze the behavior and tension of the ropes in response to currents and waves to aid in the selection of individual cage sizes for a given total volume. The numerical calculation results indicate that depending on the flow rate, the most resistant system is System 1, which has eight strays, and System 2 and System 3 have 69.4% and 54.8% of the resistance of System 1. Further, total resistance increased as the number of cages increased for all flow rates.

• Journal of Korea Game Society
• /
• v.18 no.1
• /
• pp.125-134
• /
• 2018

This paper presents a practical method to efficiently extract the rotational part of a $3{\times}3$ matrix that changes continuously in time. This is the key technique in the corotational FEM and the shape matching deformation popular in physics-based dynamic deformation. Recently, in contrast to the traditional polar decomposition methods independent of time, an iterative method was proposed that formulates the rotation extraction in a physics-based way and exploits an incremental representation of rotation. We develop an optimization method that reduces the number of iterations under the assumption that the maximum magnitude of the incremental rotation vector is limited within ${\pi}/2$. Realistic simulation of dynamic deformation employs a sufficiently small time step, and thus this assumption is not problematic in practice. We demonstrate the efficiency and practicality of our method in various experiments.

• Wind and Structures
• /
• v.30 no.4
• /
• pp.393-403
• /
• 2020

One of the growing concerns of the wind energy production is wind ramp events. To improve the wind ramp event forecasts, the nonlinear Kalman filter bias correction method was applied to 24-h wind speed forecasts issued from the WRF model at 70-m height in Zhangbei wind farm, Hebei Province, China for a two-year period. The Kalman filter shows the remarkable ability of improving forecast skill for real-time wind speed forecasts by decreasing RMSE by 32% from 3.26 m s_-1 to 2.21 m s_-1, reducing BIAS almost to zero, and improving correlation from 0.58 to 0.82. The bias correction improves the forecast skill especially in wind speed intervals sensitive to wind power prediction. The fact shows that the Kalman filter is especially suitable for wind power prediction. Moreover, the bias correction method performs well under abrupt weather transition. As to the overall performance for improving the forecast skill of ramp events, the Kalman filter shows noticeable improvements based on POD and TSS. The bias correction increases the POD score of up-ramps from 0.27 to 0.39 and from 0.26 to 0.38 for down-ramps. After bias correction, the TSS score is significantly promoted from 0.12 to 0.26 for up-ramps and from 0.13 to 0.25 for down-ramps.

• Journal of the Korean Magnetics Society
• /
• v.15 no.6
• /
• pp.307-310
• /
• 2005

To get the spatial feature of arterial pulse, we designed spatial pulse diagnostic apparatus (SPDA) using a 2-dimensional magnetoresistive sensor array. The magnetic field distribution fur magnet may was simulated using finite element method. We recognized that the field distribution of parallel magnet mays was more sensitive and uniformed than that of perpendicular one. Also the spatial displacements of magnet array were agreed with the output signal of magnetic tunnel junction (MTJ) sensor array.

• Geomechanics and Engineering
• /
• v.30 no.2
• /
• pp.187-199
• /
• 2022

Kriging metamodel, as a flexible machine learning method for approximating deterministic analysis models of an engineering system, has been widely used for efficiently estimating slope reliability in recent years. However, the autocorrelation function (ACF), a key input to Kriging that affects the accuracy of reliability estimation, is usually selected based on empiricism. This paper proposes an adaption of the Kriging method, named as Genetic Algorithm optimized Whittle-Matérn Kriging (GAWMK), for addressing this issue. The non-classical two-parameter Whittle-Matérn (WM) function, which can represent different ACFs in the Matérn family by controlling a smoothness parameter, is adopted in GAWMK to avoid subjectively selecting ACFs. The genetic algorithm is used to optimize the WM model to adaptively select the optimal autocorrelation structure of the GAWMK model. Monte Carlo simulation is then performed based on GAWMK for a subsequent slope reliability analysis. Applications to one explicit analytical example and two slope examples are presented to illustrate and validate the proposed method. It is found that reliability results estimated by the Kriging models using randomly chosen ACFs might be biased. The proposed method performs reasonably well in slope reliability estimation.

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3317-3323
• /
• 2022

Half Value Layer calculations theoretically need prior specification of linear attenuation calculations, since the HVL value is derived by dividing ln(2) by the linear attenuation coefficient. The purpose of this study was to establish a direct computational model for determining HVL, a vital parameter in nuclear radiation safety studies and shielding material design. Accordingly, a typical gamma-ray transmission setup has been modeled using MCNPX (version 2.4.0) general-purpose Monte Carlo code. The MCNPX code's INPUT file was designed with two detection locations for primary and secondary gamma-rays, as well as attenuator material between those detectors. Next, Half Value Layer values of some well-known gamma-ray shielding materials such as lead and ordinary concrete have been calculated throughout a broad gamma-ray energy range. The outcomes were then compared to data from the National Institute of Standards and Technology. The Half Value Layer values obtained from MCNPX were reported to be highly compatible with the HVL values obtained from the NIST standard database. Our results indicate that the developed INPUT file may be utilized for direct computations of Half Value Layer values for nuclear safety assessments as well as medical radiation applications. In conclusion, advanced simulation methods such as the Monte Carlo code are very powerful and useful instruments that should be considered for daily radiation safety measures. The modeled MCNPX input file will be provided to the scientific community upon reasonable request.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3450-3463
• /
• 2023

In the two-step analysis of Pebble-Bed type High-Temperature Gas-Cooled Reactor (PB-HTGR), the lattice physics calculation for the generation of homogenized cross-sections is based on the fuel pebble. However, the randomly-dispersed fuel particles in the fuel pebble introduce double heterogeneity and randomness. Compared to the deterministic method, the Monte Carlo method which is flexible in geometry modeling provides a high-fidelity treatment. Therefore, the Monte Carlo code NECP-MCX is extended in this study to perform the lattice physics calculation of the PB-HTGR. Firstly, the capability for the simulation of randomly-dispersed media, using the explicit modeling approach, is developed in NECP-MCX. Secondly, the capability for the generation of the homogenized cross-section is also developed in NECP-MCX. Finally, simplified PB-HTGR problems are calculated by a two-step neutronics analysis tool based on Monte Carlo homogenization. For the pebble beds mixed by fuel pebble and graphite pebble, the bias is less than 100 pcm when compared to the high-fidelity model, and the bias is increased to 269 pcm for pebble bed mixed by depleted fuel pebble. Numerical results show that the Monte Carlo lattice physics calculation for the two-step analysis of PB-HTGR is feasible.

• Journal of Radiation Protection and Research
• /
• v.48 no.4
• /
• pp.204-212
• /
• 2023

Background: Isotopes of the projectile may be produced along the beam path during the irradiation of a target by a heavy ion due to inelastic interactions with the media. This study analyzed the production cross-section of carbon (C) and Helium (He) projectile's isotopes resulting from the interactions of these beams with different materials along the beam path. Materials and Methods: In this study, we transport C and He ion beams through different materials. This transportation was made by the Monte Carlo simulation. Particle and Heavy Ion Transport code System (PHITS) has been used for this calculation. Results and Discussion: It has been found that ¹⁰C, ¹¹C, and ¹³C from the ¹²C ion beam and ³He from the ⁴He ion beam are significant projectile's isotopes that have higher flux than other isotopes of these projectiles. The ⁴He ion beam has a higher projectile's isotope production cross-section along the beam path, which adds more impurities to the beam than the ¹²C ion beam. These projectile's isotopes from both the ¹²C and ⁴He ion beams have higher production cross-sections in hydrogenous materials like water or polyethylene. Conclusion: It is important to distinguish these projectile's isotopes from the primary beam particles to obtain a precise and accurate cross-section result by minimizing the error during measurement with a nuclear track detector. This study will show the trend of the production probability of projectile's isotopes for these ion beams.