[논문 리뷰] Computational multiscale methods for parabolic wave approximations in heterogeneous media

Eric T. Chung, Yalchin Efendiev외 2명|arXiv (Cornell University)|2021. 04. 06.

Advanced Mathematical Modeling in Engineering참고 문헌 47인용 수 9

한 줄 요약

이 논문은 후진 오일러 시간 이산화, 공간 다중스케일 축소를 위한 제약 에너지 최소화 일반 다중스케일 유한요소법(CEM-GMsFEM), 그리고 진동 방향에 따른 차원 축소를 위한 적절한 직교 분해(POD)를 조합하여, 극도로 이방성 매질에서 편평한 파동 근사해를 구하는 계산 다중스케일 방법을 제안한다. 수치 실험에서 상대 오차가 약 10⁻³에서 10⁻⁴ 수준을 유지하며 안정적이고 정확한 해를 도출함으로써, 변수성과 층상 매질에 대해 효율성과 강건성을 입증한다.

ABSTRACT

In this paper, we develop a computational multiscale to solve the parabolic wave approximation with heterogeneous and variable media. Parabolic wave approximation is a technique to approximate the full wave equation. One benefit of the method is that: one wave propagation direction can be taken as an evolution direction, and we then can discretize it using a classical scheme like Backward Euler. Consequently, we obtain a set of quasi-gas-dynamic (QGD) models with different heterogeneous permeability fields. Then, we employ constraint energy minimization generalized multiscale finite element method (CEM-GMsFEM) to perform spatial discretization for the problem. The resulting system can be solved by combining the central difference in time evolution. Due to the variable media, we apply the technique of proper orthogonal decomposition (POD) to further the dimension of the problem and solve the corresponding model problem in the POD space instead of in the whole multiscale space spanned by all possible multiscale basis functions. We prove the stability of the full discretization scheme and give the convergence analysis of the proposed approximation scheme. Numerical results verify the effectiveness of the proposed method.

연구 동기 및 목표

• 극도로 이방성이고 변수적인 매질에서 편평한 파동 근사해를 효율적으로 풀 수 있는 계산 다중스케일 프레임워크를 개발하는 것.
• 공간적 이방성과 다중스케일 특성으로 인해 각 z-레벨에서 준기체역학적(QGD) 모델을 풀이하는 데 발생하는 계산적 과제를 해결하는 것.
• 공간과 z-진동 방향 양쪽에서 CEM-GMsFEM와 POD를 각각 활용하여 문제의 차원을 축소하는 것.
• 편평한 파동 모델에 대한 전체 이산화 스킴의 안정성과 수렴성을 확립하는 것.
• 층상 및 Marmousi 유형의 이방성 매질에서 방법의 정확성과 효율성을 수치적으로 검증하는 것.

제안 방법

• z-미분을 후진 오일러 스킴으로 이산화하여, 각 z-레벨에서 2차원 QGD 모델의 연속적 시퀀스로 문제를 변환한다.
• 거친 메esh에서 공간 이산화를 위해 CEM-GMsFEM를 사용하여 국소적 이방성을 반영하는 다중스케일 기저 함수를 구성한다.
• 유도된 QGD 모델에서 시간 미분 ∂ₜₜv 와 ∂ₜv 를 중심 차분 스킴으로 이산화한다.
• 에너지 노름에서 투영 오차를 최소화하기 위해 결합된 CEM 기저 함수 위에서 적절한 직교 분해(POD)를 적용하여 저차원 모델을 구성한다.
• POD 공간에서 축소된 시스템을 풀어 자유도를 크게 감소시키면서도 정확성을 유지한다.
• 후진 오일러 스킴의 무조건적 안정성을 증명하고, CEM-GMsFEM 거친 메쉬 수렴 결과를 바탕으로 전체 스킴의 수렴 분석을 제공한다.

실험 결과

연구 질문

• RQ1 다중스케일 유한요소법이 변수성 매질에서 편평한 파동 근사해의 공간 이방성을 효과적으로 다룰 수 있는가?
• RQ2 시간에 따라 변화하는 다중스케일 문제에서 정확도를 훼손하지 않고 z-진동 방향의 차원을 어떻게 축소할 수 있는가?
• RQ3 CEM-GMsFEM와 POD를 조합했을 때 전체 이산화 스킴의 안정성과 수렴성에 어떤 영향을 미치는가?
• RQ4 층상 및 복잡한 이방성 매질(예: Marmousi 필드)에서 상대 오차는 어떻게 척도화되는가?
• RQ5 제안된 방법은 미세 메쉬 해석기법에 비해 상당히 낮은 계산 비용으로 고정밀 해를 달성할 수 있는가?

주요 결과

• 층상 이방성 매질을 가진 첫 번째 실험에서 제안된 방법은 상대 오차 약 10⁻³을 달성한다.
• Marmousi 실험에서는 최종 시간 단계에서 상대 오차가 약 3.9×10⁻⁴로 나타나, 복잡하고 현실적인 매질에서도 높은 정확도를 입증한다.
• 첫 번째 실험에서 300개의 CEM 기저 함수 중 100개의 POD 기저 함수를 사용함으로써 주요 모드를 최소 투영 오차로 잘 포괄한다.
• Marmousi 케이스에서는 300개의 CEM 기저 함수 중 50개의 POD 기저 함수를 선택하여 효과적인 차원 축소를 입증한다.
• 후진 오일러 스킴은 무조건적 안정성을 보장하며, 제안된 프레임워크 하에서 전체 스킴이 안정적이고 수렴함을 증명하였다.
• 수치 결과는 구조적 및 비구조적 이방성 매질, 특히 현실적인 지질학적 구성까지 포함하여 방법의 효율성과 강건성을 확인한다.