최종 수정: 2026-11-04 읽기 약 9분

균일 자유류 U 속에 놓인 얇은 평판 위에서 발달하는 층류 경계층입니다. 압력구배가 없으면 속도 분포가 상사 변수 η 하나의 함수가 되고 (Blasius 1908), 벽 마찰 계수와 배제 두께가 정확히 정해집니다. 벽 no-slip, 벽 근처의 확산, 그리고 이류가 한꺼번에 들어가는 외부 유동이라 이 케이스는 Boreas의 이류 스킴 (semi-Lagrangian · upwind) 이 경계층을 얼마나 정확히 만드는가를 검증합니다. 이 문서의 Boreas 결과는 v0.10.0 GPU 계산이며, 가장 고운 격자가 430만 셀이라 CPU 대조는 하지 않았습니다.

이론해

상사 변수 η 와 무차원 유동함수 f(η) 에 대해 Blasius 방정식과 경계조건은 다음과 같습니다.

η=yUνx,f′′′+12ff′′=0,f(0)=f′(0)=0,  f′(∞)=1\eta = y\sqrt{\frac{U}{\nu x}}, \qquad f''' + \tfrac{1}{2} f f'' = 0, \qquad f(0) = f'(0) = 0,\; f'(\infty) = 1

속도 분포는 u/U = f′(η) 이고, 벽 마찰 계수와 배제 두께는 다음과 같습니다.

Cf=2τwρU2=0.664Rex,δ∗=1.7208νxUC_f = \frac{2\tau_w}{\rho U^2} = \frac{0.664}{\sqrt{\mathrm{Re}_x}}, \qquad \delta^* = 1.7208\sqrt{\frac{\nu x}{U}}

이 방정식을 사격법으로 풀어 f″(0) = 0.332057, δ* 계수 1.7208, η₉₉ = 4.91 을 얻었고 문헌값과 일치합니다. Blasius 해는 Re_x → ∞ 의 점근해라 Re_x 500 ~ 1500 에서는 앞전 · 배제 유동의 고차 보정이 Re_x^(−1/2) ≈ 3 % 차수로 남습니다.

계산 조건

항목Boreas 설정
영역2 차원 slab x ∈ [0, 1.2] m, y ∈ [0, 0.3] m
경계x = 0 속도 입구 U = 1 m/s, 바닥 x < 0.2 m 는 slip (symmetry), 평판 x ∈ [0.2, 1.2] m 는 no-slip, 출구와 위 경계는 pressure_outlet
유체ν = 5e-4 m²/s (평판 끝 Re_L = 2000) — 셀 Reynolds 수 Uh/ν 를 4 · 2 · 1 로 두어 upwind 수치 점성 (≈ Uh/2) 을 물리 점성 이하로 유지
격자균일 h = 2 · 1 · 0.5 mm (0.18 · 1.08 · 4.3 M 셀)
시간dt = 0.004 s, 6 s (5 통과 시간)
이류semi-Lagrangian 과 upwind 둘 다
비교앞전에서 0.25 · 0.5 · 0.75 m (Re_x ≈ 500 · 1000 · 1500) 의 셀 열. 국소 경계층 바깥 속도 u_e (열의 최대 u) 로 정규화해 프로파일 u/u_e vs f′(η) 의 RMS (η < 8), 첫 셀의 C_f, δ* = Σ(1 − u/u_e)h

경계층 바깥 속도 u_e 는 U 보다 3 ~ 4 % 큽니다. 속도 입구가 앞전에서 0.2 m 밖에 떨어져 있지 않아 앞전의 정체 압력이 입구까지 닿고, 바닥 쪽 유선이 더 큰 전수두를 갖고 들어오기 때문입니다 (아래 유의 사항). 측정된 경계층을 다룰 때와 같이 비교는 모두 국소 u_e 로 정규화했습니다.

속도 분포

세 지점의 Boreas 속도 분포와 Blasius 상사해 비교

h = 0.5 mm 격자에서 Re_e ≈ 500 · 1000 · 1500 세 지점의 속도 분포 (기호, 세 셀에 하나) 와 Blasius 해 (실선). 위: semi-Lagrangian, 아래: upwind. 세 지점이 한 곡선 위에 겹쳐 상사성이 재현됩니다.

벽 마찰 계수

Re_e 에 대한 Boreas 벽 마찰 계수와 Blasius 비교

세 격자 × 두 이류 스킴의 C_f (첫 셀의 벽 전단에서) 와 Blasius 0.664/√Re_e. upwind 는 2 mm 부터 격자에 무관하고, semi-Lagrangian 은 격자를 세분하면 upwind 에 다가갑니다.

격자 수렴

격자 세분에 대한 프로파일 RMS 오차와 C_f 오차의 수렴

Re_e ≈ 1000 지점의 프로파일 RMS 오차 (위, 점선은 1차 기울기) 와 C_f 오차 (아래).

이류h [mm]u_e/U프로파일 RMS / u_e [%]C_f 오차 [%]δ* 오차 [%]
upwind21.0370.27−2.0+0.8
upwind11.0370.19−2.1+0.3
upwind0.51.0380.19−2.1+0.1
semi21.0370.90−5.1+2.9
semi11.0370.43−3.4+1.2
semi0.51.0360.25−2.6+0.4

(Re_e ≈ 1000, 앞전에서 0.5 m.) upwind 는 셀 Reynolds 수 ≤ 4 인 이 조건에서 2 mm 격자부터 이미 수렴해 있고, 남는 C_f −2 % 는 격자 오차가 아니라 약한 역압력구배와 저 Re 고차 보정의 몫입니다. semi-Lagrangian 은 격자 1차로 수렴해 0.5 mm 에서 upwind 와 같은 값에 이릅니다.

수렴 이력

h = 0.5 mm 의 세 지점 프로브 속도 시간 이력

세 지점 위 y = 0.02 m 프로브의 속도 이력 (upwind 실선, semi-Lagrangian 점선). 앞전에 가까운 지점부터 차례로 정상에 닿으며, t = 5 → 6 s 의 변화는 0.06 · 0.3 · 1 % 입니다 — 가장 하류 지점은 6 s 에서 아직 미세하게 움직입니다.

정량 결과

항목값
프로파일 RMS / u_e, upwind h = 0.5 mm, Re_e ≈ 1000 [%]0.19
C_f 오차, upwind, Re_e ≈ 1000 [%]−2.1
δ* 오차, upwind, Re_e ≈ 1000 [%]+0.1
C_f 오차, upwind, Re_e ≈ 500 · 1500 [%]+0.8 · −3.6
semi-Lagrangian 격자 수렴 차수 (h 2 → 0.5 mm)1
Blasius f″(0) (사격법, 문헌 0.33206)0.332057

유의 사항

  • 경계층 바깥 속도 u_e > U. 속도 입구가 앞전에 가까운 (0.2 m) 이 덱에서 u_e/U 는 1.03 ~ 1.04 이며 격자와 무관합니다. 영역 높이를 0.3 → 1.2 m 로 늘리면 오히려 커지고 (1.056), 상류 slip 구간을 0.2 → 0.6 → 1.2 m 로 늘리면 1.05 → 1.02 → 1.007 로 1 에 다가갑니다 — 위 경계가 아니라 입구 근처에서 깨진 균일 자유류 가정이 원인입니다. u_e 로 정규화하면 프로파일 · δ* · C_f 의 결론은 이 설정에 무관합니다.
  • C_f 는 하류로 갈수록 낮게 나옵니다 (Re_e 500 · 1000 · 1500 에서 +0.8 · −2.1 · −3.6 %, upwind). u_e 가 1.041 → 1.032 로 줄어드는 약한 역압력구배와 Blasius 의 저 Re 보정이 겹친 범위입니다.
  • GPU 만 계산했습니다. 이 케이스는 CPU 결과와 대조하지 않았습니다. 같은 v0.10.0 빌드의 GPU · CPU 대조는 Couette 유동과 정사각 덕트 페이지에 있습니다.