Multispecies Bhatnagar-Gross-Krook models and the Onsager reciprocal relations

본 논문은 다성분 계를 위한 기존 속도 무관 Bhatnagar-Gross-Krook 모델 대부분이 온저가 상호역관계식을 만족하지 못함으로써 특정 유체와 일치하도록 그 수송 특성을 보정하는 작업을 복잡하게 만든다는 것을 보여준다.

핵심

두 가지 다른 유형의 구슬 (빨간색과 파란색이라고 가정해 봅시다) 이 상자를 흔들 때 어떻게 섞이는지 예측한다고 상상해 보세요. 과학자들은 이러한 구슬의 거동에 대한 "규칙집"을 가지고 있는데, 이를 볼츠만 방정식이라고 합니다. 이는 놀라울 정도로 정확하지만, 동시에 너무 복잡하여 해를 구하는 것은 허리케인이 불고 있는 동안 해변의 모래알 하나하나를 세어 보려는 것과 같습니다.

문제를 쉽게 만들기 위해 과학자들은 이 규칙집의 단순화된 버전인 BGK 모델 (바트나가르, 그로스, 크룩의 이름을 따서 명명됨) 을 만들었습니다. BGK 모델을 생각할 때, 복잡한 구슬의 거동을 모든 무거운 수학을 수행하지 않고 근사화하는 "치트시트"나 "단순화된 방법"으로 간주할 수 있습니다. 이 모델은 날개 위를 흐르는 공기부터 별 속의 플라즈마에 이르기까지 모든 것을 시뮬레이션하는 데 수십 년간 사용되어 왔습니다.

문제: 고장난 나침반

E. S. 베닐로프가 쓴 이 논문은 서로 다른 기체 (예: 산소와 질소, 또는 수증기와 공기) 의 혼합물을 다룰 때 이 "치트시트"의 가장 일반적인 버전에서 발견된 치명적인 결함을 지적합니다.

저자는 이러한 인기 있는 BGK 모델들이 온사거 상호관계라는 물리학의 근본 법칙을 위반한다는 사실을 발견했습니다.

온사거 관계에 대한 간단한 비유:
두 도시 사이를 오가는 교통이 있는 양방향 도로를 상상해 보세요.

• 규칙 A: A 도시에 언덕을 만들면 B 도시의 차량 속도에 영향을 미칩니다.
• 규칙 B: B 도시에 언덕을 만들면 A 도시의 차량 속도에 영향을 미칩니다.

온사거 관계는 이 두 가지 효과가 서로 완벽하게 "조정"되어야 한다고 말합니다. A 도시의 언덕이 B 도시의 교통을 10% 늦춘다면, B 도시의 언덕은 수학적으로 연결된 양만큼 A 도시의 교통을 늦추어야 합니다. 이는 대칭성의 규칙입니다. 우주는 이러한 상호작용이 균형을 이루도록 요구합니다.

논문이 발견한 것

베닐로프는 이 규칙에 대해 표준 BGK "치트시트"를 테스트했습니다. 그는 다음과 같은 사실을 발견했습니다.

1. 모델의 비대칭성: BGK 모델에서 A 도시의 "언덕" (온도 변화) 은 B 도시의 "교통" (질량 흐름) 에 영향을 전혀 주지 않습니다. 반면, B 도시의 "언덕" (밀도 변화) 은 A 도시의 "교통" (열 흐름) 에 영향을 미칩니다.
2. 불일치: 방정식의 한쪽은 0 이고 다른 한쪽은 0 이 아니기 때문에 대칭성이 깨집니다. 이 모델은 영구적으로 한쪽으로 기울어진 저울과 같습니다.
3. 결과: 이 모델이 이 근본적인 규칙을 위반하기 때문에 이를 "보정"하는 것은 불가능합니다. 보정은 명확한 신호를 받기 위해 라디오를 튜닝하는 것과 같습니다. 특정 유체에 대한 실제 데이터와 일치하도록 BGK 모델의 노브 (매개변수) 를 조정해 보려 해도 불가능합니다. 이 모델은 어떻게 조정하든 결코 완벽하게 정확해질 수 없도록 근본적으로 고장 난 상태입니다.

"수증기" 예외 (그리고 왜 이것이 상황을 구하지 못하는지)

"음, 아마도 이는 이상한 기체에만 중요할지도 모릅니다. 수증기와 공기 같은 일반적인 혼합물은 어떨까요?"라고 생각할 수 있습니다.

논문은 이것도 확인했습니다. 온도가 질량 흐름에 미치는 영향이 미미하다 하더라도 (수증기와 공기의 경우 그렇습니다), 모델은 여전히 실패합니다. 이 특정 경우에 모델을 작동시키려면 다이얼을 무한대로 돌려야 하는데, 이는 모든 운동을 정지시켜 모델을 완전히 무효화하는 효과를 냅니다. 따라서 이 모델은 복잡한 혼합물뿐만 아니라 단순한 혼합물에서도 실패합니다.

좋은 모델은 존재할까?

논문은 규칙을 따르는 몇 가지 더 복잡한 BGK 모델이 존재한다고 지적하지만, 이들 역시 자체적인 문제점 (예: 엔트로피가 항상 증가하도록 보장하는 "H-정리"와 같은 다른 물리 법칙을 위반하는 것) 을 가지고 있습니다.

저자는 결론적으로, 현재로서는 존재하는 BGK 모델 중 어느 것도 완벽하지 않다고 결론 내립니다. 완벽한 모델은 다음을 충족해야 합니다.

• 질량, 운동량, 에너지를 보존해야 합니다.
• 열역학 법칙 (엔트로피) 을 따라야 합니다.
• 동일한 입자를 공정하게 대우해야 합니다.
• 온도와 농도를 양수로 유지해야 합니다.
• 과학자들이 어떤 실제 유체와도 일치하도록 조정할 수 있어야 합니다.
• 그리고 온사거 상호관계 (대칭성 규칙) 를 준수해야 합니다.

현재 우리가 가진 모든 모델은 적어도 이 중 하나의 테스트에서 실패합니다.

결론

이 논문은 이러한 모델을 사용하는 과학자들에게 경고하는 것입니다. 가스 혼합물을 시뮬레이션하기 위해 표준 BGK 모델을 사용한다면, 당신은 물리 법칙과 근본적으로 "불일치"하는 도구를 사용하고 있는 것입니다. 이는 일어나는 일에 대한 대략적인 아이디어를 제공할 수는 있지만, 우주의 핵심 대칭 규칙을 위반하기 때문에 실제 유체에 대한 정밀하고 보정된 결과를 신뢰할 수 없습니다. 저자는 미래에 이러한 모든 문제를 해결하는 "완벽한" 모델을 누군가가 만들기를 희망합니다.

기술 요약

기술적 요약: 다종 Bhatnagar–Gross–Krook 모델과 Onsager 상호역 관계

문제 제기
Bhatnagar–Gross–Krook (BGK) 모델은 특히 기체 혼합물에 대해 볼츠만 운동론 방정식의 현상론적 근사로 널리 사용됩니다. 원래의 단일 종 모델과 혼합물에 대한 그 확장 (예: Gross 와 Krook [2]) 은 정성적 정확성과 적응성으로 가치 있게 평가받지만, Onsager 상호역 관계 준수를 엄격하게 검증받은 바는 없습니다. 이러한 관계는 수송 계수에 근본적인 열역학적 제약을 부과하며, 구체적으로 질량 플럭스 내의 온도 구배 계수 (열확산 또는 소레 효과) 와 열 플럭스 내의 밀도 구배 계수 (듀푸르 효과) 를 서로 연결합니다. 첫 원리에서 유도된 모델은 이러한 관계를 자동으로 만족하는 반면, 현상론적 모델은 그렇지 않을 수 있습니다. 모델이 이러한 관계를 준수하지 않는다는 사실은 그 모델의 물리적 관련성에 의문을 제기할 뿐만 아니라, 결정적으로 특정 실제 유체의 수송 특성과 일치하도록 모델 매개변수를 보정하는 것을 불가능하게 만듭니다.

방법론
본 논문은 단원자 기체의 이원 혼합물에 대한 가장 일반적인 다종 BGK 모델을 수립합니다. 이 수립은 충돌 빈도 ($\nu_{ij}$) 와 종 의존 속도와 온도를 갖는 혼합 맥스웰 분포 ($M_{ij}$) 로의 이완을 정의하는 매개변수 ($\alpha_i, \beta_i, \gamma_i$) 를 포함합니다. 저자는 이러한 매개변수에 대한 특정 제약을 부과함으로써 모델이 질량, 운동량, 에너지를 보존하도록 보장합니다.

수송 특성을 분석하기 위해, 본 논문은 표준 유체역학적 근사 대신 확산 근사를 사용합니다. 이는 다음을 포함합니다:

1. 스케일링: 느리고 약한 확산 흐름에 적합한 작은 매개변수 $\epsilon$을 도입하여 공간 좌표 ($\nabla \to \epsilon \nabla$) 와 시간 ($\partial_t \to \epsilon^2 \partial_t$) 을 확장합니다.
2. 전개: 분포 함수, 거시적 속도, 온도, 밀도를 $\epsilon$의 거듭제곱으로 전개합니다.
3. 유도: 이러한 전개를 재스케일링된 BGK 방정식에 대입하여 분포 함수의 1 차 보정 ($f^{(1)}_i$) 을 풉니다.
4. 플럭스 계산: $f^{(1)}_i$를 적분하여 밀도와 온도 구배에 대한 질량 플럭스 ($J_i$) 와 열 플럭스 ($Q$) 의 명시적 표현식을 유도합니다.
5. 비교: 유도된 BGK 플럭스를 Onsager 관계를 만족하는 표준 유체역학적 표현식 (Enskog–Chapman 이론) 과 비교합니다.

주요 기여 및 결과
본 논문의 중심 결과는 분자 속도와 무관한 충돌 빈도 및 모델 계수를 갖는 가장 일반적인 다종 BGK 모델 버전이 Onsager 상호역 관계를 만족하지 않는다는 것을 입증하는 것입니다.

구체적으로, 분석은 근본적인 불일치를 드러냅니다:

• Onsager 관계에 따르면, 온도 구배를 질량 플럭스와 연결하는 계수 (열확산율, $B$) 는 밀도 구배를 열 플럭스와 연결하는 계수 (듀푸르 효과, $C$) 가 0 이 아닐 때 0 이 아니어야 하며, 두 계수는 $C = pB$로 연결됩니다.
• BGK 모델 유도는 열확산 항이 엄격히 0 ($B=0$) 인 질량 플럭스를 산출하며, 이는 특정 비물리적 극한을 취하지 않는 한 성립합니다.
• 그러나 결과적으로 산출된 열 플럭스에는 밀도 구배에 비례하는 항 (듀푸르 효과) 이 포함되어 있으며, 이는 0 이 아니며 종의 질량 차이에 비례합니다.
• 결과적으로 BGK 모델은 0 인 소레 효과와 함께 0 이 아닌 듀푸르 효과를 예측하여 Onsager 관계 $C = pB$를 위반합니다.

본 논문은 또한 이 불일치가 모델 매개변수 (상대 속도 이완을 지배하는 매개변수 $\beta$ 등) 를 보정함으로써 해결될 수 없음을 보여줍니다. 일반적인 혼합물에 대해 표준 수송 방정식과 Onsager 관계를 동시에 만족하는 $\beta$ 값은 존재하지 않습니다. 열확산율이 무시할 수 있는 유체의 특정 경우 ($B \approx 0$) 에조차, BGK 모델은 전체 확산 플럭스를 0 으로 강제하지 않는 한 올바른 열 플럭스 거동을 재현하지 못합니다.

의의 및 주장
본 논문은 이러한 불준수가 다종 BGK 모델의 실용적 적용에 심각한 문제를 제기한다고 주장합니다. 모델에서 유도된 수송 계수가 필요한 열역학적 상호역성을 따르지 않기 때문에, 모델을 "보정"하는 것, 즉 모델의 수송 특성이 특정 실제 유체 (예: 수증기와 공기) 의 것과 일치하도록 매개변수 값을 선택하는 것이 불가능합니다.

저자는 모델이 Onsager 관계를 정확히 만족하지는 않지만, 두 가지 특정 극한에서 점근적으로 만족한다고 지적합니다:

1. 질량비 $m_1/m_2 \to 0$일 때 (전자 질량이 이온에 비해 무시할 수 있는 이온화 플라즈마에 관련됨).
2. 질량비 $m_1/m_2 \to 1$일 때 (열 플럭스 표현식의 첫 번째 항이 소멸됨).

본 논문은 보존 법칙, H-정리, 무차별성 원리, 온도/농도의 양수성, 임의의 프란틀/슈미트 수를 표현할 수 있는 능력, 그리고 Onsager 관계 준수를 포함하는 "완벽한" 다종 운동론 모델의 요구 사항을 나열하며 결론을 맺습니다. 저자는 현재까지 모든 이러한 요구 사항을 동시에 만족하는 것으로 입증된 BGK 유형 모델은 존재하지 않는다고 명시합니다. 이러한 발견은 현재 BGK 모델이 유용한 근사치이지만, Onsager 관계를 만족하지 못한다는 사실이 물리적 충실도와 일반 유체 역학 응용 분야에서 정밀한 보정을 위한 유용성을 제한한다는 것을 시사합니다.