• Journal of Korean Society on Water Environment
• /
• v.23 no.1
• /
• pp.38-45
• /
• 2007

A one-dimensional kinematic wave model was used to calculate temporal and spatial changes of the highway runoff. Infiltration into pavement was considered using Darcy's law, as a function of flow depth and pavement hydraulic conductivity ($K_p$). The model equation was calculated using the method of characteristics (MOC), which provided stable solutions for the model equation. 22 storm events monitored in a highway runoff monitoring site in west Los Angeles in the U.S. were used for the model calculation and evaluation. Using three different values of $K_p$ ($5{\times}10^{-6}$, $10^{-5}$, and $2{\times}10^{-5}cm/sec$), total runoff volume and peak flow rate were calculated and then compared with the measured data for each storm event. According to the calculation results, $10^{-5}cm/sec$ was considered a site representative value of $K_p$. The study suggested a one-dimensional method to predict hydrodynamic behavior of highway runoff, which is required for the water quality prediction.

• Journal of Navigation and Port Research
• /
• v.27 no.6
• /
• pp.677-682
• /
• 2003

The long-period gravity wave, the wave period from some ten seconds to some minutes, induces not only the big sway of a ship moored and berthed in the harbor due to the horizontal long-distance motion of a water but also strong exfoliated flow and vortices near the harbor entrance. They muse serious problems on the safety navigation of vessels entering and leaving the harbor, but this gravity wave has not been searched sufficiently yet. Then it is quite important to reveal the characteristics of this long-period gravity wave ana to solve various problems induced by this wave. The long-period gravity wave measurement system with arrayed buoys installed the kinematic GPS was already proposed, which provides the precise propagating direction of the long-period gravity wave. In this paper, the observation results of the wave measurement system are shown by the MUSIC method And the propagating wave direction was estimated precisely enough in comparison with other results used other method.

• Journal of Korea Water Resources Association
• /
• v.52 no.11
• /
• pp.887-894
• /
• 2019

We proposed a GPU (Grapic Processing Unit) accelerated kinematic wave model for rainfall runoff simulation and tested the accuracy and speed up performance of the proposed model. The governing equations are the kinematic wave equation for surface flow and the Green-Ampt model for infiltration. The kinematic wave equations were discretized using a finite volume method and CUDA fortran was used to implement the rainfall runoff model. Several numerical tests were conducted. The computed results of the GPU accelerated kinematic wave model were compared with several measured and other numerical results and reasonable agreements were observed from the comparisons. The speed up performance of the GPU accelerated model increased as the number of grids increased, achieving a maximum speed up of approximately 450 times compared to a CPU (Central Processing Unit) version, at least for the tested computing resources.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.35 no.2
• /
• pp.41-48
• /
• 2023

In this study, simplified linear numerical method that can simulate wave generation and transformation by a moving bottom is introduced. Numerical analysis is conducted in wave number domain after continuity equation, linear dynamic and kinematic free surface boundary conditions and linear kinematic bottom boundary condition are Fourier transformed, and the results are expressed in space domain by an inverse Fourier transform. In the wavenumber domain, the dynamic free water surface boundary condition and the kinematic free water surface boundary condition are numerically calculated, and the velocity potential in the mean water level (z = 0) satisfies the continuity equation and the kinematic bottom boundary condition. Wave generation and transformation are investigated when the triangular and rectangular shape of bottoms move periodically. The results of the simplified numerical method are compared with the results of previous analytical solutions and agree well with them. Stability of numerical results according to the calculation time interval (Δt) and the calculation wave number interval (Δk) was also investigated. It was found that the numerical results were appropriate when Δt ≤ T(period)/1000 and Δk ≤ π/100.

• Journal of Korea Water Resources Association
• /
• v.32 no.2
• /
• pp.185-195
• /
• 1999

The finite difference method and the method of characteristics are frequently used for the numerical analysis of kinematic wave model. Truncation errors cause the peak discharge dissipated in the solution from the finite difference method. The peak discharge is conserved in the solution from the finite difference method. The peak discharge is conserved in the solution from the method of characteristics, however, the shock may deteriorates the numerical solution. In this paper, distinctive features of each scheme are investigated for the numerical analysis of kinematic wave model, and applicability of shock fitting algorithm such as Propagating Shock Fitting and Approximated Shock Fitting methods are studied. Propagating Shock Fitting method appears to treat shock properly, however, it failed to fit the shock appropriately when applied to a sudden inflow change in a long river. Approximate Shock Sitting method, which uses finer elements, is found to be more proper shock-fitting than the Propagating Shock Fitting method. Comparisons are made between two solution from the kinematic wave theory with shock fitting and full dynamic wave theory, and the results are discussed.

• Journal of Korea Water Resources Association
• /
• v.52 no.12
• /
• pp.1057-1066
• /
• 2019

The storage function method is used as a flood prediction model for four flood control offices in Korea as a method to analyze the actual rainfall-runoff relationship with non-linearity. It is essential to accurately estimate the parameters of the storage function method for accurate runoff analysis. However, the parameters of the storage function method currently in use are estimated by the empirical formula developed by the limited hydrological analysis in 2012; therefore, they are somewhat inaccurate. The kinematic wave method is a method based on physical variables of watershed and channel and is widely used for rainfall-runoff analysis. By adopting the two-term storage function method by the conversion of the kinematic wave method, parameters can be estimated based on physical variables, which can increase the accuracy of runoff calculation. In this research, the reproducibility of the kinematic wave method by the two-term storage function method was investigated. It is very easy to estimate the parameters because equivalent roughness, which is an important physical variable in watershed runoff, can be easily obtained by using land use and land cover, and the physical variable of channel runoff can be easily obtained from the basic river planning report or topographic map. In addition, this research examined the applicability of the two-term storage function method to runoff simulation of Naechon Stream, a tributary of the Hongcheon River in the Han River basin. As a result, it is considered that more accurate runoff calculation results could be obtained than the existing one-term storage function method. It is expected that the utilization of the storage function method can be increased because the parameters can be easily estimated using physical variables even in unmeasured watersheds and channels.

• Journal of Ocean Engineering and Technology
• /
• v.10 no.2
• /
• pp.136-145
• /
• 1996

The theory is based on two thermodynamic equations for the air mass in the air column and bydrodynamic equation for the relation between the response of the air in the water column and the incident wave. The numerical model is experimented in a two dimensional water tank and the caisson model with sloped front wall is tested in the large towing tank.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2000.04a
• /
• pp.140-145
• /
• 2000

The spectral method which is composed of an eigenfunction expansion of free modes in the wave number domain is used to produce two dimensional unsteady inviscid wave simulation such as progressive waves in a numerical pneumatic wave tank. A spatial and time dependent free surface elevation and the potential are calculated by integrating ODE derived from fully nonlinear kinematic and dynamic free surface boundary condition at each time step. The nonlinear characteristics in the waves by this method were notable as increasing wave steepness. This method is very useful and powerful in terms of saving computational time caused by rapid convergence exponentially with increasing number of nodes, even preserving accurate numerical results. Moreover, it will given us many possibilities to apply to naval and ocean engineering fields.

• Water for future
• /
• v.29 no.5
• /
• pp.139-150
• /
• 1996

This paper presents the effect of numerical parameters, such as grid size and grid ratio, on the outflow hydrograph of a unit-width plane in the linear Muskingum-Cunge method. The numerical results depend on Courant number C and cell Reynolds number D, two physically and numerically meaningful parameters. As C approache 1 and D increases, the numerical dispersion-relating oscillations are difficult to occur. The numerical oscillations occur in the front of a propagating wave for C < 1, while smaller oscillations occur behind the wave for C > 1 due to the numerical diffusion effect. For a plane with a small value of characteristic reach length L (e.g., a steep plane), the numerical solution of the Muskingum-Cunge method is similar to that of the kinematic wave method, which shows no wave attenuation. However, for a plane with a large value of L (e.g., a mild plane), the Muskingum-Cunge method leads to the diffusion waves which are essentially independent of the Courant number. Accordingly, the Muskingum-Cunge method will be suited for the routing of the catchment with relatively mild slopes.

• Journal of Mechanical Science and Technology
• /
• v.19 no.12
• /
• pp.2263-2269
• /
• 2005

The ship wave phenomena in the restricted waterway were investigated by a numerical analysis. The Euler and continuity equations were employed for the present study. The boundary fitted and moving grid system was adopted to enhance the computational efficiency. The convective terms in the governing equations and the kinematic free surface boundary condition were solved by the Constrained Interpolated Profile (CIP) algorithm in order to solve accurately wave heights in far field as well as near field. The advantage of the CIP method was verified by the comparison of the computed results by the CIP and the Maker and Cell (MAC) method. The free surface flow simulation around Wigley hull was performed and compared with the experiment for the sake of the validation of the numerical method. The present numerical scheme was applied to the free surface simulation for various canal sections in order to understand the effect of the sectional shape of waterways on the ship waves. The wave heights on the side wall and the shape of the wave patterns with their characteristics of flow are discussed.