[논문 리뷰] Perfectly matched layer for second-order time-domain elastic wave equation: formulation and stability

Hisham Assi, R.S.C. Cobbold|arXiv (Cornell University)|2013. 12. 13.

Electromagnetic Simulation and Numerical Methods참고 문헌 32인용 수 3

한 줄 요약

이 논문은 복소 좌표 스트레칭과 이중 매개변수 스트레칭 함수를 사용하여 2차원 이방성 고체에서 탄성파 시뮬레이션을 위한 흡수경계조건(PML)에 대해 2차 시간영역 형식을 제안한다. 이는 불안정한 이방성 매질에서 스케일링 매개변수를 증가시킴으로써 우수한 안정성을 확보하고, 두 개의 2차 편미분방정식과 네 개의 보조방정식만을 사용하여 총 여섯 개의 방정식로 구성되어 계산 비용을 감소시키며, 수치 반사가 최소화된 안정적인 장기 시뮬레이션을 가능하게 한다.

ABSTRACT

A time domain system of equations is proposed to model elastic wave propagation in an unbounded two-dimensional anisotropic solid using perfectly matched layer (PML). Starting from a system of first-order frequency domain stress-velocity equations and using complex coordinate stretching approach with a two-parameter stretch function, a second-order formulation is obtained. The final system, which consists of just two second order equations along with four auxiliary equations, is smaller than existing formulations, thereby simplifying the problem and reducing the computational cost. The discrete stability of the solutions for a given mesh size is examined with the help of a plane-wave analysis of the corresponding continuous problem. It is shown that increasing the scaling parameter of the stretch function leads to significant stability improvements for certain anisotropic media that have known issues. Numerical computations for different isotropic and anisotropic media are used to illustrate the results.

연구 동기 및 목표

• 비유한 탄성파 시뮬레이션에서 발생하는 비합리적인 반사 현상을 안정적인 PML 방법으로 해결하기 위해.
• 시간영역 PML 형식에서 방정식 수를 최소화하여 계산 비용을 감소시키고 구현을 단순화하기 위해.
• 기존의 PML 형식에서 불안정한 것으로 알려진 이방성 매질에 대해 이산 안정성을 향상시키기 위해.
• 복잡한 다중 매개변수 조정 없이 스케일링 매개변수 조정만으로도 안정성을 향상시킬 수 있는 간단하고 효과적인 방법을 제공하기 위해.

제안 방법

• 복소 좌표 스트레칭과 이중 매개변수 스트레칭 함수를 사용하여 주파수 영역의 응력-속도 1차 방정식에서 시간영역 2차 편미분방정식계를 유도한다.
• 2차 구조를 유지하고 시간영역 구현을 가능하게 하기 위해 네 개의 보조방정식을 도입한다.
• 이산 안정성을 평가하고 PML 영역 내에서 불안정한 모드를 식별하기 위해 평면파 분석을 적용한다.
• PML 전역에서 1에서 더 높은 값(예: 20, 90)으로 증가하는 가변 스케일링 매개변수(α̃j)를 사용하여 불안정한 파동 모드를 억제한다.
• COMSOL Multiphysics를 사용한 유한요소법(FEM)을 활용하여 수치적 검증 및 등방성 및 이방성 매질에서의 파동 전파 시뮬레이션을 수행한다.
• 물리적 영역 경계에서 수치 반사를 줄이기 위해 PML 계수에 다항식 전이를 적용하여 부드러운 변화를 보장한다.

실험 결과

연구 질문

• RQ1 기존 방법에 비해 방정식 수를 줄일 수 있는 2차원 이방성 고체에서 탄성파 전파를 위한 2차 시간영역 PML 형식을 유도할 수 있는가?
• RQ2 스트레칭 함수 매개변수의 선택이 이방성 매질에서 PML의 이산 안정성에 어떤 영향을 미치는가?
• RQ3 스케일링 매개변수 α̃j를 증가시킴으로써 기존에 불안정한 것으로 알려진 이방성 매질을 PML 내에서 안정화시킬 수 있는가? 이 경우 복잡한 다중 매개변수 형식이 필요하지 않은가?
• RQ4 장기 안정성과 계산 효율성 측면에서 제안된 형식이 기존의 PML에 비해 어느 정도 뛰어나게 성능을 발휘하는가?
• RQ5 메쉬 크기와 파장 입사각이 2차 형식에서 PML의 안정성에 어떤 영향을 미치는가?

주요 결과

• 제안된 형식은 총 여섯 개의 방정식(두 개의 2차 편미분방정식과 네 개의 보조방정식)만을 사용하여 탄성파에 대한 시간영역 PML 형식 중에서 보고된 바에서 가장 작은 크기를 갖는다.
• 물질 IV의 경우 최적화된 스케일링 매개변수를 사용하여 20ms가 넘는 시뮬레이션 동안에도 PML가 안정성을 유지함을 확인하였다. 이는 장기 안정성을 입증한다.
• 물질 III에서 스케일링 매개변수 α̃j를 1에서 20 및 90으로 증가시키면 불안정성이 크게 억제되며, 특히 x₁ 방향에서 두드러진다. 그러나 x₂ 방향에서는 잔여 불안정성이 존재한다.
• 물질 V의 경우 α̃j를 10 및 1로 증가시킴으로써 모든 방향에서 안정성이 향상되었으며, 이는 이전에 낮은 값에서 관측된 불안정성은 사라졌다.
• 평면파 분석을 통해 높은 α̃j 값이 불안정한 근을 이산적으로 해석 가능한 공간 주파수 범위 외부로 이동시켜 안정성 향상을 설명할 수 있었다.
• 에너지 변화도 그래프(그림 11)는 높은 α̃j 값을 사용할수록 물리적 영역 내 반사 에너지가 뚜렷이 감소하는 것을 보여주며, 이는 향상된 에너지 흡수 및 안정성 확보를 확인한다.