[논문 리뷰] A discontinuous Galerkin multiscale method for convection-diffusion problems

Daniel Elfverson|arXiv (Cornell University)|2015. 09. 11.

Advanced Mathematical Modeling in Engineering참고 문헌 23인용 수 3

한 줄 요약

이 논문은 극도로 이질적이고 거친 계수를 가진 대류-확산 문제를 위한 비연속 갈레르킨 국소 수직 분해 다중스케일 방법을 제안한다. 크기 $\mathcal{O}(H\log(H^{-1}))$ 의 국소 패치에서 보정 기저 함수 계산에 대류 항을 통합함으로써, 계수 변동에 영향을 받지 않는 수렴 속도를 확보하며, 특히 대류 지배 영역에서도 성능이 우수함을 수치 실험으로 검증하였다.

ABSTRACT

We propose an discontinuous Galerkin local orthogonal decomposition multiscale method for convection-diffusion problems with rough, heterogeneous, and highly varying coefficients. The properties of the multiscale method and the discontinuous Galerkin method allows us to better cope with multiscale features as well as interior/boundary layers in the solution. In the proposed method the trail and test spaces are spanned by a corrected basis computed on localized patches of size $\mathcal{O}(H\log(H^{-1}))$, where $H$ is the mesh size. We prove convergence rates independent of the variation in the coefficients and present numerical experiments which verify the analytical findings.

연구 동기 및 목표

• 극도로 이질적이고 거친, 급격히 변하는 계수를 가진 대류-확산 문제를 해결하는 데 도전하는 것.
• 다중스케일 특성과 경계/내부 계면을 해석하는 데 어려움을 겪는 전통적인 유한요소 방법의 한계를 극복하는 것.
• 척도 분리나 주기성 조건 없이 강한 대류 조건에서도 안정적이고 수렴 가능한 다중스케일 방법을 개발하는 것.
• 대류 항을 기저 함수 보정 과정에 통합하여 대류 지배 영역에서의 안정성을 향상시키는 것.
• 계수 변동에 영향을 받지 않는 수렴 속도를 확립하여 다양한 다중스케일 환경에서의 안정적인 성능을 보장하는 것.

제안 방법

• 거친 확산 계수와 발산이 없는 대류 계수를 가진 대류-확산 문제에 대해 비연속 갈레르킨(DG) 변분 형식을 수립한다.
• 해결 공간을 거시적 척도와 미시적 척도로 분해하기 위해 국소 수직 분해(LOD) 프레임워크를 적용한다.
• 확산 및 대류 항을 모두 포함한 보정 기저 함수를 크기 $\mathcal{O}(H\log(H^{-1}))$ 의 국소 패치에서 구성한다.
• 동일한 보정 기저 함수에서 유도된 시험 함수를 사용하는 페트로프-갈레르킨 접근법을 통해 방법의 안정성을 확보한다.
• 고유값 해법이나 큰 패치 지지 영역을 요구하지 않도록 하며, 선형 계산 복잡도를 유지한다.
• 대류 크기의 가정 하에, 계수 변동에 영향을 받지 않는 에너지 노름에서의 수렴을 증명한다.

실험 결과

연구 질문

• RQ1 거친, 비주기적 계수를 가진 대류-확산 문제에 대해 안정적이고 수렴 가능한 비연속 갈레르킨 다중스케일 방법을 구성할 수 있는가?
• RQ2 기본 LOD와 비교해 볼 때, 기저 함수 보정 과정에 대류 항을 포함시키는 것이 대류 지배 영역에서의 안정성에 어떻게 기여하는가?
• RQ3 보정 기저 함수의 최적 패치 크기는 무엇이어야 하며, 계수 변동에 영향을 받지 않는 안정적인 수렴을 보장할 수 있는가?
• RQ4 계수의 스펙트럼 조건수 $C_A = (\beta/\alpha)^{1/2}$ 에 영향을 받지 않는 수렴 속도를 달성할 수 있는가?
• RQ5 다중스케일 대류-확산 문제에 대해 기존의 유한요소 및 LOD 방법과 비교해 수치적으로 어떻게 성능을 발휘하는가?

주요 결과

• 제안된 DG-LOD 방법은 확산 및 대류 계수의 변동에 영향을 받지 않는 에너지 노름에서의 수렴 속도를 확보한다.
• 기저 함수 보정 과정에 대류 항을 포함시킴으로써, 기존 LOD가 실패하는 강한 대류 조건에서도 안정성과 수렴성을 확보한다.
• 크기 $\mathcal{O}(H\log(H^{-1}))$ 의 국소 패치를 사용하여 척도 분리나 주기성 조건 없이도 선형 계산 복잡도를 유지한다.
• 수치 실험을 통해 이론적 수렴 속도가 확인되었으며, 다양한 다중스케일 계수 설정에서 안정적인 성능을 보였다.
• 오차 추정치는 잔여항의 미세 척도 성분이 지수적으로 감쇠됨을 보여주며, 이 감쇠 속도는 패치 크기와 대류 크기에 의존한다.
• 잔여항의 반복적인 감쇠 추정을 통해 증명된 lin-정규화된 에너지 노름에서 최적의 수렴 순서를 유지하며, 계수 변동에 영향을 받지 않는 상수를 가진다.