불규칙 지형을 위한 새로운 침투 경계법과 음향 파동 방정식 적용

초록

본 논문은 불규칙한 지형을 정확히 표현할 수 있는 일반화된 침투 경계(Immersed Boundary) 기법을 제안한다. N차원 테일러 급수를 이용해 경계면에서의 파동 변수와 그 미분을 외삽하고, 이를 하이퍼구형 지원 영역과 수정된 유도 스텐실에 적용한다. 자유면·강체면 모두에 대해 1차·2차 음향 파동 방정식에 적용 가능하며, 2D·3D 실험과 수렴 검증을 통해 높은 정확도와 자동화 가능성을 입증한다.

상세 분석

이 연구는 구조화된 유한차분(FD) 격자와 불규칙 지형 사이의 불일치를 해결하기 위해, 방정식에 독립적인 침투 경계(Immersed Boundary, IB) 프레임워크를 설계하였다. 핵심 아이디어는 경계면 근처의 격자점들에 대해 N차원 테일러 급수를 이용해 물리량을 외삽하고, 외삽식에 경계조건(예: 자유면의 트랙션 제로, 강체면의 속도 고정)을 직접 대입함으로써 미분 연산자를 수정하는 것이다. 이를 위해 ‘하이퍼구형 지원 영역(hyperspherical support region)’을 정의하여, 경계면에서 일정 거리 이내에 있는 모든 격자점을 포함시킨다. 이렇게 하면 스텐실이 비대칭이 되더라도 기존 FD 연산과 동일한 차수와 정확도를 유지할 수 있다.

특히, 외삽 과정은 차수 N을 자유롭게 선택할 수 있어 1차·2차 음향 파동 방정식 모두에 적용 가능하다. 1차 방정식에서는 압력만을 다루지만, 2차 방정식에서는 압력과 속도 벡터를 동시에 외삽함으로써 벡터 경계조건도 근사 없이 구현한다. 경계조건을 미분 연산 전에 직접 삽입함으로써 ‘방정식 전처리’ 단계에서 근사오차가 발생하지 않는다.

코드 생성(paradigm)과의 연계도 강조한다. 자동 코드 생성 시스템에 이 IB 모듈을 플러그인 형태로 삽입하면, 사용자는 지형 데이터와 경계조건만 제공하면 복잡한 스텐실 수정 과정을 자동으로 수행할 수 있다. 이는 기존 방법들이 특정 지형 형태(예: 계단식, 평탄화된 표면)나 특정 차수에만 제한된 점을 극복한다.

수치 실험에서는 2차원 및 3차원 실제 지형 데이터를 이용해 자유면과 강체면을 각각 구현하였다. 해석해와의 비교, 격자 간격을 변화시킨 수렴 테스트에서 L2 오차가 기대 차수와 일치함을 확인했다. 특히, 경계 근처에서의 파동 반사와 전파 패턴이 물리적으로 정확히 재현되었으며, 기존 ‘가상 격자’나 ‘셀레니티드 마스크’ 방식에 비해 계산 비용이 크게 증가하지 않았다.

결과적으로, 제안된 IB 기법은 (1) 경계조건을 직접 외삽식에 포함해 근사오차를 최소화, (2) 차수와 차원을 자유롭게 선택 가능, (3) 자동 코드 생성과 연계해 실무 적용성을 높이며, (4) 자유면·강체면 모두를 동일한 프레임워크로 처리한다는 장점을 제공한다. 이러한 특성은 역시계열 전이동(RTM)·전파 전도 전산(FWI) 등 고해상도 지구물리 탐사에 바로 적용될 수 있는 실용적 기반을 마련한다.