[논문 리뷰] A fully discrete Calderon Calculus for two dimensional time harmonic waves

Sijiang L. Lu|arXiv (Cornell University)|2014. 01. 01.

Numerical methods in inverse problems참고 문헌 16인용 수 11

한 줄 요약

이 논문은 2차원 정현파 해르츠홀츠 방정식에 대해 격자 이격을 이용하고, 소스 밀도로 딜타 함수를 사용하며, 듀폴 근사에는 조각상수 함수를 사용하는 완전 이산화된 캘러존 미적분을 제안한다. 이로 인한 방법들은 연속 문제의 잘 정의됨이 보장될 경우 두 번째 차수 정확도(O(h²))를 달성하며, 수치 실험을 통해 검증된다.

ABSTRACT

In this paper, we present a fully discretized Calderon Calculus for the two dimensional Helmholtz equation. This full discretization can be understood as highly non-conforming Petrov-Galerkin methods, based on two staggered grids of mesh size h, Dirac delta distri- butions substituting acoustic charge densities and piecewise constant functions for approxi- mating acoustic dipole densities. The resulting numerical schemes from this calculus are all of order h 2 provided that the continuous equations are well posed. We finish by presenting some numerical experiments illustrating the performance of this discrete calculus.

연구 동기 및 목표

• 2차원 정현파 해르츠홀츠 방정식에 대한 캘러존 미적분의 완전 이산화된 공식화를 개발하는 것.
• 캘러존 미적분의 구조를 유지하는 방식으로 경계 적분 연산자를 이산화하는 데 도전하는 것.
• 연속 문제의 잘 정의됨이 보장될 경우 높은 순서 정확도(O(h²))를 보장하는 것.
• 2차원에서 경계 적분 방정식을 통한 산란 문제를 해결하기 위한 강력한 수치적 프레임워크를 제공하는 것.

제안 방법

• 이 방법은 메쉬 크기 h를 가진 두 개의 격자 이격 직각좌표 격자를 사용하여 문제를 이산화한다.
• 음향 전하 밀도는 한 쪽 격자 위에 놓인 델타 함수 분포를 사용하여 표현된다.
• 음향 듀폴 밀도는 두 번째 격자, 즉 이격된 격자 위에서 조각상수 함수를 사용하여 근사된다.
• 이 접근법은 안정성과 일致성을 보장하는 매우 비일치하는 페트로프-갈레르킨 방법으로 공식화된다.
• 이산 밀도와 적절한 적분 연산자를 조합하여 이산 캘러존 미적분이 구성된다.
• 결과로 얻어진 방법은 연속 캘러존 미적분의 구조를 유지하면서도 완전 이산화를 가능하게 한다.

실험 결과

연구 질문

• RQ1 연속 문제의 수학적 구조를 유지하면서 2차원 헬름홀츠 방정식에 대해 완전 이산화된 캘러존 미적분을 구성할 수 있는가?
• RQ2 연속 문제의 잘 정의됨이 보장될 경우, 결과로 얻어진 이산화 방법의 수렴 속도는 어떠한가?
• RQ3 격자 이격과 비일치 함수 공간은 수치 해의 안정성과 정확도에 어떤 영향을 미치는가?
• RQ4 델타 함수 분포와 조각상수 함수를 사용하여 듀폴 근사를 효과적으로 이산화된 미적분에 구현할 수 있는가?

주요 결과

• 제안된 완전 이산화된 캘러존 미적분은 연속 헬름홀츠 문제의 잘 정의됨이 보장될 경우 두 번째 차수 정확도를 달성하며, 오차는 O(h²)로 유계이다.
• 격자 이격과 비일치 함수 공간의 사용은 경계 적분 연산자의 안정적이고 일치하는 이산화를 가능하게 한다.
• 이 방법은 이산 환경에서 연속 이론과의 일致성을 보장하는 캘러존 미적분의 수학적 구조를 유지한다.
• 수치 실험은 이론적 수렴 속도를 확인하고 산란 문제에 대한 방법의 강건성을 입증한다.