• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.6
• /
• pp.406-418
• /
• 2015

The present modified FUNWAVE-TVD model, which is a modification to its previous version 2.1, is applied to solitary wave propagation and is tested against the experiments of Vincent and Briggs(1989) and Luth et al.(1994). The eddy viscosity breaking scheme is used for comparison with the existing study in the case of breaking experiment. The symmetry of wave-induced current is maintained when the modified model is employed to Vincent and Briggs(1989) breaking experiment, but the symmetry of wave-induced current in previous model is not maintained. A better agreement with the breaking experimental data is obtained in the modified model using eddy viscosity breaking scheme than the shock capturing breaking scheme using nonlinear shallow water equation. For comparison with the schemes in the model, the fourth order MUSCL-TVD scheme by Erduran et al.(2005) and the third order MUSCL-TVD scheme using minmod limiter is applied, and the numerical solutions of solitary wave are compared.

• Atmosphere
• /
• v.33 no.4
• /
• pp.387-398
• /
• 2023

The double Fourier series (DFS) spectral dynamical core is evaluated for the two idealized test cases in comparison with the spherical harmonics (SPH) spectral dynamical core. A new approach in calculating the meridional expansion coefficients of DFS, which was recently developed to alleviate a computational error but only applied to the 2D spherical shallow water equation, is also tested. In the 3D deformational tracer transport test, the difference is not conspicuous between SPH and DFS simulations, with a slight outperformance of the new DFS approach in terms of undershooting problem. In the baroclinic wave development test, the DFS-simulated wave pattern is quantitatively similar to the SPH-simulated one at high resolutions, but with a substantially lower computational cost. The new DFS approach does not offer a salient advantage compared to the original DFS while computation cost slightly increases. This result suggests that the current DFS spectral method can be a practical and alternative dynamical core for high-resolution global modeling.

• Water for future
• /
• v.26 no.3
• /
• pp.87-101
• /
• 1993

A finite element model RIV-FEM2 for the hydrodynamic study in a river is developed based on two-dimensional shallow water wave equation and dissipative Galerkin's method. RIV-FEM2 consists of pre-processing, analysis processing and post-processing. Pre- and analysis processing is programmed with Fortran-77 and post-processing with turbo-Pascal respectively. The model is tested with two dimensional problems, including flow through bends, bridges, and symmetric contraction. The two dimensional tests shows stable and efficient results for various situations. Applicability of the model is verified by applying to natural river. The model will provide a basic contribution to the hydrodynamic analysis in a river.

• 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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.3
• /
• pp.176-183
• /
• 1991

In this paper. boundary element method is applied to the analysis of nonlinear free surface wave. A particular concern is given to the treatment of the open boundaries at the in-flow boundary and out-flow boundary, which uses the mass-flux and energy-flux considering the continuity of fluid. By assuming the fluid to be inviscid and incompressible and the flow to be irrotational. the problem is formulated mathematically as a two-dimentional nonlinear problem in terms of a velocity potential. The equation(Laplace equation) and the boundary conditions are transformed into two boundary integral equations. Due to the nonlinearity of the problem. the incremental method is used for the numerical analysis. Numerical results obtained by the present boundary element method are compared with those obtained by the finite element method and also with experimental values.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.847-857
• /
• 1994

FESWMS-2DH developed by the U.S. Department of Transportation based on two-dimensional shallow water wave equation is used in this study to simulate the flow characteristecs of the river reach between Chamsil and Shingok submerged weirs, which acts as a tidal river under low flow conditions. The model uses Galerkin F.E.M and meshes are composed of triangular or quadrangular elements. The model shows accurate and stable results concerning mass conservation as well as velocity distribution and water surface elevation. The results obtained in the present study may provide useful informations on the planning of river pollution abatement measures and artificial structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.5B
• /
• pp.429-439
• /
• 2009

When Catastrophic extreme flood occurs due to dam break, the response time for flood warning is much shorter than for natural floods. Numerical models can be powerful tools to predict behaviors in flood wave propagation and to provide the information about the flooded area, wave front arrival time and water depth and so on. But flood wave propagation due to dam break can be a process of difficult mathematical characterization since the flood wave includes discontinuous flow and dry bed propagation. Nevertheless, a lot of numerical models using finite volume method have been recently developed to simulate flood inundation due to dam break. As Finite volume methods are based on the integral form of the conservation equations, finite volume model can easily capture discontinuous flows and shock wave. In this study the numerical model using Riemann approximate solvers and finite volume method applied to the conservative form for two-dimensional shallow water equation was developed. The MUSCL scheme with surface gradient method for reconstruction of conservation variables in continuity and momentum equations is used in the predictor-corrector procedure and the scheme is second order accurate both in space and time. The developed finite volume model is applied to 2D partial dam break flows and dam break flows with triangular bump and validated by comparing numerical solution with laboratory measurements data and other researcher's data.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.11a
• /
• pp.220-221
• /
• 2005

Sub-bottom profiler(SBP) has been used extensively for the mapping of basement in the foundation design of offshore structure, for pre- and post-dredging operations within harbors and channels, for selection of pipeline routes, sitting of drilling platforms, and in the exploration for an aggregates such as sands and gravels. During the construction of Siwha embankment for irrigation water and the expansion of arable land, the breaking of an embankment unfortunately occurred so that a lot of riprap was swept away and widely dispersed by the tide and strong current. The feasibility study for the construction of the tidal-powered electric plant in Siwha embankment was performed quite recently. Therefore we made use of SBP survey to investigate the distribution of the lost riprap. We could successfully map out the distribution of the lost riprap from the reflection amplitude characteristics of the sediments in SBP data set. We demonstrated the variation of reflection amplitude versus the sediments with and/or without riprap by means of the numerical modeling of acoustic wave equation using finite difference method. Also we examined an amplitude anomaly of the ripraped area through the physical modeling using ultrasonic.

• Journal of Ocean Engineering and Technology
• /
• v.24 no.2
• /
• pp.10-17
• /
• 2010

Conventionally, because it is difficult to control a ship in shallow water and because attempting to do so creates unwanted environmental effects, maneuvering ships in the harbor area for berthing is usually done with the assistance of tugboats. In this paper, we propose a new method for berthing ships automatically by using bow and stern thrusters. Specifically, a steering motion model of a ship is considered, and parameters in the equation are evaluated by the system identification technique. An optimal controller based on observations was designed from the linearization of the non-linear ship motion in the horizontal plane. It is used to reduce the uncertainty about the ship's dynamics and reduce measurement requirements. The performance of the controller was also analyzed for its robustness relative to avoiding disturbing the environment due to winds, currents, and wave-drift forces. Experiments were conducted to estimate the potential for identifying result and the design of the controller. Specifically, in this paper, the system modeling and tracking control approach are discussed based on a two-degree-of-freedom (2DOF) servo-system design.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.3
• /
• pp.216-225
• /
• 1994

By use of laboratory experiment data set for an elliptic shoal by Berkhoff et al. (1982), both accuracy and Performance tests of numerical results between PCGM (Preconditioned Conjugate Gradient Method) and PA(Parabolic Approximation) are compared. Although both results show good agreement with the experimental data the PA model gives better reproduction of the relatively high amplitudes in the section 4-5 downwave of the shoal, in comparison with the PCGM. The PA model has been proved to be a useful tool for predicting wave transformationsin large shallow water region, but it can be applied only to the case of negligible reflection. On the other hand, there is a need to improve the computational efficiency of the PCGM model which is a finite difference scheme directly derived from the mild slope equation and can handle reflection. By taking the results of th PA model as an input data of the PCGM, the CPU time can be reduced by about 40%.