[논문 리뷰] On the numerical solution of a variable-coefficient Burgers equation arising in granular segregation

Ivan C. Christov|arXiv (Cornell University)|2017. 06. 30.

Granular flow and fluidized beds참고 문헌 2인용 수 3

한 줄 요약

이 논문은 건조한 이성분 혼합물에서의 입자 분리 현상을 모델링하는 변수 계수를 가진 버거스 방정식을 해결하기 위해 두 번째 차수 정확도를 가지며 강하게 암시적인 크랭크-니콜슨 수치적 방법을 제시한다. 이 방법은 입자 확률의 비선형 경계 조건을 처리할 수 있으며, 정확한 kink 해와의 비교를 통해 두 번째 차수 수렴성과 양호한 보존 성질을 입증한다. 또한 시뮬레이션에서 유도된 운동 에너지 스트레스 프로파일을 바탕으로 일정계수 및 변수 계수 케이스에 모두 적용된다.

ABSTRACT

We study a variable-coefficient Burgers equation arising in the modelling of segregation of dry bidisperse granular mixtures. The equation is subject to nonlinear boundary conditions for the particle flux. We construct a strongly implicit Crank--Nicolson type of numerical scheme for the latter equation. The scheme is benchmarked against a standard exact solution of kink type, showing second-order of accuracy and good discrete conservation properties. Two segregation problems considered in the literature are then solved and discussed. The first is the case of a linear kinetic stress profile, which renders the governing equation of constant-coefficient type, while the second is the case of a variable kinetic stress profile based on an expression fit to particle dynamics simulation data.

연구 동기 및 목표

• 이성분 혼합물의 입자 분리 현상에서 발생하는 변수 계수 버거스 방정식을 위한 강력한 수치적 방법을 개발하는 것.
• 기본 방정식에서 입자 확률에 대한 비선형 경계 조건을 다루는 것.
• 수치적 해에서 두 번째 차수 정확도와 이산 보존 성질을 확보하는 것.
• 검증을 위해 정확한 kink 유형 해와의 비교를 통해 방법을 평가하는 것.
• 두 가지 물리적 사례에 적용: 입자 운동 시뮬레이션에서 유도된 선형 및 변수 운동 에너지 스트레스 프로파일 기반

제안 방법

• 변수 계수 버거스 방정식에 대해 강하게 암시적인 크랭크-니콜슨 유형의 유한차분 방법을 구성한다.
• 이 방법은 입자 확률에 대한 비선형 경계 조건을 처리하도록 설계된다.
• 정확한 kink 유형 해석적 해를 사용하여 방법을 검증하여 두 번째 차수 수렴성을 확인한다.
• 이산 보존 성질을 유지하여 물리적 일관성에 필수적인 특성을 확보한다.
• 두 가지 물리적 시나리오를 시뮬레이션한다: 선형 운동 에너지 스트레스 프로파일(일정계수 케이스)과 입자 운동 시뮬레이션에서 유도된 변수 프로파일.
• 수치 결과를 기존 문헌과 비교하여 정확도와 물리적 관련성을 평가한다.

실험 결과

연구 질문

• RQ1 입자 분리 현상에서 비선형 확률 경계 조건을 가진 변수 계수 버거스 방정식을 정확하게 해결하기 위해 어떤 수치적 방법을 설계할 수 있는가?
• RQ2 제안된 수치적 방법의 수렴 차수와 보존 행동은 어떠한가?
• RQ3 선형 운동 에너지 스트레스 프로파일을 가진 일정계수 케이스에 적용했을 때 모델의 성능은 어떠한가?
• RQ4 입자 운동 시뮬레이션에서 유도된 변수 운동 에너지 스트레스 프로파일을 적용했을 때 모델의 거동는 어떠한가?
• RQ5 입자 운동 시뮬레이션 데이터에 기반한 변수 계수 방정식을 수치적으로 해결함으로써 분리 역학에 대한 어떤 통찰을 얻을 수 있는가?

주요 결과

• 제안된 크랭크-니콜슨 방법은 정확한 kink 해와의 수렴 테스트를 통해 공간 및 시간에서 두 번째 차수 정확도를 확보한다.
• 이 방법은 양호한 이산 보존 성질을 보이며, 입자 확률 모델링에 필수적인 물리적 불변량을 유지한다.
• 선형 운동 에너지 스트레스 케이스에서는 기존의 일정계수 버거스 방정식의 행동과 일치하는 해를 보인다.
• 입자 운동 시뮬레이션 데이터에 기반한 변수 계수 케이스는 일정계수 모델이 포착하지 못하는 특별한 분리 역학을 드러낸다.
• 수치 결과는 기존 문헌과 일치하여 모델의 물리적 관련성과 방법의 신뢰성을 검증한다.