• Journal of Korea Water Resources Association
• /
• v.44 no.3
• /
• pp.249-262
• /
• 2011

In this study, in order to reduce the numerical oscillation due to the unbalance between source and flux terms as the HLLC scheme is applied to the flow analysis on the irregular bed topography, a unstructured finite volume model based on the well-balanced HLLC scheme and the shallow water equations is developed and applied to problems of dam-break waves. The well-balanced HLLC scheme considers directly the gradient of bed topography as the flux terms is calculated. This scheme provides the good numerical balance between the source and flux terms in the case of the application to the steady-state transcritical flow. To verify the numerical model developed in this study, it is applied to three cases of hydraulic model experiments and a field case study of Mapasset dam failure (France). As a result of the verification, the predicted numerical results agree relatively well with available laboratory and field measurements. The model provides slightly more accurate results compared with the existing models.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.15 no.3
• /
• pp.159-166
• /
• 2003

In this study, we develop a finite element model to directly solve the Laplace equation while keeping the same computational efficiency as the models based on the extended mild-slope equation which has been widely used for calculation of wave transformation in shallow water. For this, the computational domain is discretized into finite elements with a single layer in the vertical direction. The velocity potential in the element is then expressed in terms of the potentials at the nodes located at water surface, and the Galerkin method is used to construct the numerical model. A common shape function is adopted in horizontal direction, and the cosine hyperbolic function in vertical direction, which describes the vertical behavior of progressive waves. The model was developed for vertical two-dimensional problems. In order to verify the developed model, it is applied to vertical two-dimensional problems of wave reflection and transmission. It is shown that the present finite element model is comparable to the models based on extended mild-slope equations in both computational efficiency and accuracy.

• Water for future
• /
• v.18 no.1
• /
• pp.67-73
• /
• 1985

Tidal propagation in the Keum River has been routinely handled by numerical integration of the long fravity wave equation by Dronkers. The dynamic equations include non-linear terms thereby reproducing the shallow water tides. The model was used to compute tidal distribution of the Kum River for aveage spring, mean, neap tidal conditions and further utilised to investigate the waterlevel response within tidal reaches by combined tide and flood discharge effects. The objective of this initial study is to investigate the tidal dynamics of the lower reaches of the Keum River under the condition of before-cross-channel barrage construction.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.3
• /
• pp.619-633
• /
• 1994

A wave deformation model for general purpose combined refraction, diffraction, and breaking is developed in the shallow water. A parabolic approximation equation considered a higher order diffraction term is derived from the previous mild slope equation. A wave energy dissipation term due to bottom friction and breaking is introduced from the turbulence model. The Crank-Nicoloson implicit scheme is used in the numerical calculation, then the solutions are compared with the various hydraulic experiment data in the circular, the elliptic shoal, and the surf zone. The wave height decay in the surf zone is sensitively affected by the incident wave steepness, and the wave height variation around the elliptic shoal is well explained by the non-linear dispersion relation and the wave energy dissipation term. The model is also applied to a field coastal area and reasonable results are obtained.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.389-396
• /
• 1994

We analyze existing time-dependent mild-slope equations, which were developed by Smith and Sprinks (1975) (or, equivalently, Radder and Dingemans (1985)) and Kubo et al. (1992), in terms of the dispersion relation and energy transport. One-dimensionally in the horizontal direction, we compare the modulation of wave amplitudes for the time-dependent mild-slope equations against the linear Scrodinger equation. In view of the dispersion relation and modulation of wave amplitudes, Smith and Sprinks' model is more accurate in shallower water (kh$\leq$0.2$\pi$) and satisfies the linear Scrodinger equation in very shallow water (kh>0.2$\pi$) and satisfies the linear Scrodinger equation at a point of intermediate water depth (kh=0.3$\pi$). In view of the energy transport, Kubo et al.'s model is more accurate but yields singular solutions at some higher frequency range.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.197-203
• /
• 1998

An approach of decomposing the reflecting components is proposed by using the mild-slope equation of hyperbolic type which has the similar form to the shallow water equations. The approach is verified on Booij's problem and sinusoidally varying ripples. Inclusion of higher-order bottom effect given by chamberlain and Porter(1995) yields even more satisfactory results than the Berkhoff's mild-slope equation when compared with finite element solution or experiments.

• 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 Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.279-290
• /
• 1998

The height and speed of the shock wave are critical data in flood-control operations or in the design of channel walls and bridges along rivers with high flow velocities. Therefore, a numerical model is needed for simulating flow discontinuity over a wide range of conditions. In this study, a governing equation. As a Riemann solver Roe(1981)'s one is used. The model employs the modified MUSCL for handling the unstructured grids in this research. this model that adopts the explicit tradditional twl dimmensional dam break problems, two hydraulic dam break model is simulations, and a steady state simulation in a curved channel. Conclusions of this research are as follows : 1) the finite volume method can be combined with the Godonov-type method that is useful for modeling shocks. Hence, the finite volume method is suitable for modeling shocks. 2) The finite volume model combined with the modified MUSCL is successful in modeling shock. Therefore, modified MUSCL is proved to be valid.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• Water for future
• /
• v.29 no.5
• /
• pp.173-184
• /
• 1996

In this study, a two-dimensional shallow-water equation was used to develop the mathematical model for computing water levels and flow distribution. In the discretization equations, based on the finite volume method (FVM), the third order Runge-Kutta method and the third order upwind scheme were introduced to handle the unsteady and vconvective terms in the governing equations. To determine the accuracy of the developed model, it was applied to the rectangular horizontal channel in a frictionless flow. The water depth and velocity obtained by the numerical model were found to agree closely with the exact solution. The model was also applied to the rectangular channel with both the symmetric and the non symmetric constriction. The velocity distribution of the flow and the propagation of the flood wave were simulated and the results well described the flow characteristics.