Kinetic derivation of thermal viscous models for nematic liquid crystal dynamics

기술 요약: 네마틱 액정 역학을 위한 운동론적 유도 기반 열적 점성 모델

문제 정의
최근의 연구는 운동론으로부터 유도된 압축성 비점성 변형의 레슬리-에릭센(Leslie–Ericksen) 방정식을 확립하였다. 그러나 해당 접근법은 이미 정렬된 상태를 중심으로 모멘트 폐쇄(moment closure)를 수행함으로써, 정렬이 역학으로부터 발현되도록 하는 대신 집단적 배향 행동을 사전에 부과하는 방식에 의존했다. 또한, 결과적으로 도출된 모델은 엄격하게 비점성 모델이었기에, 실제 액정 흐름에서 내재적인 소산, 열적 효과 및 스핀 확산 효과를 포착하지 못했다. 본 논문은 이러한 한계를 해결하기 위해, 정렬의 발현을 고려하고 열적 및 점성 효과를 포함하는 압축성 점성 변형의 레슬리-에릭센 방정식을 유도한다.

방법론
저자들은 정렬된 유체의 운동론으로부터 출발하여 거시적 열역학 이론을 전개한다. 방법론은 다음 단계들을 통해 진행된다:

1. 운동론적 프레임워크: 출발점은 Vlasov 포텐셜을 통해 평균장 정렬 상호작용을 포함하는 [7]의 정렬된 유체의 운동론이다. 이 시스템은 배향 완화가 병진 운동량 완화보다 빠르게 일어나는 시간 척도의 분리를 가정한다. 분포 함수 $f(x, p, \nu, \varsigma, t)$는 위치, 선운동량, 미세구조 구성(배향 $\nu$), 그리고 공액 운동량(conjugate momentum) $\varsigma$에 의존한다.
2. 충돌 연산자: 분석을 용이하게 하기 위해, 전체 볼츠만 충돌 연산자를 단일 완화 BGK(Bhatnagar–Gross–Krook) 연산자로 대체한다. 이 연산자는 전체 운동론의 보존 법칙과 엔트로피 생성 구조를 보존한다.
3. 평형 법칙: 저자들은 운동량의 모멘트를 취함으로써 질량, 선운동량, 각운동량, 에너지 및 엔트로피에 대한 거시적 평형 법칙을 유도한다. 핵심적인 특징은 각운동량 평형에서 정렬을 유도하는 평균장 토크 항을 식별하는 것이다.
4. 열역학적 구조: 헬름홀츠 자유 에너지를 열역학적 포텐셜로 식별한다. 저자들은 평형 분포(Maxwellian)가 이 자유 에너지를 최소화함을 증명하고, 이와 관련된 깁스 분포를 배향 프로파일에 대해 유도한다. 특정 상호작용 포텐셜(Kuramoto 및 Maier–Saupe 유형) 하에서, 시스템은 온도가 낮아짐에 따라 등방성에서 네마틱 상태로 상전이를 일으키며, 정렬된 평형 상태에서 Vlasov 토크가 소멸함을 입증한다.
5. 채프먼-엔스코그(Chapman–Enskog) 전개: 국소 맥스웰 평형 주변에서 1차 채프먼-엔스코그 전개를 수행한다. 표준적인 유도와 달리, 전개에서의 시간 미분은 표준 압축성 오일러 방정식이 아닌, 제0차 평형 법칙(Ericksen 응력 및 커플 응력을 포함)을 사용하여 대체된다. 이는 배향 탄성력이 1차 보정항으로 올바르게 전파되도록 보장한다.
6. 구성 폐쇄: 저자들은 Rajagopal과 Srinivasa의 제약 최대화 원리를 채택한다. 운동론적 전개로부터 유도된 엔트로피 생성율과 헬름홀츠 자유 에너지로부터 추론된 열역학적 구조를 결합한다. 평형 법칙에 대해 엔트로피 생성율을 최대화함으로써, 코시 응력(Cauchy stress), 커플 응력(couple-stress) 및 열 유속(heat flux)에 대한 구성 관계식을 유도한다.

주요 기여 및 결과

• 발현적 정렬(Emergent Alignment): 정렬을 가정한 기존의 비점성 모델들과 달리, 본 유도는 Vlasov 토크를 통해 등방성-네마틱 전이가 자연스럽게 발현되도록 한다. 분석 결과, 정렬된 깁스 평형에서 평균장 토크가 소멸함을 확인하였으며, 이는 거시적 평형 법칙과 일관성을 유지한다.
• 점성 및 열적 효과: 저자들은 압축성, 점성, 열전도 및 스핀 확산 특성을 가진 네마틱 유체를 위한 완전한 구성 방정식 세트를 유도한다. 결과적인 응력 텐서는 다음을 포함한다:
  • 뉴턴식 전단 점성.
  • 병진-회전 에너지 불균형에서 기인하는 등방성 부피 점성.
  • 에릭센 텐서 및 속도 구배와 결합된 비등방성 점성 응력.
  • 스핀 확산 효과(Eringen 유형의 마이크로폴라 항) 및 디렉터 구배와의 결합을 포함하는 커플 응력 텐서.
• 수송 계수: 전단 점성, 열전도도, 부피 점성, 비등방성 점성 및 스핀 점성에 대한 수송 계수의 명시적 표현이 BGK 완화 시간 $\tau$와 프랑크 계수(Frank modulus) $K$에 대해 유도된다. 열전도도는 이원자 가스에 대한 Eucken 관계를 만족한다.
• 비압축성 극한: 이 프레임워크는 최대화 과정 중에 라그랑주 승서(Lagrange multiplier)를 통해 발산 없는 제약 조건을 강제함으로써, 별도의 유도 없이도 자연스럽게 비압축성 대응물을 산출한다.
• 구조적 참신함: 유도된 점성 응력 텐서는 단일 점성 매개변수(등방성 BGK 근사로 인해)에 의해 지배되며, 다섯 개의 계수를 가진 완전한 레슬리 점성 텐서를 회복하지 못한다. 저자들은 완전한 레슬리 텐서를 회복하려면 여러 완화 시간 척수를 가진 비등방성 BGK 연산자가 필요하며, 이는 향후 연구 과제임을 언급한다.

의의 및 주장
본 논문은 최근 Farrell, Russo, Zerbinati [1]의 비점성 운동론을 일반화하여 점성, 열적 및 스핀 확산 효과를 설명한다고 주장한다. 주요 의의는 이러한 소산 효과를 운동론으로부터 직접적으로 유도함으로써, 미시적 입자 상호작용과 거시적 연속체 역학 사이의 간극을 메우고 열역학적으로 일관된 유도를 제공한다는 점에 있다.

저자들은 정렬을 사전에 부과하지 않고, 대신 자유 에너지 최소화와 평형에서의 Vlasov 토크 소멸을 통해 정렬이 발현되도록 함으로써 기존 모델의 한계를 극복했음을 강조한다. 결과적으로 도출된 시스템(식 8.21)은 레슬리-에릭센 방정식의 압축성, 점성, 열전도 및 스핀 확산 변형을 나타낸다. 비점성 극한($\tau \to 0$)에서 이 모델은 [1]에서 유도된 비점성 시스템으로 환원된다.

이 연구는 액정의 더 완전한 운동론을 향한 단계로 제시되며, 저자들은 다섯 개의 독립적인 점성을 가진 완전한 레슬리 점성 텐서의 회복은 더 복잡한 비등방성 충돌 연산자를 요구하는 미해결 과제임을 겸허히 밝히고 있다.