SIMWIKI DOC // 전산유체역학

경계 조건

전산유체역학 CONVERGED last update 2026-07-01 10:26:53

1. 개요[편집]

경계 조건(boundary condition, BC)은 해석 도메인의 경계에서 유동 변수 또는 그 구배에 부여하는 구속 조건으로, 편미분방정식 문제를 유일한 해를 갖는 적정(well-posed) 문제로 만들어 주는 필수 요소이다. 전산유체역학에서 지배 방정식과 격자가 아무리 훌륭해도 경계 조건이 물리적으로 틀렸다면 솔버는 “틀린 문제의 정확한 해”를 성실하게 계산해 줄 뿐이다.

경계 조건이 유독 중요한 이유는, 실제 유동은 무한한 세계 속에 있지만 계산 도메인은 어딘가에서 반드시 잘려야 하기 때문이다. 그 절단면에서 “바깥 세계가 이랬을 것”이라고 선언하는 행위가 곧 경계 조건 설정이며, 이 선언이 곧 해석의 가정이 된다. 잔차가 수렴하지 않을 때 엔지니어가 격자 다음으로 의심해야 할 용의자이기도 하다.¹

2. 수학적 배경: 적정성과 세 가지 유형[편집]

아다마르(Hadamard)의 정의에 따르면 적정 문제는 해가 존재하고, 유일하며, 입력에 연속적으로 의존해야 한다. 경계 조건의 개수와 유형은 방정식의 성격(타원형·포물형·쌍곡형)이 결정한다. 필요한 것보다 적게 주면 해가 유일하지 않고, 많이 주면 해가 존재하지 않으며, 유형이 틀리면 경계에서 비물리적 반사가 발생한다. 수학적으로 경계 조건은 세 가지로 분류된다.

디리클레(Dirichlet) 조건: 변수 값 자체를 지정. $\phi = \phi_b$
노이만(Neumann) 조건: 경계 법선 방향 구배를 지정. $\partial \phi / \partial n = g$
로빈(Robin) 조건: 값과 구배의 선형 결합을 지정. $a\phi + b \, \partial\phi/\partial n = c$

CFD 코드의 경계 조건들은 결국 이 셋의 조합이다. 속도 입구는 속도에 대한 디리클레, 출구의 zero-gradient는 노이만, 대류 열전달 벽면 $-k\,\partial T/\partial n = h(T - T_\infty)$ 은 로빈 조건의 대표적 예다.

3. 입구 경계 조건[편집]

속도 입구(velocity inlet): 속도 벡터를 직접 지정하는 디리클레 조건. 비압축성 유동의 표준 입구다. 이때 압력은 지정하지 않는다(도메인 내부에서 결정되도록 남겨둔다). 속도와 압력을 동시에 고정하면 과잉 구속으로 문제가 부적정해진다.
질량 유량 입구(mass flow inlet): 질량 유량을 지정하고 속도 분포는 코드가 역산. 압축성 유동처럼 밀도가 변해 속도를 미리 알 수 없을 때 쓴다.
압력 입구(pressure inlet): 전압력(total pressure)을 지정. 팬 흡입구나 대기 개방면처럼 유량이 결과로 나와야 하는 경우에 적합하다.

난류 모델링을 쓴다면 입구에서 난류량도 지정해야 한다. $k$ 와 $\varepsilon$ (또는 $\omega$ )을 직접 줄 수도 있지만, 실무에서는 난류 강도 $I$ 와 길이 스케일 $\ell$ (또는 수력직경, 점성비)로 환산하는 방식이 일반적이다.

k = \frac{3}{2}(U I)^2, \qquad \varepsilon = C_\mu^{3/4} \frac{k^{3/2}}{\ell}

문제는 입구 난류량의 실측값을 아는 경우가 거의 없다는 것. 그래서 태초부터 내려오는 기본값 5%가 온 세상의 입구에 복사-붙여넣기되고 있다.²

4. 출구 경계 조건[편집]

압력 출구(pressure outlet): 정압을 지정(디리클레)하고 나머지 변수는 zero-gradient(노이만)로 처리. 가장 무난하고 강건한 출구다. 출구에서 유입이 발생할 경우를 대비해 역류 시의 온도·난류량(backflow condition)을 지정하게 되어 있는데, 이 값을 대충 넣으면 역류가 생기는 순간 그 대충이 도메인 안으로 침투한다.
아웃플로(outflow): 완전발달 유동을 가정하고 모든 변수에 zero-gradient를 적용. 출구 압력조차 모를 때 쓰지만, 출구 근처에 재순환이 있거나 출구가 여러 개면서 유량 배분을 모르는 경우에는 쓰면 안 된다. 가정이 깨진 채로도 그럴듯한 그림이 나온다는 점이 오히려 위험하다.

공통 원칙: 출구는 유동이 충분히 정리된 곳, 즉 zero-gradient 가정이 실제로 성립할 만한 위치에 두어야 한다. 재순환 영역을 출구가 가로지르는 순간 수렴성과 정확도가 동시에 무너진다.

5. 벽면 경계 조건[편집]

점착(no-slip) 조건: 벽면에서 유체 속도 = 벽 속도. 점성 유동의 물리적 기본값이며, 벽 근처 경계층이 형성되는 원인이다. 난류 해석이라면 벽함수를 쓸지 벽까지 해상할지에 따라 첫 격자 높이( $y^+$ ) 요구가 완전히 달라지므로, 벽면 조건은 격자 설계와 한 몸이다.
미끄럼(slip) 조건: 법선 속도 0, 접선 응력 0. 점성 효과를 무시하고 싶은 원방 경계나 자유표면의 근사, 비점성 해석에 쓴다.
이동/회전 벽: 벽에 접선 속도를 부여한다. 회전 기계, 컨베이어, 쿠에트 유동 등. 벽이 유체를 끌고 가는 효과가 그대로 전단으로 들어간다.
거칠기(roughness): 등가 모래알 거칠기 등으로 벽함수를 수정해 표면 거칠기가 마찰과 열전달에 미치는 영향을 반영한다. 거칠기 높이가 첫 격자 높이보다 크면 모델 가정이 깨진다는 점은 의외로 자주 무시된다.

6. 대칭·주기·원방 경계 조건[편집]

대칭(symmetry): 법선 속도 0, 모든 변수의 법선 구배 0. 형상과 유동이 모두 대칭일 때 도메인을 반으로 줄여주는 고마운 조건이지만, “형상의 대칭”이 “유동의 대칭”을 보장하지 않는다는 함정이 있다. 원기둥 후류의 카르만 와열처럼 대칭 형상에서 비대칭 유동이 나오는 경우, 대칭 조건은 물리 현상 자체를 금지해 버린다.
주기(periodic): 마주 보는 두 경계를 이어 붙여 무한 반복 구조를 표현한다. 병렬 유로, 터보기계의 단일 블레이드 통로(회전 주기), 완전발달 유동 등에 쓰며, 유동 방향 주기 조건에는 구동력으로서 압력 강하(pressure jump)나 질량 유량을 지정한다.
원방(far-field): 외부 유동에서 무한히 먼 곳의 상태(마하수, 정압, 온도, 받음각)를 지정한다. 압축성 해석에서는 리만 불변량 기반의 특성 경계 조건으로 파동이 경계에서 반사되지 않고 빠져나가게 처리한다. 경계가 물체에 너무 가까우면 지정한 원방 값이 물체 주위 유동을 왜곡하므로, 익형 해석에서 코드 길이의 수십 배 밖에 두는 것이 관례다.

7. 흔한 실수와 실무 팁[편집]

경계 조건 설정 실수는 CFD 오류 원인 중 가장 흔한 축에 속하며, 대부분 다음 목록의 재방송이다.

출구를 너무 가까이 자르기: 물체 바로 뒤, 재순환이 살아 있는 곳에서 도메인을 자르면 역류 경고가 뜨고 수렴이 망가진다. 후류 방향으로 도메인을 넉넉히 연장하는 것이 정석. “reversed flow in N faces” 경고가 반복되면 그것은 솔버의 비명이다.
입구 난류 강도를 무조건 5%로: 정온한 풍동은 1% 미만, 배관 완전발달 난류는 5~10%, 터보기계 후단은 그 이상이다. 열전달 계수나 박리 위치가 중요한 해석이라면 입구 난류량 민감도 검토가 필수다.
입구와 출구의 과잉 구속: 모든 경계에 속도를 지정하면 질량 보존과 충돌해 압력장이 표류한다. 비압축성 해석에서 압력 지정 경계가 하나도 없다면 기준 압력점(reference pressure)을 잡아 줘야 한다는 것도 같은 맥락.
대칭 조건의 남용: 정상 해석으로 절반 도메인 + 대칭 조건을 쓰다가, 비정상 해석으로 넘어가면서도 그대로 두면 와류 이탈이 억제된 반쪽짜리 물리를 얻는다.
입구-출구를 물체에 너무 붙이기: 입구가 가까우면 물체가 만드는 상류 영향(upstream effect)이 입구의 균일 유속 가정과 충돌한다. 입구 쪽도 여유가 필요하다.
경계 조건 민감도 검토 생략: 격자 수렴성 검토처럼, 도메인 크기와 경계 값에 대한 민감도 검토도 보고서의 방어력을 결정한다.³

8. 관련 문서[편집]

첫 번째 용의자는 언제나 격자, 두 번째가 경계 조건, 세 번째가 완화 계수다. 솔버 버그를 의심하는 것은 앞의 셋을 모두 무죄 방면한 뒤에도 늦지 않다. 대부분 그 전에 범인이 잡힌다. ↩
상용 코드의 기본값이 5%이기 때문이다. 전 세계 수많은 해석 보고서의 입구 난류 강도가 5%로 통일되어 있는 것은 유체역학적 보편 상수가 발견되어서가 아니다. ↩
“도메인을 두 배로 키웠더니 양력계수가 3% 변했습니다”라는 문장은 검토 회의에서 당신을 구원할 수 있다. 반대로 그 검토를 안 했다는 사실은 반드시 가장 곤란한 순간에 발각된다. ↩