열 흐름의 와류와 역류: 유체역학적 열전달의 새로운 해석

초록

**
본 논문은 점탄성 열방정식(VHE)을 두 개의 변형된 바이헐름 방정식으로 분리하고, 속도 퍼텐셜과 스트림 함수로 해석한다. 이를 통해 온도 변화를 와류와 압축성 두 기여로 나누어 설명하고, 2D 흑연 스트립에서 열 와류와 음의 열저항(역류)을 예측한다.

**

상세 분석

**
본 연구는 최근 실험적으로 확인된 포논 수소역학 현상을 이론적으로 정밀히 다루며, 점탄성 열방정식(VHE)을 기존의 나비에-스토크스(NSE)와 유사한 형태로 제시한다. 핵심은 VHE를 두 개의 변형된 바이헐름 방정식으로 완전히 분리한다는 점이다. 첫 번째 방정식은 압축성(∇·u)과 관련된 열압축성 Φ에 대한 헬름홀츠형식이며, 두 번째는 회전성(∇×u)인 와류 W에 대한 변형된 헬름홀츠 방정식이다.

이 두 스칼라 양은 각각 속도 퍼텐셜 ϕ와 스트림 함수 ψ를 통해 표현된다. u = −∇ϕ + ∇×ψ 라는 헬름홀츠 분해를 이용해, ϕ와 ψ는 4차 미분 연산자를 포함한 변형 바이헐름 방정식(식 5)을 만족한다. 이는 기존의 비압축성 Stokes 흐름에서 ∇⁴ψ = 0인 경우를 γ→0 한계로 복원한다.

압축성 항은 포논 시스템에서 일반적으로 무시할 수 없는데, 전자 수소역학과 달리 포논은 에너지-결정운동량 관계에 의해 열압축성을 갖는다(α, β 계수). 따라서 Φ와 W는 온도 프로파일 T(x,y) = T_ϕ + T_ψ 로 각각 기여한다. T_ϕ는 압축성에 의해 발생하고, T_ψ는 와류에 의해 발생한다는 점은 온도 분포를 직접 측정함으로써 두 현상을 구분할 수 있는 새로운 실험적 지표를 제공한다.

논문은 2D 흑연 스트립(폭 h)에서 점 전류와 온도 구배를 주입하는 점 접촉을 고려한다. 경계조건은 u_y = U δ(x) (주입/추출), u_x = 0 (무미끄럼), 그리고 T = \bar T ± ΔT δ(x) 로 설정한다. Fourier 변환을 이용해 ψ(k,y)와 ϕ(k,y)를 일반 해로 구하고, q_ψ² = k² + γ/η, q_ϕ² = k² + γ/η + (μ+αβκ)/(η+μ) 로 정의한다.

온도 프로파일은 두 항으로 나뉘며, 식 12에서 정의된 Fourier Deviation Number(FDN) = 1/(ε+ξ) 로 비탄성 효과를 정량화한다. ε = γh²/η는 Umklapp에 의한 모멘텀 소산, ξ = καΔT h/U는 드리프트 속도 비중을 나타낸다. ε와 ξ가 작을수록(즉, 약한 소산·강한 드리프트) FDN이 커져 Fourier 법칙에서 크게 벗어나며, 온도 분포는 역방향(음의 열저항)으로 전환된다.

특히 ξ→0 한계에서 온도는 1/x² 형태의 장거리 감소와 y = ±x 선을 따라 노달 라인이 형성되는 것을 보인다. 이는 스트림 라인이 중앙 채널과 양측 와류 영역으로 분리되는 모습과 일치한다. 실험적으로는 온도 차이 ΔT ≈ 1 K, 드리프트 속도 U ≈ 2×10⁴ m/s 이상이 필요하며, 동위 원소 정제된 흑연(¹²C 99.95%)을 사용하면 백플로우가 0.4 K까지 증폭된다.

결과적으로, 압축성·와류가 동시에 존재하는 포논 유체는 전자 유체와 달리 복합적인 복사와 전도 메커니즘을 보이며, 복소 퍼텐셜 χ = iψ − ϕ 로 정의된 스트림 라인은 압축성(ϕ)와 회전성(ψ) 모두에 의존한다. 이는 기존의 비압축성 가정이 불가능함을 의미하고, 열 와류와 역류 현상을 설계·제어하기 위한 새로운 이론적 도구를 제공한다.

**