• Nuclear Engineering and Technology
• /
• v.53 no.8
• /
• pp.2767-2773
• /
• 2021

We used the GEANT4 Monte Carlo MC Toolkit to simulate carbon ion beams incident on water, tissue, and bone, taking into account nuclear fragmentation reactions. Upon increasing the energy of the primary beam, the position of the Bragg-Peak transfers to a location deeper inside the phantom. For different materials, the peak is located at a shallower depth along the beam direction and becomes sharper with increasing electron density NZ. Subsequently, the generated depth dose of the Bragg curve is then benchmarked with experimental data from GSI in Germany. The results exhibit a reasonable correlation with GSI experimental data with an accuracy of between 0.02 and 0.08 cm, thus establishing the basis to adopt MC in heavy-ion treatment planning. The Kolmogorov-Smirnov K-S test further ascertained from a statistical point of view that the simulation data matched the experimentally measured data very well. The two-dimensional isodose contours at the entrance were compared to those around the peak position and in the tail region beyond the peak, showing that bone produces more dose, in comparison to both water and tissue, due to secondary doses. In the water, the results show that the maximum energy deposited per fragment is mainly attributed to secondary carbon ions, followed by secondary boron and beryllium. Furthermore, the number of protons produced is the highest, thus making the maximum contribution to the total dose deposition in the tail region. Finally, the associated spectra of neutrons and photons were analyzed. The mean neutron energy value was found to be 16.29 MeV, and 1.03 MeV for the secondary gamma. However, the neutron dose was found to be negligible as compared to the total dose due to their longer range.

• Nuclear Engineering and Technology
• /
• v.53 no.6
• /
• pp.1769-1785
• /
• 2021

The Integral Pressurized Water Reactors (IPWRs) as the innovative advanced and generation-III + reactors are under study and developments in a lot of countries. This paper is aimed at the thermal hydraulic study of the hot and average fuel sub-channel in a Generation III + IPWR by loose external coupling to the neutronic simulation. The power produced in fuel pins is calculated by the neutronic simulation via MCNPX2.6 then fuel and coolant temperature changes along fuel sub-channels evaluated by computational fluid dynamic thermal hydraulic calculation through an iterative coupling. The relative power densities along the fuel pin in hot and average fuel sub-channel are calculated in sixteen equal divisions. The highest centerline temperature of the hottest and the average fuel pin are calculated as 633 K (359.85 ℃) and 596 K (322.85 ℃), respectively. The coolant enters the sub-channel with a temperature of 557.15 K (284 ℃) and leaves the hot sub-channel and the average sub-channel with a temperature of 596 K (322.85 ℃) and 579 K (305.85 ℃), respectively. It is shown that the spacer grids result in the enhancement of turbulence kinetic energy, convection heat transfer coefficient along the fuel sub-channels so that there is an increase in heat transfer coefficient about 40%. The local fuel pin temperature reduction in the place and downstream the space grids due to heat transfer coefficient enhancement is depicted via a graph through six iterations of neutronic and thermal hydraulic coupling calculations. Working in a low fuel temperature and keeping a significant gap below the melting point of fuel, make the IPWR as a safe type of generation -III + nuclear reactor.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1540-1553
• /
• 2019

A nuclear power plant (NPP) is a highly complex system-of-systems as manifested through its internal systems interdependence. The negative impact of such interdependence was demonstrated through the 2011 Fukushima Daiichi nuclear disaster. As such, there is a critical need for new strategies to overcome the limitations of current risk assessment techniques (e.g. the use of static event and fault tree schemes), particularly through simulation of the nonlinear dynamic feedback mechanisms between the different NPP systems/components. As the first and key step towards developing an integrated NPP dynamic probabilistic risk assessment platform that can account for such feedback mechanisms, the current study adopts a system dynamics simulation approach to model the thermal dynamic processes in: the reactor core; the secondary coolant system; and the pressurized water reactor. The reactor core and secondary coolant system parameters used to develop system dynamics models are based on those of the Palo Verde Nuclear Generating Station. These three system dynamics models are subsequently validated, using results from published work, under different system perturbations including the change in reactivity, the steam valve coefficient, the primary coolant flow, and others. Moving forward, the developed system dynamics models can be integrated with other interacting processes within a NPP to form the basis of a dynamic system-level (systemic) risk assessment tool.

• Atmosphere
• /
• v.31 no.5
• /
• pp.563-576
• /
• 2021

The performance of East Asian summer monsoon is assessed for GC2 and GC3.1, which are climate change models of the current and next climate prediction system in the Korea Meteorological Administration (KMA), GloSea5 and GloSea6. The most pronounced characteristics of GC models are strong monsoon trough and the weakening of the Western North Pacific Subtropical High (WNPSH). These are related to the weakening of the southwesterly wind and resulting weak monsoon band toward the Korean Peninsula. The GC3.1 is known to have improved the model configuration version compared to GC2, such as cloud physics and ocean parameters. We can confirm that the overall improvements of GC3.1 against GC2, especially in pressure, 850 hPa wind fields, and vertical wind shear. Also, the precipitation band stagnant in the south of 30°N in late spring is improved, therefore the biases of rainy onset and withdrawal on the Korean Peninsula are reduced by 2~4 pentad. We also investigate the impact of initialization in comparison with GloSea5 hindcast. Compared with GCs, hindcast results show better simulation within 1 month lead time, especially in pressure and 850 hPa wind fields, which can be expected to the improvement of WNPSH. Therefore, it is expected that the simulation performance of WNPSH will be improved in the result of applying the initialization of GloSea6.

• Nuclear Engineering and Technology
• /
• v.55 no.10
• /
• pp.3648-3658
• /
• 2023

The early localization of a fuel subassembly with a failed (wet rupture) fuel pin is very important in reactors to limit the associated radiological and operational consequences. This requires a fast and reliable system for failure detection and their localization in the core. In the Prototype Fast Breeder Reactor, the system specially designed for this purpose is Failed Fuel Location Modules (FFLM) housed in the control plug region. It identifies a failed sub-assembly by detecting the presence of delayed neutrons in the sodium from a failed sub-assembly. During the commissioning phase of PFBR, it is mandatory to demonstrate the FFLM effectiveness. The paper highlights the engineering and physics design aspects of FFLM and the integrated simulation towards its function demonstration with a source assembly containing a perforated metallic fuel pin. This test pin mimics a MOX pin of 1 cm² of geometrical defect area. At 10% power and 20% sodium flow rate, the counts rate in the BCCs of FFLM system range from 75 cps to 145 cps depending upon the position of DN source assembly. The model developed for the counts simulation is applicable to both metal and MOX pins with proper values of k-factor and escape coefficient.

• Progress in Medical Physics
• /
• v.20 no.3
• /
• pp.139-144
• /
• 2009

Recently, the mathematical analysis of stent simulation has been improved, with the help of development of various tool which measure mechanical property and location of stent in artery. The most crucial part of the stent modeling is how to design ideal stent and to evaluate the interaction between stent and artery. While there has been great deal of researches on the evaluation of the expansion, stress distribution, deformation of the stent in terms of the various parameters, few verification through computer simulation has been performed about deformation and stress distribution of the stent. In this study, we have produced the corresponding results between experimental test using Universal Testing Machine and computer simulation for the ideal model of stent. Also, we have analyzed and compared stress distribution of stent in the cases of that with membrane and that without membrane. The results of this study would provide minimum change of plan and good quality for ideal stent replacing damaged artery through the analysis using computer simulation in the early stage of stent design.

• Radiation Oncology Journal
• /
• v.20 no.2
• /
• pp.165-171
• /
• 2002

Purpose : In order to perform craniospinal irradiation (CSI) in the supine position on patients who are unable to lie in the prone position, a new simulation technique using a CT simulator was developed and its availability was evaluated. Materials and Method : A CT simulator and a 3-D conformal treatment planning system were used to develop CSI in the supine position. The head and neck were immobilized with a thermoplastic mask in the supine position and the entire body was immobilized with a Vac-Loc. A volumetrie image was then obtained using the CT simulator. In order to improve the reproducibility of the patients' setup, datum lines and points were marked on the head and the body. Virtual fluoroscopy was peformed with the removal of visual obstacles such as the treatment table or the immobilization devices. After the virtual simulation, the treatment isocenters of each field were marked on the body and the immobilization devices at the conventional simulation room. Each treatment field was confirmed by comparing the fluoroscopy images with the digitally reconstructed radiography (DRR)/digitally composite radiography (DCR) images from the virtual simulation. The port verification films from the first treatment were also compared with the DRR/DCR images for a geometrical verification. Results : CSI in the supine position was successfully peformed in 9 patients. It required less than 20 minutes to construct the immobilization device and to obtain the whole body volumetric images. This made it possible to not only reduce the patients' inconvenience, but also to eliminate the position change variables during the long conventional simulation process. In addition, by obtaining the CT volumetric image, critical organs, such as the eyeballs and spinal cord, were better defined, and the accuracy of the port designs and shielding was improved. The differences between the DRRs and the portal films were less than 3 mm in the vertebral contour. Conclusion : CSI in the supine position is feasible in patients who cannot lie on prone position, such as pediatric patienta under the age of 4 years, patients with a poor general condition, or patients with a tracheostomy.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.4
• /
• pp.308-318
• /
• 2002

The requirements to the numerical model of wind-generated waves in shallow water are discussed in the framework of the Korteweg-de Vries equation. The weakness of nonlinearity and dispersion required for the Korteweg-de Vries equation applicability is considered for fully developed sea, non-stationary wind waves and swell, including some experimental data. We note for sufficient evaluation of the freak wave statistics it is necessary to consider more than about 10,000 waves in the wave record, and this leads to the limitation of the numerical domain and number of realizations. The numerical modelling of irregular water waves is made to demonstrate the possibility of effective evaluation of the statistical properties of freak waves with heights equal to 2-2.3 significant wave height.

• Atmosphere
• /
• v.25 no.1
• /
• pp.1-18
• /
• 2015

The most objective way to overcome the limitation of numerical weather prediction model is to represent the uncertainty of prediction by introducing probabilistic forecast. The uncertainty of the numerical weather prediction system developed due to the parameterization of unresolved scale motions and the energy losses from the sub-scale physical processes. In this study, we focused on the growth of model errors. We performed ensemble forecast to represent model uncertainty. By employing the multi-physics scheme (PHYS) and the stochastic kinetic energy backscatter scheme (SKEBS) in simulating typhoon Rusa (2002), we assessed the performance level of the two schemes. The both schemes produced better results than the control run did in the ensemble mean forecast of the track. The results using PHYS improved by 28% and those based on SKEBS did by 7%. Both of the ensemble mean errors of the both schemes increased rapidly at the forecast time 84 hrs. The both ensemble spreads increased gradually during integration. The results based on SKEBS represented model errors very well during the forecast time of 96 hrs. After the period, it produced an under-dispersive pattern. The simulation based on PHYS overestimated the ensemble mean error during integration and represented the real situation well at the forecast time of 120 hrs. The displacement speed of the typhoon based on PHYS was closest to the best track, especially after landfall. In the sensitivity tests of the model uncertainty of SKEBS, ensemble mean forecast was sensitive to the physics parameterization. By adjusting the forcing parameter of SKEBS, the default experiment improved in the ensemble spread, ensemble mean errors, and moving speed.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2004.11a
• /
• pp.115-118
• /
• 2004

Respiratory motion in the thorax and abdomen is an important limiting factor in high-precision radiation therapy. The lung tumor and tumor(pancreas, stomach) in abdomen therefore are internal motion due to breathing. We will perform to measurement of dose distributions for these moving tumors. In preliminary study, we investigated displacement of moving tumor such as liver, lung tumor in abdomen with previously reported papers. With reference data, internal movements of tumor are displayed with phantom and moving control device(MCD), which appear three dimension (3-D) motion such as x, y and z axis. These devices are used to access dose delivered in tumor with and without internal motion. The MCD and phantom were used to evaluate a delivered dose under similar condition, although there are not same internal tumor motion. In future, we will obtain the exact evaluation of dose if improved in programed software of moving control device and measure precise internal motion using image modality such as fluoroscopy, simulator in based on this study.