• Proceedings of the Korea Water Resources Association Conference
• /
• 1989.07a
• /
• pp.123-132
• /
• 1989

자연하천에서의 흐름특성을 해석하기 위해 Boundary - fitted(BF) coordinate system을 응용하여 유한차분법에 의한 2차원모델을 개발하였다. BF coordinates는 경계면의 형상에 관계없이 적용할 수 있으며 모든 계산은 기본식의 좌표변환을 통해 직각좌표계에서 행해지므로 경계조건의 입력에 용이하다. Physical domain(X - Y 좌표계)에서 하천의 형상을 입력하면 Grid generation에 의해 모든 계산은 Computational domain($\varepsilon$ - n 좌표계)에서 행해진다. Computational domain에서의 유한차분법은 half - time step으로 ADI 방법을 이용했고, 한 방향의 유속과 수위를 Double sweep으로 풀었다. 유속, 수위 및 하상의 격점망은 Staggered grid system을 사용했으며 geometric elements는 각 격점별로 계산하였다. 본 모델을 이용하므로써 불규칙한 수로나 하천의 흐름상태를 해석할 수 있으므로 흐름의 종단, 횡단방향의 유속분포, Superelevation을 구할 수 있고 하천의 계획, 관리, 제방의 호안이나 구조물의 설치등으로 일어나는 수리학적 영향등을 예측할 수 있다.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.73-82
• /
• 2011

In this paper, the vortex-induced vibration of circular cylinders is studied using the immersed boundary method on the Cartesian mesh. The Reynolds numbers considered is from 100 to 200. Using the configuration of tendemly arranged multiple circular cylinders, the vortex shedding behind of the cylinders and their flow-induced motion are investigated. The staggered MAC grid arrangement, which is the typical grid system for the incompressible flow on the Cartesian meshes, is utilized. Pressure correction method is applied for solving the divergence-free incompressible velocity field. The body motion is described by immersed boundary technique that has advantages for moving object on the fixed computational domain. It is also discussed for the computational noise in hydrodynamic forces when body motion is represented by the immersed boundary method. The Predictor/Corrector method is used for simulating the nonlinear response of the elastically mounted cylinder excited by vortex-shedding.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.243-248
• /
• 2007

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.103-106
• /
• 2002

This paper presents an effective numerical method for analyzing three-dimensional unsteady conjugate heat transfer problems of a curved pipe subjected to infernally thermal stratification. In the present numerical analyses, the thermally stratified flows in the pipe are simulated using the standard $k-{\varepsilon}$turbulent model and the unsteady conjugate heat transfer is treated numerically with a simple and convenient numerical technique. The unsteady conjugate heat transfer analysis method is implemented in a finite volume thermal-hydraulic computer code based on a non-staggered grid arrangement, SIMPLEC algorithm and higher-order bounded convection scheme. Numerical calculations have been performed far the two cases of thermally stratified pipe flows where the surging directions are opposite each other i.e. In-surge and out-surge. The results show that the present numerical analysis method is effective to solve the unsteady flow and conjugate heat transfer in a curved pipe subjected to infernally thermal stratification.

• Journal of Advanced Marine Engineering and Technology
• /
• v.22 no.6
• /
• pp.935-941
• /
• 1998

Experiments were carried out for a cubic cavity flow. Contrinuous shear stress is supplied by driven flow for high Reynolds number and three kinds of aspect ratios. Velocity vectors are obtained by PIV and they are used as velocity components for Poisson equation for pressure, Related boundary conditions and no-slip condition at solid wall and the linear velocity extrapolation on the upper side of cavity are well examined for the present study. For calculation of pressure resolution of grid is basically $40{\times}40$ and 2-dimensional uniform mesh using MSC staggered grid is adopted. The flow field within the cavity maintains a forced-vortex formation and almost of the shear stress from the driving inflow is transformed into rotating flow energy and the size of the distorted forced-vortex increases with increment of Reynolds number

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.2 no.1
• /
• pp.63-73
• /
• 1990

A numerical analysis of initial unsteady state flow and heat transfer in an enclosed cavity has been performed by the Modified QUICK Scheme. The stable QUICK Scheme which modified the coefficient always to be positive is included in this numerical analysis. The implicit method is applied to solve the unsteady state flow; between iterations the PISO (Pressure - Implicit with Splitting of Operators) algorithm is employed to correct and update the velocity and pressure fields on a staggered grid. The accuracy of the Modified QUICK Scheme is proved by applying fewer grid systems than those which Ghia et al. and Davis applied. The initial unsteady mixed convection in an enclosed cavity is analyzed using the above numerical procedure. This study focuses on the development of the large main vortex and secondary vortex in forced convection, the effects of the Rayleigh Number in natural convection and the relative direction of the forced and natural convection.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.11 no.6
• /
• pp.766-773
• /
• 1999

A numerical analysis has been performed solidification problem using the fixed grid-enthalpy method with enthalpy-porosity relation. A modified standard $k-\varepsilon$ model was applied to describe the influence of turbulent flow. Computational procedures are based on the finite volume method and the non-staggered grid system. Comparisons with the different three experimental results show that applying a modified standard $k-\varepsilon$model in mushyzone is better than the previous computation results. This paper includes another EMBR's influences such as change of velocity field, Increasement of temperature and dispersion of flow out of nozzle into the flow field. These EMBR's influences are compared to case without EMBR.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.229-232
• /
• 2008

The Central East Sea Tsunami caused huge damage to the eastern coast of the Korean Peninsula, eapecially Imwon port was damaged relatively strongly beacause of water depth variation which makes the wave concentration on this port. there are many virtual tsunami in east sea which has a possibility of happening. So, it is very important to expect the region which may be damaged by vritual tsunmis. In this study, modified dispersion-correction terms are used. The modified scheme has the advantage of using the constant spatial grid size and time step size even in real topography. Dynamic linking technique and staggered grid system are used. Using this model some cases of virtual tsunami was simulated and check the region which is occured maximum wave heights on the eastern sea of the Korean peninsula.

• 한국방재학회:학술대회논문집
• /
• 2008.02a
• /
• pp.431-434
• /
• 2008

Pratical dispersion-correction scheme is applicated to simulate the distant propagation of tsunami. This scheme is based on the leap-frog finite difference scheme for the linear shallow-water equations. The new scheme has the advantage of using the constant spatial grid size and time step size even in area of variable depths. And this new model constructed by using the 2nd upwind scheme, dynamic linking method, and staggered grid system. This model is simulated to near Sokcho harbor about The Central East Sea Tsunami in 1983. And this result is compared to tide gage and result of former model.

• Journal of Ship and Ocean Technology
• /
• v.5 no.1
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.