[논문 리뷰] On Unsteady Heat Conduction in a Harmonic Crystal

А. М. Кривцов|arXiv (Cornell University)|2015. 09. 07.

Advanced Mathematical Modeling in Engineering참고 문헌 2인용 수 11

한 줄 요약

이 논문은 1차원 조화 결정에서 비정상 열전도를 모델링하기 위해 비국소 온도 형식을 제안하며, 유한한 열 전 fronts 속도와 길이에 독립적인 열류를 예측하는 시간 반전 가능 열파 방정식을 도출한다. 이는 푸아송의 법칙과는 근본적으로 대조된다. 이 방법은 분자역학 시뮬레이션과 해석적 추정치와 일치하는 구성 법칙을 도출하며, 저차원 시스템에서의 이방성 열전도를 드러낸다.

ABSTRACT

An analytical model of unsteady heat transfer in a one-dimensional harmonic crystal is presented. A nonlocal temperature is introduced as a generalization of the kinetic temperature. A closed equation determining unsteady thermal processes in terms of the nonlocal temperature is derived. For an instantaneous heat perturbation a time-reversible equation for the kinetic temperature is derived and solved. The resulting constitutive law for the heat flux in the considered system is obtained. This law significantly differs from Fourier's law and it predicts a finite velocity of the heat front and independence of the heat flux on the crystal length. The analytical results are confirmed by computer simulations.

연구 동기 및 목표

• 1차원 조화 결정에서 비정상 열전도의 엄밀한 거시적 기술을 개발하여 푸아송의 법칙을 넘어서는 것.
• 특히 나노구조에서 이방성 열전도가 관측되는 저차원 시스템에서 푸아송의 법칙의 실패를 다루는 것.
• 비국소 온도로 일반화된 운동온도의 일반화를 통해 열과정에 대한 폐쇄형 방정식을 유도하는 것.
• 유한 속도로 열이 전파되고 결정 길이에 영향을 받지 않는 열류의 구성 법칙을 수립하는 것.
• 분자역학 시뮬레이션과 이전의 해석 결과와의 비교를 통해 분석 모델을 검증하는 것.

제안 방법

• 입자 속도의 공분산으로 정의된 비국소 온도 θₙ(x)를 도입하여 표준 운동온도 T(x) = m⟨ẋ²⟩/k_B를 일반화한다.
• 비틀림 근사와 공간적 매끄러움 가정을 사용하여 θₙ(x,t)에 대한 폐쇄형 미분-차분 방정식을 유도한다: ÿₙ + ½c²(θₙ₋₁ − 2θₙ + θₙ₊₁)′′ = 0.
• 점근적 분석을 적용하여 미분-차분 방정식을 국소 운동온도 T(x,t)에 대한 편미분 방정식으로 간소화하여 시간 반전 가능 열파 방정식을 도출한다: ÿ = c²T′′.
• 순간적 열변동에 대해 해를 해석적으로 구해 자가유사 해 T = T(x/(ct))를 얻는다.
• 시스템 크기가 증가함에 따라 수렴하는 것을 확인하기 위해 최대 10⁸ 개 입자를 사용한 분자역학 시뮬레이션과 분석 해를 비교한다.
• 분석 해를 고전적 모델들(열방정식, 열파 방정식, 파동 방정식)과 비교하여 물리적 행동의 차이를 부각한다. 특히 x = ct에서 부드럽고 수직 탄젠트를 갖는 열 전 fronts 가 특징이다.

실험 결과

연구 질문

• RQ1 1차원 조화 결정에서 비정상 열전도는 비국소 온도 형식을 통해 어떻게 기술할 수 있는가?
• RQ2 조화 결정에서 열류의 구성 법칙은 어떤 형태이며, 푸아송의 법칙과 어떻게 다를까?
• RQ3 조화 결정에서 열 전 fronts 는 유한한 속도로 전파되는가? 그리고 열류는 결정 길이에 독립적인가?
• RQ4 분석 예측은 온도 프로파일과 열 전 fronts 다이내믹스 측면에서 분자역학 시뮬레이션과 어떻게 비교되는가?
• RQ5 조화 시스템에서 미세한 역학에서 엄밀하게 시간 반전 가능한 거시적 열전도 방정식을 유도할 수 있는가?

주요 결과

• 비국소 온도 형식은 운동온도를 성공적으로 일반화하며, 조화 결정에서 비정상 열과정에 대한 폐쇄형 방정식 유도를 가능하게 한다.
• 최종적으로 도출된 운동온도에 대한 거시적 방정식은 시간 반전 가능 열파 방정식 ÿ = c²T′′이며, 이는 유한한 열 전 fronts 속도 c를 예측한다.
• 해석적 해는 자가유사성 T = T(x/(ct))를 보이며, x = ct에서 수직 탄젠트를 갖는 부드러운 열 전 fronts 를 나타내어 시뮬레이션 결과와 일치한다.
• 열류의 구성 법칙은 열류가 결정 길이에 영향을 받지 않는다는 것을 예측하며, 1차원 조화 시스템에서의 열 초전도성과 일치한다.
• 모델은 고전적 열방정식의 무한한 전 fronts 속도와 열파 방정식의 단계적 전 fronts 와는 대조적으로, 유한한 열 전 fronts 속도를 예측한다.
• 최대 10⁸ 개 입자를 사용한 분자역학 시뮬레이션은 분석 해와 뛰어난 일치를 보이며, 시스템 크기가 증가함에 따라 수렴함을 확인한다.