[논문 리뷰] A hybrid Eulerian-Lagrangian flow solver

Artur Palha, Lento Manickathan외 2명|arXiv (Cornell University)|2015. 05. 13.

Fluid Dynamics Simulations and Interactions참고 문헌 50인용 수 12

한 줄 요약

이 논문은 고체 경계 근처에서 유한요소법을 사용하는 오일러형 유한요소 방법과 난류 영역에서 라그랑주형 비선형 입자 방법을 조합한 하이브리드 오일러-라그랑주형 유동 해석기법을 제시한다. 이는 점성 경계층의 정확한 해석과 낮은 확산을 갖는 난류 영역의 역학을 가능하게 하며, 기준 솔루션과 양호한 일치를 보이는 벤치마크 유동(다이폴, 실린더, 정지한 공기역학체)에서 높은 정확도를 달성한다. 이는 풍력단지와 같은 복잡한 유동의 대규모이고 효율적인 시뮬레이션 잠재력을 보여준다.

ABSTRACT

Currently, Eulerian flow solvers are very efficient in accurately resolving flow structures near solid boundaries. On the other hand, they tend to be diffusive and to dampen high-intensity vortical structures after a short distance away from solid boundaries. The use of high order methods and fine grids, although alleviating this problem, gives rise to large systems of equations that are expensive to solve. Lagrangian solvers, as the regularized vortex particle method, have shown to eliminate (in practice) the diffusion in the wake. As a drawback, the modelling of solid boundaries is less accurate, more complex and costly than with Eulerian solvers (due to the isotropy of its computational elements). Given the drawbacks and advantages of both Eulerian and Lagrangian solvers the combination of both methods, giving rise to a hybrid solver, is advantageous. The main idea behind the hybrid solver presented is the following. In a region close to solid boundaries the flow is solved with an Eulerian solver, where the full Navier-Stokes equations are solved (possibly with an arbitrary turbulence model or DNS, the limitations being the computational power and the physical properties of the flow), outside of that region the flow is solved with a vortex particle method. In this work we present this hybrid scheme and verify it numerically on known 2D benchmark cases: dipole flow, flow around a cylinder and flow around a stalled airfoil. The success in modelling these flow conditions presents this hybrid approach as a promising alternative, bridging the gap between highly resolved and computationally intensive Eulerian CFD simulations and fast but less resolved Lagrangian simulations.

연구 동기 및 목표

• 순수 오일러형 해석기법은 난류에서 수치적 확산 문제를 야기하고, 순수 라그랑주형 해석기법은 복잡한 기하구조에서 경계 조건을 정확히 구현하는 데 어려움을 겪는다 능력의 한계를 해결하기 위해.
• 고체 경계 근처에서는 오일러형 방법의 정확성을, 난류 영역에서는 라그랑주형 비선형 입자 방법의 낮은 확산 특성을 활용하는 하이브리드 해석기법을 개발하기 위해.
• 각 기구체에 대해 오일러형 하위영역을 분리하고 난류 영역에는 메쉬 기반의 라그랑주형 해석기법을 사용함으로써, 복수의 이동하는 물체(예: 풍력단지)를 포함한 복잡한 유동의 효율적이고 확장 가능한 시뮬레이션을 가능하게 하기 위해.
• 기본적인 2차원 시험 케이스에서 하이브리드 접근법을 검증하고, 높은 레이놀즈 수의 유동에서의 능력을 입증하기 위해.
• 향후 3차원 확장, 이동 또는 변형 가능한 물체, 상하위영역 간 난류 결합을 고려한 난류 흐름으로의 확장 기반을 마련하기 위해.

제안 방법

• 계산 영역은 두 개의 하위영역으로 분할된다: 고체 경계 근처에서 나비에-스토크스 방정식을 유한요소법으로 해석하는 오일러형 영역과, 정규화된 비선형 입자 방법을 사용해 소용돌이 강도를 입자로 추적하는 라그랑주형 영역.
• 오일러형 및 라그랑주형 하위영역 간의 인터페이스 결합은 상호 보간을 통해 이루어지며, 보존성과 일致성을 확보하기 위해 정의된 겹침 영역이 포함된다.
• 전체 순환량을 정확히 유지하여, 반복적 샤우츠 기반 결합 기법에서 흔히 발생하는 순환량 오류를 방지한다.
• 라그랑주형 하위영역 내에서 속도장을 계산하기 위해 빠른 다중체 방법(FMM)이 사용되어 대규모 입자 수에 대한 확장성을 보장한다.
• 시간 적분은 각 하위영역에서 독립적으로 수행되며, 안정성과 정확도를 확보하기 위해 오일러형 시간 간격은 라그랑주형 시간 간격보다 작다.
• 벽면 경계 조건은 표준 기법(예: 벽 함수 포함)을 사용해 오일러형 해석기법에서 구현되며, 라그랑주형 해석기법은 자연스럽게 외부 유동 및 원거리 경계 조건을 만족한다.

실험 결과

연구 질문

• RQ1 순수 오일러형 해석기법에 비해 수치적 확산이 감소한 채 고체 경계 근처 및 난류 영역의 유동 구조를 정확하게 시뮬레이션할 수 있는가?
• RQ2 다이폴, 실린더, 정지한 공기역학체 등의 벤치마크 유동 솔루션에서 비선형 입자 방법의 정확도가 소용돌이 강도 분포 및 기계적 힘(양력, 저항력) 측면에서 얼마나 잘 재현되는가?
• RQ3 특히 소용돌이 분리 및 오일러형 영역으로의 재진입 상황에서, 고레이놀즈 수의 흐름에서 정확도와 안정성이 어느 정도 유지되는가?
• RQ4 각 기구체의 오일러형 하위영역를 별도로 처리함으로써, 풍력단지와 같은 복잡한 다체 유동에 대해 이 하이브리드 기법이 효과적으로 적용될 수 있는가?
• RQ5 단일 오일러형 해석기법에 비해 이 하이브리드 기법이 확장성과 병렬 처리 측면에서 어떤 계산적 이점이 있는가?

주요 결과

• 다이폴 흐름, 실린더 주위의 흐름, 정지한 공기역학체 주위의 흐름에 대해 하이브리드 해석기법은 소용돌이장과 기계적 힘(양력, 저항력)을 기준 유한요소 솔루션과 강한 일치를 보이며, t = 2s까지 정확하게 재현한다.
• 레이놀즈 수 Re = 5000인 정지한 타원체에 대해, 하이브리드 해석기법은 t = 2s까지 기준 솔루션과 밀도 높은 일치를 보이며, 소용돌이 구조와 힘 계수 모두 동일한 범위 내에서 유지된다. 소용돌이 분리 이후 점차 증가하는 편차가 발생하지만, 이는 예상 가능한 범위 내이다.
• 이 방법은 소용돌이 분리 및 오일러형 영역으로의 재진입과 같은 복잡한 유동 역학을 성공적으로 포착하여, 도전적인 흐름 영역에서의 강건성을 입증한다.
• 반복적이지 않고 순환량을 유지하는 결합 기법을 사용함으로써, 이전의 하이브리드 접근법에 비해 순환량 오류를 제거하고 안정성을 향상시켰다.
• 각 기구체에 대해 오일러형 하위영역를 분리함으로써 효율적이고 확장 가능한 시뮬레이션을 가능하게 하여, 이동하거나 변형 가능한 물체의 처리를 크게 단순화시켰다.
• 결과는 이 하이브리드 접근법이 풍력단지 및 기타 복잡한 외부 유동에서 난류 역학과 경계층 해석이 핵심적인 대규모 시뮬레이션에 유망한 대안이 될 수 있음을 시사한다.