• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.11a
• /
• pp.177-184
• /
• 2003

Numerical simulation on marine wind and sea surface elevation was carried out using both three-dimensional hydrostatic and non-hydrostatic models and a simple oceanic model from 0900 LST, August 13 to 0900 LST, August 15, 1995. As daytime easterly meso-scale sea-breeze from the eastern sea penetrates Kangnung city in the center part as basin and goes up along the slope of Mt. Taegullyang in the west, it confronts synoptic-scale westerly wind blowing over the top of the mountain at the mid of the eastern slope and then the resultant wind produces an upper level westerly return flow toward the East Sea. In a narrow band of weak surface wind within 10km of the coastal sea, wind stress is generally small, less than l${\times}$10E-2 Pa and it reaches 2 ${\times}$ 10E-2 Pa to the 35 km. Positive wind stress curl of 15 $\times$ 10E-5Pa $m^{-1}$ still exists in the same band and corresponds to the ascent of 70 em from the sea level. This is due to the generation of northerly wind driven current with a speed of 11 m $S^{-1}$ along the coast under the influence of south-easterly wind and makes an intrusion of warm waters from the southern sea into the northern coast, such as the East Korea Warm Current. On the other hand, even if nighttime downslope windstorm of 14m/s associated with both mountain wind and land-breeze produces the development of internal gravity waves with a hydraulic jump motion of air near the coastal inland surface, the surface wind in the coastal sea is relatively moderate south-westerly wind, resulting in moderate wind stress. Negative wind stress curl in the coast causes the subsidence of the sea surface of 15 em along the coast and south-westerly coastal surface wind drives alongshore south-easterly wind driven current, opposite to the daytime one. Then, it causes the intrusion of cold waters like the North Korea Cold Current in the northern coastal sea into the narrow band of the southern coastal sea. However, the band of positive wind stress curl at the distance of 30km away from the coast toward further offshore area can also cause the uprising of sea waters and the intrusion of warm waters from the southern sea toward the northern sea (northerly wind driven current), resulting in a counter-clockwise wind driven current. These clockwise and counter-clockwise currents much induce the formation of low clouds containing fog and drizzle in the coastal region.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.4B
• /
• pp.343-352
• /
• 2011

The water quantity by intake station as well as the tributary flow discharge acting as sink or source were added to the main flow rate in the present study and RMA-2 and RMA-4 models were applied to the reach from Pal-dang dam to Jam-sil submerged weir to investigate the influence of water supply withdrawal on the river flow and water quality. The numerical results revealed that the water supply withdrawal from 5 intake stations located upstream of Jam-sil submerged weir changed the total flow rate and therby induced different hydraulic characteristics in terms of water surface elevation and velocity. The changed flow field by the inclusion of water intake quantity led to the variation of water quality. By the consideration of the water supply withdrawal, the velocity structure was significantly disturbed by the outflowing flow condition nearby Gu-ui, Ja-yang, and Pung-nap intake stations. Furthermore, the mean velocity was lowered by 25% and the stage upstream of Gu-ui station rose upto 1.5 cm compared with the result by exclusion of water intake. In case of no water withdrawal, the distribution of BOD concentration was parallel throughout the domain. However, when the water withdrawal is considered, the distribution of BOD concentration nearby the Gu-ui, Am-sa, and Ja-yang station was signifiantly changed. In addition, the BOD concentration including the intake stations showed higher value at the downstream of the reach due to the loss of the discharge by water withdrawal effect. It is concluded that both the inflow and outflow discharges from tributaries and water intake stations should be included in the numerical simulation to analyze the hydrodynamic behaviors and mixing characteristics more accurately.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.11
• /
• pp.69-78
• /
• 2012

Conventional numerical analysis for rainfall-induced slope stability has been estimated by separating seepage and stress-strain behavior, respectively. Many researchers' models from commercial softwares and literatures define that partially saturated permeability is the only function of degree of saturation (or matrix suction) and then they do not consider hydraulic-mechanical characteristics for the analysis. However, in practice, the water flow processes in a deformable soil are influenced by soil skeleton movement and the pore water pressure changed due to seepage will lead to changes in stress and to deformation of a soil. The relationship between seepage and soil behavior causes a change of partially saturated permeability as well as saturated permeability with the lapse of time. Instability of partially saturated soil slopes due to infiltration would be analyzed from reduction of negative pore water pressure calculating the process of water flow based on predicted partially saturated permeability. Therefore, partially saturated permeability should be defined by the function of degree of saturation (or matric suction) and porosity. The paper presents the comparison between staggered and monolithic coupled analysis regarding seepage and stress deformation problems. As a result, the decrease in matric suction on soil slope from monolithic analysis is slower than that from staggered analysis.

• Journal of the Korean GEO-environmental Society
• /
• v.11 no.5
• /
• pp.35-43
• /
• 2010

Partially saturated permeability should be defined by the function of suction (or degree of saturation) and porosity. However, commercial software and most researchers' model often describe as the function of suction. The stability of a soil slope can be affected by both hydraulic and shear strength properties of partially saturated soils. For both studies, we generally use an uncoupled seepage analysis program Seep/W(Geo-Slope, 2007) and a series stress-deformation analysis program Sigma/W, or slope stability analysis program Slope/W. Seep/W is performed for simulations of partially saturated flow problems in non-deformable soil media. However, under real situations, the water flow processes in a deformable soil are influenced by soil skeleton movement and the pore water pressure changed due to seepage will lead to changes in stresses and to deformation of a soil. Many researchers are currently developing their models for solving coupled hydro-mechanical problems to simulate slope stability during a rainstorm. For a proper implementation in the field, the developed model should be still needed in order to achieve appropriate accuracy of the solution for coupled hydro-mechanical problems in soil slope stability. Thus, the paper presents the comparison of slope stability between uncoupled and coupled analyses of seepage and stress deformation problems.

• The Journal of the Acoustical Society of Korea
• /
• v.37 no.3
• /
• pp.147-155
• /
• 2018

In this paper, we tried to find out sound quality metrics to express discomfort of overload excavator noise and to develop sound quality indexes through multiple regression analysis by using them. For this purpose, the interior noise of cabin under overload condition was recorded for six excavator models with different noise properties and Jury test was carried out by PCM (Paired Comparison Method) and MEM (Magnitude Estimation Method). Jury test result with low consistency was classified into two groups with different preference tendencies by cluster analysis and multiple regression analysis was conducted in order to find out which sound quality metrics have significant effects on discomfort(low preference). As a result, we figured out that the sound quality metrics to express the discomfort were the partial loudness (= $PN_{10Bark}$) between 0 and 10 Bark in case of group1 and the difference between engine noise(= $dB_{EG}$) and hydraulic system noise ($dB_1$) in case of group2. Using the results of preference ranking and tendency analysis of PCM followed by the correlation analysis between PCM and MEM, the more reliable results were adopted by excluding the data with low consistency obtained from Jury test via MEM.

• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.138-146
• /
• 2007

This paper describes 2.5D induced polarization (IP) modeling and inversion algorithms using complex resistivity. The complex resistivity method has merits for acquiring more valuable information about hydraulic parameters and pore fluid than the conventional IP methods. The IP modeling and inversion algorithms are developed by allowing complex arithmetic in existing DC modeling and inversion algorithms. The IP modeling and inversion algorithms use a 2.5D DC finite-element algorithm and a damped least-squares method with smoothness constraints, respectively. The accuracy of the IP modeling algorithm is verified by comparing its responses of two synthetic models with two different approaches: linear filtering for a three-layer model and an integral equation method for a 3D model. Results from these methods are well matched to each other. The inversion algorithm is validated by a synthetic example which has two anomalous bodies, one is more conductive but non-polarizable than the background, and the other is polarizable but has the same resistivity as the background. From the inverted section, we can cleary identify each anomalous body with different locations. Furthermore, in order to verify its efficiency to the real filed example, we apply the inversion algorithm to another three-layer model which includes phase anomaly in the second layer.

• Water for future
• /
• v.14 no.4
• /
• pp.13-18
• /
• 1981

The rainfall pattern analysis on time distribution characteristics of rainfall rates in important in determination of design flow for hydraulic structures, particularly in urban area drainage network system design. The historical data from about 400 storm samples during 31 years in Seoul have been used to investigate the time distribution of 5-minute rainfall in the warm season. Time distribution relations have been deveolped for heavy stroms over 20mm in total rainfall and represented by relation percentage of total storm rainfall to percentage of total storm time and grouping the data according to the quartile in which rainfall was heaviest. And also time distribution presented in probability terms to provide quantitative information on inter-strom variability. The resulted time distribution relations are applicable to construction of rainfall hyetograph of design storm for determination of design flow hydrograph and identification of rainfall pattern at given watershed area. They can be used in conjuction with informations on spatstorm models for hydrologic applications. It was found that second-quartile storms occurred most frequently and fourth-quartile storms most infrequently. The time distribution characteristics resulted in this study have been presented in graphic forms such as time distribution curves with probability in cumulative percent of storm-time and precipitation, and selected histograms for first, second, third, and fourth quartile stroms.

• Journal of Korea Water Resources Association
• /
• v.53 no.11
• /
• pp.985-997
• /
• 2020

SWAT (Soil and Water Assessment Tool)-MODFLOW (Modular Groundwater Flow) is a coupled model that linking semi-distributed watershed hydrology with fully-distributed groundwater behavior. In this study, the groundwater simulation results of SWAT and SWAT-MODFLOW were compared for Bokhacheon watershed in Namhan river basin. The models were calibrated and validated with 9 years (2009~2017) daily streamflow (Q) data of Heungcheon (HC) water level gauge station and the daily groundwater level observation data of Yulheon (YH). For SWAT, the groundwater parameters of GW_DELAY, GWQMN, and ALPHA_BF affecting baseflow and recession phase were treated. The SWAT results showed the coefficient of determination (R²) of 0.7 and Nash-Sutcliffe model efficiencies (NES_Q, NSE_inQ) for Q and 1/Q with 0.73 and -0.1 respectively. For SWAT-MODFLOW, the spatio-temporal aquifer hydraulic conductivity (K, m/day), specific storage (S_s, 1/m), and specific yield (S_y) were applied. The SWAT-MODFLOW showed R², NSE_Q, and NSE_inQ of 0.69, 0.74, and 0.51 respectively. The SWAT-MODFLOW considerably enhanced the low flow simulation with the help of aquifer physical information. The total streamflow of SWAT and SWAT-MODFLOW were 718.6 mm and 854.9 mm occupying baseflow of 342.9 mm and 423.5 mm respectively.

• Nuclear Engineering and Technology
• /
• v.27 no.2
• /
• pp.176-188
• /
• 1995

A particle tracking scheme was developed in order to model radionuclide transport through a tortuous flow Held in a rock fracture. The particle tacking method may be used effectively in a heterogeneous flow field such as rock fracture. The parallel plate representation of the single fracture fails to recognize the spatial heterogeneity in the fracture aperture and thus seems inadequate in describing fluid movement through a real fracture. The heterogeneous flow field une modeled by a variable aperture channel model after characterizing aperture distribution by a hydraulic test. To support the validation of radionuclide transport models, a radionuclide migration experiment was performed in a natural fracture of granite. $^3$$H_2O$ and $^{131}$ I are used as tracers. Simulated results were in agreement with experimental result and therefore support the validity of the transport model. Residence time distributions display multipeak curves caused by the fast arrival of solutes traveling along preferential fracture channels and by the much slower arrival of solutes following tortous routes through the fracture. Results from the modelling of the transport of nonsorbing tracer through the fracture show that diffusion into the interconnected pore space in the rock mass has a significant effect on retardation.