[논문 리뷰] Onsager reciprocity principle for kinetic models and kinetic schemes

A.K. Mahendra, Ramesh Singh|arXiv (Cornell University)|2013. 08. 19.

Gas Dynamics and Kinetic Theory참고 문헌 47인용 수 8

한 줄 요약

이 논문은 비평형 열역학과 온제거의 상호성 원리를 기반으로 한 새로운 운동론적 모델을 제안하여 운동론적 시뮬레이션에서 정확한 운반 계수와 엔트로피 생성을 보장한다. 충돌 확률 함수 $\mathbb{P}(\bar{y})$ 를 통합함으로써 희박유체 및 연속체 유체를 위한 단일 통합 솔버를 가능하게 하며, 비용이 많이 드는 다중 척도 결합 없이도 스케일을 초월해 슬립 유동을 정확히 포착하는 운동론적 벽 경계 조건을 제공한다.

ABSTRACT

Boltzmann equation requires some alternative simpler kinetic model like BGK to replace the collision term. Such a kinetic model which replaces the Boltzmann collision integral should preserve the basic properties and characteristics of the Boltzmann equation and comply with the requirements of non equilibrium thermodynamics. Most of the research in development of kinetic theory based methods have focused more on entropy conditions, stability and ignored the crucial aspect of non equilibrium thermodynamics. The paper presents a new kinetic model formulated based on the principles of non equilibrium thermodynamics. The new kinetic model yields correct transport coefficients and satisfies Onsager's reciprocity relationship. The present work also describes a novel kinetic particle method and gas kinetic scheme based on this linkage of non-equilibrium thermodynamics and kinetic theory. The work also presents derivation of kinetic theory based wall boundary condition which complies with the principles of non-equilibrium thermodynamics, and can simulate both continuum and rarefied slip flow in order to avoid extremely costly multi-scale simulation.

연구 동기 및 목표

• 볼츠만 방정식의 기본 성질을 유지하면서 비평형 열역학, 특히 온제거의 상호성과 엔트로피 생성 원칙을 준수하는 운동론적 모델을 개발하기 위해.
• 기존 운동론적 방법에서 비정상적인 열역학 일관성이 떨어지는 문제를 해결하기 위해, 일반적으로 엔트로피 생성을 우선시하면서도 공액 열역학력 간의 엔트로피 분포를 정확히 반영하지 못하는 경향이 있음을 해결하기 위해.
• 비평형 열역학을 충돌 모델에 통합하여 희박유체에서 연속체에 이르기까지 전체 유동 영역에서 원활하게 작동하는 통합 운동론 입자 방법과 기체 운동론적 스킴을 개발하기 위해.
• 비평형 열역학을 충족하고 다중 척도 모델링 없이도 연속체 및 희박유체 영역의 슬립 유동을 정확히 시뮬레이션할 수 있는 기반 운동론적 벽 경계 조건을 유도하기 위해.
• 연속체 및 희박유체 영역을 아우르는 시험 케이스를 통해 제안된 프레임워크의 정확성과 강건성을 검증하기 위해.

제안 방법

• 온제거의 변분 원리를 사용하여 운반 계수의 대칭성과 온제거 상호성 관계 준수를 보장하는 새로운 운동론적 모델을 수립한다.
• 정규 거리 $\bar{y} = y/\lambda_e$ 를 함수로 하는 충돌 확률 함수 $\mathbb{P}(\bar{y})$ 를 도입하여, 이는 이완 시간 $\tilde{t}_{R(j)} = \mathbb{P}(\bar{y})t_{R(j)}$ 를 조절한다.
• 새로운 모델에 기반한 운동론적 플럭스 벡터 분할 및 입자 방법을 유도하며, 점성 플럭스는 비평형 효과를 반영하기 위해 $\mathbb{P}(\bar{y})$ 로 스케일링된다.
• 맥스웰 기체 상호작용 모델을 사용하여 벽 경계 조건을 구성하며, 입사, 반사, 확산 분포를 조정 계수 $\sigma$ 를 통해 가중치를 부여한다. 비평형 부분은 벽에서 사라지며 ($\mathbb{P}(0) = 0$), 이는 비평형 효과의 감쇠를 의미한다.
• 확장된 비가역 열역학(EIT) 프레임워크를 사용하여 쿤들렌 층 내의 비선형 플럭스-력 관계를 모델링하며, EIT 플럭스는 $\boldsymbol{J}_i \approx \mathbb{P}(\bar{y}) \sum_j \boldsymbol{L}_{ij} \odot \boldsymbol{X}_j$ 로 근사한다.
• 쿤들렌 층 내에서 채플만-엔스코그 분포를 사용하며, 비평형 부분은 $f_1(\bar{y}) = f_0 - \mathbb{P}(\bar{y}) \sum_j \boldsymbol{\Upsilon}_j \odot \boldsymbol{X}_j$ 로 표현되어 열역학 일관성이 보장된다.

실험 결과

연구 질문

• RQ1 어떻게 하면 온제거의 상호성 원리를 만족하면서도 프란틀 수와 같은 정확한 운반 계수를 유지하는 운동론적 모델을 구축할 수 있는가?
• RQ2 비평형 열역학을 충돌 모델에 통합함으로써 하나의 운동론적 스킴이 희박유체와 연속체 유체를 정확히 시뮬레이션할 수 있는가?
• RQ3 충돌 확률 함수 $\mathbb{P}(\bar{y})$ 는 쿤들렌 층 전체에 걸쳐 점성 플럭스를 통합적으로 처리하는 데 어떤 역할을 하는가?
• RQ4 비평형 열역학을 충족하고 다중 척도 모델링 없이도 희박유체 및 연속체 영역의 슬립 유동을 정확히 모델링할 수 있는 기반 운동론적 벽 경계 조건을 어떻게 도출할 수 있는가?
• RQ5 제안된 모델은 응력과 열기울기와 같은 다양한 열역학력 간의 엔트로피 생성 분포를 얼마나 정확히 보장하는가?

주요 결과

• 제안된 운동론적 모델은 온제거의 상호성 관계를 만족하며, 정확한 프란틀 수를 생성하여 운반 성질의 열역학 일관성을 확보한다.
• 운동론적 스킴 내 점성 플럭스 성분은 충돌 확률 함수 $\mathbb{P}(\bar{y})$ 로 스케일링되어 쿤들렌 층 내 비평형 효과를 정확히 해석할 수 있다.
• 맥스웰 기체 상호작용 모델에서 유도된 벽 경계 조건은 정반사 및 확산 반사를 모두 정확히 포착하며, $\mathbb{P}(0) = 0$ 로 인해 비평형 부분이 벽에서 사라진다.
• 쿤들렌 층 내 상태 갱신 방정식은 $\mathbb{P}(\bar{y})$ 의 함수로 된 이완 시간을 포함하여 벽 근처에서 비평형에서 평형 행동으로의 부드러운 전이를 가능하게 한다.
• 새로운 모델에 기반한 운동론적 입자 방법과 플럭스 벡터 분할 스킴은 연속체 및 희박유체 슬립 유동 시험 케이스를 성공적으로 시뮬레이션하여 프레임워크의 강건성과 정확성을 입증한다.
• 비평형 열역학을 직접 운동론 프레임워크에 통합함으로써 다중 척도 시뮬레이션의 고비용을 피할 수 있으며, 모든 유동 영역에서 단일 솔버를 구현할 수 있다.