• Radiation Oncology Journal
• /
• v.1 no.1
• /
• pp.21-24
• /
• 1983

The thickness of the part being irradiated is finite. Percent depth dose tables being used routinely are generally obtained from dosimetry in a phantom much thickner than usual patient. At or close to exit surface, the dose should be less than that obtained from the percent depth dose tables, because of insufficient volume for backscattering. To know the difference between the true absorbed dose and the dose obtained from percent depth dose table, the doses at or close to the exit surface were measured with plate type ionization chamber with volume of 0.5ml. The results are as follows; 1. In the case of $^{60}Co$, percent depth dose at a given depth increases with underlying phantom thickness up to the 5cm. 2. In the case of $^{60}Co$, the dose correction factor at exit surface which is less than 1, increases with part thickness and decreases with field size. 3. Exposure time may not be corrected when the part above 10cm in thickness is treated by $^{60}Co$. 4. In the case of 10MV x-ray, the dose correction factor is nearly 1 and constant for the underlying phantom thickness and field size, so the correction of monitor unit is not necessary for part thickness.

• The Journal of Engineering Geology
• /
• v.5 no.3
• /
• pp.289-300
• /
• 1995

At the initial stage of the underground reservoir design one should thoroughly consider surface and subsurface hydrology, hydrogeologic characteristics of aquifer system, and the function of cut - off wall because it is linked to the effective management. In this study, three dimensional finite difference model was applied to analyse the function of Ian underground reservoir at Kyungbuk Province. The steady and unsteady state conditions after construction of the underground dam were simulated through the model, and from these results the groundwater budget and the safe yield were determined. The model simulation indicates the infiltration of irrigation water to be one of the major factors of seasonal fluctuation of groundwater level. The recharge rates of irrigation water were estimated as 4.3mm/d during May and June, and 1.7mm/d during July and Agust. Groundwater recharge from the watershed area estimated to about $0.04m^3/s$, almost consistent through the year. In 1984, groundwater discharge through the transverse section of the dam was $0.002m^3/s$ and the optimum yield for two momths(July and Aguest)was $254000m^3$, however, the discharge became $0.013m^3/s$ in1993, implying the failure of cut -off function. without appropaiate of the cut - off wall, optiumum yield during the irrigaton period would be $93, 000m^3$.

• Journal of the Korean GEO-environmental Society
• /
• v.2 no.3
• /
• pp.47-56
• /
• 2001

In the present study, mainly to determine the embarkation time of equipment in dredged clay fills, an analytical approach is performed to predict a variation of the undrained shear strength in the outermost layer. In this approach, Gibson's non-dimensional linear constant defining the relationship between the void ratio and the effective stress is employed. Also in this approach, void ratios and settlements associated with the volume change due to the self-consolidation and the desiccation shrinkage are evaluated at various elapsed times based on the finite difference solution technique proposed by the authors(1999) and the developed computer program named as DSCON. Predicted results(water content ratio, unit weight and undrained shear strength) are compared with those of laboratory model tests conducted with soil samples obtained from the Koheung site. Based on the predicted undrained shear strengths, possible embarkation time of a equipment is also evaluated. In addition, further analyses are made to indirectly verify the efficiency of the analytical approach proposed in the present study using the PSDDF computer program which can consider the drainage efficiency.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.179-192
• /
• 1993

Three-dimensional cohesive sediment transport model, COSETM-3, is develpoed using a finite difference method. The model results are compared with the physical experimental results for the relative concentration with time at the mid-depth of the recirculating flume and are found to be in good agreement. This model is applied to Suyoung Bay in Pusan of Korea to verify the field applicability of the model and to investigate on the SS (suspended solids) diffusion phenomena at the bay. Behaviors of discharging SS from Suyoung River at normal river flow and flood river flow are predicted. The numerical results appear to be reasonable and qualitative agreement with field data. The influence of settling velocity on the concentration distribution of SS is also investigated. In case of not considering settling velocity, SS concentration at surface layer is higher than that at lower layer, but in case of considering settling velocity, SS concentration at lower layer is higher than that at surface layer. The fluctuation of SS concentration at surface layer is large due to the strong mixing, but the fluctuation of the concentration at lower layer is small due to the weak mixing. SS diffusion patterns at flood river flow are similar to those at normal river flow, while the concentration at that flow is so much higher than that at this flow. SS concentration increases with time until the peak discharge occurs, but the concentration decreases with time with decreasing river flow after the peak discharge.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.3
• /
• pp.279-285
• /
• 2011

We perform a numerical study to determine the time of onset of natural convection in a transient hot wire (THW) device for measuring the thermal conductivity of nanofluids. The samples used in this simulation are water-based $Al_2O_3$ nanofluids with volume fractions of 1%, 4%, and 10%, and the properties are calculated by theoretical models and experimental correlations. The THW apparatus using coated wire is modeled by the control-volume-based finite difference method, and the start of natural convection is determined by observing the temperature rise of the wire under a gravity field. The onset time is 11.5 s for water and 41.6 s for water-based $Al_2O_3$ nanofluids predicted by Maxwell thermal conductivity model with a 10% volume fraction. We confirm that the onset time of natural convection of nanofluids in the cylinder increases with the nanoparticle volume fraction. We suggest a correlation for predicting the onset time on the basis of the numerical results. Finally, it is shown that the measurement error due to natural convection is negligible if the measurement using the transient hot wire method is completed before the onset of natural convection in the base fluid.

• Economic and Environmental Geology
• /
• v.22 no.4
• /
• pp.381-393
• /
• 1989

Hydrogeological modelling was performed to evaluate groundwater flow system in Igsan Area. The study area extends over $790km^2$. The geology consists of Jurassic Daebo granite and gneissose granite and Precambrian metamorphic rocks. The capability of pumping yield is the highest in gneissose granite region among them due to comparatively thick weathered zone with thickness ranging from 10m to 25m. The Colorado State University Finite Difference Model was used for the model simulation. The model was divided into 28 rows and 31 columns with variable grid spacing. The model was calibrated under steady-state and unsteady-state conditions. In the steady-state simulation, the model results were compared with measured water table contours in September 1985 with determining hydraulic conductivities and net recharge rates during rainy season. Unsteady state simulation was done to know the aquifer response due to groundwater abstraction. The non- steady state calibration was conducted to determine the distribution and magnitudes of specific yields and discharge/recharge rates during dry season as matching water level altitudes in May 1986. The calibrated model was used to simulate water level vaiation caused by groundwater withdrawal and natural recharge from 1 October, 1985 until 30 September, 1995. The calibrated model can be used to groundwater development schemes on regional groundwater levels, but it cannot be used to simulate local groundwater level change at a specific site.

• Tunnel and Underground Space
• /
• v.19 no.3
• /
• pp.203-212
• /
• 2009

The prediction of near field behavior around an underground high-level radioactive waste repository is important for the repository design as well as the safety assessment. In this study, a sensitivity analysis for seven parameters consisted of design parameters and material properties was carried out using a three-dimensional finite difference code. From the sensitivity analysis, it was found that the effects of borehole spacing, tunnel spacing, cooling time and rock thermal conductivity were more significant than the other parameters. For getting a statistical distribution of buffer and rock temperatures around the repository, an artificial neural network, backpropagation, was applied. The reliability of the trained neural network was tested with the cases with randomly chosen input parameters. When the parameter variation is within ${\pm}10%$, the prediction from the network was found to be reliable with about a 1% error. It was possible to calculate the temperature distribution for many cases quickly with the trained neural network. The buffer and rock temperatures showed a normal distribution with means of $98^{\circ}C$ and $83.9^{\circ}C$ standard deviations of $3.82^{\circ}C$ and $3.67^{\circ}C$, respectively. Using the neural network, it was also possible to estimate the required change in design parameters for reducing the buffer and rock temperatures for $1^{\circ}C$.

• Journal of Soil and Groundwater Environment
• /
• v.6 no.2
• /
• pp.69-81
• /
• 2001

The constitutive relation among capillary pressure, saturation and relative permeability should be predetermined in order to simulate immiscible water-gas flow in porous media. The relation between saturation and relative permeability becomes more important when the capillary force can be disregarded and viscous friction force governs the flow. In this study, a 2-dimensional finite difference numerical model was developed, in which the variation of viscosity with pressure and that of relative permeability with water saturation can be treated. Seven cases of parallel plate tests were performed in order to obtain the characteristic equation of relative permeability which would be used in. the developed numerical model. It was not possible, however, to match the curves of relative permeability from the plate tests with the existing emperical models. Consequently a logistic equation was proposed as a new emperical model. As this model was composed of the parameter involving aperture size, any aperture size of fracture can be applied to the model. For the purpose of verification, the characteristic equation of relative permeability was applied to the developed numerical model and the computed results were compared with those of plate test. As a result of application of numerical model, in order to check the field applicability, to single fracture surrounding an underground storage cavern, the simultaneous flow of water and propane gas was able to be simulated properly by the model.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.