식도 내 비뉴턴 식품 볼루스의 파동식 연동 운송 메커니즘

본 논문은 식도 내 음식 볼루스의 연동 운송을 물리학적으로 해석하기 위해, 비뉴턴 유체인 오스트발드‑드‑와엘 파워‑로우 모델을 채택하였다. 연구의 주요 목표는 (1) 임의 파형과 유한 길이(250~300 mm)의 관에서 전파되는 단일 파동과 파동열의 압력·유속·입자 궤적을 비교 분석하고, (2) 흐름 지수 n이 압력 구배와 유량에 미치는 영향을 규명하는 것이다. **모델 설정** 식도는 원통형 관으로 가정하고, 방사형 좌표 (R,θ,Z)와 시간 t를 사용한다. 관의 반경 a는 최소 반경 ε와 평균 반경 H(Z,t) 사이에서 변동한다. 파동 길이 λ와 파동 진폭 φ를 정의하고, 파동 번호 δ = a/λ ≪ 1 로 설정해 윤활 이론을 적용한다. 유체는 일관성 계수 α와 흐름 지수 n을 갖는 파워‑로우이며, 전단 얇음(n<1)·뉴턴(n=1)·전단 두꺼움(n>1) 세 경우를 모두 고려한다. **지배 방정식 및 비차원화** 연속 방정식과 운동량 방정식을 원통 좌표계에서 전개하고, 레이놀즈 수가 매우 낮아 관성 항을 무시한다. 비차원화 변수( Z̄=Z/λ, R̄=R/a, Ū=U/c 등)를 도입해 차원 없는 형태로 정리하면, 축방향 압력 구배 p=∂P/∂Z와 전단 응력 |∂U/∂R|ⁿ⁻¹ 사이의 비선형 관계식(8)이 도출된다. 경계 조건은 관 중심에서 대칭성(∂U/∂R=0, V=0)과 벽면에서 무슬립( U=0, V=∂H/∂t )을 적용한다. **해석적 해** 식(8)을 적분하고 경계 조건을 적용하면, 축방향 속도 U(R,Z,t)와 방사형 속도 V(R,Z,t)가 (12)와 (13) 형태로 얻어진다. 여기서 k=1/n, p₁=k ∂p/∂Z / p 등 추가 파라미터가 등장한다. 압력 구배 p는 (15)식에 의해 H(Z,t)와 시간 미분 ∂H/∂t에 의존한다. 압력 차 ΔP는 (22)식으로 전체 관 길이에 걸쳐 적분된다. **수치 구현** 임의 파형을 사인·코사인 조합으로 표현하고, 시간에 따라 변하는 H(Z,t)를 전산적으로 계산한다. 비선형 방정식(15)·(18)에서 나타나는 적분 상수 c₁은 직접 해석이 불가능하므로, 반복적 뉴턴-라프슨 방법으로 수치적으로 구한다. 파동열(η=1)과 단일 파동(η=L/λ) 두 경우에 대해 압력 분포, 유속 프로파일, 입자 궤적을 시뮬레이션한다. **주요 결과** 1. **압력 민감도**: 전단 얇음 유체(n=0.5)에서는 압력 피크가 파동 진행 중반에 집중돼 역류 위험이 높다. 전단 두꺼움 유체(n=1.5)에서는 압력 피크가 파동 말기에 나타나며, 전체 압력 구배가 부드럽다. 뉴턴 유체(n=1)는 중간 특성을 보인다. 2. **유량·압력 차**: 파동열에서는 ΔP가 거의 일정하게 유지돼 평균 유량 Q̅가 크게 증가한다. 반면, 단일 파동에서는 파동 간 간격이 넓어 압력 회복이 충분히 이루어지지 않아 Q̅가 감소하고 역류가 발생한다. 3. **입자 궤적**: 전단 얇음 경우 입자들이 파동 전후로 급격히 뒤로 이동하는 현상이 관찰되며, 전단 두꺼움 경우는 전방 이동이 지배적이다. 파동열에서는 입자들이 연속적으로 전진해 전체 볼루스가 효율적으로 이동한다. 4. **파라미터 영향**: 파동 진폭 φ가 클수록 압력 구배가 커져 운송 효율이 상승하지만, 동시에 역류 위험도 증가한다. 파동 길이 λ가 짧을수록(δ 증가) 윤활 가정이 약해져 모델 정확도가 떨어진다. **임상·공학적 함의** - 식도 근육의 수축 파형을 조절하거나, 음식의 레오시티(전단 지수)를 조절함으로써 역류성 식도염 등 병리적 상황을 완화할 수 있다. - 인공 연동 펌프 설계 시, 파동열을 구현하면 뉴턴·비뉴턴 유체 모두에 대해 높은 운송 효율을 달성할 수 있다. - 파동 간격을 최적화하면 약물 전달 시스템에서 목표 부위로의 정확한 전달이 가능하다. **결론** 본 연구는 비뉴턴 유체의 전단 특성이 식도 내 연동 운송에 미치는 영향을 정량적으로 밝히고, 파동열이 단일 파동에 비해 압력·유량·역류 측면에서 현저히 우수함을 입증하였다. 향후 연구에서는 삼차원 비선형 변형, 점성·탄성 복합 모델, 그리고 실제 식도 근육 전기생리학적 데이터와의 연계가 필요하다.