탄성파 PML 영역에서 변환 탄성 텐서의 대칭 공식

초록

본 논문은 주파수 영역에서 탄성파 전파를 다루는 완전 매치 레이어(PML) 구현 시, 기존 변환 운동량 방정식으로부터 도출된 탄성 텐서가 비대칭성을 띠는 문제를 지적한다. 저자는 약한 형태(weak form)를 이용해 대칭성을 보존하는 새로운 변환 탄성 텐서 공식을 제시함으로써, 전통적인 탄성 이론을 PML 영역에서도 그대로 적용할 수 있음을 증명한다.

상세 요약

PML은 무한 매질을 유한 영역으로 잘라내어 인위적인 경계 반사를 최소화하는 기법으로, 전자기·음향·탄성 파동 해석에 널리 사용된다. 탄성 파동 문제는 2계 텐서인 탄성계수 C_{ijkl}와 변위 u_i가 결합된 2차 편미분 방정식으로 기술되며, 이때 C_{ijkl}=C_{klij}=C_{jikl}=C_{ijlk}와 같은 완전한 대칭성이 물리적 의미(에너지 보존, 작용-반작용 법칙)를 보장한다. 기존 연구에서는 복소 좌표 스트레칭을 적용한 변환 운동량 방정식으로부터 새로운 탄성 텐서를 도출했지만, 변환 행렬 A와 그 역행렬 A^{-1}의 곱셈 순서가 비대칭적으로 처리되어 C’{ijkl}=A{im}A_{jn}C_{mnop}A^{-1}{ok}A^{-1}{pl} 형태가 된다. 이 식은 일반적인 경우 C’{ijkl}≠C’{klij}를 만족하지 않아, 변환된 매질이 비탄성(비대칭) 특성을 갖게 된다.

비대칭 텐서는 유한 요소 해석에서 행렬의 강정성(positive‑definiteness)을 손상시켜 수치적 불안정성을 초래하고, 물리적으로는 에너지 흐름이 비보존적이라는 모순을 내포한다. 저자는 이러한 이론적 모순을 해결하기 위해, 변환 전후의 약한 형태를 동일하게 유지하는 접근법을 제안한다. 구체적으로, 원래의 변위 시험함수 v_i와 변환된 시험함수 \tilde v_i를 복소 좌표 변환에 따라 매핑하고, 변환된 변위와 시험함수 사이의 내적을 보존하도록 가중치를 도입한다. 그 결과 얻어지는 변환 탄성 텐서는

\