레이놀즈 응력 모델

난류 전산유체역학 마지막 수정: 2026-07-20 04:13:52

1. 개요[편집]

레이놀즈 응력 모델(Reynolds Stress Model, RSM)은 난류의 레이놀즈 응력 텐서 6개 성분 각각에 대해 별도의 수송방정식을 풀어 난류를 닫는(closure) 2차 폐합(second-moment closure) 난류 모델이다. k-엡실론 모델이나 k-오메가 SST 모델 같은 와점성(eddy-viscosity) 모델이 “난류 응력은 평균 변형률에 비례한다”는 등방성 가정을 깔고 가는 반면, RSM은 그 가정을 아예 버리고 응력의 이방성(anisotropy)을 직접 추적한다. 한마디로 난류 모델링의 정공법이자, 그만큼 값비싼 사치다.¹

RANS 계열 난류 모델링 중에서는 가장 정교한 축에 속한다. 곡률, 스월, 2차 유동처럼 등방성 가정이 무너지는 상황에서 강점을 보이지만, 방정식이 많고 수렴이 까다로워서 “쓸 줄 알면 좋은데 잘 안 쓰게 되는” 모델의 대표주자이기도 하다.

2. 왜 와점성을 버리는가[편집]

나비에-스토크스 방정식을 시간 평균(레이놀즈 평균)하면, 비선형 대류항에서 $-\rho\,\overline{u_i' u_j'}$ 라는 항이 남는다. 이것이 레이놀즈 응력으로, 대칭 텐서라 미지수가 6개다. 방정식은 4개(연속 1 + 운동량 3)인데 미지수가 늘었으니 계가 닫히지 않는다. 이 폐합 문제(closure problem)를 어떻게 때우느냐가 난류 모델의 정체성이다.

와점성 모델은 부시네스크 가정(Boussinesq hypothesis)으로 때운다.

-\rho\,\overline{u_i' u_j'} = \mu_t\!\left(\frac{\partial \overline{u_i}}{\partial x_j} + \frac{\partial \overline{u_j}}{\partial x_i}\right) - \frac{2}{3}\rho k\,\delta_{ij}

즉 6개 응력을 스칼라 하나( $\mu_t$ )로 뭉뚱그린다. 계산은 싸지만, 이 식은 난류가 방향에 무관하게 똑같이 거칠다고 가정하는 것과 같다. 현실의 난류는 벽 근처나 회전 유동에서 성분마다 크기가 다른 강한 이방성을 띠므로, 이 가정은 근본적으로 틀렸다. 대개는 그럭저럭 쓸 만한 근사일 뿐이다.

3. 응력 수송방정식[편집]

RSM은 각 $\overline{u_i' u_j'}$ 에 대해 정확한 수송방정식을 세운다. 개념적으로 다음과 같은 꼴이다.

\frac{D\,\overline{u_i' u_j'}}{Dt} = P_{ij} + \Phi_{ij} - \varepsilon_{ij} + D_{ij}

우변의 네 항은 각각 다음을 뜻한다.

생성항 $P_{ij}$ : 평균 유동의 변형률에서 난류 에너지가 만들어지는 항. 미지수 없이 정확히 계산되는 유일한 항이라 RSM의 최대 강점의 원천이다. 곡률·스월의 효과가 자연스럽게 들어온다.
압력-변형률 재분배항 $\Phi_{ij}$ : 성분 간에 에너지를 재분배해 이방성을 완화하는 항. 정확한 형태를 모르므로 모델링이 필요하며, RSM의 정확도와 성패를 좌우하는 가장 어려운 조각이다.²
소산항 $\varepsilon_{ij}$ : 점성에 의한 소산. 보통 등방성으로 가정하고 별도의 $\varepsilon$ (또는 $\omega$ ) 방정식으로 스케일을 공급한다.
확산항 $D_{ij}$ : 난류·점성 확산.

결국 6개 응력 방정식 + 소산율 방정식 하나로, 3차원에서 최소 7개의 추가 수송방정식을 푼다. 2방정식 모델(k와 ε 둘)의 3.5배 부담이다.

4. 강점: 이방성과 곡률[편집]

RSM이 진가를 발휘하는 곳은 와점성 모델이 원리적으로 실패하는 유동이다.

강한 스월·회전: 사이클론 분리기, 스월 버너, 회전 좌표계 기법으로 푸는 회전기계. 회전이 응력 성분을 재분배하는 효과를 생성항이 직접 담는다.
곡률이 큰 유동: 곡관, U-벤드, 터빈 블레이드 통로. 볼록/오목 벽에서 난류가 억제·증폭되는 방향성을 잡아낸다.
2차 유동: 사각 덕트 코너의 프란틀 2종 2차 유동(Prandtl’s second kind)은 응력 이방성이 유일한 구동력이라, 와점성 모델은 아예 이 현상을 재현하지 못한다. RSM은 예측한다.

이런 문제에서 k-엡실론 모델은 정성적으로도 틀린 답을 내는 반면, RSM은 정성·정량 모두 개선된 결과를 준다.

여기서 중요한 통찰 하나. RSM의 강점은 “더 많은 방정식을 풀어서”가 아니라 “생성항을 근사 없이 정확히 계산해서” 나온다. 와점성 모델은 생성항조차 부시네스크 가정을 통과시키며 곡률·회전 정보를 흘려버리는데, RSM은 응력마다 생성 메커니즘을 있는 그대로 담는다. 방정식 개수는 결과이지 원인이 아니라는 점을 이해하면, RSM이 왜 특정 유동에서만 확 좋아지는지가 명확해진다.

5. 단점: 비용과 수렴성[편집]

그런데도 RSM이 산업계 기본값이 못 되는 이유는 명확하다.

계산비용: 방정식이 7개로 늘어 자유도당 메모리·연산이 늘고, 수렴까지 반복 횟수도 많다. 셀당 부담이 2방정식 모델의 배 이상.
수렴성·강건성: 응력 방정식들이 서로 강하게 연성되어 있고 압력-변형률 모델이 뻣뻣해서, 초기값에 민감하고 잘 발산한다. 실무 국룰은 “k-엡실론 모델로 수렴시킨 장을 초기값으로 얹고 RSM으로 갈아탄다”이다.
경계 근처 민감성: 벽 근처 처리와 벽함수 선택에 결과가 예민하다.
모델 상수의 불확실성: 압력-변형률 모델 계수들이 표준화가 덜 되어 있어, 문제마다 최적 조합이 다르다.

정리하면 RSM은 “정확도를 위해 강건성과 비용을 갈아 넣는” 트레이드오프의 극단이다.

6. 대비: k-ε, k-ω와의 위치[편집]

모델	추가 방정식	이방성	비용	강건성
k-엡실론 모델	2	없음(등방성)	낮음	높음
k-오메가 SST 모델	2	없음(등방성)	낮음	높음
RSM	5~7	직접 해석	높음	낮음
대와류 모사	(여과)	대와류 직접	매우 높음	중간

RSM은 RANS의 정확도 상한선이자, 대와류 모사로 넘어가기 직전의 마지노선쯤에 위치한다. 스케일 분해 없이 이방성을 담고 싶다면 RSM, 스케일까지 분해할 예산이 있으면 LES로 가는 것이 요즘의 판단 기준이다. 다만 시간 정확 대와류 해석의 비용을 감당 못 하는 정상상태 산업 문제에서, 곡률·스월이 핵심이라면 RSM은 여전히 대체 불가능한 선택지다.

7. 관련 문서[편집]

RSM은 넓게 보면 2차 폐합의 대명사지만, 정확히는 응력 방정식 미분형(differential RSM, DRSM)을 뜻한다. 대수식으로 근사한 대수 응력 모델(ASM/EARSM)도 사촌뻘로 존재하는데, 이쪽은 이방성 일부를 담으면서 비용은 2방정식 수준이라 절충안으로 인기다. ↩
압력-변형률 항 모델링의 역사가 곧 RSM의 역사다. 로타(Rotta)의 선형 복귀 모델에서 시작해 LRR, SSG 모델로 발전했다. 이 항을 얼마나 잘 모델링하느냐가 “RSM인데 왜 k-ε보다 못하냐”와 “역시 RSM”을 가른다. ↩