[논문 리뷰] An embedded corrector problem to approximate the homogenized coefficients of an elliptic equation

Éric Cancès, Virginie Ehrlacher외 2명|arXiv (Cornell University)|2014. 12. 19.

Advanced Mathematical Modeling in Engineering참고 문헌 16인용 수 18

한 줄 요약

이 논문은 고주기 성분을 가진 타원형 PDE에서 동질화 계수를 근사하기 위한 임bedded 보정자 문제를 제안한다. 전체 공간에서 일정한 확산 행렬을 가진 확산 방정식을 구해, 반경 R의 구 외부에서 R → ∞일 때 동질화 행렬 $A^\star$의 세 가지 수렴 근사값을 정의함으로써, 기존의 초세포 방법에 비해 수치적으로 효율적인 대안을 제공한다.

ABSTRACT

We consider a diffusion equation with highly oscillatory coefficients that admits a homogenized limit. As an alternative to standard corrector problems, we introduce here an embedded corrector problem, written as a diffusion equation in the whole space in which the diffusion matrix is uniform outside some ball of radius $R$. Using that problem, we next introduce three approximations of the homogenized coefficients. These approximations, which are variants of the standard approximations obtained using truncated (supercell) corrector problems, are shown to converge when $R \ o \\infty$. We also discuss efficient numerical methods to solve the embedded corrector problem.

연구 동기 및 목표

• 고주기 성분을 가진 타원형 PDE에서 동질화 계수 $A^\star$를 계산하는 데 도전하는 문제를 다루며, 특히 비주기적 또는 무작위 정적 설정에서의 적용을 고려한다.
• 전체 공간에서 보정자 문제를 해결하는 데 어려움이 있음을 고려해, 새로운 임bedded 보정자 설정을 제안한다.
• 동질화된 행렬 $A^\star$의 세 가지 새로운 근사값을 개발하여, 도메인 반경 $R \to \infty$일 때 진정한 동질화된 행렬로 수렴하도록 한다.
• 특히 경계 적분 형식을 통해 기존의 초세포(자르기된 보정자) 방법에 비해 수치적으로 효율적인 대안을 제공한다.
• 제안된 근사값의 이론적 수렴성을 확립하고, 등방성 동질화와 같은 특수 케이스에서 해의 존재성을 증명한다.

제안 방법

• 전체 공간 $\mathbb{R}^d$에서의 확산 방정식으로서 임bedded 보정자 문제를 설정하며, 구 $B_R$ 내부에서는 비균일한 확산 행렬 $\mathbb{A}$를, 외부에서는 일정한 행렬 $A^R_1$을 사용한다.
• 첫 번째 근사값 $A^R_1$은 임bedded 보정자 문제의 해를 이용해 $B_R$ 위에서의 유량 평균으로 정의한다.
• 두 번째 근사값 $A^R_2$는 비균일 시스템의 에너지와 동질성 시스템의 에너지를 동일시하는 형식적 에너지 등식을 통해 정의한다. 이때 동질성 시스템의 행렬은 $A^R_2$이다.
• 세 번째 근사값 $A^R_3$은 자기일관성 방정식 $A^R_3 = G^{R,\mathbb{A}}(A^R_3)$를 통해 도입한다. 여기서 $G^{R,\mathbb{A}}$는 행렬을 임bedded 보정자 해로부터 유량 평균으로 매핑하는 함수이다.
• 임bedded 보정자 문제의 효율적 수치적 해법을 위해 경계 적분 형식을 활용한다.
• 가정 2.1 하에, $R \to \infty$일 때 세 근사값 모두 진짜 동질화된 행렬 $A^\star$로 수렴함을 증명한다.

실험 결과

연구 질문

• RQ1 고주기 성분을 가진 타원형 PDE에서 동질화 계수를 근사하기 위해 새로운 임bedded 보정자 문제를 설정할 수 있는가?
• RQ2 제안된 세 근사값인 $A^R_1$, $A^R_2$, $A^R_3$이 도메인 반경 $R \to \infty$일 때 진짜 동질화된 행렬 $A^\star$로 수렴하는가?
• RQ3 경계 적분 형식과 같은 수치적 방법을 통해 임bedded 보정자 문제를 효율적으로 풀 수 있는가?
• RQ4 자기일관성 방정식 $A^R_3 = G^{R,\mathbb{A}}(A^R_3)$가 어떤 조건에서 해를 가지며, 이로 인해 수렴하는 근사값을 얻을 수 있는가?
• RQ5 제안된 근사값들은 정규 초세포 방법과 비교해 정확도와 계산 효율성 측면에서 어떻게 다른가?

주요 결과

• 구 $B_R$ 위에서의 유량 평균으로 정의된 근사값 $A^R_1$은 가정 2.1 하에 $R \to \infty$일 때 $A^\star$로 수렴한다.
• 에너지 등식을 통해 정의된 근사값 $A^R_2$ 역시 $R \to \infty$일 때 $A^\star$로 수렴한다.
• 자기일관성 근사값 $A^R_3$은 $A^R_3 = G^{R,\mathbb{A}}(A^R_3)$를 만족하며, 수열 $R_k \to \infty$이면서 해가 존재할 경우 $A^\star$로 수렴한다.
• 등방성 케이스에서 $A^\star = a^\star I_d$일 경우, 모든 $R > 0$에 대해 스칼라 해 $a^R_3$이 존재하며 $R \to \infty$일 때 $a^R_3 \to a^\star$를 만족한다.
• 일차원 문제에서는 $A^R_1$, $A^R_2$, $A^R_3$의 명시적 표현을 유도할 수 있으며, 이들이 모두 $A^\star$로 수렴함을 보일 수 있다.
• 만약 $\mathbb{A}$가 $B_R$ 내부에서 일정하다면, 세 근사값은 모두 동일한 상수 행렬이 되며: $A^R_1 = A^R_2 = A^R_3 = A$.