[논문 리뷰] A Linear, Decoupled and Energy stable scheme for smectic-A Liquid Crystal Flows

Xiaofeng Yang, A. V. Brylev|arXiv (Cornell University)|2017. 01. 25.

Fluid Dynamics and Turbulent Flows참고 문헌 5인용 수 3

한 줄 요약

이 논문은 비압축성 나비에-스토크스 방정식과 두 개의 2차 선형 타원형 방정식(배향 및 층 배열 순서 매개변수에 대한)을 결합한 비선형 시스템에 의해 지배되는 스메크틱-A 액정 유체의 흐름을 시뮬레이션하기 위한 선형적이고 분리된, 그리고 조건부 에너지 안정성의 시간 이산화 방법을 제시한다. 비선형 항에 대해 명시적-암시적 처리를 적용함으로써, 이 방법은 조건부 에너지 안정성을 확보하고 각 시간 단계에서 선형적이고 분리된 방정식을 풀어, 시간이 지남에 따라 에너지가 소산되는 것을 보장하는 효율적이고 견고한 장기 시뮬레이션을 가능하게 한다.

ABSTRACT

In this paper, we consider numerical approximations for the model of smectic-A liquid crystal flows. The model equation, that is derived from the variational approach of the de Gennes free energy, is a highly nonlinear system that couples the incompressible Navier-Stokes equations, and two nonlinear coupled second-order elliptic equations. Based on some subtle explicit--implicit treatments for nonlinear terms, we develop a unconditionally energy stable, linear and decoupled time marching numerical scheme. We also rigorously prove that the proposed scheme obeys the energy dissipation law at the discrete level. Various numerical simulations are presented to demonstrate the accuracy and the stability thereafter.

연구 동기 및 목표

• 스메크틱-A 액정 유체를 모델링하는 복잡하고 비선형적인 시스템을 위한 수치적으로 효율적이고 열역학적으로 일관된 방법을 개발하기 위해.
• 유속, 압력, 배향장, 층 함수 사이의 강한 결합을 대류 및 응력 항을 통해 다루기 위해.
• 시간 이산화가 조건부 안정성이며 이산화 수준에서 에너지 소산 법칙을 유지하도록 설계하기 위해.
• 각 시간 단계에서 비선형 시스템을 풀지 않고도 선형적이고 분리된 방정식을 얻기 위해.
• 이 방법의 정확도, 안정성 및 자석장과 비틀림 흐름 하에서의 물리적 현상(예: 케브론 패턴)을 잘 포착할 수 있는지 검증하기 위해.

제안 방법

• 비선형 항에 대해 명시적 처리, 선형 항에 대해 암시적 처리를 적용한 안정화된 시간 이산화 방법을 설계하였다.
• 이 방법은 각 시간 단계에서 유속, 압력, 배향장, 층 함수에 대해 별개의 선형 방정식으로 시스템을 분리한다.
• 비선형 대류 및 응력 항을 다루기 위해 미묘한 분할 전략을 적용하여 에너지 안정성을 확보하였다.
• 이산 에너지 법칙을 엄밀히 증명하여, 시간 간격의 크기와 관계없이 조건부 에너지 소산이 보장됨을 보였다.
• 적절한 경계 조건을 갖춘 유한 차분 또는 유한 요소 공간 이산화 방법을 사용하여 방법을 구현하였다.
• 수치 실험을 통해 시간에 대한 일계 정밀도와 에너지 안정성이 수렴성 및 장기 시뮬레이션을 통해 검증되었다.

실험 결과

연구 질문

• RQ1 강한 비선형성을 띤 스메크틱-A 액정 모델에 대해, 에너지 안정성을 유지하면서도 선형적이고 분리된 시간 이산화를 구성할 수 있는가?
• RQ2 유속, 압력, 배향장, 층 함수 사이의 강한 결합을 각 시간 단계에서 비선형 시스템을 풀지 않고 어떻게 다룰 수 있는가?
• RQ3 이산 설정에서 조건부 에너지 안정성을 보장하기 위해 비선형 항에 대해 어떤 명시적-암시적 처리가 필요한가?
• RQ4 이 방법은 자석장과 비틀림 흐름에 의해 유도되는 물리적 현상(예: 케브론 패턴)을 정확히 포착할 수 있는가?
• RQ5 장기 시뮬레이션에서 이 방법은 일계 시간 정밀도와 에너지 소산을 유지하는가?

주요 결과

• 제안된 방법은 조건부 에너지 안정성을 확보하여, 시간 간격의 크기와 관계없이 매 시간 단계에서 이산 에너지가 감소함을 보였다.
• 이 방법은 선형적이고 분리되어 있어, 각 시간 단계에서 유속, 압력, 배향장, 층 함수에 대해 별개의 선형 시스템을 풀어 계산 비용을 크게 감소시켰다.
• 수치적 수렴 테스트를 통해 모든 변수에 대해 시간에 대해 일계 정밀도를 확보하였으며, 오차 기울기가 이론적 기대와 일치함을 확인하였다.
• 헬프리치-후라우트 효과의 시뮬레이션은 자석장 하에서 케브론 패턴의 형성을 보여주었고, 에너지 소산 곡선은 안정성을 확인하였다.
• 비틀림 흐름을 포함함으로써 불규칙한 파동의 대칭성이 깨지며, 복잡한 유체역학적 효과를 포착할 수 있음을 보여주었다.
• 장기 시뮬레이션 동안 에너지 소산 곡선은 단조 감소하는 경향을 유지하여, 이 방법이 장기 동역학에 대해 견고함을 검증하였다.