• 한국전산유체공학회지
• /
• 제12권2호
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• 한국에너지공학회:학술대회논문집
• /
• 한국에너지공학회 1999년도 춘계 학술발표회 논문집
• /
• pp.127-132
• /
• 1999

우리나라의 도시폐기물 소각로는 화격자 위에 폐기물을 공급하고 화격자 밑에서 공기를 공급하는 스토커식 소각로를 대부분 채택하고 있다. 이러한 스토커 소각로 연소실내에서는 매우 복잡한 연소현상이 발생하는데, 연소실로 투입된 쓰레기는 먼저 건조부에서 수분의 건조가 일어나고, 화격자의 구동에 의해 쓰레기가 혼합 및 이송되면서 열분해, 가스화, 가연성분의 탈휘발화 및 연소, 일부 고정탄소의 표면연소 등의 반응이 일어난다. 그리고 1,2차 연소실에서는 휘발분 및 비산된 고체의 연소가 일어나는데, 이때 대류 및 복사열전달 등의 복잡한 현상을 수반하는 유동장이 형성된다. 더욱이 불균질한 특성을 갖는 쓰레기층 내에서의 복잡한 현상으로 인하여 발생하는 경계조건 설정의 불확실성으로 연소실내의 연소 현상을 전산해석하는 데에는 상당한 어려움이 있다.(중략)

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.27-33
• /
• 2009

In this paper, the flow past a rotary oscillating circular cylinder is simulated. The high-order and high-resolution numerical schemes with the characteristic boundary conditions are used for the compressible Navier-Stokes equation. The frequencies of rotating oscillation are $0.19\;{\leq}\;S_f\;{\leq}\;0.25$ for the maximum angular $\theta_{max}=10^{\circ}$ and $17^{\circ}$. The flow conditions are Mach number of 0.3 and Reynolds number of 1000. At Lock-on and Non-lock-on region which are defined by the relation between the vortex shedding frequency and the oscillating frequency, the drag and lift coefficient are analyzed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.51-59
• /
• 2001

An evaluation of one zero-equation and two one-equation eddy viscosity-transport turbulence closure models as implemented CFDS(Characteristic Flux Difference Splitting ) code is presented herein. Comparisons of Baldwin-Lomax model as zero-equation and Baldwin-Barth and Spalart-Allmaras model as one-equation are presented for three test cases, first inlvolving the 3 dimensional supersonic flow at M=1.98 over tangent ogive cylinder, second involving the 2 dimensional transonic flow at M=0.79 over RAE 2822 airfoil, third involving the 3 dimensional transonic flow at M=0.84 over ONERA M6 wing. The numerical results of CFDS code will also examined through direct comparison with experimental data.

• 해양환경안전학회:학술대회논문집
• /
• 해양환경안전학회 2018년도 공동국제학술발표회
• /
• pp.161-161
• /
• 2018

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2010년 춘계학술대회논문집
• /
• pp.380-388
• /
• 2010

A thermal-hydraulic code, named CUPID, has been developed for the analysis of transient two-phase flows in nuclear reactor components. A two-fluid three-field model was used for steam-water two-phase flows. To obtain numerical solutions, the finite volume method was applied over unstructured cell-centered meshes. In steam-water two-phase flows, a phase change, i.e., evaporation of condensation, results in a great change in the flow field because of substantial density difference between liquid and vapor phases. Thus, two-phase flows are very sensitive to the local pressure that determines the phase change. This in turn puts emphasis on the accurate evaluation of local pressure gradient. This paper presents a new numerical scheme to evaluate the pressure gradient at cell centers on unstructured meshes. The results of the new scheme for a simple test function a gravity-driven cavity, and a wall boiling two-phase flow are compared with those of the previous schemes in the cupid code.

• 한국전산유체공학회지
• /
• 제14권3호
• /
• pp.86-94
• /
• 2009

The non-conservative form of momentum equations is often used for some two-phase flow codes instead of a conservative form because of numerical convenience. Another non-conservative form, so called, a semi-conservative form can improve the numerical solution of these codes maintaining the numerical convenience. It is close to the conservative form but still maintains the feature of the non-conservative form. A semi-conservative form of the momentum equations and a non-conservative form of the momentum equations are implemented in CUPID[1] code. The numerical results of the semi-conservative and the non-conservative forms are compared against analytical solutions and the solutions of the FLUENT code that uses the conservative form. The results clearly showed that the semi-conservative form of the momentum equations provides better solutions than the non-conservative form, especially for heterogeneous two-phase flows.

• 한국해양공학회지
• /
• 제26권1호
• /
• pp.9-16
• /
• 2012

This paper presents the numerical results of a simulation of the free surface flow around a blunt bow ship model and focuses on the validation of the proposed method with a brief investigation of the relation between the resistance and free surface behavior. A finite volume method based on the Reynolds Averaged Navier-Stokes (RANS) approach is used to solve the governing flow equations, where the free surface, including wave breaking,is captured by using a two-phase Level-Set (LS) method. For turbulence closure, a two equation k-${\varepsilon}$ model with the standard wall function technique is used. Finally, the numerical results are compared with the available experimental data, showing good agreement.

• 한국에너지공학회:학술대회논문집
• /
• 한국에너지공학회 2008년도 춘계학술 발표회
• /
• pp.179-186
• /
• 2008

현재 대부분의 산업용 열기관은 효율을 높이기 위하여 연소에 사용되는 공기를 예열하는 방법을 사용한다. 하지만, 산업용 열기관에서 평균적으로 발생되는 $1000^{\circ}C$ 이상의 배기가스는 일반 금속 열교환기에는 적합하지 않다. 이에 반해 세라믹 열교환기의 경우 고온에서 견디는 장점이 있다. 본 연구에서는 기본적인 열교환기 설계 이론을 이용하여 설계프로그램을 제작하였다. 또한 세라믹 열교환기 내 열 유체 거동을 CFD 상용코드인 FlUENT 6.2를 이용한 전산해석을 수행하여 설계결과를 비교 검증하였다. 설계 결과에서 휜의 형태 변화에 따라 열전달율과 온도구배는 무시할 수 있을 정도로 작았으나, 압력강하는 크게 변동되는 결과가 도출되었다. 제한된 모듈 크기에서 휜 간거리는 휜의 두께에 비해 약 3배 이상 클 경우가 적당하며, 판(plate)의 두께는 작을수록 압력손실이 적고, 열전달율이 상승하지만 두께가 너무 얇게 된다면 제작상의 어려움이 생긴다. 향후 연구에서는 단순한 구조에서 벗어나 off-set이나 판형구조를 고려하여 설계함으로서 열전달 면적을 넓히거나 난류유동을 발생시켜 열전달율을 높이는 연구를 진행 할 필요가 있다.

• 한국음향학회지
• /
• 제37권2호
• /
• pp.92-98
• /
• 2018

본 논문은 모형 프로펠러를 대상으로 공동수조 시험, 수중 충격시험, 유한요소해석 및 전산유체해석에 기반하여 수행한 명음 발생 메커니즘 연구이다. 선미 유동을 모사하기 위해 반류망, 프로펠러 및 방향타를 설치하고 수중청음기와 가속도계로 프로펠러 명음 현상의 발생과 소멸을 계측하였다. 유한요소해석을 통해 프로펠러 날개의 고유진동수를 예측하고 접촉 및 비접촉식 충격시험으로 이를 검증하였다. RANS(Reynolds Averaged Navier-Stokes) 방정식 기반 전산유체해석을 통하여 프로펠러 날개 각 단면의 유속과 유효 받음각을 계산하였으며, DES(Detached Eddy Simulation) 기반 고해상도 해석을 통해 명음 발생 위치에서 2-D 날개 단면 뒷전의 와류흘림주파수(vortex shedding frequency) 계산을 수행하였다. 수치적으로 예측된 와류흘림주파수는 모형시험으로 계측한 명음 발생 주파수 및 날개 고유진동수와 일치함을 확인하였다.