[논문 리뷰] On the homogenization of the Helmholtz problem with thin perforated walls of finite length

Adrien Semin, Bérangère Delourme외 1명|arXiv (Cornell University)|2016. 11. 18.

Advanced Mathematical Modeling in Engineering참고 문헌 34인용 수 14

한 줄 요약

이 논문은 두께가 얇고 주기적이며 유한 길이의 구멍이 뚫린 벽을 가진 2차원 영역에서 헬름홀츠 방정식에 대한 두 번째 차수 점근 전개를 제시한다. 주기적 표면 균질화와 일치하는 점근 전개를 결합하여 경계층과 모서리 특이성을 포착한다. 이 방법은 효과적인 인터페이스에서의 전송 조건과 발산하는 특이성을 해결하는 모서리 조건을 도출하며, 작은 주기 δ에 독립적인 계산 비용으로 고정밀도 해를 가능하게 한다.

ABSTRACT

In this work, we present a new solution representation for the Helmholtz transmission problem in a bounded domain in $\\mathbb{R}^2$ with a thin and periodic layer of finite length. The layer may consists of a periodic pertubation of the material coefficients or it is a wall modelled by boundary conditions with an periodic array of small perforations. We consider the periodicity in the layer as the small variable $\\delta$ and the thickness of the layer to be at the same order. Moreover we assume the thin layer to terminate at re-entrant corners leading to a singular behaviour in the asymptotic expansion of the solution representation. This singular behaviour becomes visible in the asymptotic expansion in powers of $\\delta$ where the powers depend on the opening angle. We construct the asymptotic expansion order by order. It consists of a macroscopic representation away from the layer, a boundary layer corrector in the vicinity of the layer, and a near field corrector in the vicinity of the end-points. The boundary layer correctors and the near field correctors are obtained by the solution of canonical problems based, respectively, on the method of periodic surface homogenization and on the method of matched asymptotic expansions. This will lead to transmission conditions for the macroscopic part of the solution on an infinitely thin interface and corner conditions to fix the unbounded singular behaviour at its end-points. Finally, theoretical justifications of the second order expansion are given and illustrated by numerical experiments. The solution representation introduced in this article can be used to compute a highly accurate approximation of the solution with a computational effort independent of the small periodicity $\\delta$.

연구 동기 및 목표

• 박막하고 주기적이며 유한 길이의 미세한 구멍 또는 물질의 비균질성으로 구성된 양면 2차원 영역에서 헬름홀츠 전파 문제를 모델링한다.
• 두 가지 별개의 특이 행동을 고려한다: 주기적 인터페이스 근처의 경계층과, 내측 모서리로 인한 끝점에서 발생하는 모서리 특이성.
• macroscopic 행동과 국소 보정을 동시에 포착하는 두 번째 차수 점근 전개를 구성한다. 이는 일치하는 점근 전개와 주기적 균질화를 활용한다.
• 무한히 얇은 인터페이스에서 효과적인 전송 조건과 끝점에서 발생하는 발산 특이성을 정규화하는 모서리 조건을 도출한다.
• 이론적으로 타당하고 계산 비용이 작은 매개변수 δ에 독립적인 해의 근사해를 제공한다.

제안 방법

• 해는 매크로스코픽 부분, 경계층 보정항(접선 방향으로 주기적이고 법선 방향으로 지수적으로 감쇠하는), 그리고 모서리 특이성을 위한 근접장 보정항으로 분해된다.
• 경계층 보정항은 주기적 표면 균질화를 통해 유도되며, 주기성 δ와 인터페이스에서의 거리 d/δ에 따라 감쇠하는 정준 문제를 푸는 방식이다.
• 근접장 보정항은 레이어의 끝점 부근에서 일치하는 점근 전개를 통해 구성되며, 진동하는 경계를 가진 무한한 원뿔 영역에서의 해를 기반으로 한다.
• 이 방법은 (R±)λ(R±)k ln(R±)m 형태의 이중 척도 점근 전개를 사용한다. 여기서 R±는 모서리에서의 극좌표 거리이며, λ는 꼭짓점의 괴이각 Θ에 의해 결정된다.
• 모서리 특이성은 λn−q ∈ π/Θ ℤ* 또는 λn−q = 0일 때 프로파일에 수직 조건을 도입하여 유일성을 확보함으로써 해결된다.
• 전체 해는 매크로스코픽 항, 경계층 보정항, 근접장 보정항의 초월합으로 표현되며, 계수는 원뿔 영역에서의 정준 문제를 풀어 결정된다.

실험 결과

연구 질문

• RQ1 두께가 얇고 주기적이며 유한 길이의 구멍이 뚫린 벽을 가진 영역에서 헬름홀츠 방정식에 대한 두 번째 차수 점근 전개는 어떻게 구성할 수 있는가?
• RQ2 해의 매크로스코픽 행동을 포착하는 효과적인 인터페이스에서의 전송 조건은 무엇인가?
• RQ3 박막하고 주기적인 레이어의 끝점에서 발생하는 모서리 특이성은 수학적으로 어떻게 특성화하고 해결할 수 있는가?
• RQ4 괴이각 Θ은 특이 행동과 점근 전개에서 δ의 거듭제곱에 어떤 영향을 미치는가?
• RQ5 점근 해는 이론적으로나 수치적으로 타당한가? 계산 비용은 작은 매개변수 δ에 독립적인가?

주요 결과

• 점근 전개는 두 번째 차수까지 구성되었으며, 매크로스코픽 항, 경계층 보정항, 근접장 보정항을 포함하고 있으며, 엄밀한 이론적 근거가 제시되었다.
• 경계층 보정항은 접선 방향으로 주기적이며, 법선 방향으로 e−d/δ 비례로 감쇠하는 지수 감쇠 함수이다. 여기서 d는 인터페이스로부터의 거리이다.
• 근접장 보정항은 진동하는 경계를 가진 무한한 원뿔 영역에서의 해로부터 유도되며, 내측 모서리에서의 특이 행동을 δ의 거듭제곱에 따라 포착한다. 이 거듭제곱은 괴이각 Θ에 따라 달라진다.
• 이 방법은 주기적 표면 균질화로부터 도출된 무한히 얇은 인터페이스에서의 효과적 전송 조건과 발산 특이성을 정규화하는 모서리 조건을 도출한다.
• 해의 표현은 수치 실험을 통해 확인되었듯이, 계산 비용이 δ에 독립적인 고정밀도 해 근사에 적합하다.
• λn−q ∈ π/Θ ℤ* 또는 λn−q = 0일 때 수직 조건을 도입함으로써 보정항의 유일성이 회복되며, 이는 점근 시스템의 잘 정의됨을 보장한다.