• Journal of Soil and Groundwater Environment
• /
• v.10 no.6
• /
• pp.20-31
• /
• 2005

Three dimensional Monte-Carlo simulations of non-ergodic transport of a lion-reactive solute plume by steady-state groundwater flow under a uniform mean velocity in isotropic heterogeneous aquifers were conducted. The log-normally distributed hydraulic conductivity, K(x), is modeled as a random field. Significant efforts are made to reduce tile simulation uncertainties. Ensemble averages of the second spatial moments of the plume and plume centroid variances were simulated with 1600 Monte Carlo runs for three variances of log K, ${\sigma}_Y^2=0.09,\;0.23$, and 0.46, and three dimensionless lengths of line plume sources normal to the mean velocity. The simulated second spatial moment and the plume centroid variance in longitudinal direction fit well to the first order theoretical results while the simulated transverse moments are generally larger than the first order results. The first order theoretical results significantly underestimated the simulated dimensionless transverse moments for the aquifers of large ${\sigma}_Y^2$ and large dimensionless time. The ergodic condition for the second spatial moments is far from reaching in all cases simulated, and transport In transverse directions may reach ergodic condition much slower than that in longitudinal direction. The evolution of the contaminant transported in a heterogeneous aquifer is not affected by the shape of the initial plume but affected mainly by the degree of the heterogeneity and the size of the initial plume.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.18 no.3
• /
• pp.111-121
• /
• 2013

Accurate prediction of sea water temperature has been emphasized to make precise local weather forecast and to understand change of ecosystem. The Yellow Sea, which has turbid water and strong tidal current, is an unique shallow marginal sea. It is essential to include the effects of the turbidity and the strong tidal mixing for the realistic simulation of temperature distribution in the Yellow Sea. Evaluation of ocean circulation model response to vertical mixing scheme and turbidity is primary objective of this study. Three-dimensional ocean circulation model(Regional Ocean Modeling System) was used to perform numerical simulations. Mellor- Yamada level 2.5 closure (M-Y) and K-Profile Parameterization (KPP) scheme were selected for vertical mixing parameterization in this study. Effect of Jerlov water type 1, 3 and 5 was also evaluated. The simulated temperature distribution was compared with the observed data by National Fisheries Research and Development Institute to estimate model's response to turbidity and vertical mixing schemes in the Yellow Sea. Simulations with M-Y vertical mixing scheme produced relatively stronger vertical mixing and warmer bottom temperature than the observation. KPP scheme produced weaker vertical mixing and did not well reproduce tidal mixing front along the coast. However, KPP scheme keeps bottom temperature closer to the observation. Consequently, numerical ocean circulation simulations with M-Y vertical mixing scheme tends to produce well mixed vertical temperature structure and that with KPP vertical mixing scheme tends to make stratified vertical temperature structure. When Jerlov water type is higher, sea surface temperature is high and sea bottom temperature is low because downward shortwave radiation is almost absorbed near the sea surface.

• The Journal of Korean Society for Radiation Therapy
• /
• v.26 no.1
• /
• pp.137-147
• /
• 2014

Purpose : The purpose of this study is to verify the accuracy of dose delivery according to the patient's breathing cycle in Gated Volumetric Modulated Arc Therapy Materials and Methods : TrueBeam STxTM(Varian Medical System, Palo Alto, CA) was used in this experiment. The Computed tomography(CT) images that were acquired with RANDO Phantom(Alderson Research Laboratories Inc. Stamford. CT, USA), using Computerized treatment planning system(Eclipse 10.0, Varian, USA), were used to create VMAT plans using 10MV FFF with 1500 cGy/fx (case 1, 2, 3) and 220 cGy/fx(case 4, 5, 6) of doserate of 1200 MU/min. The regular respiratory period of 1.5, 2.5, 3.5 and 4.5 sec and the patients respiratory period of 2.2 and 3.5 sec were reproduced with the $QUASAR^{TM}$ Respiratory Motion Phantom(Modus Medical Devices Inc), and it was set up to deliver radiation at the phase mode between the ranges of 30 to 70%. The results were measured at respective respiratory conditions by a 2-Dimensional ion chamber array detector(I'mRT Matrixx, IBA Dosimetry, Germany) and a MultiCube Phantom(IBA Dosimetry, Germany), and the Gamma pass rate(3 mm, 3%) were compared by the IMRT analysis program(OmniPro I'mRT system software Version 1.7b, IBA Dosimetry, Germany) Results : The gamma pass rates of Case 1, 2, 3, 4, 5 and 6 were the results of 100.0, 97.6, 98.1, 96.3, 93.0, 94.8% at a regular respiratory period of 1.5 sec and 98.8, 99.5, 97.5, 99.5, 98.3, 99.6% at 2.5 sec, 99.6, 96.6, 97.5, 99.2, 97.8, 99.1% at 3.5 sec and 99.4, 96.3, 97.2, 99.0, 98.0, 99.3% at 4.5 sec, respectively. When a patient's respiration was reproduced, 97.7, 95.4, 96.2, 98.9, 96.2, 98.4% at average respiratory period of 2.2 sec, and 97.3, 97.5, 96.8, 100.0, 99.3, 99.8% at 3.5 sec, respectively. Conclusion : The experiment showed clinically reliable results of a Gamma pass rate of 95% or more when 2.5 sec or more of a regular breathing period and the patient's breathing were reproduced. While it showed the results of 93.0% and 94.8% at a regular breathing period of 1.5 sec of Case 5 and 6, it could be confirmed that the accurate dose delivery could be possible on the most respiratory conditions because based on the results of 100 patients's respiratory period analysis as no one sustained a respiration of 1.5 sec. But, pretreatment dose verification should be precede because we can't exclude the possibility of error occurrence due to extremely short respiratory period, also a training at the simulation and careful monitoring are necessary for a patient to maintain stable breathing. Consequently, more reliable and accurate treatments can be administered.

• Progress in Medical Physics
• /
• v.21 no.3
• /
• pp.291-297
• /
• 2010

This study examined the dosimetric influence of implanted gold markers in proton therapy and the effects of their positions in the spread-out Bragg peak (SOBP) proton beam. The implanted cylindrical gold markers were 3 mm long and 1.2 mm in diameter. The dosimetric influence of the gold markers was determined with markers at various locations in a proton-beam field. Spatial dose distributions were measured using a three-dimensional moving water phantom and a stereotactic diode detector with an effective diameter of 0.5 mm. Also, a film dosimetry was performed using Gafchromic External Beam Treatment (EBT) film. The GEANT4 simulation toolkit was used for Monte-Carlo simulations to confirm the measurements and to construct the dose-volume histogram with implanting markers. Motion data were obtained from the portal images of 10 patients to investigate the effect of organ motions on the dosimetric influence of markers in the presence of a rectal balloon. The underdosed volume due to a single gold marker, in which the dose was less than 95% of a prescribed amount, was 0.15 cc. The underdosed volume due to the presence of a gold marker is much smaller than the target volume. However, the underdosed volume is inside the gross tumor volume and is not smeared out due to translational prostate motions. The positions of gold markers and the conditions of the proton-beam field give different impacts on the dose distribution of a target with implanted gold markers, and should be considered in all clinical proton-based therapies.

• Journal of Korea Water Resources Association
• /
• v.55 no.12
• /
• pp.1041-1052
• /
• 2022

This study is to estimate the optimal ecological flow and analysis the spatial distribution of fish habitat for Andong dam downstream reach (4,565.7 km²) using PHABSIM (Physical Habiat Simulation System) and River2D. To establish habitat models, the cross-section informations and hydraulic input data were collected uisng the Nakdong river basic plan report. The establishment range of PHABSIM was set up about 410.0 m from Gudam streamflow gauging station (GD) and about 6.0 km including GD for River2D. To select representative fish species and construct HSI (Habitat Suitability Index), the fish survey was performed at Pungji bridge where showed well the physical characteristics of target stream located downstream of GD. As a result of the fish survey, Zacco platypus was showed highly relative abundance resulting in selecting as the representative fish species, and HSI was constructed using physical habitat characteristics of the Zacco platypus. The optimal range of HSI was 0.3~0.5 m/s at the velocity suitability index, 0.4~0.6 m at the depth suitability index, and the substrate was sand to fine gravel. As a result of estimating the optimal ecological flow by applying HSI to PHABSIM, the optimal ecological flow for target stream was 20.0 m³/sec. As a result of analysis two-dimensional spatial analysis of fish habitat using River2D, WUA (Weighted Usable Area) was estimated 107,392.0 m²/1000 m under the ecological flow condition and it showed the fish habitat was secured throughout the target stream compared with Q₃₅₅ condition.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.2
• /
• pp.63-75
• /
• 2013

The Inner Port Phase 2 area of the Pyeongtaek-Dangjin Port is enclosed by a total of three permeable sea-walls, and the disposal site to the east of the Inner Port Phase 2 is also enclosed by two permeable sea-walls. The maximum tidal range measured in the Inner Port Phase 2 and in the disposal site in May 2010 is 4.70 and 2.32 m, respectively. It reaches up to 54 and 27%, respectively of 8.74 m measured simultaneously in the exterior. Regression formulas between the difference of hydraulic head and the rate of interior water volume change, are induced. A three-dimensional numerical hydrodynamic model for the Asan Bay is constructed incorporating a module to compute water discharge through the permeable sea-walls at each computation time step by employing the formulas. Hydrodynamics for the period from 13th to 27th May, 2010 is simulated by driving forces of real-time reconstructed tide with major five constituents($M_2$, $S_2$, $K_1$, $O_1$ and $N_2$) and freshwater discharges from Asan, Sapkyo, Namyang and Seokmoon Sea dikes. The skill scores of modeled mean high waters, mean sea levels and mean low waters are excellent to be 96 to 100% in the interior of permeable sea-walls. Compared with the results of simulation to obstruct the flow through the permeable sea-walls, the maximum current speed increases by 0.05 to 0.10 m/s along the main channel and by 0.1 to 0.2 m/s locally in the exterior of the Outer Sea-wall of Inner Port. The maximum bottom shear stress is also intensified by 0.1 to 0.4 $N/m^2$ in the main channel and by more than 0.4 $N/m^2$ locally around the arched Outer Sea-wall. The module developed to compute the flow through impermeable seawalls can be practically applied to simulate and predict the advection and dispersion of materials, the erosion or deposion of sediments, and the local scouring around coastal structures where large-scale permeable sea-walls are maintained.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.11 no.4
• /
• pp.192-205
• /
• 2009

The unexpected wind over the Mt. Hwawang on 9 February 2009 was deadly when many spectators were watching a traditional event to burn dried grasses and the fire went out of control due to the wind. We analyzed the fatal wind based on wind flow simulations over a digitized complex terrain of the mountain with a localized heating area using a three dimensional computational fluid dynamics model, CFD_NIMR_SNU (Computational Fluid Dynamics_National Institute of Meteorological Research_Seoul National University). Three levels of fire intensity were simulated: no fire, $300^{\circ}C$ and $600^{\circ}C$ of surface temperature at the site on fire. The surface heat accelerated vertical wind speed by as much as $0.7\;m\;s^{-1}$ (for $300^{\circ}C$) and $1.1\;m\;s^{-1}$ (for $600^{\circ}C$) at the center of the fire. Turbulent kinetic energy was increased by the heat itself and by the increased mechanical force, which in turn was generated by the thermal convection. The heating together with the complex terrain and strong boundary wind induced the unexpected high wind conditions with turbulence at the mountain. The CFD_NIMR_SNU model provided valuable analysis data to understand the consequences of the fatal mountain fire. It is suggested that the place of fire was calm at the time of the fire setting due to the elevated terrain of the windward side. The suppression of wind was easily reversed when there was fire, which caused updraft of hot air by the fire and the strong boundary wind. The strong boundary wind in conjunction with the fire event caused the strong turbulence, resulting in many fire casualties. The model can be utilized in turbulence forecasting over a small area due to surface fire in conjunction with a mesoscale weather model to help fire prevention at the field.

• Investigative Magnetic Resonance Imaging
• /
• v.8 no.1
• /
• pp.42-49
• /
• 2004

Purpose : The quantization noise in magnetic resonance imaging (MRI) systems is analyzed. The signal-to-quantization noise ratio (SQNR) in the reconstructed image is derived from the level of quantization in the signal in spatial frequency domain. Based on the derived formula, the SQNRs in various main magnetic fields with different receiver systems are evaluated. From the evaluation, the quantization noise could be a major noise source determining overall system signal-to-noise ratio (SNR) in high field MRI system. A few methods to reduce the quantization noise are suggested. Materials and methods : In Fourier imaging methods, spin density distribution is encoded by phase and frequency encoding gradients in such a way that it becomes a distribution in the spatial frequency domain. Thus the quantization noise in the spatial frequency domain is expressed in terms of the SQNR in the reconstructed image. The validity of the derived formula is confirmed by experiments and computer simulation. Results : Using the derived formula, the SQNRs in various main magnetic fields with various receiver systems are evaluated. Since the quantization noise is proportional to the signal amplitude, yet it cannot be reduced by simple signal averaging, it could be a serious problem in high field imaging. In many receiver systems employing analog-to-digital converters (ADC) of 16 bits/sample, the quantization noise could be a major noise source limiting overall system SNR, especially in a high field imaging. Conclusion : The field strength of MRI system keeps going higher for functional imaging and spectroscopy. In high field MRI system, signal amplitude becomes larger with more susceptibility effect and wider spectral separation. Since the quantization noise is proportional to the signal amplitude, if the conversion bits of the ADCs in the receiver system are not large enough, the increase of signal amplitude may not be fully utilized for the SNR enhancement due to the increase of the quantization noise. Evaluation of the SQNR for various systems using the formula shows that the quantization noise could be a major noise source limiting overall system SNR, especially in three dimensional imaging in a high field imaging. Oversampling and off-center sampling would be an alternative solution to reduce the quantization noise without replacement of the receiver system.

• Radiation Oncology Journal
• /
• v.19 no.1
• /
• pp.66-73
• /
• 2001

Purpose : For the research of Boron Neutron Capture Therapy (BNCT), fast neutrons generated from the MC-50 cyclotron with maximum energy of 34.4 MeV in Korea Cancer Center Hospital were moderated by 70 cm paraffin and then the dose characteristics were investigated. Using these results, we hope to establish the protocol about dose measurement of epi-thermal neutron, to make a basis of dose characteristic of epi-thermal neutron emitted from nuclear reactor, and to find feasibility about accelerator-based BNCT. Method and Materials : For measuring the absorbed dose and dose distribution of fast neutron beams, we used Unidos 10005 (PTW, Germany) electrometer and IC-17 (Far West, USA), IC-18, ElC-1 ion chambers manufactured by A-150 plastic and used IC-l7M ion chamber manufactured by magnesium for gamma dose. There chambers were flushed with tissue equivalent gas and argon gas and then the flow rate was S co per minute. Using Monte Carlo N-Particle (MCNP) code, transport program in mixed field with neutron, photon, electron, two dimensional dose and energy fluence distribution was calculated and there results were compared with measured results. Results : The absorbed dose of fast neutron beams was $6.47\times10^{-3}$ cGy per 1 MU at the 4 cm depth of the water phantom, which is assumed to be effective depth for BNCT. The magnitude of gamma contamination intermingled with fast neutron beams was $65.2{\pm}0.9\%$ at the same depth. In the dose distribution according to the depth of water, the neutron dose decreased linearly and the gamma dose decreased exponentially as the depth was deepened. The factor expressed energy level, $D_{20}/D_{10}$, of the total dose was 0.718. Conclusion : Through the direct measurement using the two ion chambers, which is made different wall materials, and computer calculation of isodose distribution using MCNP simulation method, we have found the dose characteristics of low fluence fast neutron beams. If the power supply and the target material, which generate high voltage and current, will be developed and gamma contamination was reduced by lead or bismuth, we think, it may be possible to accelerator-based BNCT.