[논문 리뷰] Simulation of strong nonlinear waves with vectorial lattice Boltzmann schemes

François Dubois|arXiv (Cornell University)|2014. 01. 02.

Lattice Boltzmann Simulation StudiesEngineering참고 문헌 20인용 수 22

한 줄 요약

이 논문은 비선형 수류 방정식의 초순수 보존법칙에서 강한 비선형 파동을 시뮬레이션하기 위해 이중 엔트로피의 운동학적 표현을 기반으로 한 벡터형 격자 볼츠만 방법을 제안한다. 특히 1차원 및 2차원에서 얕은 수면 방정식에 적용한다. 이 방법은 충격파, 반사 및 비정상 흐름을 성공적으로 포착하며, 고두노프 유한부피 방법과 유사한 결과를 보이지만, 더 높은 수치 점성도를 갖는다.

ABSTRACT

We show that an hyperbolic system with a mathematical entropy can be discretized with vectorial lattice Boltzmann schemes with the methodology of kinetic representation of the dual entropy. We test this approach for the shallow water equations in one and two space dimensions. We obtain interesting results for a shock tube, reflection of a shock wave and unstationary two-dimensional propagation. This contribution shows the ability of vectorial lattice Boltzmann schemes to simulate strong nonlinear waves in unstationary situations.

연구 동기 및 목표

• 1D에서 다중차원 초순수 보존법칙 시스템으로 이중 엔트로피의 운동학적 표현을 확장하기 위해.
• 양수 방정식과 같은 시스템에서 충격 형성 및 상호작용을 포함한 강한 비선형 파동을 시뮬레이션하는 데 도전하는 데에.
• 엔트로피 구조를 유지하고 불연속성을 다룰 수 있는 안정적이고 일致한 격자 볼츠만 프레임워크를 개발하기 위해.
• 정적 및 비정상 충격 역학을 포함한 벤치마크 시험 케이스에서의 성능 평가를 위해.

제안 방법

• 이 방법은 N개의 보존 변수를 위해 N개의 입자 분포 f_j^k와 이산 속도 v_j를 사용하여 모멘트 기반 방정식을 통해 시스템을 표현한다.
• 퍼스테임-부샤 라이프를 적용하여 이중 엔트로피를 이산 속도 집합과 관련된 볼록 함수로 분해한다.
• 평형 분포는 이중 엔트로피와 엔트로피 변수를 사용하여 구성되며, 이는 기초가 되는 초순수 시스템과의 일致성을 보장한다.
• 안정성을 확보하기 위해 동일한 회복 파rameter s_j = 1.8을 사용하는 BGK 유형의 충돌 모델을 적용하고, 비평형 모멘트에 대해 다중 회복 시간을 구현한다.
• 벽, 유입 및 유출 조건에 대해 반사 및 반대 반사 방법을 사용하여 경계 조건을 적용한다.
• 구조적 격자에서 정확한 해와 고두노프 유한부피 결과와의 비교를 통해 방법을 검증한다.

실험 결과

연구 질문

• RQ1 이중 엔트로피 분해를 기반으로 한 벡터형 격자 볼츠만 방법은 초순수 시스템에서 강한 비선형 파동을 정확하게 시뮬레이션할 수 있는가?
• RQ2 이 방법은 얕은 수면 유동에서 정적 충격 반사 및 비정상 충격 상호작용을 어떻게 포착하는가?
• RQ3 고해상도 유한부피 방법과 비교할 때 수치 점성의 영향이 이 방법의 정확도에 미치는 영향은 무엇인가?
• RQ4 다중차원 시스템과 다수의 보존 변수를 갖는 시스템으로 이중 엔트로피의 운동학적 표현을 일반화할 수 있는가?
• RQ5 메쉬 정밀도 향상과 파rameter 선택은 비정상 시뮬레이션에서 안정성과 수렴성에 어떤 영향을 미치는가?

주요 결과

• 벡터형 격자 볼츠만 방법은 2차원 얕은 수면 흐름에서 정적 충격 반사를 성공적으로 포착하며, 140×80 격자에서 정확한 해와 고두노프 유한부피 방법 결과와 일치한다.
• 프루드 수 3인 비정상적인 에머리 시험 케이스에서, 480×120 격자에서 t=1/2 및 t=4 시점에서 복잡한 충격파 상호작용과 비선형 파동 전파를 재현한다.
• 이 방법은 높은 수치 점성도를 보이며, 이는 고두노프 방법보다 비정상 상태에서 더 낮은 정확도를 초래한다. 고두노프 방법은 1차 정확도를 갖지만 일시적 영역에서 더 정밀하다.
• 이 방법은 기초가 되는 엔트로피 구조와 일치하여 시스템의 열역학적으로 일致한 진화를 보장한다.
• 메쉬 정밀도 향상은 수렴성을 보이며, 가장 정밀한 격자에서 격자 볼츠만 방법과 고두노프 방법의 밀도 프로파일이 매우 유사하게 일치한다.
• λ=80 및 s_j=1.8을 사용하는 D2Q5Q4Q4 방법 구성은 수치적 안정성을 보장하지만, 높은 λ 값으로 인해 작은 시간 간격과 느린 계산이 발생한다.