수송·침식 모델링 혁신: 새로운 혼합길이와 혼합속도 프로파일

초록

본 논문은 von Kármán 유사성 가설을 확장한 새로운 혼합길이(lₘ)와 혼합속도(uₘ) 프로파일을 제안한다. 네 가지 기존 난류 점성 모델과 비교·검증을 위해 균일 흐름에서의 모래바닥 침식, 2차원 흐름 내 농도 확산, 유럽 SANDPIT 프로젝트의 채취구 채우기 세 가지 사례를 수치해석하였다. 유한요소 기반 h‑s 스펙트럼 방법과 Lax‑Wendroff‑Richtmeyer·암시적 Euler 스킴을 사용했으며, 실험 데이터와의 비교에서 제안 모델이 가장 높은 정확도를 보였다.

상세 분석

논문은 먼저 기존의 난류 확산 모델이 실제 현장에서 관측되는 혼합속도의 지수적 감소를 충분히 반영하지 못한다는 점을 지적한다. 이를 보완하기 위해 von Kármán의 스케일 불변성 가정을 기반으로 혼합길이 lₘ(z)를 일반 함수 f(z) 형태로 도입하고, 이를 통해 혼합속도 uₘ(z)=k·C·vₘ(z)와 난류 점성 νₜ(z)=lₘ(z)·uₘ(z)를 유도한다. 네 가지 모델을 정리한 표에서 기존 모델 1~3은 각각 선형·쌍곡선·지수형 lₘ을 사용하지만, 제안 모델(모델 4)은 지수형 uₘ와 하이브리드 형태의 lₘ을 결합해 로그형 평균속도 프로파일을 재현한다. 특히 모델 4는 식(22)·(23)에서 제시된 z₀(거칠기)와 보정계수 C₁을 포함해 실험적 경계조건을 만족한다.

수치해석에서는 2‑D 유한요소(h‑s) 모델을 채택해 수평 방향은 전통적인 2‑D 요소, 수직 방향은 레전드르 다항식 기반의 직교함수 전개를 사용하였다. 흐름 해석은 암시적 Euler 스킴으로, 침식·침전은 Lax‑Wendroff‑Richtmeyer 스킴으로 처리하였다. 첫 번째 테스트(균일 흐름에서 모래바닥 침식)에서는 초기 농도 C=0, 경계조건 C(zₐ)=Cₐ, 상부 무플럭스 조건을 적용했으며, 모델 4가 실험 데이터(Van Rijn, 1985)와 거의 일치함을 확인했다. 보정계수 C₁을 0.7으로 조정하면 평균 농도 오차가 최소화된다.

두 번째 테스트에서는 수직 확산계수 εₛ를 상수·이차형·모델 4(지수형)로 바꾸어 농도 프로파일의 민감도를 조사하였다. 결과는 모델 4가 확산 경계층을 보다 정확히 재현하고, 특히 x=0.5 m, 1 m, 2 m 지점에서 실험적 등농도선과 일치한다는 점을 보여준다. 마지막으로 SANDPIT 프로젝트의 채취구 채우기 사례에서는 실제 현장 측정값과 비교했을 때, 모델 4가 침전량과 시간 변화를 가장 잘 예측한다. 전체적으로 제안된 혼합길이·속도 프로파일은 난류 점성 모델링의 물리적 근거를 강화하고, 실무 엔지니어가 현장 데이터를 기반으로 모델을 보정할 수 있는 실용적인 프레임워크를 제공한다.