• Magazine of the Korean Society of Agricultural Engineers
• /
• v.19 no.1
• /
• pp.4279-4295
• /
• 1977

Two types of model were established in order to product the soil moisture content by which information on irrigation could be obtained. Model-I was to represent the soil moisture depletion and was established based on the concept of water balance in a given soil profile. Model-II was a mathematical model derived from the analysis of soil moisture variation curves which were drawn from the observed data. In establishing the Model-I, the method and procedure to estimate parameters for the determination of the variables such as evapotranspirations, effective rainfalls, and drainage amounts were discussed. Empirical equations representing soil moisture variation curves were derived from the observed data as the Model-II. The procedure for forecasting timing and amounts of irrigation under the given soil moisture content was discussed. The established models were checked by comparing the observed data with those predicted by the model. Obtained results are summarized as follows: 1. As a water balance model of a given soil profile, the soil moisture depletion D, could be represented as the equation(2). 2. Among the various empirical formulae for potential evapotranspiration (Etp), Penman's formula was best fit to the data observed with the evaporation pans and tanks in Suweon area. High degree of positive correlation between Penman's predicted data and observed data with a large evaporation pan was confirmed. and the regression enquation was Y=0.7436X+17.2918, where Y represents evaporation rate from large evaporation pan, in mm/10days, and X represents potential evapotranspiration rate estimated by use of Penman's formula. 3. Evapotranspiration, Et, could be estimated from the potential evapotranspiration, Etp, by introducing the consumptive use coefficient, Kc, which was repre sensed by the following relationship: Kc=Kco$.$Ka＋Ks‥‥‥(Eq. 6) where Kco : crop coefficient Ka : coefficient depending on the soil moisture content Ks : correction coefficient a. Crop coefficient. Kco. Crop coefficients of barley, bean, and wheat for each growth stage were found to be dependent on the crop. b. Coefficient depending on the soil moisture content, Ka. The values of Ka for clay loam, sandy loam, and loamy sand revealed a similar tendency to those of Pierce type. c. Correction coefficent, Ks. Following relationships were established to estimate Ks values: Ks=Kc-Kco$.$Ka, where Ks=0 if Kc,=Kco$.$K0$\geq$1.0, otherwise Ks=1-Kco$.$Ka 4. Effective rainfall, Re, was estimated by using following relationships : Re=D, if R-D$\geq$0, otherwise, Re=R 5. The difference between rainfall, R, and the soil moisture depletion D, was taken as drainage amount, Wd. {{{{D= SUM from { {i }=1} to n (Et-Re-I+Wd)}}}} if Wd=0, otherwise, {{{{D= SUM from { {i }=tf} to n (Et-Re-I+Wd)}}}} where tf=2∼3 days. 6. The curves and their corresponding empirical equations for the variation of soil moisture depending on the soil types, soil depths are shown on Fig. 8 (a,b.c,d). The general mathematical model on soil moisture variation depending on seasons, weather, and soil types were as follow: {{{{SMC= SUM ( { C}_{i }Exp( { - lambda }_{i } { t}_{i } )+ { Re}_{i } - { Excess}_{i } )}}}} where SMC : soil moisture content C : constant depending on an initial soil moisture content $\lambda$ : constant depending on season t : time Re : effective rainfall Excess : drainage and excess soil moisture other than drainage. The values of $\lambda$ are shown on Table 1. 7. The timing and amount of irrigation could be predicted by the equation (9-a) and (9-b,c), respectively. 8. Under the given conditions, the model for scheduling irrigation was completed. Fig. 9 show computer flow charts of the model. a. To estimate a potential evapotranspiration, Penman's equation was used if a complete observed meteorological data were available, and Jensen-Haise's equation was used if a forecasted meteorological data were available, However none of the observed or forecasted data were available, the equation (15) was used. b. As an input time data, a crop carlender was used, which was made based on the time when the growth stage of the crop shows it's maximum effective leaf coverage. 9. For the purpose of validation of the models, observed data of soil moiture content under various conditions from May, 1975 to July, 1975 were compared to the data predicted by Model-I and Model-II. Model-I shows the relative error of 4.6 to 14.3 percent which is an acceptable range of error in view of engineering purpose. Model-II shows 3 to 16.7 percent of relative error which is a little larger than the one from the Model-I. 10. Comparing two models, the followings are concluded: Model-I established on the theoretical background can predict with a satisfiable reliability far practical use provided that forecasted meteorological data are available. On the other hand, Model-II was superior to Model-I in it's simplicity, but it needs long period and wide scope of observed data to predict acceptable soil moisture content. Further studies are needed on the Model-II to make it acceptable in practical use.

• Progress in Medical Physics
• /
• v.20 no.4
• /
• pp.290-297
• /
• 2009

In case of radiation treatment using small field high-energy photon beams, an accurate dosimetry is a challenging task because of dosimetrically unfavorable phenomena such as dramatic changes of the dose at the field boundaries, dis-equilibrium of the electrons, and non-uniformity between the detector and the phantom materials. In this study, the absorbed dose in the phantom was measured by using an ion chamber and a diode detector widely used in clinics. $GAFCHROMIC^{(R)}$ EBT films composed of water equivalent materials was also evaluated as a small field detector and compared with ionchamber and diode detectors. The output factors at 10 cm depth of a solid phantom located 100 cm from the 6 MV linear accelerator (Varian, 6 EX) source were measured for 6 field sizes ($5{\times}5\;cm^2$, $2{\times}2\;cm^2$, $1.5{\times}1.5\;cm^2$, $1{\times}1\;cm^2$, $0.7{\times}0.7\;cm^2$ and $0.5{\times}0.5\;cm^2$). As a result, from $5{\times}5\;cm^2$ to $1.5{\times}1.5\;cm^2$ field sizes, absorbed doses from three detectors were accurately identified within 1%. Wheres, the ion chamber underestimated dose compared to other detectors in the field sizes less than $1{\times}1\;cm^2$. In order to correct the observed underestimation, a convolution method was employed to eliminate the volume averaging effect of an ion chamber. Finally, in $1{\times}1\;cm^2$ field the absorbed dose with a diode detector was about 3% higher than that with the EBT film while the dose with the ion chamber after volume correction was 1% lower. For $0.5{\times}0.5\;cm^2$ field, the dose with the diode detector was 1% larger than that with the EBT film while dose with volume corrected ionization chamber was 7% lower. In conclusion, the possibility of $GAFCHROMIC^{(R)}$ EBT film as an small field dosimeter was tested and further investigation will be proceed using Monte Calro simulation.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.6 no.2
• /
• pp.81-92
• /
• 2001

According to analyses of high-resolution seismic profiles (air gun, sparker, and SBP) and a deep-drill core(YSDP 105) in the mid-eastern Yellow Sea, stratigraphic and geoacoustic models have been established and seismo-acoustic modeling has been fulfilled using ray tracing of finite element method. Stratigraphic model reflects seismo-, litho-, and chrono-stratigraphic sequences formed under a significant influence of Quaternary glacio-eustatic sea-level fluctuations. Each sequence consists of terrestrial to very-shallow-marine coarse-grained lowstand systems tract and tidal fine-grained transgressive to highstand systems tract. Based on mean grain-size data (121 samples) of the drill core, bulk density and P-wave velocity of depositional units have been inferred and extrapolated down to a depth of the recovery using the Hamilton's regression equations. As goo-acoustic parameters, the 121 pairs of bulk density and P-wave velocity have been averaged on each unit of the stratigraphic model. As a result of computer ray-tracing simulation of the subsurface strata, we have found that there are complex ray paths and many acoustic-shadow zones owing to the presence of irregular layer boundaries and low-velocity layers.

• 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.

• The Korean Journal of Nuclear Medicine
• /
• v.32 no.1
• /
• pp.32-49
• /
• 1998

For estimation of regional myocardial blood flow with O-15 water PET, a few modifications considering partial volume effect based on single compartment model have been proposed. In this study, we attempted to quantify the degree of heterogeneity and to show the effect of tissue flow heterogeneity on partition coefficient(${\lambda}$) and to find the relation between perfusable tissue index(PTI) and ${\lambda}$ by computer simulation using two modified models. We simulated tissue curves for the regions with homogeneous and heterogeneous blood flow over a various flow range(0.2-4.0ml/g/min). Simulated heterogeneous tissue composed of 4 subregions of the same or different size of block which have different homogeneous flow and different degree of slope of distribution of blood flow. We measured the index representing heterogeneity of distribution of blood flow for each heterogeneous tissue by the constitution heterogeneity(CH). For model I, we assumed that tissue recovery coefficient ($F_{MME}$) was the product of partial volume effect($F_{MMF}$) and PTI. Using model I, PTI, flow, and $F_{MM}$ were estimated. For model II, we assumed that partition coefficient was another variable which could represent tissue characteristics of heterogeneity of flow distribution. Using model II, PTI, flow and ${\lambda}$ were estimated. For the simulated tissue with homogeneous flow, both models gave exactly the same estimates, of three parameters. For the simulated tissue with heterogeneous flow distribution, in model I, flow and $F_{MM}$ were correctly estimated as CH was increased moderately. In model II, flow and ${\lambda}$ were decreased curvi-linearly as CH was increased. The degree of underestimation of ${\lambda}$ obtained using model II, was correlated with CH. The degree of underestimation of flow was dependent on the degree of underestimation of ${\lambda}$. PTI was somewhat overestimated and did not change according to CH. We conclude that estimated ${\lambda}$ reflect the degree of tissue heterogeneity of flow distribution. We could use the degree of underestimation of ${\lambda}$ to find the characteristic heterogeneity of tissue flow and use ${\lambda}$ to recover the underestimated flow.

• Korean Circulation Journal
• /
• v.52 no.9
• /
• pp.699-711
• /
• 2022

Background and Objectives: We investigated whether extra-pulmonary vein (PV) ablation targeting a high maximal slope of the action potential duration restitution curve (Smax) improves the rhythm outcome of persistent atrial fibrillation (PeAF) ablation. Methods: In this open-label, multi-center, randomized, and controlled trial, 178 PeAF patients were randomized with 1:1 ratio to computational modeling-guided virtual Smax ablation (V-Smax) or empirical ablation (E-ABL) groups. Smax maps were generated by computational modeling based on atrial substrate maps acquired during clinical procedures in sinus rhythm. Smax maps were generated during the clinical PV isolation (PVI). The V-Smax group underwent an additional extra-PV ablation after PVI targeting the virtual high Smax sites. Results: After a mean follow-up period of 12.3±5.2 months, the clinical recurrence rates (25.6% vs. 23.9% in the V-Smax and the E-ABL group, p=0.880) or recurrence appearing as atrial tachycardia (11.1% vs. 5.7%, p=0.169) did not differ between the 2 groups. The post-ablation cardioversion rate was higher in the V-Smax group than E-ABL group (14.4% vs. 5.7%, p=0.027). Among antiarrhythmic drug-free patients (n=129), the AF freedom rate was 78.7% in the V-Smax group and 80.9% in the E-ABL group (p=0.776). The total procedure time was longer in the V-Smax group (p=0.008), but no significant difference was found in the major complication rates (p=0.497) between the groups. Conclusions: Unlike a dominant frequency ablation, the computational modeling-guided V-Smax ablation did not improve the rhythm outcome of the PeAF ablation and had a longer procedure time.

• 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.