• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.1B
• /
• pp.11-20
• /
• 2008

The spillway type of small and midsize dams in Korea is almost overflow weir. To examine flood control capacity of overflow spillway, FLOW-3D was applied to Daesuho dam and analysis was focused on the discharge of dam spillway by changing weir shape. Overflow phases and discharges of linear labyrinth weir and curved labyrinth weir were compared with those of existing linear ogee weir. Hydraulic model experiment was performed to verify numerical result. Verification results showed that overflow behaviors and flow characteristics in the side channel by hydraulic model experiment and numerical simulation are well matched, and water surface elevation at side wall coincides with each other. When the reservoir elevation was increased up to design flood level, in case of the linear ogee weir the flow over the crest ran through smoothly in the side channel, whereas in cases of linear labyrinth weir and curved labyrinth weirs, the flow discharge was increased by 40 cms, and the flow over the weir crest, rotating counter-clockwise, was submerged in the side channel. The results of the water level-discharge curve revealed that labyrinth weir can increase discharge by 71% compared to the discharge of linear ogee weir at low reservoir elevation since it can have longer effective length. But as water surface elevation rises, the slope of water level-discharge curve of labyrinth weir becomes milder by submergence and nappe interference in the side channel.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.2B
• /
• pp.215-224
• /
• 2008

The effect of streamflow considering future land use change and vegetation index information by climate change scenario was assessed using SLURP (Semi-distributed Land-Use Runoff Process) model. The model was calibrated and verified using 4 years (1999-2002) daily observed streamflow data for the upstream watershed ($260.4km^2$) of Gyeongan water level gauging station. By applying CA-Markov technique, the future land uses (2030, 2060, 2090) were predicted after test the comparison of 2004 Landsat land use and 2004 CA-Markov land use by 1996 and 2000 land use data. The future land use showed a tendency that the forest and paddy decreased while urban, grassland and bareground increased. The future vegetation indices (2030, 2060, 2090) were estimated by the equation of linear regression between monthly NDVI of NOAA AVHRR images and monthly mean temperature of 5 years (1998-2002). Using CCCma CGCM2 simulation result based on SRES A2 and B2 scenario (2030s, 2060s, 2090s) of IPCC and data were downscaled by Stochastic Spatio-Temporal Random Cascade Model (SST-RCM) technique, the model showed that the future runoff ratio was predicted from 13% to 34% while the runoff ratio of 1999-2002 was 59%. On the other hand, the impact on runoff ratio by land use change showed about 0.1% to 1% increase.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.4D
• /
• pp.671-678
• /
• 2006

In this study, to describe vehicle-structure dynamic interaction phenomena with 1/4 vehicle model, nonlinear Hertzian contact spring and nonlinear contact damper are adopted. The external loads acting on 1/4 vehicle model are selfweight of vehicle and geometry information of running surface. The constraint equation on contact surface is implemented by the Penalty method with stabilization and the reaction from constraint violation. To describe pitching motion of various vehicles two types of the displacement constraint equations are exerted to connect between car bodies and between bogie frames, i.e., the rigid body connection and the rigid body connection with pin, respectively. For the time integration of dynamic equations of vehicles and structure Newmark time integration scheme is adopted. To reduce the error caused by inadequate time step size, adaptive time-stepping technique is also adopted. Thus, it is expected that more versatile dynamic interaction phenomena can be described by this approach and it can be applied to various railway dynamic problems with low computational cost.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5B
• /
• pp.447-457
• /
• 2006

This paper presents a numerical model for the simulations of 2D depth-averaged flows. The shallow water equations are solved numerically by the Characteristic Dissipative Galerkin (CDG) finite element method. For validation, the developed model is applied to the hydraulic jump. The computed results are compared with the analytical solution, revealing good agreement. In addition, flow in a contracting channel showing standing waves is simulated. The calculated water surface profile appears to be qualitatively consistent with the observed data. The foregoing results indicate that the model is capable of simulating the abrupt change in flow field. Next, the model is applied to the flow in a $180^{\circ}$ curved channel. The simulated results show that the velocity near the inner bank is faster than that near the outer bank and the water depth near the inner bank is shallower than that near the outer bank. However, the simulated results show that the velocity distribution across the channel is almost uniform in the bend except the reach close to the end of the bend. This is due to the limitation of the governing equations in which the transverse convection of momentum by the secondary flows along a channel bend is not taken into account.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.2C
• /
• pp.133-142
• /
• 2006

Many methods have been tried to more correctly measure rock joint roughness. However, true roughness may be distorted and underestimated due to the sampling interval and measurement method. Thus, currently used measurement methods produce a dead zone and distort roughness profiles. The purpose of this study is to suggest new roughness measurement method by a camera-type 3D scanner as an alternative of currently used methods. First, the underestimation of artificial roughness is analyzed by using the current measurement method such as laser profilometry. Second, we replicate eight specimens from two rock joint surfaces, and digitize by a 3D scanner. Then, the roughness coefficient values obtained from eight numbers of 3D surface data and from three hundred twenty numbers of 2D profiles data are analyzed by using current and new measurement methods. The artificial simulation confirms that the sampling interval is one of main factors for the distortion of roughness and shows that inclination of waviness may not be considered any current methods. The experimental results show that the camera-type 3D scanner produces 10% larger roughness values than current methods. As the proposed new method is a fast, high precision and more accurate method for the roughness measurement, it should be a promising technique in this area.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.4B
• /
• pp.339-346
• /
• 2009

In this study, the flow characteristics of the sea region, where tidal current and shock waves are combined, are identified using a three-dimensional numerical model (Princeton Ocean Model, POM). The model is adopted and applied for simulating the flows of the sea region near the open sections during the seadike closure work of Sihwa Seadike which was closed in 1994. The simulation results show that the shock waves with high velocities propagate through the sections toward the inside and outside of the seadike during the periods of the spring and ebb tides, respectively. It is found that the phenomena of flow separation occur near the shock waves; as the shock waves extend to wider zones after passing the sections, their effects on the tidal current become weak. In addition, the longitudinal velocity profiles of the flows are revealed to be affected by the shock waves. For all the simulations, at the ebb tide, the drawdown of the water levels occurs in front of the open section, respectively, especially, hydraulic jump occurs when simulating the case of maximum difference in water level between the inside and outside of the seadike. As a result, it is thought that the flow characteristics of the sea region dominated by shock waves need to be identified employing three-dimensional analysis approach, which is expected to provide the information for ocean engineering works and facility management.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.559-573
• /
• 2008

In this study, we newly proposed 3-D nonlinear wave driver utilizing the Navier-Stokes Eq. the numerical integration of which is carried out using SPH (Smoothed Particle Hydrodynamics), an internal wave generation with the source function of Gaussian distribution and an energy absorbing layer. For the verification of new 3-D nonlinear wave driver, we numerically simulate the sloshing problem within a parabolic water basin triggered by a Gaussian hump and uniformly inclined water surface by Thacker (1981). It turns out that the qualitative behavior of sloshing caused by relaxing the external force which makes a free surface convex or uniformly inclined is successfully simulated even though phase error is visible and an inundation height shrinks as numerical simulation more proceeds. For the more severe test, we also simulate the nonlinear shoaling and refraction over uniform beach of wedge shape. It is shown that numerically simulated waves are less refracted than the linear counterpart by Hamiltonian ray theory due to nonlinearity, energy dissipation at the bottom and side walls, energy loss induced by breaking, and the hydraulic jump occurring when breaking waves encounter a down-rush by the preceding wave.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.697-707
• /
• 2008

The grid-based KIneMatic wave STOrm Runoff Model (KIMSTORM) by Kim (1998) predicts the temporal variation and spatial distribution of overland flow, subsurface flow and stream flow in a watershed. The model programmed with C++ language on Unix operating system adopts single flowpath algorithm for water balance simulation of flow at each grid element. In this study, we attempted to improve the model by converting the code into FORTRAN 90 on MS Windows operating system and named as ModKIMSTORM. The improved functions are the addition of GAML (Green-Ampt & Mein-Larson) infiltration model, control of paddy runoff rate by flow depth and Manning's roughness coefficient, addition of baseflow layer, treatment of both spatial and point rainfall data, development of the pre- and post-processor, and development of automatic model evaluation function using five evaluation criteria (Pearson's coefficient of determination, Nash and Sutcliffe model efficiency, the deviation of runoff volume, relative error of the peak runoff rate, and absolute error of the time to peak runoff). The modified model adopts Shell Sort algorithm to enhance the computational performance. Input data formats are accepted as raster and MS Excel, and model outputs viz. soil moisture, discharge, flow depth and velocity are generated as BSQ, ASCII grid, binary grid and raster formats.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3D
• /
• pp.423-429
• /
• 2008

The object of the present study is to estimate the potential effects of climate change and land use on soil erosion in the mid-east Korea. Simulated precipitation by CCCma climate model during 2030-2050 is used to model predicted soil erosion, and results are compared to observation. Simulation results allow relative comparison of the impact of climate change on soil erosion between current and predicted future condition. Expected land use changes driven by socio-economic change and plant growth driven by the increase of temperature and are taken into accounts in a comprehensive way. Mean precipitation increases by 17.7% (24.5%) for A2 (B2) during 2030-2050 compared to the observation period (1966-1998). In general predicted soil erosion for the B2 scenario is larger than that for the A2 scenario. Predicted soil erosion increases by 48%~90% under climate change except the scenario 1 and 2. Predicted soil erosion under the influence of temperature-induced fast plant growth, higher evapotranspiration rate, and fertilization effect (scenario 5 and 6) is approximately 25% less than that in the scenario 3 and 4. On the basis of the results it is said that precipitation and the corresponding soil erosion is likely to increase in the future and care needs to be taken in the study area.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6D
• /
• pp.809-817
• /
• 2008

This paper deals with the development of ViscoElastic Continuum Damage Finite Element Program (VECD-FEP++) and its verification with the results from both field and laboratory accelerated pavement tests. Damage characteristics of asphalt concrete mixture have been defined by Schapery's work potential theory, and uniaxial constant crosshead rate tests were carried out to be used for damage model implementation. VECD-FEP++ predictions were compared with strain responses (longitudinal and transverse strains) under moving wheel loads running at different constant speeds. To this end, an asphalt pavement section (A5) of Korea Expressway Corporation Test Road (KECTR) instrumented with strain gauges were loaded with a dump truck. Also, a series of accelerated pavement fatigue tests have been conducted at pavement sections surfaced with four asphalt concrete mixtures (Dense-graded, SBS, Terpolymer, CR-TB). Planar strain responses were in good agreement with field measurements at base layers, whereas strains at both surface and intermediate layers were found different from simulation results due to the complexity of tire-road contact pressures. Finally, fatigue characteristics of four asphalt mixtures were reasonably described with VECD-FEP++.