[논문 리뷰] Generalized Multiscale Finite Element Method for Elasticity Equations

Eric T. Chung, Yalchin Efendiev외 1명|arXiv (Cornell University)|2014. 08. 25.

Advanced Mathematical Modeling in Engineering참고 문헌 13인용 수 8

한 줄 요약

이 논문은 고이종성, 고대비 매질에서 탄성 방정식을 해결하기 위한 일반화된 다스케일 유한요소법(GMsFEM)을 제안한다. 이는 지하 심층파 모델링에서 흔한 매질에서의 응용을 대상으로 하며, 국소 스펙트럼 분해를 통한 스냅샷 및 오프라인 공간 구축과 연속 및 비연속 갈레르킨 결합 기법을 활용하여, 특히 오버샘플링과 조화 확장 기반 스냅샷 공간을 사용할 경우 자유도 수를 크게 감소시키면서도 정확한 해를 달성한다.

ABSTRACT

In this paper, we discuss the application of Generalized Multiscale Finite Element Method (GMsFEM) to elasticity equation in heterogeneous media. Our applications are motivated by elastic wave propagation in subsurface where the subsurface properties can be highly heterogeneous and have high contrast. We present the construction of main ingredients for GMsFEM such as the snapshot space and offline spaces. The latter is constructed using local spectral decomposition in the snapshot space. The spectral decomposition is based on the analysis which is provided in the paper. We consider both continuous Galerkin and discontinuous Galerkin coupling of basis functions. Both approaches have their cons and pros. Continuous Galerkin methods allow avoiding penalty parameters though they involve partition of unity functions which can alter the properties of multiscale basis functions. On the other hand, discontinuous Galerkin techniques allow gluing multiscale basis functions without any modifications. Because basis functions are constructed independently from each other, this approach provides an advantage. We discuss the use of oversampling techniques that use snapshots in larger regions to construct the offline space. We provide numerical results to show that one can accurately approximate the solution using reduced number of degrees of freedom.

연구 동기 및 목표

• 스케일 분리가 없는 경우 고전적 균질화 기법이 실패하는 다기질, 고대비 매질에서 탄성 방정식을 해결하는 데 도전하는 것.
• 균열, 용굴, 공동 등 복잡한 지하 형성에서 정적인 탄성 문제를 위한 모델 축소 기법을 개발하는 것.
• 자유도 수를 극적으로 감소시켜 계산 비용을 줄이면서도 정확도를 유지하는 다스케일 방법을 설계하는 것.
• 연속 갈레르킨(CG) 및 비연속 갈레르킨(DG) 방법을 활용한 다스케일 기저 함수에 대한 두 가지 글로벌 결합 전략을 비교하고 분석하는 것.
• 오버샘플링 및 다양한 스냅샷 공간 구성 방식이 수렴성과 정확도에 미치는 영향을 조사하는 것.

제안 방법

• 각 거시격자 블록에서 모든 미세격자 함수 또는 조화 확장을 활용하여 국소 해 행동을 캡처하는 스냅샷 공간을 구축한다.
• 스냅샷 공간 내에서 국소 스펙트럼 분해를 수행하여 주요 변동 모드를 나타내는 다스케일 기저 함수를 추출한다.
• 스펙트럼 분해에서 가장 중요한 모드만 유지하여 차원을 감소시키는 오프라인 공간을 구성한다.
• 두 가지 글로벌 결합 방법을 구현한다: 분할 단위 함수를 사용하는 연속 갈레르킨(CG) 및 대칭 내부 페널티 방법을 사용하는 비연속 갈레르킨(DG).
• 고해상도 요소를 초과하는 영역으로 고유값 문제의 정의역을 확장하여 오버샘플링을 적용함으로써 다스케일 기저 함수의 품질을 향상시키고 차원을 감소시킨다.
• 결과적으로 도출된 다스케일 기저 함수를 사용하여 자유도 수를 최소화하는 감소된 차원 공간에서 전역 탄성 문제를 해결한다.

실험 결과

연구 질문

• RQ1 고대비 및 스케일 분리가 없는 매질에서 탄성 방정식에 대한 다스케일 기저 함수를 효율적으로 구성하는 방법은 무엇인가?
• RQ2 GMsFEM의 탄성 문제에 적용할 때, 연속 갈레르킨 대비 비연속 갈레르킨 결합의 정확도 및 안정성에 대한 상대적 성능은 어떠한가?
• RQ3 오버샘플링은 탄성 문제에 대한 GMsFEM의 수렴성과 강건성에 어느 정도 향상시키는가?
• RQ4 다양한 스냅샷 공간 선택(예: 조화 확장 대비 모든 미세격자 함수)은 다스케일 해의 정확도에 어떤 영향을 미치는가?
• RQ5 제안된 GMsFEM는 고대비 탄성 문제에서 자유도 수를 줄였을 때 최적의 수렴 속도를 달성할 수 있는가?

주요 결과

• 표준 미세격자 스냅샷 공간에 비해 조화 확장 기반 스냅샷 공간을 사용할 경우 더 나은 수렴 행동을 보였다.
• 오버샘플링 기법은 다스케일 기저 함수의 품질을 향상시켜 GMsFEM의 정확도와 수렴 속도를 크게 향상시켰다.
• 비연속 갈레르킨 접근법은 기저 함수를 수정 없이 독립적으로 구성할 수 있어 연속 갈레르킨 대비 계산 및 구현 측면에서 유리한 점을 가졌다.
• 연속 갈레르킨 방법은 분할 단위 함수가 필요하며, 이는 다스케일 기저 함수의 성질를 변경할 수 있지만 페널티 파rameter를 피한다.
• 수치 결과는 고대비 매질에서도 자유도 수를 줄였음에도 불구하고 GMsFEM가 미세스케일 해를 정확하게 근사함을 확인하였다.
• 이론적 오차 한계는 스펙트럼 갭과 경계 유량 항에 의존하며, 오프라인 공간이 주요 모드를 포괄할 경우 수렴 속도가 향상됨을 보여주었다.