• The Journal of Engineering Geology
• /
• v.15 no.3
• /
• pp.245-255
• /
• 2005

Three dimensional Monte-Carlo simulations of non-ergodic transport of a non-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 the simulation uncertainties. Ensemble averages of the second spatial moments of the plume, $$lt;S_{ij}'(t',l')$gt;$ and plume centroid variances, $$lt;R_{ij}'(t',l')$gt;$ were simulated with 3200 Monte Carlo runs for three variances of log K, $\omega^2_y1.0,,2.5,$ and 5.0, and three dimensionless lengths of line plume sources ( l=,5 and 10) 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 not fit well with the first order results. The first order theoretical results definitely underestimated the simulated transverse second spatial moments for the aquifers of large u: and small initial plume sources. The ergodic condition for the second spatial moments is far from reaching, and the first order theoretical results of the transverse second spatial moment of the ergodic plume slightly underestimated the simulated moments.

• Journal of Korean Society of Environmental Engineers
• /
• v.31 no.5
• /
• pp.358-363
• /
• 2009

The HSPF (Hydrological Simulation Program-Fortran) model was applied to the Kyoungan stream watershed to analyze effects of watershed management plans on receiving water qualities. Utilizing BASINS 3.1 GIS program, the Kyoungan stream watershed was divided into 57 sub-basins and model input parameters were obtained, from DEM (Digital Elevation Model), land use type, stream map, and wastewater treatment facilities, etc.. The hydrologic module of the model was validated based on the measured meteorological data and stream flow data. Then the model was calibrated and verified against the field measurements of water qualities, including temperature, DO, BOD, $NO_3-N$, $NH_3-N$, Org-N, TN and TP. In most cases, there were reasonable agreements between measurements and predictions. The validated model was used to analyze the water quality improvements in the main stream of Kyoungan stream according to the watershed management plans in sub-basins, which are three different scenarios: water quality improvement in tributaries through watershed management activities, expansion and up-grade of wastewater treatment plants, and application of first and second scenarios together. It was concluded that expansion and upgrade of wastewater treatment plants would be more effective than watershed management activities. In order to improve water qualities to the satisfactory level, both watershed management and point source control must be required in the Kyoungan stream.

• Journal of Korea Water Resources Association
• /
• v.48 no.12
• /
• pp.1023-1035
• /
• 2015

This study applies 3-D N-S solver based on PBM (Porous Body Model), LED-WASS-3D ver 2.0 to directly analyze non-linear interaction of seawater-freshwater-coastal aquifer in order to simulate the seawater infiltration into coastal aquifer. This numerical simulation is the first trial in Korea, as well as unusual and new numerical analysis abroad. Firstly, to validate the applied numerical model, the validity and effectiveness was verified for the numerical model by comparing and considering it with the result of laboratory experiment for seawater-freshwater interface in coastal aquifer. And then it simulated the seawater infiltration into coastal aquifer considering the changed levels of seawater and groundwater in order to analyze the distribution characteristics of flow field and seawater-freshwater interface of coastal aquifer as the level difference between seawater and groundwater and rate of seawater level (${\Delta}h/h$) increased. In addition, the characteristics of seawater infiltration were analyzed from the vertical salinity in the coastal aquifer by ${\Delta}h/h$, which cannot be obtained from existing non-diffusion numerical models. Finally, it analyzed the effect of ${\Delta}h/h$ on the seawater infiltration distance in coastal aquifer, which was indexed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.2
• /
• pp.92-108
• /
• 2017

Present work shows a numerical method to analysis of interaction analysis between solitary wave and onshore bridge. Numerical simulation is carried out by TWOPM-3D (three-dimensional one-field model for immiscible two-phase flows), which is based on Navier-Stokes solver. To do this, the solitary wave is generated numerically in numerical wave channel, and numerical results and experimental results were compared and analyzed in order to verify the applicability of force acting on an onshore bridge. From this, we discussed precisely the characteristics of horizontal and vertical forces (uplift and downward forces) changes including water level and velocity changes due to the variation of solitary wave height, water depth, onshore bridge's location and type, and number of girder. Furthermore, It is revealed that the maximum horizontal and vertical forces acting on the girder bridge show different varying properties according to the number of girder, although each maximum force acting on the girder bridge is proportional to the increasement of incident solitary wave height, and the entrained air in the fluid flow affects the vertical force highly.

• Journal of Korea Water Resources Association
• /
• v.38 no.8 s.157
• /
• pp.645-653
• /
• 2005

Water cycle analysis in the Cheonggyecheon watershed(river length: 13.75 km, area: $50.96\;km^2$) was performed using WEP model, a physically based distributed rainfall-runoff model. As the application results of the model, the hydrological characteristics of the Cheonggyecheon watershed are significantly consistent with those of a typical urbanized watershed. The direct runoff from the watershed was larger and the evapotranspiration. was lower, and the response of runoff to rainfall was occurred very fast, as compared to forest watersheds. The river channel routing simulation results are similar to the change pattern and scale of the field data. The possible supply period of instream flow from Cheonggyecheoon watershed itself was estimated using WEP. According to the WEP simulation results for the annual water balance of the Cheonggyecheon watershed in 2002, the amount of direct runoff, infiltration and evapotranspiration were 830 mm, 388 mm and 397 mm respectively for an annual precipitation of 1,388 mm. The runoff to rivers was 1,288 mm. And the proportion of direct runoff, intermediate runoff and groundwater runoff were $67.6\%,\;12.7\%$ and $19.7\%$ respectively.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.49 no.2
• /
• pp.49-60
• /
• 2007

The alpine agricultural activities are usually performed at higher and steep areas in nature. Thus, significant amounts of soil erosion are occurring compared with those from other areas. Thus, the soil erosion induced environmental impacts in these areas are getting greater. The Doam watershed is located at alpine areas and it has been well known that the agricultural activities in the watershed are causing accelerated soil erosion and water quality degradations. Many modeling approaches were employed to solve soil erosion and water quality issues. In this study, the Soil and Water Assessment Tool (SWAT) model was utilized to simulate the hydrologic and sediment behaviors in the Doam watershed. In many previous modeling studies, the digital soil map and its corresponding soil properties were used without modification to reflect soil conditioning at many agricultural fields of the Doam watershed. Thus, the soil sample was taken at the agricultural field within the Doam watershed and analyzed for its physical properties. In this study, the digital topsoil properties in the agricultural fields within the Doam watershed were replaced with the soil properties for reconditioned soil analyzed in this study to simulate the impacts of using soil properties for reconditioned soil in hydrologic and sediment modeling at the Doam watershed using the SWAT model. The hydrologic component of the SWAT model was calibrated and validated for measured flow data from 2002 to 2003. The $R^2$ value was 0.79 and the EI value was 0.53 for weekly simulated data. The calibrated model parameters were used for hydrologic component validation and the $R^2$ value was 0.86 and the EI value was 0.74 for weekly data. For sediment comparison, the $R^2$ value was 0.67 and the EI value was 0.59. These statistics improved with the use of soil properties of the reconditioned soil in the field compared with the results obtained without considering soil reconditioning. The simulated sediment amounts with and without considering the soil properties of the reconditioned soil were 284,813 ton and 158,369 ton, respectively. This result indicates that there could be approximately 79% of errors in estimated sediment yield at the Doam watershed, although the model comparison with the measured data gave similar satisfactory statistics with and without considering soil properties from the reconditioned soil.

• Wind and Structures
• /
• v.27 no.1
• /
• pp.11-27
• /
• 2018

The straight-cone steel cooling tower is a novel type of structure, which has a distinct aerodynamic distribution on the internal surface of the tower cylinder compared with conventional hyperbolic concrete cooling towers. Especially in the extreme weather conditions of strong wind and heavy rain, heavy rain also has a direct impact on aerodynamic force on the internal surface and changes the turbulence effect of pulsating wind, but existing studies mainly focus on the impact effect brought by wind-driven rain to structure surface. In addition, for the indirect air cooled cooling tower, different additional ventilation rate of shutters produces a considerable interference to air movement inside the tower and also to the action mechanism of loads. To solve the problem, a straight-cone steel cooling towerstanding 189 m high and currently being constructed is taken as the research object in this study. The algorithm for two-way coupling between wind and rain is adopted. Simulation of wind field and raindrops is performed with continuous phase and discrete phase models, respectively, under the general principles of computational fluid dynamics (CFD). Firstly, the rule of influence of 9 combinations of wind sped and rainfall intensity on flow field mechanism, the volume of wind-driven rain, additional action force of raindrops and equivalent internal pressure coefficient of the tower cylinder is analyzed. On this basis, the internal pressures of the cooling tower under the most unfavorable working condition are compared between four ventilation rates of shutters (0%, 15%, 30% and 100%). The results show that the 3D effect of equivalent internal pressure coefficient is the most significant when considering two-way coupling between wind and rain. Additional load imposed by raindrops on the internal surface of the tower accounts for an extremely small proportion of total wind load, the maximum being only 0.245%. This occurs under the combination of 20 m/s wind velocity and 200 mm/h rainfall intensity. Ventilation rate of shutters not only changes the air movement inside the tower, but also affects the accumulated amount and distribution of raindrops on the internal surface.

• Radiation Oncology Journal
• /
• v.20 no.1
• /
• pp.81-90
• /
• 2002

Purpose : To establish and verify the proper and the practical IMRT (Intensity--modulated radiation therapy) patient QA (Quality Assurance). Materials and Methods : An IMRT QA which consists of 3 steps and 16 items were designed and examined the validity of the program by applying to 9 patients, 12 IMRT cases of various sites. The three step OA program consists of RTP related QA, treatment information flow QA, and a treatment delivery QA procedure. The evaluation of organ constraints, the validity of the point dose, and the dose distribution are major issues in the RTP related QA procedure. The leaf sequence file generation, the evaluation of the MLC control file, the comparison of the dry run film, and the IMRT field simulate image were included in the treatment information flow procedure QA. The patient setup QA, the verification of the IMRT treatment fields to the patients, and the examination of the data in the Record & Verify system make up the treatment delivery QA procedure. Results : The point dose measurement results of 10 cases showed good agreement with the RTP calculation within $3\%$. One case showed more than a $3\%$ difference and the other case showed more than $5\%$, which was out side the tolerance level. We could not find any differences of more than 2 mm between the RTP leaf sequence and the dry run film. Film dosimetry and the dose distribution from the phantom plan showed the same tendency, but quantitative analysis was not possible because of the film dosimetry nature. No error had been found from the MLC control file and one mis-registration case was found before treatment. Conclusion : This study shows the usefulness and the necessity of the IMRT patient QA program. The whole procedure of this program should be peformed, especially by institutions that have just started to accumulate experience. But, the program is too complex and time consuming. Therefore, we propose practical and essential QA items for institutions in which the IMRT is performed as a routine procedure.

• Economic and Environmental Geology
• /
• v.47 no.6
• /
• pp.589-599
• /
• 2014

This study selected optimal sites in Myeongsasimni located in west coast of Korea for stratigraphic research containing extreme climate event during quaternary period by spatio-temporal analyses of changes in sedimentary environment and land use employing 1918 topographic map, 2000 digital terrain map, 1976 and 2012 air photographies. The study area shows no significant changes in topographic characteristics that hilly areas with relatively large variations in elevation are distributed over north and south part of the study area, and sand dues are developed along the coast line. Moreover, flat low lying areas are located at the back side of the sand dues. The movement of surface run off and sediment loads shows two major trends of inland direction flow from back sides of sand dunes and outland direction flow from high terrains inland, and the two flows merge into the stream located in the center of the study area. Two sink with individual area of $0.2km^2$ are observed in Yongjeong-ri and Jaryong-ri which are located in south central part and south part of the study area, respectively. In addition, sea level change simulation reveals that $3.4km^2$ and $3.64km^2$ are inundated with 3 m of sea level rise in 1918 and 2000, respectively, and it would contribute to chase sea level change records preserved in stratigraphy. The inundated areas overlaps well with sink areas where it indicates the low lying areas located in south cental and south part of the study area are identical for sediment accumulation. The areas with minimal human impact on sediment records over last 100 years are $3.51km^2$ distributed over central and south part of the study area with the land use changes of mud and rice field in 1918 to rice field in 2012. The candidate sites of $0.15km^2$ in central part and $0.09km^2$ in south part are identified for preferable locations of geologic record of extreme climate events during quaternary period based on the overlay analysis of optimal sedimentary environment and land use changes.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.174-183
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (II).