• Progress in Medical Physics
• /
• v.16 no.4
• /
• pp.192-201
• /
• 2005

Experiments and simulation were done to study the impact of contrast agent when CT scan was used to attenuation correction for PET Images in PET/CT system. Whole body phantom was imaged with various concentration of iodine-based contrast agent using CT. Mathematical emission and transmission density map with liver were made to simulate for whole body FDG Imaging. A variety of factors were estimated, including non-uniform enhancement of contrast agent, concentration and distribution size of contrast agent, noise level, image resolution, reconstruction algorithm, hypo-attenuation of contrast agent, and different time phases for contrast agent. Experimental studies showed that Hounsfield unit depends on the concentration of contrast agent and tube voltage. From the simulation data, contrast agents Introduced artifacts and degraded image quality on the attenuation-corrected PET images. The severity of these effects depends on a variety of factors, including the concentration and distribution size of contrast agent, the noise levels, and the Image resolution. These results Indicated that the impact of contrast agents should be considered with a full understanding of their potential problems in clinical PET/CT images.

• Progress in Medical Physics
• /
• v.21 no.4
• /
• pp.311-322
• /
• 2010

In proton therapy, the analysis of secondary particles is important due to delivered dose outside the target volume and thus increased potential risk for the development of secondary cancer. The purpose of this study is to analyze the influence of secondary particles from proton beams on fluence and energy deposition in the presence of inhomogeneous material by using Geant4 simulation toolkit. The inhomogeneity was modeled with the condition that the adipose tissue, bone and lung equivalent slab with thickness of 2 cm were inserted at 30% (Plateau region) and 80% (Bragg peak region) dose points of maximum dose in Bragg curve. The energy of proton was varied with 100, 130, 160 and 190 MeV for energy dependency. The results for secondary particles were presented for the fluence and deposited energy of secondary particles at inhomogeneous condition. Our study demonstrates that the fluence of secondary particles is neither influenced insertion of inhomogeneties nor the energy of initial proton, while there is a little effect by material density. The deposited energy of secondary particles has a difference in the position placed inhomogeneous materials. In the Plateau region, deposited energy of secondary particles mostly depends on the density of inserted materials. Deposited energy in the Bragg region, in otherwise, is influenced by both density of inserted material and initial energy of proton beams. Our results suggest a possibility of prediction about the distribution of secondary particles within complex heterogeneity.

• Progress in Medical Physics
• /
• v.15 no.4
• /
• pp.215-219
• /
• 2004

Due to their excellence for the high-energy therapy range of photon beams, researchers show increasing interest in applying MOSFET dosimeters to low- and medium-energy applications. In this energy range, however, MOSFET dosimeter is complicated by the fact that the interaction probability of photons shows significant dependence on the atomic number, Z, due to photoelectric effect. The objective of this study is to develop a very detailed 3-dimensional Monte Carlo simulation model of a MOSFET dosimeter for radiological characterizations and calibrations. The sensitive volume of the High-Sensitivity MOSFET dosimeter is very thin (1 ${\mu}{\textrm}{m}$) and the standard MCNP tallies do not accurately determine absorbed dose to the sensitive volume. Therefore, we need to score the energy deposition directly from electrons. The developed model was then used to study various radiological characteristics of the MOSFET dosimeter. the energy dependence was quantified for the energy range 15 keV to 6 MeV; finding maximum dependence of 6.6 at about 40 keV. A commercial computer code, Sabrina, was used to read the particle track information from an MCNP simulation and count the tracks of simulated electrons. The MOSFET dosimeter estimated the calibration factor by 1.16 when the dosimeter was at 15 cm depth in tissue phantom for 662 keV incident photons. Our results showed that the MOSFET dosimeter estimated by 1.11 for 1.25 MeV photons for the same condition.

• Progress in Medical Physics
• /
• v.13 no.1
• /
• pp.37-43
• /
• 2002

Effect of scatter media materials and thickness, located between radioactivity and small gamma camera, on imaging characteristics was evaluated. The small gamma camera developed for breast imaging was consisted of collimator, NaI(T1) crystal(60$\times$60$\times$6 ㎣), PSPMT(position sensitive photomultimplier tube), NIMs and personal computer. Monte Carlo simulation was performed to evaluate the system sensitivity with different scatter media thickness(0~8 cm) and materials(air and acrylic) with parallel hole collimator and diverging collimator. The sensitivity and spatial resolution was measured using the small gamma camera with the same condition applied to the simulation. Counts was decreased by 10%(air) and 54%(acryl) with the parallel hole collimator and by 35%(air) and 63%(acryl) with the diverging collimator. Spatial resolution was decreased as increasing the thickness of scatter media. This study substantiate the importance of a gamma camera positioning and the minimization of the distance between detector and target lesion in the clinical application of a gamma camera.

• Progress in Medical Physics
• /
• v.24 no.4
• /
• pp.303-314
• /
• 2013

Image quality of computed tomography (CT) is very vulnerable to metal artifacts. Recently, the thickness and background normalization techniques have been introduced. Since they provide flat sinograms, it is easy to determine metal traces and a simple linear interpolation would be enough to describe the missing data in sinograms. In this study, we have developed a theory describing two normalization methods and compared two methods with respect to various sizes and numbers of metal inserts by using simple numerical simulations. The developed theory showed that the background normalization provide flatter sinograms than the thickness normalization, which was validated with the simulation results. Numerical simulation results with respect to various sizes and numbers of metal inserts showed that the background normalization was better than the thickness normalization for metal artifact corrections. Although the residual artifacts still existed, we have showed that the background normalization without the segmentation procedure was better than the thickness normalization for metal artifact corrections. Since the background normalization without the segmentation procedure is simple and it does not require any users' intervention, it can be readily installed in conventional CT systems.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2013.08a
• /
• pp.221-221
• /
• 2013

Vertical LED (VLED) has been recognized as a way to obtain the high-power LED due to their advantages [1]. However, approximately 4% of the light generated from the active region is extracted, if the light extraction from side walls and back side is neglected because of Fresnel reflection (FR) and total internal reflection (TIR) [2,3]. In this study, the optical simulation of the VLED with the various microstructures was performed. Among them, the microlens having the diameter of 3 ${\mu}m$ and the height of 1.5 ${\mu}m$ shown the best result was chosen, and then, optimized microlens was formed on a GaN template using conventional semiconductor process. Various microstructures were proposed to improve the light extraction efficiency (LEE) of the VLED for the simulation. The LEE was simulated using LightTools based on a Monte Carlo ray tracing. The microstructures with hemisphere, cone, truncated and cylinder pattern having diameter of 3 ${\mu}m$ were employed on the top layer of the VLED respectively. The improvement of the LEE by using the microstructure is 87% for the hemisphere, 77% for the cone, 53% for the truncated, 21% for the cylinder, compared with the LEE of the flat surface at the reflectance of 85%. The LEE was increased by 88% at the height of 1.5 ${\mu}m$, compared with the LEE of the flat surface. We found that the microlens on the top layer is the most suitable for increasing the LEE. In order to apply the proposed microlens on n-GaN surface, we fabricated microlens on a GaN template. A photoresist array having hexagonal-closed packed microlens was fabricated on the GaN template. Then, optimization of etching the GaN template was performed using a dry etching process with ICP-RIE. The dry etching carried out using a gas mixture of Cl2 and Ar, each having a flow rate of 16 sccm and 10 sccm, respectively with RF power of 50 W, ICP power of 900 W and chamber pressure of 2 mTorr was the optimum etching condition as shown in Fig. 2(a).

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.43 no.1
• /
• pp.49-61
• /
• 2007

Otter boards in the trawl are the one of essential equipments for the net mouth to be spread to the horizontal direction. Its performance should be considered in the light of the spreading force to the drag and the stability of towing in the water. Up to the present, studies of the otter boards have focused mainly on the drag and lift force, but not on the stability of otter boards movement in 3 dimensional space. In this study, the otter board is regarded as a rigid body, which has six degrees of freedom motion in three dimensional coordinate system. The forces acting on the otter boards are the underwater weight, the resistance of drag and spread forces and the tension on the warps and otter pendants. The equations of forces were derived and substituted into the governing equations of 6 degrees of freedom motion, then the second order of differential equations to the otter boards were established. For the stable numerical integration of this system, Backward Euler one of implicit methods was used. From the results of the numerical calculation, graphic simulation was carried out. The simulations were conducted for 3 types of otter boards having same area with different aspect ratio(${\lambda}=0.5,\;1.0,\;1.5$). The tested gear was mid-water trawl and the towing speed was 4k't. The length of warp was 350m and all conditions were same to each otter board. The results of this study are like this; First, the otter boards of ${\lambda}=1.0$ showed the longest spread distance, and the ${\lambda}=0.5$ showed the shorted spread distance. Second, the otter boards of ${\lambda}=1.0$ and 1.5 showed the upright at the towing speed of 4k't, but the one of ${\lambda}=0.5$ heeled outside. Third, the yawing angles of three otter boards were similar after 100 seconds with the small oscillation. Fourth, it was revealed that the net height and width are affected by the characteristics of otter boards such as the lift coefficient.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2003.09a
• /
• pp.44-44
• /
• 2003

As a basic measurement tool in the areas of animal models of human disease, gene expression and therapy, and drug discovery and development, small animal PET imaging is being used increasingly. An ideal small animal PET should have high sensitivity and high and uniform resolution across the field of view to achieve high image quality. However, the combination of long narrow pixellated crystal array and small ring diameter of small animal PET leads to the degradation of spatial resolution for the source located at off center. This degradation of resolution can be improved by determining the depth of interaction (DOI) in the crystal and by taking into account the information in sorting the coincident events. Among a number of 001 identification schemes, dual layer phsowich detector has been widely investigated by many research groups due to its practicability and effectiveness on extracting DOI information. However, the effects of each crystal length composing dual layer phoswich detector on DOI measurements and image qualities were not fully characterized. In order to minimize the DOI effect, the length of each layer of phoswich detector should be optimized. The aim of this study was to perform simulations using a simulation tool, GATE to design the optimum lengths of crystals composing a dual layer phoswich detector. The simulated small PET system employed LSO front layer LuYAP back layer phoswich detector modules and the module consisted of 8${\times}$8 arrays of dual layer crystals with 2 mm ${\times}$ 2 mm sensitive area coupled to a Hamamatsu R7600 00 M64 PSPMT. Sensitivities and variation of radial resolutions were simulated by varying the length of LSO front layer from 0 to 10 mm while the total length (LSO + LuYAP) was fixed to 20 mm for 10 cm diameter ring scanner. The radial resolution uniformity was markedly improved by using DOI information. There existed the optimal lengths of crystal layers to minimize the variation of radial resolutions. In 10 cm ring scanner configuration, the radial resolution was kept below 3.4 mm over 8 cm FOV while the sensitivity was higher than 7.4% for LSO 5 mm : LuYAP 15 mm phoswich detector. In this study, the optimal length of dual layer phoswich detector was derived to achieve high and uniform radial resolution.

• Journal of the Korean Magnetics Society
• /
• v.17 no.3
• /
• pp.124-127
• /
• 2007

Magnetic tunnel junctions (MTJs), which consisted of amorphous CoFeSiB layers, were investigated. The CoFeSiB layers were used to substitute for the traditionally used CoFe and/or NiFe layers with an emphasis given on understanding the effect of the amorphous free layer on the switching characteristics of the MTJs. CoFeSiB has a lower saturation magnetization ($M_s\;:\;560\;emu/cm^3$) and a higher anisotropy constant ($K_u\;:\;2800\;erg/cm^3$) than CoFe and NiFe, respectively. An exchange coupling energy ($J_{ex}$) of $-0.003\;erg/cm^2$ was observed by inserting a 1.0 nm Ru layer in between CoFeSiB layers. In the Si/$SiO_2$/Ta 45/Ru 9.5/IrMn 10/CoFe 7/$AlO_x$/CoFeSiB 7 or CoFeSiB (t)/Ru 1.0/CoFeSiB (7-t)/Ru 60 (in nm) MTJs structure, it was found that the size dependence of the switching field originated in the lower $J_{ex}$ using the experimental and simulation results. The CoFeSiB synthetic antiferromagnet structures were proved to be beneficial for the switching characteristics such as reducing the coercivity ($H_c$) and increasing the sensitivity in micrometer size, even in submicrometer sized elements.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.5
• /
• pp.320-333
• /
• 2019

In this study, a physics-based 3D morphology model for the estimation of an erosion rate of nourished beach is newly proposed. As a hydrodynamic module, IHFOAM toolbox having its roots on the OpenFoam is used. On the other hand, the morphology model comprised a transport equation for suspended sediment, and Exner type equation derived from the viewpoint of sediment budget with the bed load being taken to accounted. In doing so, the incipient motion of sediment is determined based on the Shields Diagram, while the bottom suspended sediment concentration, the bed load transport rate is figured out using the bottom shearing stress directly calculated from the numerically simulated flow field rather than the conventional quadratic law and frictional coefficient. In order to verify the proposed morphology model, we numerically simulate the nonlinear shoaling, breaking over the uniform beach of 1/m slope, and its ensuing morphology change. Numerical results show that the partially skewed, and asymmetric bottom shearing stresses can be successfully simulated. It was shown that sediments suspended and eroded at the foreshore by wave breaking are gradually drifted toward a shore and accumulated in the process of up-rush, which eventually leads to the formation of swash bar. It is also worth mentioning that the breaker bar formed by the sediments dragged by the back-wash flow which commences at the pinnacle of up-rush as the back-wash flow gets weakened due to the increased depth was successfully duplicated in the numerical simulation.