스팔라트-알마라스 모델

편집 역사 토론
난류 전산유체역학 마지막 수정: 2026-07-11 04:12:22

1. 개요[편집]

스팔라트-알마라스 모델(Spalart-Allmaras, 흔히 SA 모델)은 난류 와점성(eddy viscosity)과 관련된 변수 하나에 대한 수송방정식만 푸는 1방정식 난류 모델링 기법이다. 1992년 필립 스팔라트(Philippe Spalart)와 스티븐 알마라스(Steven Allmaras)가 항공우주 응용, 특히 벽에 붙어 흐르는 경계층 유동을 겨냥해 만들었다.1 방정식이 하나뿐이라 k-엡실론 모델이나 k-오메가 SST 모델 같은 2방정식 모델보다 가볍고 튼튼해서, 외부 공력해석의 국룰로 오래 사랑받아 왔다.

“난류를 방정식 하나로 퉁친다”는 발상이 무모해 보이지만, 대상을 항공기 날개·동체 주위의 부착 경계층으로 좁히면 이야기가 달라진다. 그 좁은 영역에서는 SA가 놀랄 만큼 잘 맞고, 무엇보다 수렴이 얌전하다. 리뷰어가 “왜 SA 썼냐”고 물으면 “안정적이고 빠르니까”라고 당당히 답할 수 있는 몇 안 되는 모델이다.

2. 왜 1방정식인가[편집]

난류 모델의 방정식 개수는 대개 “얼마나 많은 물리를 수송방정식으로 실어 나르느냐”를 뜻한다. 2방정식 모델은 난류운동에너지와 그 소산 관련 척도를 각각 따로 수송하지만, SA는 그 과정을 통째로 생략하고 와점성에 직결된 변수 ν~\tilde{\nu} 하나만 수송한다.

이렇게 하면 얻는 게 세 가지다.

  • 계산 비용 절감: 수송방정식이 하나 줄면 그만큼 반복당 연산과 메모리가 준다. 대형 격자·다케이스 최적화 반복에서 이 차이가 누적된다.
  • 강건한 수렴: 방정식 간 커플링이 단순해 잔차가 얌전히 떨어진다. 2방정식 모델의 ε\varepsilon이나 ω\omega가 벽 근처에서 날뛰며 발산시키는 사고가 거의 없다.
  • 경계값 둔감: 자유류 경계값에 대한 민감성이 원조 k-엡실론 모델·k-ω보다 낮다.

대가는 물론 일반성이다. 여러 물리를 하나로 압축했으니, SA가 조율된 영역(부착 경계층, 가벼운 박리)을 벗어나면 정확도가 떨어진다.

3. 수송방정식[편집]

SA는 수정 난류점성 ν~\tilde{\nu}에 대한 하나의 수송방정식을 푼다. 개념적 형태는 다음과 같다.

∂ν~∂t+u⋅∇ν~=cb1S~ν~⏟생성−cw1fw(ν~d)2⏟소멸+1σ[∇⋅((ν+ν~)∇ν~)+cb2∣∇ν~∣2]⏟확산\frac{\partial \tilde{\nu}}{\partial t} + \mathbf{u}\cdot\nabla\tilde{\nu} = \underbrace{c_{b1}\tilde{S}\tilde{\nu}}_{\text{생성}} - \underbrace{c_{w1}f_w\left(\frac{\tilde{\nu}}{d}\right)^2}_{\text{소멸}} + \underbrace{\frac{1}{\sigma}\left[\nabla\cdot\left((\nu+\tilde{\nu})\nabla\tilde{\nu}\right) + c_{b2}|\nabla\tilde{\nu}|^2\right]}_{\text{확산}}

세 항의 역할은 직관적이다. 생성항은 와도 크기 S~\tilde{S}에 비례해 난류점성을 키우고, 소멸항은 벽까지의 거리 dd에 반비례해 벽에 가까울수록 난류점성을 억누른다. 확산항은 이를 공간적으로 퍼뜨린다. 실제 나비에-스토크스 방정식에 들어가는 와점성 μt=ρν~fv1\mu_t = \rho\tilde{\nu}f_{v1}은 감쇠함수 fv1f_{v1}을 통해 점성저층에서 부드럽게 0으로 꺼진다.2

여기서 핵심은 벽까지의 거리 dd가 방정식에 직접 등장한다는 점이다. 덕분에 SA는 벽함수 없이 벽까지 직접 적분하는 것을 전제로 설계됐고, 반대로 복잡한 형상에서 벽 거리 계산이 부실하면 성능이 흔들린다.

4. 격자와 벽 처리[편집]

SA는 원래 저-레이놀즈수 모델, 즉 벽 첫 격자점을 점성저층 안(y+≈1y^+ \approx 1)에 넣고 벽까지 적분하는 방식으로 만들어졌다. 그래서 정확한 항력·박리 예측을 원하면 벽 법선방향 격자를 촘촘히 갈아 넣어야 한다.

다만 현대 상용 구현은 대부분 SST처럼 성긴 격자에서 벽함수로 자동 전환하는 처리를 얹어, y+y^+가 어중간해도 결과가 급격히 무너지지 않게 완충한다. 그래도 “격자가 다 한다”는 원칙은 여기서도 예외 없이 적용된다.

5. 강점과 한계[편집]

강점:

  • 부착 경계층과 가벼운 역압력구배 유동에서 항력·표면마찰 예측이 정확하다.
  • 항공기 외부 유동, 익형, 미사일·발사체 등 외부 공력 문제에서 검증이 두텁다.
  • 수렴이 빠르고 튼튼해 형상 최적화의 수천 회 반복 계산에 적합하다.

한계:

  • 자유전단류(제트·후류), 강한 박리, 복잡한 재순환 영역에서는 부정확하다. 애초에 그런 유동을 겨냥해 만든 모델이 아니다.
  • 등방성 와점성 가정을 깔고 있어 강한 곡률·회전·2차 유동을 못 잡는다.
  • 열전달이나 화학반응처럼 난류 스칼라 수송이 중요한 문제에는 정보가 부족하다.

이런 이유로 큰 박리가 지배하는 문제에서는 SA를 RANS 부분으로 삼는 DES(Detached Eddy Simulation) 계열로 확장해 대와류 모사와 결합하는 방식이 널리 쓰인다. 사실 DES라는 개념 자체가 1997년 스팔라트가 SA를 토대로 제안한 것이라3, SA는 하이브리드 난류 해석의 어머니 격이기도 하다.

6. 언제 SA를 고르나[편집]

실무 감각으로 정리하면 이렇다. 벽에 얌전히 붙어 흐르는 외부 유동이고, 항력이 관심사이며, 최적화 반복을 빨리 돌려야 한다면 SA가 좋은 첫 선택이다. 반대로 강한 박리·재순환·자유전단류가 지배하거나 열유동·다물리가 얽히면 k-오메가 SST 모델이나 그 이상의 모델로 넘어가는 편이 낫다. 방정식 하나짜리 모델에 만능을 기대하는 것은 처음부터 무리다.4

7. 관련 문서[편집]

8. Footnotes[편집]

  1. Spalart, P. R. & Allmaras, S. R. (1992). “A One-Equation Turbulence Model for Aerodynamic Flows”. AIAA Paper 92-0439. 두 사람 이름을 딴 모델이 이렇게 널리 쓰일 줄, 정작 본인들도 몰랐을 것이다. ↩

  2. 감쇠함수 fv1,fv2,fwf_{v1}, f_{v2}, f_w 같은 보정함수가 잔뜩 붙어 있어서, “1방정식이라 단순하다”는 말은 방정식 개수 이야기일 뿐 실제 구현은 결코 단순하지 않다. 상수 표를 처음 보면 대부분 조용히 창을 닫는다. ↩

  3. Spalart et al. (1997)의 DES. 격자 간격이 벽 거리보다 작아지는 영역에서 소멸항의 길이척도를 격자 크기로 바꿔치기해, 벽 근처는 RANS·박리 후류는 LES처럼 굴러가게 만든 아이디어다. ↩

  4. 그럼에도 “일단 SA로 돌려보고 안 맞으면 SST”라는 워크플로가 현장에 넓게 퍼져 있다. 빠른 모델로 대충 훑고 무거운 모델로 확인하는 건 언제나 합리적인 전략이다. ↩