• 한국전산유체공학회지
• /
• 제19권3호
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2004년도 학술대회지
• /
• pp.309-313
• /
• 2004

Prediction of green water loads acting an the bow deck is an essential part for the design of bow structures against the green water impact. Proper technique of the green water simulation is highly required for the prediction of green water loads. In this paper, a new numerical method for green water simulation, which is based an predictor-corrector-upwind finite difference scheme of the 2nd kind, is introduced. Through the comparisons between computed' results and experimental measurements, it is verified that the present numerical tool is adequate as a practical calculation tool for the green water problem.

• 한국환경과학회지
• /
• 제13권7호
• /
• pp.655-660
• /
• 2004

This paper demonstrates the numerical simulation of three dimensional flow pattern for vehicular exhaust dispersion in the street canyons. The wind flow around buildings in urban is computed by the SIMPLEST method. The convection-diffusion equation was used to compute the $NO_X$ concentration level near buildings. Details are given of important boundary conditions and turbulence quantities variations. The simple turbulence model was used for unisotropic viscous effect. A control-volume based finite-difference method with the upwind scheme is employed for discretization equation. The simple turbulence model applied in this study has been verified through comparison between predicted and measured data near buildings. By the predictive results, the updraft induced by the presence of high-rise buildings is important in the transport of street level pollutant out from the street canyons. Our suggestion for reducing ground level pollution is to have high-rise buildings constructed or to reduce the channelling effect of street canyons.

• International Journal of Aeronautical and Space Sciences
• /
• 제3권1호
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• 한국수자원학회논문집
• /
• 제44권6호
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 에너지공학
• /
• 제5권1호
• /
• pp.50-55
• /
• 1996

상용 code인 Fluent를 이용하여 도자기의 1차 소성로에 대해 온도장과 속도장을 계산하고 열효율을 산출하였다. 수치해석의 변수로는 배출구 및 도자기의 위치를 변화시켰으며, 수치해석 방법은 검사체적에 기초한 유한차분방법 및 Upwind scheme과 SIMPLEC Algorithm을 사용하였고 난류모델로는 표준 k-$\varepsilon$ 모델을 사용하였다. 계산결과 출구 위치가 전체 소성로내 벡터유동의 양상을 크게 좌우하는 것으로 나타났으며, 전체 온도장에 대해서는 복사의 영향으로 큰 차이를 보이지 않았으나 예열대 상부에 출구가 있는 경우와 비교할 때 예열대 또는 냉각대의 측면에 출구가 있는 경우에 그 영역의 온도가 다소 높게 나타났다. 소성품의 위치는 로내 유동장 및 온도 분포에 크게 영향을 끼치지는 않으나 소성품 내 온도는 그 위치하는 영역의 온도 분포에 크게 영향을 받는 것으로 나타났다. 예열대 상부에 출구가 있는 경우 열효율은 17%로 매우 저조하였고 출구에서의 가스온도는 약 1000 K로 매우 높았다.

• 대한조선학회논문집
• /
• 제33권3호
• /
• pp.35-47
• /
• 1996

포텐셜을 기저로 하는 판요소법을 사용하여 자유 표면이 존재하는 유동장에서 일정 속도로 전진하는 3차원 물체의 형상을 설계하였다. 설계 방법으로는 원하는 압력 분포를 경계 조건으로 부여하고 이를 만족하는 물체 형상을 찾아내는 역해석법(inverse method)을 사용하였다. 즉, 주어진 압력으로부터 물체 표면에 분포된 법선 다이폴의 세기인 포텐셜 값을 결정하게 되며, 이는 물체 표면에 대한 Dirichlet형태의 경계 조건으로서 Green의 정리로부터 유도된 적분 방정식을 해석하게 된다. 전체 속도 포텐셜은 기본 유동인 선속에 대한 성분과 선제에 의하여 교란되는 성분으로 구성되어진다고 가정하였으며, 교란 포텐셜을 사용하여 선형화된 자유 표면 경계 조건을 적용하였다. 적분 방정식에 대한 수치 해석을 위해 물체 표면에 법선 다이폴과 Rankine 쏘오스를 분포하였으며, 자유 표면에는 Rankine 쏘오스를 분포하고 4점 유한 차분법을 사용하여 자유 표면 경계 조건이 만족되도록 하였다. 해로서 얻어지는 각 판요소에서의 Rankine 쏘오스의 세기는 가상의 유동 출입량으로서 형상 수정항으로 사용되었다. 몰수 회전 타원체의 형상 설계에 대하여 본 설계법을 적용한 결과 무한 수심에서나 조파 상태에서 $4{\sim}6$회의 반복 계산으로 충분히 수렴된 해를 얻을 수 있었다. 또한 자유 표면을 가르고 전진하는 Wigley 수학적 선형에 대한 형상 설계를 수행하여 만족스러운 결과를 얻어내었으며, 얻어진 수치해는 매우 안정적이고 빠른 수렴성을 보였다. 선형의 우열 비교를 통해 조파 저항을 감소시킬 수 있는 압력 분포의 형태를 파악하였으며, 이를 바탕으로 조파 저항의 관점에서의 5500TEU급 콘테이너 운반선의 설계를 수행하였다. 설계되어진 새로운 선형은 조파 저항의 관점에서 기존의 선형보다 계산과 실험에서 모두 우수하게 개량된 것으로 나타났다.