• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.469-473
• /
• 2005

수질모형의 지배방정식에서 정의되는 대표적인 수질매개변수는 유역 및 대기로부터의 오염부하량 퇴적물로부터의 오염물질 용출부하량, 확산계수. 반응계수 등으로 직접적인 관측이 곤란할 뿐만 아니라 많은 관측비용을 필요로 한다. 본 연구에서는 매개변수를 포함한 오염물질 수지방정식을 구성하고, 구성된 선형 연립방정식을 이용함으로써 계산된 농도분포자료와 관측된 시계열 농도분포자료를 이용하여 계산한 질량변화량의 차이를 최소화하는 역산문제를 구성하여 모형의 매개변수를 추정하는 방법을 제시하였다. 이 방법을 이용하여 천수만, 울산만(울산항) 해역에서 관측된 연직방향 농도분포 자료를 이용하여 확산계수 및 대기로부터의 오염부하량, 퇴적물로부터의 오염물질 용출율, 확산$\cdot$반응에 의한 오염물질 변화량 등을 추정하였으며, 추정 매개변수는 시기적으로 변동이 크게 나타났다. 반면, 추정매개변수를 이용한 관측자료와 계산결과를 비교한 결과, RMS 오차는 관측자료 범위의 $5.0\% 이하, 일치지수는 0.95 이상으로 본 방법을 이용한 매개변수 추정결과의 신뢰성은 우수한 것으로 파악되었다.

• Journal of Korea Water Resources Association
• /
• v.36 no.3 s.134
• /
• pp.413-422
• /
• 2003

Numerical analysis for the Bragg reflection due to sinusoidally varying seabeds tying on a sloping beach was performed by using a couple of ordinary differential equations derived from the Boussinesq equations. Incident waves were wane groups generated by two short waves with slightly different phases. Effects of the slope of a seabed to the reflection were investigated in detail. It is shown that the reflection of long waves enhanced by increasing the slope of a seabed. This phenomenon caused by increase of wave amplitude due to increase of nonlinearity and shoaling.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.2
• /
• pp.77-89
• /
• 1991

Based on Boussinesq equation, the parabolic approximation equation is used to analyse the propagation of shallow water waves with currents over slowly varying depth. Rip currents (jet-like) occur mainly in shallow waters where the Ursell parameter significatly exceeds the range of application of Stokes wave theory. We employ the nonlinear parabolic approximation equation which is valid for waves of large Ursell parameters and small scale currents. Two types of currents are considered; relatively strong and relatively weak currents. The wave propagating over rip currents on a sloping bottom experiences a shoaling due to the variations of depth and current velocity as well as refraction and diffraction due to the vorticity of currents. Numerical analyses for a nonlinear theory are valid before the breaking point.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.47-58
• /
• 1994

Theoretical calculations are carried out for the estimation of linear maneuvering coefficients of a ship moving in shallow water region. Hydrodynamic forces and moments acting on a maneuvering ship are modelled based on a slender body theory, from which integro-differential equation for the unknown inner stream velocity is derived. Numerical algorithms fur solving this equation are described in detail. By considering water depth effects in the mathematical model, variations of maneuvering coefficients with water depth are studied. Programs are developed according to this method and calculations are done for Mariner, Series 60 and Wigley hull forms. For the verification of the programs, calculated results are compared with some analytic solutions and with published experimental results, which show good agreements in spite of many assumptions included in the mathematical model. It is expected that this method can be used as a preliminary tool for the estimation of maneuverability coefficients of a ship in shallow water region at its initial design stage.

• Journal of Ocean Engineering and Technology
• /
• v.22 no.5
• /
• pp.44-51
• /
• 2008

The planning and design of coastal structures against wave attack is required to accurately predict wave transformation, wave run-up, and fluid. particlevelocities an a slope. On tire other hand, in tire swash and surf zones of a natural beach, where coastal erosion and accretion occur at tire land-sea boundary, hydrodynamic analysis is essential. In this study, a RBREAK2 numerical model was created based on the nonlinear shallow water equation and laboratory measurements were carried out in terms of tire free surface elevations and velocities for tire cases of regular and irregular waves on 1 : 10 and 1 : 5 impermeable slopes. The data were used to evaluate tire applicability and limitations of tire RBREAK2 numerical model. The numerical mode1 could predict tire cross-shore variation of the wave profile reasonably well, but showed more accurate results for slopes that were steeper than 1 : 10. Except near tire wave crest, tire computed depth averaged velocities could represent tire measured profile below tire trough level fairly well.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.11 no.2
• /
• pp.95-106
• /
• 1999

A numerical model is represented to calculate the wave fields such as the reflected waves, the transmitted waves and the depth-averaged velocities over submerged breakwaters for the normally incident wave trains of nonlinear mono-chromatic wave and solitary wave. The finite amplitude shallow water equations with the effects of bottom friction are solved numerically in time domain using an explicit dissipative Lax-Wendroff finite difference method. The numerical model is verified by comparisons with the other numerical results and the measured data. It is found that the submerged breakwater may be more useful for protecting the energies of monochromatic waves rather than solitary waves. Finally, the armor stability on submerged breakwater is indirectly analyzed using the hydrodynamic characteristics of flow fields.

• Journal of Korea Water Resources Association
• /
• v.44 no.2
• /
• pp.145-156
• /
• 2011

Two dimensional numerical model of high-order accuracy is developed to analyze complex flow including transition flow, discontinuous flow, and wave propagation to dry bed emerging at natural river flow. The bed slope term of two dimensional shallow water equation consisting of integral conservation law is treated efficiently by applying quasi-steady wave propagation scheme. In order to apply Finite Volume Method using Fractional Step Method, MUSCL scheme is applied based on HLL Riemann solver, which is second-order accurate in time and space. The TVD method is applied to prevent numerical oscillations in the second-order accurate scheme. The developed model is verified by comparing observed data of two dimenstional levee breach experiment and dam breach experiment containing structure at lower section of channel. Also effect of the source term is verified by applying to dam breach experiment considering the adverse slope channel.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.1
• /
• pp.41-55
• /
• 2003

A dynamically-interfaced multi-grid finite difference model for simulation of tsunamis in the East Sea(Choi et al.) was established and further applied to produce detailed feature of coastal inundations along the whole eastern coast of Korea. The computational domain is composed of several sub-regions with different grid sizes connected in parallel of inclined directions with 16 innermost nested models. The innermost sub-region represents the coastal alignment reasonably well and has a grid size of about 30 meters. Numerical simulations have been performed in the framework of shallow-water equations(linear, as well as nonlinear) over the plane or spherical coordinate system, depending on the dimensions of the sub-region. Results of simulations show the general agreements with the observed data of run-up height for both tsunamis. The evolution of the distribution function of tsunami heights is studied numerically and it is shown that it tends to the log-normal curve for long distance from the source.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.1
• /
• pp.1-7
• /
• 2002

In the past, safety assessment on the site of Ulchin Nuclear Power Plants against tsunamis was carried out with probable maximum earthquake magnitude and related tsunamigenic fault parameters. Recently, however, based on the seismic gap theory, some seismologists warned about earthquakes of larger magnitudes than had been expected. In this study, we revaluated tsunami risk with a finite difference model based on linear and nonlinear shallow water equations. Firstly, we simulated the\`83 tsunami and compared the calculated water surface profile with the observed wave heights. Secondly, we evaluated the rise and drop of sea water level at the site of Ulchin Nuclear Power Plant with fault parameters of the past '83, '93 tsunamis and some dangerous faults. Finally, we showed that the cooling water intake facility of Ulchin Nuclear Power Plants would be safely operated in disastrous tsunamis.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.355-365
• /
• 1994

An analytical model of return flow is presented outside the surf zone. The governing equation is derived from the Navier-Stokes equation and the continuity. Each term of the governing equation is evaluated by the ordering analysis. Then the infinitesimal terms, i.e. the turbulent normal stress, the squared vertical velocity of water particle and the streaming velocity, are neglected. The driving forces of return flow are calculated using the linear wave theory for the shallow water approximation. Especially, the space derivative of local wave heights is described considering a shoaling coefficient. The vertical distribution of eddy viscosity is discussed to the customary types which are the constant, the linear function and the exponential function. Each coefficient of the eddy viscosities which sensitively affect the precision of solutions is uniquely decided from the additional boundary condition which the velocity becomes zero at the wave trough level. Also the boundary conditions at the bottom and the continuity relation are used in the integration of the governing equation. The theoretical solutions of present model are compared with the various experimental results. The solutions show a good agreement with the experimental results in the case of constant or exponential function type eddy viscosity.