[논문 리뷰] Flexible Macroscopic Models for Dense-Fluid Shockwaves: Partitioning Heat and Work; Delaying Stress and Heat Flux; Two-Temperature Thermal Relaxation

Wm. G. Hoover, Carol G. Hoover외 1명|arXiv (Cornell University)|2010. 05. 10.

Gas Dynamics and Kinetic Theory참고 문헌 21인용 수 7

한 줄 요약

이 논문은 밀도 높은 유체 충격파를 위한 유연한 거시적 모델을 제안하며, 종방향 및 횡방향 온도를 구분하고, 응력, 열류, 열에너지에 대한 세 가지 회복 시간을 포함하며, 일과 열 기여를 분할한다. 이 모델은 초순수한 특성—예를 들어 응력 및 열류 지연, 텐서 온도, 이온도 온도 열화학적 회복—을 맥스웰-카타네오 및 크룩-유형의 회복을 포함하는 쌍곡형 편미분 방정식계를 통해 재현하며, 수치 시뮬레이션에서 안정적이고 물리적으로 일관된 충격파 프로파일을 달성한다.

ABSTRACT

Macroscopic models which distinguish the longitudinal and transverse temperatures can provide improved descriptions of the microscopic shock structures as revealed by molecular dynamics simulations. Additionally, we can include three relaxation times in the models, two based on Maxwell's viscoelasticity and its Cattaneo-equation analog for heat flow, and a third thermal, based on the Krook-Boltzmann equation. This approach can replicate the observed lags of stress (which lags behind the strain rate) and heat flux (which lags behind the temperature gradient), as well as the eventual equilibration of the two temperatures. For profile stability the time lags cannot be too large. By partitioning the longitudinal and transverse contributions of work and heat and including a tensor heat conductivity and bulk viscosity, all the qualitative microscopic features of strong simple-fluid shockwave structures can be reproduced.

연구 동기 및 목표

• 분자역학 시뮬레이션에서 관찰된 밀도 높은 유체 충격파의 비평형 거동을 포괄하는 거시적 모델을 개발하는 것.
• 클래식한 뉴턴-스토크스 및 버너 방정식의 한계를 보완하기 위해 응력과 열류에 대한 국소적이지 않고 지연된 반응을 포함하는 것.
• 강력한 충격파에서 종방향 및 횡방향 온도가 별개로 존재하는 물리적 현실을 모델링하여, 미세구조적 관측과 일치시키는 것.
• 응력과 열류에 대한 맥스웰-카타네오 유형의 모델에서 유래한 두 개의 회복 시간과 열화학적 회복에 대해 크룩-볼츠만 방정식에서 기인한 하나의 회복 시간을 도입함으로써 수치적 안정성과 물리적 일관성을 확보하는 것.

제안 방법

• 질량, 운동량, 에너지에 대한 다섯 개의 쌍곡형 편미분 방정식계를 구성하며, 응력과 열류를 독립 변수로 포함한다.
• 세 가지 회복 시간을 도입한다: 응력에 대한 τσ (맥스웰 유형의 점탄성), 열류에 대한 τq (카타네오 방정식 유사), 열화학적 회복에 대한 τT (크룩-볼츠만 모델).
• 텐서 형태의 열전도도와 체적 점성도를 사용하여 일과 열 기여를 종방향 및 횡방향 성분으로 분할한다.
• 이차 공간 정확도와 사차 수준 룬게-쿠타 시간 적분을 사용하는 스타게드 유한차분 스킴을 적용하여 시스템을 해결한다.
• 입자 기반 평균화에서의 비국소성을 다루기 위해 루시 유형의 가중치 함수를 사용한 공간 평균화를 적용하여 원자적 시뮬레이션과의 일관성을 확보한다.
• 경계 조건은 후고니용 관계에 기반하며, 양끝에 냉각 및 난방 유체 상태를 설정하고, 경계에서 열류 및 비틀림 응력이 0이 되도록 보장한다.

실험 결과

연구 질문

• RQ1 밀도 높은 유체 충격파에서 변형률률과 온도 기울기의 변화에 비해 응력과 열류의 관찰된 지연을 거시적 모델이 어떻게 포괄할 수 있는가?
• RQ2 스칼라 온도 근사 이상의 충격파 구조를 정확히 기술하는 데 있어 텐서 형태의 온도(종방향 대 횡방향)는 어떤 역할을 하는가?
• RQ3 세 가지 서로 다른 회복 시간을 가진 통합 모델이 강력한 충격파의 분자역학 시뮬레이션에서 관찰된 전반적인 정성적 특징을 재현할 수 있는가?
• RQ4 일과 열을 종방향 및 횡방향 성분으로 분할하는 것이 충격파 프로파일의 안정성과 형태에 어떤 영향을 미치는가?
• RQ5 회복 시간에 대한 안정성 제약 조건은 무엇이며, 이는 물리적으로 타당한 해의 허용 가능한 매개변수 공간을 어떻게 제한하는가?

주요 결과

• 모델은 분자역학 시뮬레이션과 일치하는 수준의 응력 및 열류 지연을 성공적으로 재현하였으며, 이는 입자 충돌 시간과 같은 주기의 지연을 보였다.
• 응력과 열류가 변형률률 및 온도 기울기의 구동 기여에 뒤따르는 지연을 보이며, 강력한 충격파에서 인과관계 기반의 회복 현상을 확인하였다.
• 충격파 전면에서 종방향 및 횡방향 온도가 서로 다른 방식으로 변화하며, 횡방향 온도가 먼저 상승하고, 회복 후에 두 온도가 통합된다.
• 회복 시간이 너무 크지 않은 경우에만 수치적 해가 안정되며, 이는 비평형 상태에서 유체의 기억 능력에 물리적 한계가 있음을 시사한다.
• 텐서 열전도도와 체적 점성도의 포함으로 비평형 온도 이방성까지 포함한 충격파의 전반적인 정성적 구조를 정확히 재현할 수 있었다.
• 모델은 광범위한 매개변수 범위에서 안정적이고 수렴 가능한 해를 달성하였으며, 적응형 시간 간격과 회복 매개변수의 인위적 시간 의존성으로 인해 수치적 안정성이 향상되었다.