• Journal of Ocean Engineering and Technology
• /
• v.23 no.1
• /
• pp.16-23
• /
• 2009

The three-dimensional ship motion with forward speed was solved by a finite element method in the time domain. A boundary value problem was described in the frame of a fixed-body reference, and the problem was formulated according to Double-Body and Neumann-Kelvin linearizations. Laplace's equation with boundary conditions was solved by a classical finite element method based on the weak formulation. Chebyshev filtering was used to get rid of an unwanted saw-tooth wave and a wave damping zone was adopted to impose a numerical radiation condition. The time marching of the free surface was performed by the 4th order Adams-Bashforth-Moulton method. Wigley I and Wigely III models were considered for numerical validation. The hydrodynamic coefficients and wave exciting forces were validated by a comparison with experimental data and the numerical results of the Wigley I. The effects of the linearization are also discussed. The motion RAO was also checked with a Wigley III model through mono-chromatic and multi-chromatic regular waves.

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

This paper presents a simultaneous solution algorithm for one-dimensional unsteady flow routing through a dendritic channel system. This simulations solution algorithm is based on the double-sweep method and utilizes separate recursion equations for continuity, momentum and energy equations for each of the individual components of a dendritic channel system. Through separate recursion equations for each of the components. the new algorithm converts a dendritic channel network problem into a single-channel problem. The new algorithm is utilized in conjunction with a linearized unsteady flow model using full dynamic flow equations. The required computer storage for the coefficient matrix of the whole system is reduced significantly from the $2N{\times}2N$ matrix to a $2N{\times}4$ matrix, where N is the number of cross sections used in the computation of flow variables in a dendritic channel system. The algorithm presented in this paper provides an efficient and accurate modeling of unsteady flow events through a dendritic channel system.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.4
• /
• pp.48-59
• /
• 1996

The numerical step in the unsteady viscous flow analysis can be divided in the space analysis step satisfying continuity equation and the time marching step. In this study the spectral method is applied to solve the pressure Poisson equation in the space analysis step. If the highest order differential term of the pressure Poisson equation is transformed by Fourier series, pressure arid its first derivatives can be expressed by the integral form of Fourier series. So Gibb's phenomena can be eliminated and the spectral method can be applied to non-periodic problems. The numerical analysis of unsteady viscous flow around 2-dimensional circular cylinder and wing is carried out and compared for verification.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.39 no.1
• /
• pp.27-32
• /
• 2003

Under the assumption of potential flow, free-surface flow around a hydrofoil is calculated by the high-order spectra1!boundary-integral method, This method is one of the most efficient numerical methods by which the nonlinear interactions between hydrofoil and free-surface can be simulated in time-domain. In this method. the wave potential which represents the nonlinear evolution of free-surface is solved by the high-order spectral method and the body potential which provides the effects of hydrofoil and shed vortex is solved by the boundary-integral method. The calculated free-surface profiles which are generated by a uniformly translating hydrofoil are compared with other experimental results. And they show relatively good agreements each other. As another example, free-surface flow generated by a heaving and translating hydrofoil is calculated and discussed.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.2
• /
• pp.127-138
• /
• 1996

The strip method, unified theory and 3-D panel method are commonly used methods for the seakeeping analysis of high-speed vessels. The strip method which is basically 2-dimensional method is known to give incorrect hydrodynamic coefficients and motion responses for the cases of high speed and low frequency region. And the unified theory which uses two dimensional approach in inner domain and slender body theory in outer domain is very complicate in computational modelling. Though the 3-D panel method requires comparatively long computation time, it is believed that the method gives good results without any limitation in ship speed and range of frequency for computation. In the 3-D panel method the source singularity representing translating and pulsating Green function is used and Hoff's method is adopted for the numerical calculation of the Green function. The computation time can be reduced by using the symmetry relationship with respect to longitudinal axis. In this paper the strip method and the 3-D panel method are compared for the seakeeping analysis of a high-speed catamaran. The Compared items are the hydrodynamic coefficients, wave exciting forces, frequency response functions and short-term responses in irregular waves.

• Journal of Korea Water Resources Association
• /
• v.55 no.4
• /
• pp.257-266
• /
• 2022

Inundation damage is increasing every year due to localized heavy rain and an increase of rainfall exceeding the design frequency. Accordingly, the importance of hydraulic structures for flood control and defense is also increasing. The hydraulic structures are designed according to its purpose and performance, and the amount of flood is an important calculation factor. However, in Korea, design rainfall is used as input data for hydrological analysis for the design of hydraulic structures due to the lack of sufficient data and the lack of reliability of observation data. Accurate probability rainfall and its temporal distribution are important factors to estimate the design rainfall. In practice, the regression equation of temporal distribution for the design rainfall is calculated using the cumulative rainfall percentage of Huff's quartile method. In addition, the 6^th order polynomial regression equation which shows high overall accuracy, is uniformly used. In this study, the optimized regression equation of temporal distribution is derived using the variable selection method according to the principle of parsimony in statistical modeling. The derived regression equation of temporal distribution is verified through the significance test. As a result of this study, it is most appropriate to derive the regression equation of temporal distribution using the stepwise selection method, which has the advantages of both forward selection and backward elimination.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.7-12
• /
• 1998

Optimized high-order compact(OHOC) schemes were proposed, which have high spatial order of truncation and resolution to simulate the aeroacoustic problems due to unsteady compressible flows. Generally, numerical schemes are categorized explicit or implicit by time-marching method. In this research, OHOC differences which were developed with explicit time-marching method is used to have implicit formulation and the implicit OHOC differences result in block hepta-diagonal matrix. This paper presents the comparisons between the explicit and implicit OHOC schemes with a second order accuracy of time in the 1-d linear wave convection problem, and between the explicit OHOC scheme of 4th-order accuracy in time and the implicit OHOC scheme of 1st-order accuracy in tine for the 1-d nonlinear wave convection problem. With these comparisons, the characteristics of implicit OHOC scheme are shown in the point of CFL number.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.3
• /
• pp.229-235
• /
• 2019

In this study, a two-dimensional oscillating foil with forward speed in a propagating wave flow field was considered. The time-mean power to maintain the heaving and pitching motions of the foil was analyzed using the perturbation theory in an ideal fluid. The power, which was a non-linear quantity of the second-order, was expressed in terms of the quadratic transfer functions related to the mutual product of the heaving and pitching motions and incoming vertical flow. The effects of the pivot point and phase difference among the disturbances were studied. The negative power, which indicates energy extraction from the fluid, is shown as an example calculation.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.77-94
• /
• 1991

A finite difference method based on MAC method is used to simulate free-surface waves around a ship. Euler equations and continuity equation are differentiated using the forward time and central space, and solved by time marching scheme. By the employment of variable mesh system in horizontal and vertical direction, the numerical accuracy of wave simulation results is grossly improved. To verify the improvement of numerical accuracy, some numerical simulations are accomplished for Wigley, Series 60($C_{b}$=0.6) and a bulk carrier model. The computational results are compared to the various experimental data and show good agreements.