본 연구에서는 비정렬격자계와 체적포착법을 사용하여 표면장력이 지배적인 다상유동의 수치해석 방법을 제시하였다. 먼저 표면장력에 대한 CSF(Continuum Surface Force) 모델을 비정렬격자계에 적용할 수 있도록 수치해석 방법을 확립시켜 Myong(2009)이 개발한 비정렬격자계와 체적포착법을 사용한 수치 해석코드에 삽입하였다. 테스트 문제로 오직 표면장력만이 존재하는 평형상태의 정적(static) 액적 및 비평형상태의 동적(dynamic) 액적 문제에 적용하여, 이 해석방법의 유용성과 정확도를 평가하였다. 연구결과, 매끄러운 곡률 계산을 위해 필요한 필터로 본 연구에서 제안한 Laplacian 필터와 함께 CSF 모델로는 밀도보정(density-scaled)한 CSF 모델이 예측성능이 우수한 것으로 나타났다. 또한 표면장력 계산을 위한 이 모델을 채용한 본 수치해석방법은 표면장력이 지배적인 다상유동인 평형상태의 정적 액적 및 비평형상태의 동적 액적 문제 모두에 대해 정확성과 유용성이 입증되었다.
A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.
