[논문 리뷰] Constraint Energy Minimizing Generalized Multiscale Finite Element Method for high-contrast linear elasticity problem

Shubin Fu, Eric T. Chung|arXiv (Cornell University)|2018. 09. 11.

Advanced Mathematical Modeling in Engineering참고 문헌 19인용 수 7

한 줄 요약

이 논문은 고대비 선형 탄성 문제를 위한 오프라인 및 온라인 제약 에너지 최소화 일반화 다중스케일 유한요소법(CEM-GMsFEM)을 제안한다. 다중스케일 기저 함수는 더 큰 오버스ampled 도메인에서 에너지 최소화를 통해 구성되어 정확도와 견고성을 향상시키며, 잔차 기반 기저 함수를 사용한 적응형 온라인 보완을 도입하여 2D 및 3D 수치 실험에서 빠른 수렴과 높은 정확도를 달성하고 오차를 크게 감소시킨다.

ABSTRACT

In this paper, we consider the offline and online Constraint Energy Minimizing Generalized Mul- tiscale Finite Element Method (CEM-GMsFEM) for high-contrast linear elasticity problem. Offline basis construction starts with an auxiliary multiscale space by solving local spectral problems. We select eigenfunctions that correspond to a few small eigenvalues to form the auxiliary space. Using the auxiliary space, we solve a constraint energy minimization problem to construct offline multiscale spaces. The minimization problem is defined in the oversampling domain, which is larger than the target coarse block. To get a good approximation space, the oversampling domain should be large enough. We also propose a relaxed minimization problem to construct multiscale basis functions, which will yield more accurate and robust solution. To take into account the influence of input pa- rameters, such as source terms, we propose the construction of online multiscale basis and an adaptive enrichment algorithm. We provide extensive numerical experiments on 2D and 3D models to show the performance of the proposed method.

연구 동기 및 목표

• 미세 스케일 이질성이 존재해 계산 비용이 높은 고대비 선형 탄성 문제를 해결하는 데 도전하는 데 목적을 두며.
• 오버스ampled 도메인에서 제약 에너지 최소화를 통합하여 다중스케일 방법의 정확도와 견고성을 향상시키는 데 목적을 두며.
• 잔차 정보를 활용해 다중스케일 공간을 동적으로 보완하는 적응형 온라인 보완 전략을 개발하는 데 목적을 두며.
• 다양한 대비와 거친 메esh 크기를 가진 광범위한 2D 및 3D 수치 실험을 통해 방법의 효과성을 입증하는 데 목적을 두며.

제안 방법

• 최소 고유값에 해당하는 고유함수를 선택하여 국소 스펙트럼 문제를 풀어 보조 다중스케일 공간을 구성한다.
• 거친 블록보다 큰 오버스ampled 도메인에서 제약 에너지 최소화 문제를 풀어 오프라인 다중스케일 기저 함수를 형성한다.
• 제약 조건 없이 정확도와 안정성을 향상시키기 위해 최소화 문제의 완화된 버전을 도입한다.
• 거친 스케일 해의 잔차를 오버스ampled 도메인에서 국소 에너지 최소화 문제의 소스 항으로 사용하여 온라인 다중스케일 기저 함수를 구성한다.
• 오차 지표에 기반해 기저 보완 영역을 선택하는 적응형 보완 알고리즘을 구현하며, 잔차가 높은 영역을 우선순위로 삼는다.
• 균일 보완 및 적응형 보완 전략을 모두 사용하여 수렴 특성과 오차 감소 효율성을 비교한다.

실험 결과

연구 질문

• RQ1 오버스ampled 도메인에서 제약 에너지 최소화가 고대비 선형 탄성 문제에 대한 다중스케일 근사의 정확도와 견고성 향상에 기여하는가?
• RQ2 완화된 에너지 최소화 공식화는 제약 조건이 있는 버전에 비해 해의 정확도와 안정성 측면에서 어떻게 비교되는가?
• RQ3 특히 해의 기울기가 높은 영역에서, 적응형 온라인 보완은 균일 보완에 비해 오차 감소에 얼마나 효과적인가?
• RQ4 오버스ampling 층 수와 기저 함수의 수는 매질 대비에 대해 수렴 속도와 견고성에 어떤 영향을 미치는가?

주요 결과

• 완화된 CEM-GMsFEM는 특히 고대비 매질에서 표준 제약 버전보다 더 빠른 수렴과 높은 정확도를 달성한다.
• 적응형 온라인 보완은 균일 보완보다 L² 오차를 더 효과적으로 감소시키며, 기저 함수를 추가할수록 오차가 급격히 감소한다.
• 모델 1에서 오프라인 기저 함수 4개와 오버스ampling 층 4개를 사용할 경우, L² 오차는 1.55e-02에서 3.54e-07로 감소하여 높은 정확도를 입증한다.
• 모델 2에서는 오버스ampling 층 수를 3에서 4로 증가시켜 H¹ 오차를 1.10e-01에서 2.66e-02로 감소시켜 대비 변화에 대한 견고성 향상을 보였다.
• 다양한 거친 메쉬 크기에서 좋은 성능을 유지하며, H가 1/8에서 1/16로 감소함에 따라 수렴이 관찰되었다.
• 적응형 선택을 통한 온라인 보완은 27,798개의 자유도 이후 L² 오차를 3.10e-05로 감소시켜 오차 감소 속도 측면에서 균일 보완을 능가했다.