채널 교차점에서의 침식·퇴적 거동과 형상 파라미터 영향 분석

초록

본 연구는 유한체적법(FVM) 기반 FORTRAN 코드를 개발하여 두 채널 교차점에서의 흐름과 퇴적 현상을 모사한다. 주 채널·보조 채널 폭비, 교차각, 유량비, 입구 프루드 수 등 형상·수리 파라미터가 퇴적 높이와 유속 분포에 미치는 영향을 실험 데이터와 비교 검증하였다. 결과는 보조 채널 유량비가 0.51에 이를 때 최대 퇴적 높이가 32 % 증가하고, 프루드 수 비가 감소하면 퇴적 높이가 감소함을 보여준다.

상세 분석

본 논문은 하천공학에서 가장 빈번히 발생하는 침식·퇴적·형상 변화 문제를 수치해석적으로 접근한다. 저자는 유한체적법(Finite Volume Method, FVM)을 기반으로 한 2차원 비정압 흐름·퇴적 모델을 FORTRAN으로 구현하였다. 모델은 연속 방정식, 운동량 방정식, 그리고 비균질 퇴적-침식 방정식을 결합하여, 물-퇴적 상호작용을 시간·공간적으로 해석한다. 격자 생성은 구조화된 직교 격자를 사용했으며, 대류항은 2차 고차 스킴, 확산항은 중앙 차분으로 처리해 수치적 안정성을 확보하였다. 시간 적분은 1차 명시적 Euler 방법을 적용했으며, Courant–Friedrichs–Lewy(CFL) 조건에 따라 시간 단계가 자동 조정된다.

실험 검증을 위해 저자는 물리 모델(스케일 1:50)에서 두 채널이 교차하는 구조를 제작하고, 입구 유량·프루드 수를 조절하여 다양한 조건을 실험하였다. 실험에서 측정된 수위·퇴적 두께 프로파일을 수치 결과와 비교했으며, 평균 오차는 5 % 이하로, 모델의 신뢰성을 입증하였다.

파라미터 스터디는 네 가지 주요 변수에 초점을 맞춘다. 첫째, 보조 채널 폭을 주 채널 폭의 0.3, 0.5, 0.7배로 변환한 경우, 보조 채널 폭이 커질수록 교차점 근처의 유속이 감소하고 퇴적량이 증가하는 경향을 보였다. 둘째, 교차각을 30°, 45°, 60°로 변화시켰을 때, 각이 클수록 흐름이 더 부드럽게 전이되어 속도 구배가 완화되고, 결과적으로 퇴적 고도가 낮아졌다. 셋째, 보조 채널 유량비(Q₂/Q₁)를 0.2, 0.35, 0.51로 변화시킨 결과, 유량비가 0.51에 도달했을 때 최대 퇴적 높이가 평균 32 % 증가하였다. 이는 보조 채널이 주 흐름을 크게 방해하면서 저류 현상이 강화되기 때문이다. 넷째, 입구 프루드 수(F₁/F₂) 비를 0.8, 1.0, 1.2로 조정한 결과, 프루드 수 비가 감소하면 흐름의 관성력이 약해져 퇴적이 억제되는 현상이 관찰되었다.

또한, 저자는 유속·퇴적 프로파일이 교차각에 따라 일정한 형태를 유지한다는 점을 강조한다. 이는 교차각이 변하더라도 흐름 재배치 메커니즘이 근본적으로 동일하다는 것을 의미한다. 모델의 한계로는 3차원 효과와 퇴적 물질 입도 분포를 고려하지 않은 점, 그리고 경계 조건을 단순화한 점을 들었다. 향후 연구에서는 3D CFD와 다상 흐름 모델을 결합하고, 입도별 침식·퇴적 모델을 도입하여 정확도를 높일 필요가 있다.