• 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 the Korea Institute of Military Science and Technology
• /
• v.9 no.2 s.25
• /
• pp.152-160
• /
• 2006

Time marching methods are well-suited for high speed compressible flow computations. However, it is well known that the time marching methods suffer a slow down in convergence due to disparity in Eigenvalues. A local preconditioning method is one of numerical methods to enhance convergence characteristics of low mach number flows by modifying Eigenvalues of the governing equations. In this paper, the local preconditioning method of Weiss is applied to a 2 dimensional Navier-Stokes code and the efficiency of the preconditioning method is shown through a number of computational examples.

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

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.11a
• /
• pp.938-941
• /
• 2004

Sound generation and radiation from the duct-rotor system are calculated numerically. The wake geometries of a two-bladed rotor are calculated by using a time-marching fiee-wake method without a non-physical model of the far wake. Acoustic free field due to a rotating rotor is obtained by Lowson's equation. Using Kirchhoff source, rotating sources are modeled as stationary ones and can be inserted in the thin body boundary element method. The Kirchhoff source is validated through calculation of acoustic pressure due to a rotating point force. The thin body boundary element method (thin body BEM) is validated through calculation of acoustic radiation of ducted dipole. Using Kirchhoff source and thin body BEM, acoustic radiation of a ducted rotating source is calculated. Acoustic shielding is observed by inserting a duct and shows different phenomena at each major frequency. Acoustic radiation of a real duct-rotor system is also calculated using this method and the ducted acoustic field is significantly different from rotor only.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.16 no.5
• /
• pp.37-46
• /
• 2012

The calculation model for heat transfer analysis of rocket engine combustion chamber considering film-cooling has been established. Convective, radiative heat transfers and film-cooling effect in combustion chamber were evaluated using empirical equations especially for rocket engine combustors, and for heat transfer outward from chamber wall general convective and radiative equations were applied. Structural grid has been generated inside chamber wall for FVM calculations, and transient thermal analyses were carried out by time-marching techniques. LOx/kerosene rocket engine with chamber pressure of 50 bar has been analysed, and it is shown that, in that case, the film-cooling less than 4% remarkably contributes to reduce wall temperature, but the effect of the effect of film-cooling more than about 4% is not significantly increased.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2002.04a
• /
• pp.59-59
• /
• 2002

램제트 추진기관은 압축과정을 별도의 부품 없이 형상에 의해서 감속하여 연소 압력비를 얻는다. 따라서 구동 마하수와 형상에 의해 흡입과정의 압축 효율이 결정된다. 설계점은 충분한 유량을 확보 할 수 있는 유량과 충격파 각을 조절하여 전압력 손실을 줄이도록 고려되어야 한다. 또한 연소가 일어나면 연소실 압력이 배압으로 작용하고 비행시에 받음각은 변하므로 이에 따른 성능 분석도 고려 되어야 할 사항이다. 본 연구는 국내에서 실험한 형상에 대해 수치계산을 수행하여 코드의 검증과 아울러 램제트 유동장의 수치적 시뮬레이션도 설계단계에서 하나의 도구로 이용할 수 있음을 보여준다. 실험에서는 배압 조건을 얻기 위해 유동 블록키지를 유로 내에 두어 상응하는 배압을 얻었지만 본 계산에서는 압력 경계조건을 직접 사용하였다. 유동이 비정상 특성을 가지므로 시간 정확도를 이차로 가지도록 이중시간 전진법을 사용하였다. 사용한 압력비는 충격파가 카울 끝에 닿는 임계상태에 가까운 12, 13, 14에 대해 계산을 수행하였고 부스터모드로 흡입구 끝이 막혀 있다가 램제트 모드로 바뀌어 연소실 압력이 위의 압력비라고 가정할 때의 비정상 천이 과정을 계산해 보았다. 본 계산은 흡입구 부분만을 떼어놓고 적절한 가정 하에 수행되었지만 연소실 내부도 비정상 특성을 가지므로 흡입구와 연소실을 동시에 같이 계산해야한다. 추후에 전체적인 계산을 진행하기 위한 전 단계로 흡입구 계산만을 수행하여 실험과 잘 일치하는 계산 결과를 얻었고 전체 계산을 위한 연구는 진행 중에 있다.