[논문 리뷰] A short-memory operator splitting scheme for constant-Q viscoelastic wave equation

Yunfeng Xiong, Xu Guo|arXiv (Cornell University)|2021. 07. 11.

Fractional Differential Equations Solutions참고 문헌 42인용 수 13

한 줄 요약

이 논문은 다수의 0에 가까운 순서의 Caputo 도함수를 갖는 분수계 응력-변형률 관계를 사용하여 일정-Q 점탄성파 방정식을 해결하기 위한 단기 기억 연산자 분할(SMOS) 기법을 제안한다. 약간의 특이핵을 갖는 라그에르 함수를 통해 약한 특이핵을 국소화하는 확장 방법을 통해 문제를 재구성함으로써, 정확도를 유지하면서도 메모리와 계산 비용을 감소시킨다. 주요 기여는 β > 1인 스케일링 기법을 통해 시간에 따른 정규화 오차의 증가를 억제함으로써, 1차원 및 2차원 파동 문제에서 높은 효율성과 정확도로 안정적인 장기 시뮬레이션을 가능하게 한다.

ABSTRACT

We propose a short-memory operator splitting scheme for solving the constant-Q wave equation, where the fractional stress-strain relation contains multiple Caputo fractional derivatives with order much smaller than 1. The key is to exploit its extension problem by converting the flat singular kernels into strongly localized ones, so that the major contribution of weakly singular integrals over a semi-infinite interval can be captured by a few Laguerre functions with proper asymptotic behavior. Despite its success in reducing both memory requirement and arithmetic complexity, we show that numerical accuracy under prescribed memory variables may deteriorate in time due to the dynamical increments of projection errors. Fortunately, it can be considerably alleviated by introducing a suitable scaling factor $\beta > 1$ and pushing the collocation points closer to origin. An operator splitting scheme is introduced to solve the resulting set of equations, where the auxiliary dynamics can be solved exactly, so that it gets rid of the numerical stiffness and discretization errors. Numerical experiments on both 1-D diffusive wave equation and 2-D constant-Q $P$- and $S$-wave equations are presented to validate the accuracy and efficiency of the proposed scheme.

연구 동기 및 목표

• 매우 작은 순서의 분수도함수로 인해 발생하는 일정-Q 점탄성파 방정식 해법의 높은 메모리 및 계산 비용 문제를 해결하기 위해.
• 비감쇠 성분이 존재하는 원천으로 인해 축적되는 정규화 오차로 인해 발생하는 단기 기억 방법의 시간에 따른 정밀도 저하 문제를 극복하기 위해.
• 내재된 감쇠를 갖는 파동 모델링에서 경직성과 이산화 오차를 피하는 효율적이고 안정적이며 정확한 수치 기법을 개발하기 위해.
• 분수도함수 근사에 대한 일반화된 라그에르 보간 프레임워크 내에서 라그에르 스펙트럼 방법의 엄밀한 분석을 제공하기 위해.
• 실용적인 지질학적 응용을 고려하여 1차원 확산파 및 2차원 일정-Q P-파 및 S-파 방정식에 대해 방법을 검증하기 위해.

제안 방법

• 보조 기억 변수 y ∈ (0, ∞)를 통해 평탄한 특이핵을 강하게 국소화하는 확장 문제를 통해 일정-Q 파동 방정식을 재구성하기 위해.
• 반무한 구간 위의 약한 특이적 적분을 근사하기 위해 일반화된 라그에르-가우스 구적법을 적용하여, 소수의 콜로케이션 점으로 주요 기여를 포괄하기 위해.
• 라그에르 함수에 β > 1인 스케일링 인자를 도입하여 콜로케이션 점들을 원점에 더 가깝게 이동시켜, 시간에 따른 정규화 오차의 동적 증가를 줄이기 위해.
• 파동 역학과 보조 기억 역학을 분리하는 연산자 분할 기법을 구현하여 보조 시스템의 정확한 해를 도출하고 수치적 경직성을 제거하기 위해.
• 메인아르디 함수 평가를 위한 하이브리드 수치적 접근법을 사용: 소수의 x에 대해서는 라그에르-가우스 구적법, 큰 x에 대해서는 점근 전개를 적용하고, 정확도와 안정성을 확보하기 위해 임계값 기반 전환 기법을 도입하기 위해.
• 유안-아그라왈 방법을 이론적 분석과 결합하여 단기 기억 근사의 타당성을 입증하고, 일반화된 라그에르 보간의 맥락에서 오차 전파를 정량화하기 위해.

실험 결과

연구 질문

• RQ1 다수의 0에 가까운 순서의 Caputo 분수도함수를 갖는 일정-Q 점탄성파 방정식에 대해 단기 기억 원리가 효과적으로 적용될 수 있는가?
• RQ2 라그에르 스펙트럼 방법 기반의 단기 기억 방법에서 정규화 오차의 동적 증가가 장기 수치 정밀도에 어떤 영향을 미치는가?
• RQ3 β > 1인 스케일링 기법이 이러한 방법에서 정규화 오차의 시간에 따른 악화 행동을 현저히 억제할 수 있는가?
• RQ4 연산자 분할이 경직성이나 이산화 오차를 도입하지 않고도 안정적이고 정확한 시간 적분을 가능하게 하는가?
• RQ5 다양한 매개변수 영역에서 메인아르디 함수 평가를 위한 최적의 전략은 무엇인가, 이를 통해 강건성과 정확도를 보장할 수 있는가?

주요 결과

• 제안된 SMOS 기법은 메모리 사용량을 최소화하고 산술 복잡도를 낮추면서도 1차원 확산파 및 2차원 일정-Q P-파 및 S-파 방정식의 시뮬레이션에서 높은 정확도와 효율성을 달성한다.
• β > 1인 스케일링 인자를 도입함으로써 정규화 오차의 시간에 따른 악화 행동이 효과적으로 억제되어 장기 시뮬레이션의 안정성이 확보된다.
• 수치 실험 결과, 소수의 기억 변수(예: 10~20개)로도 높은 정확도를 유지함을 확인하였으며, 이는 3차원 문제에서 약 90GB에서까지 감소된 저장 요구량을 가능하게 한다.
• 연산자 분할 접근법은 파동 역학과 기억 역학을 성공적으로 분리하여 보조 시스템의 정확한 통합을 가능하게 하고, 경직성 관련 오차를 제거한다.
• 소수의 x에 대해서는 라그에르 구적법, 큰 x에 대해서는 점근 전개를 사용하는 하이브리드 평가 전략은 임계값 기반 전환 규칙을 통해 모든 매개변수 영역에서 높은 정확도를 보장한다.
• 이론적 분석을 통해 방법의 오차 증가가 제어 가능하며, 스케일링 기법이 표준 단기 기억 근사에서 발생하는 본질적 불안정성에 실용적인 해결책을 제공한다는 것이 확인되었다.