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

This study is to develop an advanced storm runoff analysis program which takes geomorphological characteristics of watershed into consideration in determining model parameters. Basic concept of storm runoff modelling is based upon the kinematic wave theory. And numerical solution is obtained by the characteristic curve method. The storm runoff analysis program developed by this study is composed of multiple equivalent roughness sub-basins, each of which has two equivalent catchments on both side of a stream. Because it is based upon the stream-order of the Strahler system, the equivalent catchment-stream network reflects the stochastic geomorphological characteristics in the model parameter. Applicability and reliability of the storm runoff analysis program is evidenced by model calibration and verification process utilizing geographical and hydrological data of the Bocheong-river area which is a representative watershed of IHP projects in Korea. This study will hopefully contribute to hydrological calculation essentially required to understand water quality effect caused by regional development.

• Journal of Korea Water Resources Association
• /
• v.33 no.3
• /
• pp.315-330
• /
• 2000

The GIUH (Geomorphologic Instantaneous Unit Hydrograph) is to be applied to the ungauged or insufficiently gauged basins. For tIris purpose, an accurate estimation of the charactenstlc velocity is one very important part, but any proper method for this has not been developed yet. In case that we have enough rainfall and runoff clata, the estimation of the characteristic velocity may be an easy job, but it is out of the purpose of the GIUH. Remindmg that the purpose of GIUH the characterisbc veloclty should be estimated based on the geomorpholog1c analysis and also be snnple for easy apphcation. In tIris research analysis cmd application of the GruH was given to several sub-basins in Wi-stream river basin, Gono, Donggok and Hyoryung. After deriving the characteristic velocity througn a optimizatlOn process with real data, it is compared w1th several velOCIties der1ved from geOlnoI1Jhoclimatic instantaneous unit hydrograph theory and several other concentration time formulae. The estimated charactenstic velocities using Kerby, Kim, KInematic Wave, and Brasby- Williams formulae found to g1ve the appropriate results. Hmvever, as the Kerby, and the Kinematic Wave require user's decision of the IvIanning's n value, the K1m and the Braby-Williams seem to be more applicable and recommended as characteristic velocity formula.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.4
• /
• pp.1-7
• /
• 2003

A virtual reality technology for multipurpose numerical simulation is developed to reproduce and investigate a variety of ocean environmental problems in a 3D Numerical Wave Tank(NWT). The governing equations for solving incompressible fluid motion are Navier-Stokes equation and continuity equation. The Marker-Density function technique is adopted to implement the fully nonlinear freesurface kinematic condition. The marine environmental situations, i.e., waves, currents, etc., are reproduced by use of multi-segmented wavemakers on the basis of the so-called ″snake-principle″. In this paper, some numerical reproduction techniques for regular, and irregular waves, multi-directional waves, Bull's-eye wave. wave-current, and solitary wave are presented, and a model test in motion with large amplitude of roll angle is conducted in the developed 3D-NWT, using a overlaid grid system.

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

• Structural Engineering and Mechanics
• /
• v.71 no.1
• /
• pp.99-107
• /
• 2019

Wave propagation analysis of a porous graphene platelet reinforced (GPLR) nanocomposite shell is investigated for the first time. The homogenization of the utilized material is procured by extending the Halpin-Tsai relations for the porous nanocomposite. Both symmetric and asymmetric porosity distributions are regarded in this analysis. The equations of the shell's motion are derived according to Hamilton's principle coupled with the kinematic relations of the first-order shear deformation theory of the shells. The obtained governing equations are considered to be solved via an analytical solution which includes two longitudinal and circumferential wave numbers. The accuracy of the presented formulations is examined by comparing the results of this method with those reported by former authors. The simulations reveal a stiffness decrease in the cases which porosity influences are regarded. Also, one must pay attention to the effects of longitudinal wave number on the wave dispersion curves of the nanocomposite structure.

• Journal of Korea Water Resources Association
• /
• v.38 no.7 s.156
• /
• pp.555-564
• /
• 2005

In this study, a simple rainfall-runoff model was proposed by using the hydraulic routing model that requires relatively few parameters. The parameters of this model were estimated by the watershed characteristics data, and were applied to the Soyang watershed and Ui stream watershed by using the kinematic wave for overland flow and dynamic wave routing for channel routing. In order to demonstrate validity, the proposed approach was compared to the HEC-1 model for the Soyang watershed. As the results of modeling have shown, the hydraulic model shows reasonable results similar to that of the HEC-1 model. This model also represents good results for the Ui stream watershed. Hence, even if this model is a simple rainfall-runoff model using general methodology, it is competitive to the natural watershed. However, it is still difficult to estimate the roughness coefficient and the catchment width, and therefore this model is in need of such supplements.

• 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 the Korean Society of Fisheries and Ocean Technology
• /
• v.42 no.1
• /
• pp.38-43
• /
• 2006

The 3-D modeling of ocean finite element method(OFEM) using $k-{\varepsilon}$ turbulent model and tetrahedron grids has been used to investigate the internal wave generation during the expansion of the deep water from the open sea to the shelf with a simple shape, which can be widely used in the fields of submarine development, ocean environment and meteorology, etc. In this paper, the detailed configuration of internal wave with its length and height and also the distribution of salinity and turbulent kinematic energy, etc. were derived. It is hoped that this OFEM method can be successfully applied to the numerical calculation of internal wave for and the oceanographic problems (tidal flows around underwater hill, plateau, Georges Bank, etc.) and ocean engineering problems(flow past artificial sea reefs) in future.

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

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.81-86
• /
• 2001

Interfacial pressure jump terms based on the physics of phasic interface and bubble dynamics are introduced into the momentum equations of the two-fluid model for bubbly flow. The pressure discontinuity across the phasic interface due to the surface tension force is expressed as the function of fluid bulk moduli and bubble radius. The consequence is that we obtain from the system of equations the real eigenvalues representing the void-fraction propagation speed and the pressure wave speed in terms of the bubble diameter. Inversely, we obtain an analytic closure relation for the radius of bubbles in the bubbly flow by using the kinematic wave speed given empirically in the literature. It is remarkable to see that the present mechanistic model using this practical bubble radius can indeed represent both the mathematical well-posedness and the physical wave speeds in the bubbly flow.