연소 시뮬레이션

시뮬레이션 전산유체역학 계산화학 마지막 수정: 2026-07-16 04:51:07

1. 개요[편집]

연소 시뮬레이션 Combustion Simulation
지배 물리	유동 + 화학반응 + 열전달
반응속도	아레니우스 식
화염 분류	예혼합 / 확산 / 부분예혼합
대표 모델	화염면, PDF, EDC
핵심 난제	강성 ODE, 난류-화학 상호작용

CFD만 해도 어려운데 거기에 화학을 얹었다. 그리고 화학은 CFD보다 어렵다.

연소 시뮬레이션은 화학반응을 동반하는 유동을 수치적으로 푸는 해석 분야로, 유체 수송·화학종 반응·열 방출·복사가 서로 강하게 얽힌 다물리 연성해석 문제다. 가스터빈 연소기, 내연기관, 산업용 버너, 화재 안전, 로켓 엔진, 그리고 최근에는 수소·암모니아 무탄소 연료가 주 무대다.

무엇이 그렇게 어려운가. 일반 CFD는 속도 3개와 압력 1개, 많아야 온도까지 4~~5개 변수를 푼다. 연소는 여기에 **화학종 수십~~수백 개**를 더한다. 게다가 이 종들이 반응하는 시간 스케일이 나노초부터 밀리초까지 6~9자릿수에 걸쳐 퍼져 있고, 화염면 두께는 밀리미터인데 연소기는 미터다. 시간과 공간 양쪽으로 다중 스케일이 터지는 문제 — 그래서 연소 CFD 논문의 절반은 “어떻게 화학을 싸게 만들 것인가”에 관한 것이다.

2. 지배 방정식[편집]

기본은 반응성 나비에-스토크스다. 질량·운동량 보존에 더해 화학종 $k$ 마다 화학종 수송 방정식이 붙는다.

\frac{\partial (\rho Y_k)}{\partial t} + \nabla \cdot (\rho \mathbf{u} Y_k) = \nabla \cdot (\rho D_k \nabla Y_k) + \dot{\omega}_k

$Y_k$ 는 질량분율, $D_k$ 는 확산계수, $\dot{\omega}_k$ 가 화학 생성률이다. 좌변과 우변 첫 항까지는 그냥 대류-확산 방정식이고, 모든 악의 근원은 마지막 소스항 $\dot{\omega}_k$ 다.

에너지 방정식은 엔탈피 형태로 쓰는 것이 편하다. 화학종의 생성 엔탈피를 총 엔탈피에 포함시키면 반응열이 별도 소스항 없이 자동으로 따라오기 때문이다.¹

$K$ 개 화학종 중 독립인 것은 $K-1$ 개( $\sum Y_k = 1$ )라, 상세 메탄 메커니즘(GRI-Mech 3.0, 53종)이면 방정식이 유동 5개 + 종 52개 = 57개. 항공유 상세 메커니즘은 종이 수천 개다. 그냥 못 푼다.

3. 아레니우스 반응속도[편집]

기본 반응 $\sum \nu'_{ki} \mathcal{M}_k \rightarrow \sum \nu''_{ki} \mathcal{M}_k$ 의 진행률은 질량작용 법칙을 따르고, 속도상수는 아레니우스(Arrhenius) 식으로 온도에 의존한다.

k_f = A T^{\beta} \exp\left( -\frac{E_a}{R_u T} \right)

$A$ 는 지수앞인자, $E_a$ 는 활성화 에너지. 이 지수 함수가 연소 시뮬레이션의 모든 고통의 근원이다. 전형적인 $E_a/R_u$ 가 1만~2만 K 수준이라, 1500 K 부근에서 온도가 10 K만 올라가도 반응률이 수 % ~ 십수 % 뛴다. 온도장에 작은 수치 오차가 끼면 반응률이 통째로 흔들리고, 반응률이 흔들리면 발열이 흔들리고, 발열이 온도를 다시 흔든다. 강한 양의 되먹임 — 연소 해석이 잘 발산하는 이유다.

3.1. 상세 vs 축소 메커니즘[편집]

상세(detailed) 메커니즘: 소기 반응(elementary reaction)을 다 넣은 것. 수소-공기가 ~10종/~30반응으로 가장 작고, 메탄 GRI-Mech 3.0이 53종/325반응, n-헵탄이 수백 종, 실제 항공유 대체연료는 수천 종. 0차원 반응기나 1차원 화염 계산에는 쓸 수 있지만 3차원 LES에는 언감생심.
축소(reduced) 메커니즘: 준정상상태 근사(QSSA), 부분평형, 민감도·반응경로 해석(DRG, DRGEP)으로 종과 반응을 쳐낸 것. 메탄을 4~5스텝으로 줄이는 것이 고전.
전역(global) 반응: 1~2스텝. “연료 + 산소 → CO2 + H2O” 한 줄. 화염 온도와 위치 정도만 보면 될 때는 이걸로 충분하고, 실제로 산업 현장에서 제일 많이 쓴다. 대신 CO나 NOx는 절대 못 맞춘다.

메커니즘을 고르는 것은 예산을 고르는 것과 같다. 화염 위치만 알고 싶으면 2스텝, 배출가스를 예측해야 하면 축소 메커니즘 이상. “그냥 상세로 하면 안 되나요”라는 질문에는 계산 시간이 답한다.

4. 강성 ODE 문제[편집]

유동 없이 균질 반응기 하나만 봐도 화학은 이런 상미분방정식계다.

\frac{dY_k}{dt} = \frac{\dot{\omega}_k}{\rho}

문제는 이 계가 극도로 강성(stiff)하다는 것. 라디칼(H, OH, CH3)의 특성 시간은 나노초, 주요 종의 소모는 밀리초. 자코비안 고유값의 비가 $10^6$ ~ $10^9$ 에 달한다. 명시적 룽게-쿠타법으로 풀면 가장 빠른 라디칼에 시간 스텝이 묶여, 밀리초짜리 문제에 나노초 스텝을 밟아야 한다.

그래서 암시적 강성 적분기가 필수다. CVODE/DVODE의 후진 미분 공식(BDF), DASSL, Radau5 등. 이들은 매 스텝 자코비안 행렬을 만들고 뉴턴 반복(뉴턴-랩슨법)으로 푼다. 종이 $K$ 개면 자코비안이 $K \times K$ — 종 100개짜리 메커니즘을 셀 100만 개에서 스텝마다 푼다고 상상해 보자.²

가속 기법이 발달할 수밖에 없었다.

ISAT (In Situ Adaptive Tabulation): 이미 계산한 조성 근처면 선형 보간으로 때운다. 캐시 히트율이 90%를 넘으면 극적으로 빨라진다.
동적 셀 클러스터링: 조성이 비슷한 셀들을 묶어 한 번만 적분.
화학 테이블링: 아예 미리 다 계산해 저장 (아래 화염면 모델).
GPU 오프로딩: 셀마다 독립적인 ODE라 병렬성이 완벽하다. GPU 컴퓨팅의 좋은 먹잇감.

5. 예혼합과 확산 화염[편집]

연소는 연료와 산화제가 언제 섞이느냐로 갈린다.

구분	혼합 시점	예시	특징
예혼합 (premixed)	반응 전 완전 혼합	가솔린 엔진, 분젠 버너 안쪽 파란 불꽃	화염 전파 속도 존재, 역화·소염 위험
확산 (non-premixed)	반응 지점에서 혼합	디젤, 촛불, 산업 버너	안전하지만 국소 고온 → NOx, 그을음
부분예혼합	어중간	린번 가스터빈, GDI 엔진	현실의 대부분

예혼합 화염의 핵심 물성은 층류 화염 속도 $S_L$ 이다. 화염면이 미연 혼합기 속으로 전파하는 속도로, 메탄-공기 화학양론이 대략 0.38 m/s, 수소-공기는 2~3 m/s. 수소가 10배 가까이 빠르다는 것이 수소 연소기 설계에서 역화(flashback)가 악몽인 이유다.

확산 화염은 화염 속도라는 개념 자체가 없다. 대신 혼합분율(mixture fraction) $Z$ 로 기술한다 — 유체 요소 안 연료 유래 질량의 비율로, 화학반응에 대해 보존되는 스칼라다. 이 아이디어가 강력한 이유는, 무한히 빠른 화학을 가정하면 모든 조성과 온도가 $Z$ 하나의 함수가 되기 때문이다. 종 방정식 52개가 스칼라 1개 + 룩업 테이블로 붕괴한다.

6. 난류-화학 상호작용[편집]

진짜 문제는 여기다. 난류 유동에서 반응률의 시간 평균은 평균 조성의 반응률과 전혀 같지 않다.

\overline{\dot{\omega}(Y, T)} \neq \dot{\omega}(\bar{Y}, \bar{T})

아레니우스 항이 지수 함수라 온도 요동이 평균 반응률을 극적으로 끌어올린다. 이 항의 닫힘(closure)이 난류 연소 모델링의 전부다. RANS든 LES든 반응률 소스항은 격자보다 훨씬 작은 화염면 안에서 일어나므로 절대 해상되지 않는다.³

주요 접근법은 셋으로 나뉜다.

6.1. 화염면 모델 (Flamelet)[편집]

난류 화염을 층류 화염 조각(flamelet)들의 앙상블로 본다. 화염 두께가 콜모고로프 스케일보다 얇으면(콜모고로프 스케일), 난류는 화염면을 구기고 늘일 뿐 내부 구조는 못 건드린다는 물리적 근거가 있다.

그래서 1차원 층류 대향류 화염을 미리 상세 화학으로 계산해 테이블로 저장하고, 3차원 계산에서는 혼합분율 $Z$ 와 스칼라 소산율 $\chi$ (또는 진행변수 $c$ )만 수송한다. FGM, FPI, SLFM이 이 계열. 압도적으로 싸다 — 상세 화학의 정확도를 룩업 테이블 비용에 얻는다. 단점은 화염 구조가 테이블에 갇혀 있어서, 국소 소염이나 자착화 같은 테이블 밖 현상에는 취약하다는 것.

6.2. PDF 방법[편집]

조성의 확률밀도함수를 직접 수송한다. 수송 PDF 방정식의 화학 소스항은 닫힌 형태로 정확히 나타난다 — 근사가 필요 없다. 이게 이 방법의 결정적 매력이다. 대신 모델링 부담이 분자 혼합(micromixing) 항으로 옮겨가고, 조성 공간이 고차원이라 유한차분이 불가능해 몬테카를로 방법 기반 확률 입자로 푼다. 정확하지만 비싸고, 통계 수렴을 위해 셀당 입자 수십 개를 굴려야 한다.

6.3. EDC[편집]

Eddy Dissipation Concept. 마그누센(Magnussen)이 제안한 것으로, 반응이 격자 안의 미세구조(fine structure) — 난류 소산이 집중되는 작은 영역 — 에서만 일어나고 그것을 완전혼합 반응기(PSR)로 본다. 그 미세구조의 부피 분율과 체류 시간을 난류량 $k$ , $\varepsilon$ 에서 대수적으로 추정한다.

\gamma^* = C_\gamma \left( \frac{\nu \varepsilon}{k^2} \right)^{1/4}, \qquad \tau^* = C_\tau \left( \frac{\nu}{\varepsilon} \right)^{1/2}

상세 화학을 쓸 수 있으면서 구현이 단순해 ANSYS Fluent나 OpenFOAM에서 널리 쓰인다. 단점은 셀마다 강성 ODE를 풀어야 해서 느리고, 원래 고레이놀즈수 난류를 전제한 모델이라 MILD 연소나 저난류 영역에서는 상수를 조정하지 않으면 이상한 답을 낸다. 그 상수 $C_\gamma = 2.1377$ , $C_\tau = 0.4083$ 을 만지는 것이 이 바닥의 유서 깊은 취미.

더 단순한 Eddy Break-Up이나 Eddy Dissipation Model은 “반응은 혼합 속도로 결정된다”고 밀어붙인다. 화학을 아예 안 푸는 셈이라 빠르지만, 착화·소염을 원리적으로 못 본다.

7. 배출가스 예측[편집]

연소 시뮬레이션이 돈이 되는 이유의 절반은 규제 대응이다.

NOx는 세 경로로 생긴다.

써멀(Zeldovich) NOx: 고온에서 질소 분자가 깨진다. 활성화 에너지가 커서 1800 K 부근부터 급증 — 온도가 지배한다. 린번(lean burn)과 배기가스 재순환(EGR)이 화염 온도를 낮춰 이걸 잡는 전략.
프롬프트(Fenimore) NOx: 화염면의 탄화수소 라디칼이 질소를 공격. 농후(rich) 조건에서.
연료 NOx: 연료 속 질소 화합물에서. 석탄·바이오매스.

써멀 NOx가 온도의 강한 지수 함수라, 온도 예측 오차 50 K가 NOx 예측을 2배 틀리게 만든다. 그래서 NOx는 유동 후처리(post-processing)로 계산하는 것이 관행이다 — 수렴한 온도·조성장 위에 NOx 수송 방정식만 따로 푼다. NOx 농도가 ppm 수준이라 주 유동에 영향을 안 주기 때문에 가능한 편법. 여기에 온도 요동을 반영하려고 PDF를 가정해 적분하는 보정까지 붙인다.

그을음(soot)은 더 나쁘다. 핵생성·표면성장·응집·산화가 얽힌 입자 개수밀도 문제라, 모멘트법(MOMIC)이나 섹셔널 모델을 동원해도 오차 한 자릿수(order of magnitude)면 성공으로 친다.⁴

8. 관련 문서[편집]

이걸 안 하고 감온 엔탈피(sensible enthalpy)를 쓰면 반응열 소스항 $\sum h_{f,k}^0 \dot{\omega}_k$ 를 손으로 넣어야 하는데, 부호 하나 틀리면 불꽃 대신 냉동고가 나온다. 총 엔탈피 형식이 사랑받는 이유다. ↩
그래서 상세 화학 LES의 계산 시간 중 화학 적분이 차지하는 비율이 흔히 70~90%다. 유동을 푸는 CFD 코드가 실제로는 대부분의 시간을 ODE 솔버 안에서 보낸다는 뜻. “이거 CFD 맞나요?”라는 질문이 나오는 지점. ↩
그래서 연소 LES에는 격자를 줄이면 답이 좋아진다는 CFD의 상식이 통하지 않는다. 화염 두께가 0.1 mm인데 격자가 1 mm면 뭘 해도 화염면은 안 보인다. 인위적으로 화염을 두껍게 만들어 격자에 맞추는 두꺼운 화염(thickened flame) 모델이 이 배짱에서 나왔다 — 화학을 왜곡해서 수치에 맞추는, 어떤 의미로는 반칙이지만 잘 통한다. ↩
그을음 예측에서 “실험과 팩터 2 안에 들어왔다”고 하면 학회에서 박수를 받는다. 다른 CFD 분야였으면 100% 오차로 논문 리젝인데, 여기서는 자랑거리다. 기대치는 상대적이다. ↩