[논문 리뷰] Finite element formulation of general boundary conditions for incompressible flows

Roland Becker, Daniela Capatina외 2명|arXiv (Cornell University)|2015. 05. 06.

Advanced Numerical Methods in Computational Mathematics참고 문헌 27인용 수 7

한 줄 요약

이 논문은 유체의 점성(나비에-스토퍼스 방정식)과 비점성(유레어 방정식) 영역을 하나의 약한 수식화로 통합하여 비압축성 유동에서 일반적인 경계 조건을 다루기 위한 유한요소 수식을 제시한다. 유속의 분포를 균형 잡힌 스펙트럼 분해를 통해 처리하고, 일관된 운동에너지 안정화 항을 추가함으로써, 점성 수준이 모두 포함된 경우에도 안정성과 일관성을 확보한다. 특히 쌍곡선 극한에서도 성능이 유지된다. 수치 실험을 통해 이론적 접근은 복잡한 유동 특징을 진동 없이 잘 포착하며, 강한 강제 조건보다도 더 뛰어난 안정성과 스케일 불변성을 확보한다.

ABSTRACT

We study the finite element formulation of general boundary conditions for incompressible flow problems. Distinguishing between the contributions from the inviscid and viscid parts of the equations, we use Nitsche's method to develop a discrete weighted weak formulation valid for all values of the viscosity parameter, including the limit case of the Euler equations. In order to control the discrete kinetic energy, additional consistent terms are introduced. We treat the limit case as a (degenerate) system of hyperbolic equations, using a balanced spectral decomposition of the flux Jacobian matrix, in analogy with compressible flows. Then, following the theory of Friedrich's systems, the natural characteristic boundary condition is generalized to the considered physical boundary conditions. Several numerical experiments, including standard benchmarks for viscous flows as well as inviscid flows are presented.

연구 동기 및 목표

• 모든 점성 수준, 특히 비점성인 유레어 극한까지도 안정성과 일관성을 유지하는 비압축성 유동의 일반적인 경계 조건을 위한 통합된 유한요소 수식을 개발하는 것.
• 표준 방법이 실패하는 고 페클레 수 또는 쌍곡선 영역에서 이산 운동에너지 제어 문제를 해결하는 것.
• 유동의 특성에 따라 경계 조건을 일반화하기 위해, 압축성 유동 이론과 유사하게, 유속의 분포 행렬의 균형 잡힌 스펙트럼 분해를 통해 비압축성 유동에 대한 특성 경계 조건을 일반화하는 것.
• 특히 제트 충격 및 후방 경사면 문제와 같은 복잡한 흐름에 대해 수치 시뮬레이션의 안정성과 스케일 불변성을 확보하는 것.
• 니츠에의 방법을 통한 약한 강제 조건과 강한 강제 조건을 비교하여, 경계층과 내측 모서리에서 더 뛰어난 안정성과 더 적은 진동을 보이는 것을 입증하는 것.

제안 방법

• 유동의 비점성(에일러)과 점성(스토크스) 기여를 구분하여, 점성 계수 µ = 0인 경우를 포함해 모든 점성 수준에서 니츠에의 방법을 사용해 경계 조건을 약하게 강제한다.
• 경계 행렬 B를 정의하기 위해, 유속의 분포 행렬의 균형 잡힌 스펙트럼 분해를 사용하여, 열린 영역 방정식과 차원 일치성을 확보한다. 이는 특이한 쌍곡선 극한에서도 유지된다.
• 이산 운동에너지 제어를 위해 추가적인 일관된 안정화 항을 도입하여, 수치 해에서 비물리적 에너지 증가를 방지한다.
• 나비에-스토크스 방정식에 적용하기 위해, 스토크스 흐름의 약한 수식과 유속의 분포 행렬 고유값의 절댓값에 기반한 일관된 안정화 항을 조합한다.
• 연속 유한요소를 사용해 수식을 구현하며, 공간 이산화에 집중한다. 유동성 강제 조건 대비 SUPG 안정화를 사용한 강한 강제 조건 버전과 비교한다.
• 스케일 불변성을 검증하기 위해, 경계 데이터를 적절히 변환한 스케일된 영역에서 제트 충격 문제를 해결한다. 이로써 기하학적 및 물리적 스케일링 하에서 해의 불변성을 검증한다.

실험 결과

연구 질문

• RQ1 어떻게 하나의 유한요소 수식이 점성 및 비점성 비압축성 유동, 특히 에일러 극한까지도 안정성과 일관성을 유지할 수 있는가?
• RQ2 고 페클레 수 또는 쌍곡선 특성을 띠는 비압축성 유동에서 이산 에너지 경계를 유지하는 데에 운동에너지 안정화의 역할은 무엇인가?
• RQ3 유속의 분포 행렬의 균형 잡힌 스펙트럼 분해를 사용하면 비압축성 유동에 대해 특성 경계 조건을 어떻게 자연스럽게 일반화할 수 있는가?
• RQ4 복잡한 흐름 특징(예: 소용돌이, 경계층 불안정성)을 포착할 때, 니츠에의 방법을 통한 약한 강제 조건이 강한 강제 조건보다 어떤 수치적 이점이 있는가?
• RQ5 제안된 방법이 제트 충격과 같은 벤치마크 문제에서 이산 해에 대해 어느 정도의 스케일 불변성을 유지하는가?

주요 결과

• 제안된 니츠에 기반의 약한 수식화는 점성 수준 전반, 특히 비점성인 에일러 극한까지도 일관성과 안정성 손실 없이 잘 처리한다.
• 일관된 운동에너지 안정화 항을 추가함으로써, 고 페클레 수 및 쌍곡선 영역에서 비물리적 에너지 증가를 효과적으로 억제한다.
• 후방 경사면 및 제트 충격 문제에 대한 수치 실험 결과, 약한 강제 조건 수식화는 강한 강제 조건보다 더 적은 비정상적인 진동을 보이며, 특히 내측 모서리와 벽면 근처에서 유리하다.
• 이 방법은 스케일 불변성을 확보한다: 영역과 경계 데이터를 적절히 스케일링할 경우, 이산 해는 기계 정밀도 수준까지 불변성을 유지한다. 반면 강한 강제 조건 방법은 이를 충족하지 못한다.
• 제트 충격 문제에서 압력 등압선과 속도 크기 프로파일을 분석한 결과, 약한 수식화는 진동을 경계 근처에 국한시키는 반면, 강한 수식화는 도메인 내부로 비정상적인 진동이 퍼지는 경향을 보였다.
• 유속의 분포 행렬의 스펙트럼 분해를 통해 비압축성 유동에 대해 특성 경계 조건을 물리적으로 의미 있는 방식으로 일반화할 수 있었으며, 이는 제1차 미분방정식계에 대한 프리드리히스 이론과 일치한다.