화염의 운동학 모사를 위한 G-equation solver 개발

Abstract

본 연구에서는 화염의 운동학적 특성을 모사하기 위한 G-equation solver를 개발하였다. 이 solver는 1차 업윈드 스킴과 Back and Forth Error Compensation and Correction (BFECC) 알고리즘을 결합한 수치 기법을 사용하며, 시간과 공간에 대해 2차 정확도를 제공한다. 계산 속도를 향상시키기 위해 OpenMP 병렬 프로그래밍이 적용되었다. G-equation은 두 가지 방법, 즉 전역 level-set 방법과 화염 면 주변 영역만을 계산하는 국부적 level-set 방법을 사용하여 계산하였다. Zalesak 디스크 문제를 사용한 계산 효율 비교를 통해 국부적 level-set 방법이 전역 level-set 방법보다 더 효율적인 방법임을 확인하였다. G-equation solver의 정확성은 (i) 회전 이동하는 Zalesak 디스크의 형상 비교, (ii) 진동하는 V-화염의 형상 비교, (iii) 화염전달함수의 이론 모델과의 비교를 통해 검증하였다. 화염 면의 열 방출량 계산에 있어서는 화염 면 적분을 이산 델타 함수를 활용한 유동 영역 전체의 부피 적분으로 전환하고, 이를 수치 적분하였다. 개발된 G-equation solver를 이용하여, 종 방향으로 조화 진동하는 유동장에서 원추형 분젠 화염의 화염전달함수를 계산하였다. 특히, 진폭이 증가하는 경우에 대한 이론적 모델과의 차이를 비교하였다. 또한, 횡 방향으로 조화 진동하는 유동장에서의 화염전달함수에 대해서 2차원과 3차원 시뮬레이션 결과를 비교 분석하였다. 본 연구를 통해 개발된 G-equation solver는 향후 음향 해석 툴과 결합하여 연소 불안정 예측에 활용하고자 한다.