[논문 리뷰] A coupled boundary element / finite element method for the convected Helmholtz equation with non-uniform flow in a bounded domain

Fabien Casenave, Alexandre Ern외 1명|arXiv (Cornell University)|2013. 03. 27.

Electromagnetic Simulation and Numerical Methods인용 수 1

한 줄 요약

이 논문은 비균일한 속도 프로파일을 가진 초음속 유동에서의 운반된 헬름홀츠 방정식을 위한 경계요소법과 유한요소법의 새로운 결합(BEM-FEM)을 제안한다. 전체 유동 영역에 프란틀-플로이어트 변환을 적용함으로써 문제를 외부 영역에서는 고전적 헬름홀츠 방정식, 내부 영역에서는 비대칭 항을 가진 이방성 편미분방정식으로 재구성하여 자연스러운 전달 조건을 확보하고 표준 해석기로 효율적으로 해결할 수 있도록 한다. 두 가지 제안된 형식 중 하나는 공명에 취약하고, 다른 하나는 이를 피한다.

ABSTRACT

We consider the convected Helmholtz equation modeling linear acoustic propagation at a fixed frequency in a subsonic flow around a scattering object. The flow is supposed to be uniform in the exterior domain far from the object, and potential in the interior domain close to the object. Our key idea is the reformulation of the original problem using the Prandtl--Glauert transformation on the whole flow domain, yielding (i) the classical Helmholtz equation in the exterior domain and (ii) an anisotropic diffusive PDE with skew-symmetric first-order perturbation in the interior domain such that its transmission condition at the coupling boundary naturally fits the Neumann condition from the classical Helmholtz equation. Then, efficient off-the-shelf tools can be used to perform the BEM-FEM coupling, leading to two novel variational formulations for the convected Helmholtz equation. The first formulation involves one surface unknown and can be affected by resonant frequencies, while the second formulation avoids resonant frequencies and involves two surface unknowns. Numerical simulations are presented to compare the two formulations.

연구 동기 및 목표

• 비균일한 초음속 유동에서 산란체 주위의 유한 영역에서 운반된 헬름홀츠 방정식을 해결하는 데 도전하는 것.
• 특히 산란체 근처에서 흐름이 비균일할 경우 경계요소법과 유한요소법을 결합하는 데 어려움을 해결하는 것.
• 유동 영역 간의 경계에서 자연스러운 전달 조건을 보장하는 변분형식을 개발하는 것.
• 문제를 더 다룰 수 있는 형태로 재구성함으로써 표준 오프더쇼프 BEM 및 FEM 해석기를 효과적으로 사용할 수 있도록 하는 것.

제안 방법

• 전체 흐름 영역에 프란틀-플로이어트 변환을 적용하여 운반된 헬름홀츠 방정식을 변환된 시스템으로 매핑하는 것.
• 외부 영역에서는 고전적 헬름홀츠 방정식을 만족하고, 내부 영역에서는 비대칭 항을 가진 이방성이고 확산성 편미분방정식을 만족하는 변환된 시스템을 유도하는 것.
• 내부-외부 경계에서의 전달 조건이 고전적 헬름홀츠 방정식의 노이만 조건과 자연스럽게 일치하도록 보장하는 것.
• 단일 표면 미지수(공명에 취약한)와 이중 표면 미지수(공명이 없는)를 사용하는 두 가지 변분 접근법을 제안하는 것.
• 외부 영역에는 표준 BEM, 내부 영역에는 표준 FEM을 적용하여 재구성된 문제의 구조 덕분에 기존의 해석기를 활용하는 것.
• 수치 시뮬레이션을 통해 두 형식의 성능과 안정성 비교를 수행하는 것.

실험 결과

연구 질문

• RQ1 프란틀-플로이어트 변환을 전역적으로 적용하여 운반된 헬름홀츠 방정식을 고전적 헬름홀츠 방정식과 이방성 편미분방정식으로 분리할 수 있는가?
• RQ2 유도된 결합 형식이 내부 및 외부 영역 간에 자연스러운 전달 조건을 허용하는가?
• RQ3 이러한 변환된 프레임워크에서 표준 BEM 및 FEM 해석기를 효과적으로 사용할 수 있는가?
• RQ4 단일 표면 미지수와 이중 표면 미지수를 사용하는 두 가지 제안된 형식은 안정성 및 공명 행동 측면에서 어떻게 비교되는가?
• RQ5 비균일한 초음속 유동에서 음향 산란을 모델링할 때 이 형식들의 수치 성능는 어떠한가?

주요 결과

• 프란틀-플로이어트 변환은 외부 영역에서는 고전적 헬름홀츠 방정식, 내부 영역에서는 비대칭 항을 가진 이방성 편미분방정식으로 운반된 헬름홀츠 방정식을 성공적으로 분리한다.
• 내부-외부 경계에서의 전달 조건이 고전적 헬름홀츠 방정식의 노이만 조건과 자연스럽게 일치하여 원활한 결합이 가능하다.
• 단일 표면 미지수를 사용하는 첫 번째 형식은 적분방정식 형식의 특성상 공명 주파수에 취약하다.
• 이중 표면 미지수를 포함하는 두 번째 형식은 이중 트레이스 구조를 도입함으로써 시스템을 안정화시켜 공명 문제를 피한다.
• 수치 시뮬레이션은 두 형식 모두 정확한 해를 도출함을 확인했지만, 공명이 없는 형식이 주파수 스윕에서 더 뛰어난 강건성을 보였다.
• 재구성된 문제의 구조 덕분에 표준 BEM 및 FEM 해석기의 사용이 가능하고 효과적이며, 효율적인 구현이 가능하다.