벽함수

전산유체역학 난류 마지막 수정: 2026-07-03 11:38:52

1. 개요[편집]

벽함수 Wall Function
분야	난류 모델링 벽 처리
핵심 무차원수	$y^+$ (무차원 벽거리), $u^+$ (무차원 속도)
이론적 근거	로그 법칙 (law of the wall)
목적	벽 근처 격자 비용 절감

벽까지 다 풀 자신이 없으면, 벽 근처는 공식으로 때우면 된다. 문제는 그 공식이 언제 거짓말을 하느냐다.

벽함수(wall function)는 벽 근처 난류 경계층을 격자로 직접 해상하는 대신, 실험과 이론으로 다져진 반경험식으로 벽과 첫 격자점 사이를 이어 붙이는 벽 처리 기법이다. 벽 아주 가까이에서 난류 구조는 극도로 작아지고 속도구배는 가팔라져, 이곳을 격자로 곧이곧대로 풀려면 셀을 어마어마하게 촘촘히 깔아야 한다. 벽함수는 이 비싼 영역을 계산에서 생략하고 알려진 속도 분포 공식으로 대체함으로써 격자 수와 계산 비용을 극적으로 줄인다.

RANS 기반 산업 CFD에서 벽함수는 오랫동안 사실상의 기본 설정이었다. “벽을 다 풀 것인가, 공식으로 때울 것인가”는 격자를 짜기 전에 반드시 내려야 하는 결정이며, 이 선택 하나가 격자 규모를 몇 배씩 바꿔 놓는다.

2. 무차원 벽거리와 무차원 속도[편집]

벽 근처 유동을 기술하려면 먼저 좌표와 속도를 벽 단위로 무차원화해야 한다. 기준이 되는 것은 벽면 전단응력 $\tau_w$ 에서 정의되는 마찰속도(friction velocity)다.

u_\tau = \sqrt{\frac{\tau_w}{\rho}}

이 마찰속도로 벽으로부터의 거리 $y$ 와 유동 속도 $u$ 를 각각 무차원화한 것이 무차원 벽거리 $y^+$ 와 무차원 속도 $u^+$ 다.

y^+ = \frac{\rho \, u_\tau \, y}{\mu}, \qquad u^+ = \frac{u}{u_\tau}

$y^+$ 는 본질적으로 벽 근처에 국소적으로 정의된 레이놀즈수로, “이 격자점이 벽에서 얼마나 떨어져 있는가”를 물리적으로 의미 있는 척도로 재 준다. 놀랍게도 벽 근처 난류 경계층의 평균 속도 분포는 유동 종류에 거의 무관하게 $u^+$ 가 $y^+$ 만의 보편 함수로 표현된다 — 이 보편성이 벽함수라는 편법을 가능하게 하는 근본 토대다.

3. 벽 근처의 층 구조[편집]

$y^+$ 를 기준으로 벽 근처는 성격이 뚜렷이 다른 몇 개 층으로 나뉜다.

층	대략적 범위	지배 물리	속도 분포
점성저층	$y^+ < 5$	점성 지배	$u^+ = y^+$ (선형)
완충층	$5 < y^+ < 30$	점성 ≈ 난류	어느 공식도 안 맞음
로그층	$30 < y^+ < 300$	난류 지배	로그 법칙

점성저층(viscous sublayer)에서는 점성이 완전히 지배하여 전단응력이 일정하다고 보면 속도가 벽거리에 정비례한다. 즉 $u^+ = y^+$ 라는 아름답도록 단순한 선형 관계가 성립한다.

로그층(log-law region)에서는 난류가 지배하며 속도가 벽거리의 로그에 비례한다. 이것이 그 유명한 로그 법칙(law of the wall)이다.

u^+ = \frac{1}{\kappa}\ln y^+ + B

여기서 $\kappa \approx 0.41$ 은 폰 카르만 상수, $B \approx 5.0$ 은 매끄러운 벽에 대한 상수다. 이 두 값은 특정 유동에서 유도된 것이 아니라 수많은 실험을 관통하는 보편 상수라는 점이 핵심이다.

그리고 그 사이에 낀 완충층(buffer layer)이 문제아다. 점성과 난류가 엇비슷하게 겨루는 이 지대에서는 선형 법칙도 로그 법칙도 맞지 않는다. 벽함수를 쓸 때 첫 격자점을 완충층( $5 < y^+ < 30$ )에 두는 것이 최악의 선택인 이유가 바로 이것 — 어느 공식으로도 이어 붙일 근거가 없는 사각지대이기 때문이다.¹

4. 벽함수 접근 vs 벽해상 접근[편집]

벽 처리는 크게 두 전략으로 갈린다. 첫 격자점을 어디에 두느냐가 전부를 결정한다.

4.1. 벽함수 접근 ( $y^+ \approx 30 \sim 300$ )[편집]

첫 격자점을 로그층 안에 두고, 벽과 첫 격자점 사이는 로그 법칙으로 이어 붙인다. 점성저층과 완충층을 통째로 격자에서 생략하므로 벽에 수직한 방향 셀 수가 확 줄어든다. 계산이 싸고 강건해 k-ε 모델 같은 고레이놀즈수 모델과 궁합이 좋다. 단점은 로그 법칙이 성립한다는 전제가 깔려 있다는 것 — 박리·재부착·강한 압력구배처럼 로그 법칙 자체가 깨지는 유동에서는 정확도를 보장할 수 없다.

4.2. 벽해상 접근 ( $y^+ \approx 1$ )[편집]

첫 격자점을 점성저층 안( $y^+ \approx 1$ )에 두고 벽까지 격자로 직접 적분한다. 저레이놀즈수(low-Re) 처리라고도 부른다. 열전달 계수, 마찰 항력, 박리 시작점처럼 벽 근처 물리가 결과를 좌우하는 문제에서는 이쪽이 정답이다. 대신 프리즘 레이어를 겹겹이 얇게 쌓아야 하므로 격자 수와 셀 종횡비가 크게 늘어난다.

\begin{array}{lcc} & \text{벽함수} & \text{벽해상(low-Re)} \\ \text{첫 격자 } y^+ & 30 \sim 300 & \approx 1 \\ \text{격자 비용} & \text{낮음} & \text{높음} \\ \text{박리·열전달 정확도} & \text{제한적} & \text{우수} \end{array}

현대 상용 코드들은 이 둘 사이 어중간한 $y^+$ 에서도 자동으로 적절히 전환해 주는 확장 벽함수(scalable wall function, all- $y^+$ 처리)를 제공한다. 그럼에도 원리적으로 완충층에 첫 격자점을 걸치는 것은 피하는 게 상책이다.²

5. 격자 비용 절감이라는 본질[편집]

벽함수의 존재 이유는 단 하나, 돈이다(정확히는 격자와 계산 시간). 벽해상을 하려면 벽에 수직한 방향으로 첫 셀을 $y^+ \approx 1$ 에 맞춰야 하는데, 레이놀즈수가 높아질수록 이 물리적 높이는 눈물겹게 작아진다. 자동차 외부 유동이나 항공기 익형처럼 레이놀즈수가 수백만을 넘는 문제에서 벽면 전체를 $y^+ \approx 1$ 로 감싸면 격자 수가 억 단위로 폭발한다.

벽함수는 점성저층·완충층을 격자에서 들어냄으로써 이 부담을 통째로 회피한다. 벽에 수직한 방향 셀을 여러 겹 줄일 수 있어, 전체 격자가 몇 배씩 가벼워지고 그만큼 해석이 빨라진다. 산업 CFD가 수십 년간 벽함수에 의존해 온 이유는 낭만이 아니라 마감과 예산이었다.³

물론 공짜 점심은 없다. 벽함수는 “로그 법칙이 성립하는 평형 경계층”이라는 전제를 딛고 서 있고, 그 전제가 무너지는 곳(강한 박리, 충돌 제트, 회전 유동)에서는 조용히 틀린 답을 예쁘게 수렴시켜 내놓는다. 잔차가 잘 떨어졌다고 벽 처리가 옳았다는 뜻은 결코 아니다.⁴

6. 관련 문서[편집]

완충층에 첫 격자점을 걸쳐 놓고 “왜 결과가 이상하죠?”라고 묻는 것은 CFD 입문자의 통과의례에 가깝다. $5 < y^+ < 30$ 은 유체 해석의 버뮤다 삼각지대다. ↩
확장 벽함수가 나왔다고 완충층을 마음 놓고 밟아도 된다는 뜻은 아니다. “자동으로 처리해 준다”와 “정확하다”는 전혀 다른 말이다. 자동 변속기가 있다고 아무 데나 기어를 넣어도 되는 건 아니듯이. ↩
DNS 논문 저자들은 벽함수를 미개한 근사라 여기고, 마감에 쫓기는 현업 엔지니어는 벽해상 격자를 사치라 여긴다. 둘 다 각자의 자리에서 옳다. ↩
수렴은 대수방정식을 잘 풀었다는 뜻이지, 올바른 방정식을 풀었다는 뜻이 아니다. 이 문장은 CFD를 하는 한 평생 되뇌게 된다. ↩

벽함수