• Magazine of the Korean Society of Agricultural Engineers
• /
• v.37 no.3_4
• /
• pp.34-47
• /
• 1995

It is experienced fact as a regular annual event that the structure to he designed on unreasonable flood for the agricultural structures including reservoirs have been brought not only loss of lives, but also enormous property damage. For the solution of this problem at issue, this study was conducted to develop an optimal runoff hydrograph model by comparison of the peak flows and time to peak between observed and simulated flows derived by linear time-invariant and linear time-variant models under the condition of having a short duration of heavy rainfall with uniform rainfall intensity at nine small watersheds which are within the range of 55.9 to 140.7 square kilometers in area in Han, Geum, Nagdong and Yeongsan Rivers. The results obtained through this study can be summarized as follows. 1. Storage constants and Gamma function arguments were calculated within the range of 1.2 to 6.42 and of 1.28 to 8.05 respectively by the moment method as the parameters for the analysis of runoff hydrograph based on linear time-invariant model. 2. Parameters for both linear time-invariant and linear time-variant models were calibrated with nine gaged watershed data, using a trial and error method. The resulting parameters including Gamma function argument, N and storage constant, K for linear time-invariant model were related statistically to watershed characteristic variables such as area, slope, length of main stream and the centroid length of the basin. 3. Average relative errors of the simulated peak discharge of calibrated runoff hydrographs by using linear time-variant and linear time-invariant models were shown to be 0.75 and 5.42 percent respectively to the peak of observed runoff hydrographs. Correlation coefficients for the statistical analysis in the same condition were shown to be 0.999 and 0.978 with a high significance respectively. Therefore, it can be concluded that the accuracy of a linear time-variant model is approaching more closely to the observed runoff hydrograph than that of a linear time-invariant model in the applied watersheds. 4. Average relative errors of the time to peak of calibrated runoff hydrographs by using linear time-variant and linear time-invariant models were shown to be 16.44 and 19.89 percent respectively to the time to peak of observed runoff hydrographs. Correlation coefficients in the same condition were also shown to be 0.999 and 0.886 with a high significance respectively. 5. It can be seen that the shape of simulated hydrograph based on a linear time- variant model is getting closer to the observed runoff hydrograph than that of a linear time-invariant model in the applied watersheds. 6. Two different models were verified with different rainfall-runoff events from data for the calibration by relative error and correlation analysis. Consequently, it can be generally concluded that verification results for the peak discharge and time to peak of simulated runoff hydrographs were in good agreement with those of calibrated runoff hydrographs.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.15 no.3
• /
• pp.1-9
• /
• 2011

The purpose of this study is to examine the dynamic characteristics of low, medium and high speed Maglev trains and guideways through dynamic interaction analysis. The coupled dynamic equations of motion for a vehicle of 10-dof and the associated guideway girders are developed by superposing vibration modes of the girder itself. The controller used in the UTM-01 Maglev vehicle is adopted to control the air gap between the bogie and guideway in this study. The effect of roughness, the guideway deflection-ratio and vehicle speed on the dynamic response of the maglev vehicle and guideway are then investigated using the 4th Runge-Kutta method. From the numerical simulation, it is found that the air gap increases with an increase of vehicle speed and the roughness condition. In particular, the dynamic magnification factor of the guideway girder is small at low and medium speeds, but the factor is noticeable at super-high speeds.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.6
• /
• pp.692-710
• /
• 2018

The Three-dimensional (3-D) dynamical behaviors of a fluid-conveying pipe subjected to vortex-induced vibration are investigated with different internal flow velocity ${\nu}$. The values of the internal flow velocity are considered in both subcritical and supercritical regimes. During the study, the 3-D nonlinear equations are discretized by the Galerkin method and solved by a fourth-order Runge-Kutta method. The results indicate that for a constant internal flow velocity ${\nu}$ in the subcritical regime, the peak Cross-flow (CF) amplitude increases firstly and then decrease accompanied by amplitude jumps with the increase of the external reduced velocity. While two response bands are observed in the In-line (IL) direction. For the dynamics in the lock-in condition, 3-D periodic, quasi-periodic and chaotic vibrations are observed. A variety of CF and IL responses can be detected for different modes with the increase of ${\nu}$. For the cases studied in the supercritical regime, the dynamics shows a great diversity with that in the subcritical regime. Various dynamical responses, which include 3-D periodic, quasi-periodic as well as chaotic motions, are found while both CF and IL responses are coupled while ${\nu}$ is beyond the critical value. Besides, the responses corresponding to different couples of ${\mu}_1$ and ${\mu}_2$ are obviously distinct from each other.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.406-418
• /
• 2020

In this study, a nonlinear wave simulation code is developed using a higher-order spectral (HOS) method. The HOS method is very efficient because it can determine the solution of the boundary value problem using fast Fourier transform (FFT) without matrix operation. Based on the HOS order, the vertical velocity of the free surface boundary was estimated and applied to the nonlinear free surface boundary condition. Time integration was carried out using the fourth order Runge-Kutta method, which is known to be stable for nonlinear free-surface problems. Numerical stability against the aliasing effect was guaranteed by using the zero-padding method. In addition to simulating the initial wave field distribution, a nonlinear adjusted region for wave generation and a damping region for wave absorption were introduced for wave generation simulation. To validate the developed simulation code, the adjusted simulation was carried out and its results were compared to the eighth order Stokes theory. Long-time simulations were carried out on the irregular wave field distribution, and nonlinear wave propagation characteristics were observed from the results of the simulations. Nonlinear adjusted and damping regions were introduced to implement a numerical wave tank that successfully generated nonlinear regular waves. According to the variation in the mean wave steepness, irregular wave simulations were carried out in the numerical wave tank. The simulation results indicated an increase in the nonlinear interaction between the wave components, which was numerically verified as the mean wave steepness. The results of this study demonstrate that the HOS method is an accurate and efficient method for predicting the nonlinear interaction between waves, which increases with wave steepness.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).

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

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.8
• /
• pp.531-540
• /
• 2022

In this paper, we present an algorithm that identifies artillery type and rapidly predicts the impact point based on the IMM filter. The ballistic trajectory equation is used as a system model, and three models with different ballistic coefficient values are used. Acceleration was divided into three components of gravity, air resistance, and lift. And lift acceleration was added as a new state variable. The kinematic condition that the velocity vector and lift acceleration are perpendicular was used as a pseudo-measurement value. The impact point was predicted based on the state variable estimated through the IMM filter and the ballistic coefficient of the model with the highest mode probability. Instead of the commonly used Runge-Kutta numerical integration for impact point prediction, a semi-analytic method was used to predict impact point with a small amount of calculation. Finally, a state variable initialization method using the least-square method was proposed. An integrated algorithm including artillery type identification, impact point prediction and initialization was presented, and the validity of the proposed method was verified through simulation.