• 대한토목학회논문집
• /
• 제28권5B호
• /
• pp.575-589
• /
• 2008

본 연구는 지진해일단파(tsunami bore) 혹은 조석단파(tidal bore)와 같은 단파의 동수역학적인 거동특성을 검토할 목적으로, 댐파괴류에서 단파의 형성과 동일한 방법, 즉 수조의 한쪽 끝단에 있는 고수위의 저수조(貯水槽) 게이트를 순간적으로 제거하는 방법으로 단파를 발생시킨다. 이러한 단파의 형성과 전파에 관한 수치시뮬레이션에 이상유(二相流)모델에 기초한 Navier-Stokes식을 적용하며, 이 때 비압축성 및 비혼합성의 액체와 기체흐름을 각각 고려한다. 기체와 액체의 접면을 VOF법으로 추적하고, Navier-Stokes방정식을 수치적으로 풀기 위하여 MCIP법을 적용한다. 1차원인 CIP법을 분할스텝기법을 사용하여 고차원으로 확장한 MCIP법은 수치확산이 매우 작고, 또한 안정된 스킴으로 알려져 있다. 게다가, 난류를 시뮬레이션하기 위하여 그의 유용성이 잘 알려져 있는 LES모델을 사용한다. 단파의 형성과 전파에 관한 수치해석결과를 검증하기 위하여 수리실험을 수행하였으며, 시간경과에 따른 수위변동과 평균유속변동에 대한 수치해석결과 및 실험결과를 비교하여 매우 양호한 상호대응관계를 확인할 수 있었다.

• 대한조선학회논문집
• /
• 제48권6호
• /
• pp.560-568
• /
• 2011

In SPH(Smoothed Particle Hydrodynamics) method, the fluid has been assumed that it is weakly compressible to solve the basic equations composed of Navier-Stokes equations and continuity equation. That leads to some drawbacks such as non-physical pressure fluctuations and a restriction as like small time steps in computation. In this study, to improve these problems we assume that the fluid is incompressible and the velocity-pressure coupling problem is solved by a projection method(that is, by ISPH method). The two-dimensional computation results of dam breaking and gravitational wave generation are respectively compared with the results of finite volume method and analytical method to confirm the accuracy of the present numerical computation technique. And, the agreements are comparatively acceptable. Subsequently, the green water simulations of a two-dimensional fixed barge are carried out to inspect the possibility of practical application to ship hydrodynamics, those correspond to one of the violent free surface motions with impact loads. The agreement between the experimental data and the present computational results is also comparatively good.

• 한국전산유체공학회지
• /
• 제13권2호
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• 대한토목학회논문집
• /
• 제31권3B호
• /
• pp.293-303
• /
• 2011

비압축성 점성 흐름을 수치해석하기 위한 효율적인 대각행렬화된 근사 인수분해(DAF) 알고리즘을 개발하였다. 압력에 근거한 인공압축성(AC) 기법을 이용하여 3차원 정상 비압축성 Navier-Stokes 방정식을 계산한다. AC 형태로 변형된 지배방정식은 2차 정확도의 유한차분법을 이용하여 공간에 대해서 이산화하였다. 이산화된 방정식계를 2차 정확도로 분할하기 위해서 본 연구에서 개발한 DAF 기법을 적용한다. 이 연구의 목적은 이 DAF 기법의 계산상 효율성을 검토하는 것이다. 만곡부를 갖는 사각형 덕트에서 완전히 발달한 층류 흐름과 발달하는 층류흐름 그리고 공동에서의 층류흐름에 대한 DAF 기법의 해석결과를 잘 알려진 4단계 Runge-Kutta(RK4)기법에 의한 해석해와 상대적으로 비교평가 하였다. 공간에 대해서 동일한 이산화기법을 이용하므로 동일한 격자상에서 계산된 DAF기법과 RK4기법의 해는 근본적으로 동일한 반면에, 이들 두기법의 계산상 효율성은 확연히 다른 것으로 나타났다. 본 연구에서 개발된 DAF기법은 적용한 모든 흐름 문제에 대해서 RK4기법에 비해 최소 2배 이상 적은 계산 시간만을 필요로 하는 것으로 나타났다. 이러한 DAF 기법의 계산상 효율성은 계산용량의 추가나 프로그래밍의 추가적인 복잡함이 없이 확보된다.

• 한국유체기계학회 논문집
• /
• 제1권1호
• /
• pp.72-80
• /
• 1998

Computational analysis of incompressible flows by numerically solving Navier-Stokes equations using multi-block finite volume method is conducted on a parallel computing system. Numerical algorithms adopted in this study $include^{(1)}$ QUICK upwinding scheme for convective $terms,^{(2)}$ central differencing for other terms $and^{(3)}$ the second-order Euler differencing for time-marching procedure. Structured grids are used on the body-fitted coordinate with multi-block concept which uses overlaid grids on the block-interfacing boundaries. Computational code is parallelized on the MPI environment. Numerical accuracy of the computational method is verified by solving a benchmark test case of the flow inside two-dimensional rectangular cavity. Computation in the axial compressor cascade is conducted by using 4 PE's md, as results, no numerical instabilities are observed and it is expected that the present computational method can be applied to the turbomachinery flow problems without major difficulties.

• 한국유체기계학회 논문집
• /
• 제2권2호
• /
• pp.57-63
• /
• 1999

The flow through multistage turbomachinery is affected by the interaction between a rotor and a stator. The interaction is due to the inviscid potential effect and viscous effect between closely spaced rotor and stator airfoils. Three-dimensional, unsteady, incompressible Navier-Stokes equations with a standard $k-{\epsilon}$ model are solved using a non-staggered grid system. This method is applied to the flow through a multistage compressor measured by Stauter et al. The results have shown strong interaction between the rotating and stationary flow field. The decay of rotor wake and the pressure profiles agree very well with experimental data. The wake produced by rotor causes unsteady pressure on the surface of a stator. The rotor/stator interaction produces the unsteady pressure force on the rotor and stator blades.

• 한국전산유체공학회지
• /
• 제9권3호
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• 대한기계학회논문집B
• /
• 제33권4호
• /
• pp.242-247
• /
• 2009

AILU type preconditioners for a two-dimensional combined P2P1 finite element formulation of the interaction of rigid cylinder with incompressible fluid flow have been devised and tested by solving fluid-structure interaction (FSI) problems. The FSI code simulating the interaction of a rigid cylinder with an unsteady flow is based on P2P1 mixed finite element formulation coupled with combined formulation. Four different preconditioners were devised for the two-dimensional combined P2P1 finite element formulation extending the idea of Nam et al., which was proposed for the preconditioning of a P2P1 mixed finite element formulation of the incompressible Navier-Stokes equations. It was found that PC-III or PC-IV among them perform well with respect to computational memory and convergence rate for some bench-mark problems.

• 한국유체기계학회 논문집
• /
• 제9권3호
• /
• pp.14-21
• /
• 2006

Leakage reduction through annular type seals of turbomachinery is necessary for enhancing their efficiency and the precise prediction method of seal leakage is needed. The analysis based on Bulk-flow concept has been mainly used in predicting seal leakage. However, full Navier-Stokes Equations with turbulent model derived in the seal flow passage have to be solved for improving the prediction of seal leakage. FLUENT 6 which is commercial CFD(Computational Fluid Dynamics) code based on FVM(Finite Volume Method) and SIMPLE algorism has been used to analyze leakage of various non-contact-type seals in this presentation. Comparing with the results of Bulk-flow model analysis and experiment, the result of CFD analysis shows good agreement with that of existing theoretical analysis for the incompressible grooved seal and compressive plain and staggered seal. The CFD analysis also shows improvement on the leakage prediction of the incompressible plain seal and compressive see-through-type labyrinth seal.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.160-166
• /
• 1998

In this paper, incompressible two-dimensional Navier-Stokes equations are numerically solved for the study of steady laminar flow around a body with the wall effect. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. To investigate the wall effect, numerical computations are carried out for the NACA 0012 section at the various blockage ratios. The pressure and skin friction on the foil surface, velocity pronto in its wake and drag coefficient are investigated as functions of the blockage ratio.