• Progress in Medical Physics
• /
• v.27 no.2
• /
• pp.86-92
• /
• 2016

An absorbed dose calculation method using a digital phantom is implemented in normal organs. This method cannot be employed for calculating the absorbed dose of tumor. In this study, we measure the S-value for calculating the absorbed dose of each organ and tumor. We inject a radioisotope into a torso phantom and perform Monte Carlo simulation based on the CT data. The torso phantom has lung, liver, spinal, cylinder, and tumor simulated using a spherical phantom. The radioactivity of the actual absorbed dose is measured using the injected dose of the radioisotope, which is Cu-64 73.85 MBq, and detected using a glass dosimeter in the torso phantom. To perform the Monte Carlo simulation, the information on each organ and tumor acquired using the PET/CT and CT data provides anatomical information. The anatomical information is offered above mean value and manually segmented for each organ and tumor. The residence time of the radioisotope in each organ and tumor is calculated using the time activity curve of Cu-64 radioactivity. The S-values of each organ and tumor are calculated based on the Monte Carlo simulation data using the spatial coordinate, voxel size, and density information. The absorbed dose is evaluated using that obtained through the Monte Carlo simulation and the S-value and the residence time in each organ and tumor. The absorbed dose in liver, tumor1, and tumor2 is 4.52E-02, 4.61E-02, and 5.98E-02 mGy/MBq, respectively. The difference in the absorbed dose measured using the glass dosimeter and that obtained through the Monte Carlo simulation data is within 12.3%. The result of this study is that the absorbed dose obtained using an image can evaluate each difference region and size of a region of interest.

• Journal of the Korea Society for Simulation
• /
• v.28 no.2
• /
• pp.1-14
• /
• 2019

A full scale hydrodynamic simulation that requires an accurate reproduction of shock-induced detonation was conducted for design of an energetic component system. A detailed hydrodynamic analysis SW was developed to validate the reactive flow model for predicting the shock propagation in a train configuration and to quantify the shock sensitivity of the energetic materials. The pyrotechnic device is composed of four main components, namely a donor unit (HNS+HMX), a bulkhead (STS), an acceptor explosive (RDX), and a propellant (BPN) for gas generation. The pressurized gases generated from the burning propellant were purged into a 10 cc release chamber for study of the inherent oscillatory flow induced by the interferences between shock and rarefaction waves. The pressure fluctuations measured from experiment and calculation were investigated to further validate the peculiar peak at specific characteristic frequency (${\omega}_c=8.3kHz$). In this paper, a step-by-step numerical description of detonation of high explosive components, deflagration of propellant component, and deformation of metal component is given in order to facilitate the proper implementation of the outlined formulation into a shock physics code for a full scale hydrodynamic simulation of the energetic component system.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.16 no.1
• /
• pp.176-182
• /
• 2012

Recently, a new theoretical and methodological approach for computational science is becoming more and more popular for analyzing and solving scientific problems in various scientific disciplines and applied research. Computational science is a field of study concerned with constructing mathematical models and quantitative analysis techniques and using large computing resources to solve the problems which are difficult to approach in a physical experimentally. In this paper, we present R&D of EDISON open integration platform that allows anyone like professors, researchers, industrial workers, students etc to upload their advanced research result such as simulation SW to use and share based on the cyber infrastructure of supercomputer and network. EDISON platform, which consists of 3 tiers (EDISON application framework, EDISON middleware, and EDISON infra resources) provides Web portal for education and research in 5 areas (CFD, Chemistry, Physics, Structural Dynamics, Computational Design) and user service.

• Bulletin of the Korean Chemical Society
• /
• v.29 no.8
• /
• pp.1554-1560
• /
• 2008

The probe diffusion and friction constants of methyl yellow (MY) in liquid n-alkanes of increasing chain length were calculated by equilibrium molecular dynamics (MD) simulations at temperatures of 318, 418, 518 and 618 K. Lennard-Jones particles with masses of 225 and 114 g/mol are modeled for MY. We observed that the diffusion constant of the probe molecule follows a power law dependence on the molecular weight of nalkanes, DMY${\sim}M^{-\gamma}$ well. As the molecular weight of n-alkanes increases, the exponent $\gamma$ shows sharp transitions near n-dotriacontane ($C_{32}$) for the large probe molecule (MY2) at low temperatures of 318 and 418 K. For the small probe molecule (MY1) $D_{MY1}$ in $C_{12}$ to C80 at all the temperatures are always larger than Dself of n-alkanes and longer chain n-alkanes offer a reduced friction relative to the shorter chain n-alkanes, but this reduction in the microscopic friction for MY1 is not large enough to cause a transition in the power law exponent in the log-log plot of DMY1 vs M of n-alkane. For the large probe molecule (MY2) at high temperatures, the situation is very similar to that for MY1. At low temperatures and at low molecular weights of n-alkanes, $D_{MY2}$ are smaller than $D_{self}$ of n-alkanes due to the relatively large molecular size of MY2, and MY2 experiences the full shear viscosity of the medium. As the molecular weight of n-alkane increases, $D_{self}$ of n-alkanes decreases much faster than $D_{MY2}$ and at the higher molecular weights of n-alkane, MY2 diffuses faster than the solvent fluctuations. Therefore there is a large reduction of friction in longer chains compared to the shorter chains, which enhances the diffusion of MY2. The calculated friction constants of MY1 and MY2 in liquid n-alkanes supported these observations. We deem that this is the origin of the so-called“solventoligomer”transition.

• Applied Microscopy
• /
• v.40 no.1
• /
• pp.21-28
• /
• 2010

In a simulation study of atomic resolution transmission electron microscope images of single layer and bilayer graphene, it is demonstrated that the conventional Bloch wave formulations can be used when high-order Laue zone reflections are properly taken into account in the theory. The simulated images for bilayer graphene show 3-fold rotational lattice symmetry rather than the 6-fold one under certain conditions. This result can be understood as revealed the 3-fold rotational lattice symmetry of bilayer graphene in three dimensions along [0001]. For single layer graphene the observed phase images showing 3-fold rotational lattice symmetry were particularly noted. This phenomenon has been explained by an assumption of the re-configuration of electron density on the surface of graphene. And the matching images have been obtained as simulated with up to the second order Laue zone reflections only, reflecting the re-configuration of electrons on the surface.

• Journal of Korean Society of Transportation
• /
• v.20 no.3
• /
• pp.105-113
• /
• 2002

As the number of pedestrian accident increases, the reconstruction of an accident becomes important to find the source of the fault. Generally, accidents are reconstructed by the intuition of experts or primitive physics. A reconstruction method is proposed using sophisticated optimization technology. At first, a dynamic simulation model is established for the accident environment. Occupant analysis for automobile crashworthiness is employed. The situation before an accident is identified by optimization. The impact velocity and the position of the pedestrian are utilized as design variables. The design variables are found by minimizing the difference between the simulation and the real accident. The optimization process is performed by linking an occupant analysis program MADYMO to an optimization program VisualDOC. Since the involved analysis is dynamics and highly nonlinear, response surface method is selected for the optimization process. Problems are solved for various situations.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2013.02a
• /
• pp.537-537
• /
• 2013

Cyclotron-accelerated ion beams are used for various researches, such as nuclear physics, nuclear chemistry, biotechnology, and material sciences including radio-isotope production. Recently considerable applications are asked to the cyclotron development undertaken to meet user requirements of various ions'energies, intensities, and their pulsed beams. For instance, a cocktail beam acceleration technique rapidly changing the ion species and energies was developed to irradiating integrated circuit chips. Also a chopping system in a cyclotron injection line is considered for producing a pulsed ion beam with a relatively long period compared with that generated by the resonance frequency. For the research in neutron time-of-flight measurement, a single-pulsed beam with a repetition interval of the order of mili-seconds or longer is necessary to have a good resolution and to remove background events. In this paper a feasibility of pulsed beam with an external ion source is simulated by adopting a combination system of a chopper accompanying with a bunching stage in the injection line and an additional chopper after the exit of the cyclotron in order to produce beam pulses with a range of $1{\mu}s{\sim}1ms$ periods from a resonance RF cycle. The pulseperiod will be adjusted by chopping the number of beam bunches from the injected pulses in the injection line. However, the longer pulses will have reduced number of beam pulses and sacrificed beam currents. Because the beam users need an intense single pulsed beam, a careful tuning of the acceleration phase and a high-intense external ion source are necessary to achieve an intense single-pulsed beam from the cyclotron. It is essential to strictly match the acceleration phase of injected beams in the central region of the cyclotron to improve its efficiency. An effect of space charge at each pulse from the ion source will be also considered.

• Journal of Magnetics
• /
• v.19 no.4
• /
• pp.399-403
• /
• 2014

We designed an ultra-high frequency (UHF: 300MHz to 3 GHz) ferrite circulator to investigate magnetic parameters, which are suitable for a self-biased GHz circulator design. The size of the ferrite disk was 1.58 mm in thickness and 13.5 mm in diameter. The saturation magnetization ($4{\pi}M_s$) of 3900 Gauss, internal magnetic field ($H_{in}$) of 1 kOe, and ferromagnetic linewidth (${\Delta}H$) of 354 Oe were used in circulator performance simulation. The simulation results show the isolation of 36.4 dB and insertion loss of 2.76 dB at 2.6 GHz and were compared to measured results. A Ni-Zn ferrite circulator was fabricated based on the above design parameters. An out-of-plane DC magnetic field ($H_0$) of 4.8 kOe was applied to the fabricated circulator to measure isolation, insertion loss, and bandwidth. Experimental magnetic parameters for the ferrite were $H_{in}$ of about 1.33 kOe and $4{\pi}M_s$ of 3935 Gauss. The isolation 43.9 dB and insertion loss of 5.6 dB measured at 2.5 GHz are in close agreement with the simulated results of the designed ferrite circulator. Based on the simulated and experimental results, we demonstrate that the following magnetic parameters are suitable for 2 GHz self-biased circulator design: $4{\pi}M_r$ of 3900 Gauss, $H_a$ of 4.5 kOe, $H_c$ greater than 3.4 kOe, and ${\Delta}H$ of 50 Oe.

• Korean Journal of Optics and Photonics
• /
• v.15 no.1
• /
• pp.39-45
• /
• 2004

A computer simulation about removal of the 0-th order diffraction is achieved by using numerical reconstruction in digital holography and the Fourier transform method. A light intensity distribution hologram is generated through numerical calculation of the interference pattern. Additionally a phase hologram without the 0-th order diffraction is generated. The removal function for elimination of the 0-the order diffraction is introduced and the numerical reconstructions with several conditions for the removal of the 0-th order diffraction and the production of high quality numerically reconstructed images are tested and compared. The removal function is proven to be more effective at the suppression of the 0-th order diffraction compared with the DC suppression method.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.49 no.3
• /
• pp.175-187
• /
• 2013

Fuel consumption in fisheries is a primary concern due to environmental effects and costs to fishermen. Much research has been carried out to reduce the fuel consumption related to fishing operations. The fuel consumption of fishing gear during fishing operation is generally related to hydrodynamic resistance on the gear. This research demonstrates a new approach using numerical methods to reduce fuel consumption. The results from the simulation were verified with results that mirrored the model experiments. By designing the fishing gear using drawing software, the whole and partial resistance force on the gear can be calculated as a result of simulations. The simulation results will suggest suitable materials or gear structure for reducing the hydrodynamic forces on the gear while maintaining the performance of the gear. Furthermore, the efficiency of low energy used anchovy dragnet as economic point of view will be dealt. This research will helpful to reduce the GHG emissions from fishing operations and lead to reduce fishing costs due to fuel savings.