[논문 리뷰] Computing homogenized coefficients via multiscale representation and hierarchical hybrid grids

Antti Hannukainen, Jean-Christophe Mourrat외 1명|arXiv (Cornell University)|2019. 05. 16.

Advanced Mathematical Modeling in Engineering참고 문헌 44인용 수 4

한 줄 요약

이 논문은 무작위로 정의된 조각별로 일정한 계수를 가진 발산형 연산자에 대해 효율적으로 동질화 계수를 계산하기 위해 계층적 하이브리드 메esh를 사용하는 다스케일 유한요소 방법을 제시한다. 동질화 행렬을 척도에 따라 분리된 항으로 분해하고 국소화된 계산 및 적응형 메쉬 다듬힘을 활용함으로써, 중간 정도의 탄성 대비 비율과 몇 퍼센트의 정밀도를 갖는 2차원 및 3차원 문제에서 랩탑에서 1초 이내에서 수분 이내의 수렴을 달성한다.

ABSTRACT

We present an efficient method for the computation of homogenized coefficients of divergence-form operators with random coefficients. The approach is based on a multiscale representation of the homogenized coefficients. We then implement the method numerically using a finite-element method with hierarchical hybrid grids, which is a semi-implicit method allowing for significant gains in memory usage and execution time. Finally, we demonstrate the efficiency of our approach on two- and three-dimensional examples, for piecewise-constant coefficients with corner discontinuities. For moderate ellipticity contrast and for a precision of a few percentage points, our method allows to compute the homogenized coefficients on a laptop computer in a few seconds, in two dimensions, or in a few minutes, in three dimensions.

연구 동기 및 목표

• 무작위로 정의된 조각별로 일정한 계수를 가진 발산형 타원형 연산자의 동질화 계수를 계산하기 위한 효율적인 수치적 방법을 개발하는 것.
• 고정된 탄성 대비 비율을 갖는 랜덤 매질에서 타일링 모서리에서 발생하는 특이성으로 인한 계산 도전 과제를 해결하는 것.
• 적응형 척도 분리 메쉬를 사용하여 표준 하드웨어에서 고정밀도 동질화 계수 추정을 가능하게 하는 것.
• 재계산이 필요할 경우에도 다시 계산하지 않고도 안정적이고 재정밀화 가능한 알고리즘을 제공하는 것.
• 2차원 및 3차원 예제에서 탄성 대비 비율이 다양할 경우에 대해 방법을 검증하여 랩탑에서 1초 이내에서 수분 이내의 실행 시간을 보여주는 것.

제안 방법

• 동질화 행렬을 다양한 길이 척도의 기여도로 분해하는 다스케일 분해 기법을 사용하며, 점점 더 큰 영역에서 국소화된 PDE의 재귀적 해를 통해 정의된다.
• 함수 $ v_k $의 순서열이 $ (2^{-k} - \nabla\cdot\mathbf{a}\nabla)v_k = 2^{-k}v_{k-1} $ 를 만족하며, 점점 더 거시적인 척도의 영향을 포착한다.
• 제곱형태의 평균화를 위해 $ \chi_{r_k} $ 라는 컴팩트 지지 집합을 가진 스케일링된 번프 함수를 사용하여 $ v_{k-1}v_k $ 와 $ v_k^2 $ 의 곱에 대해 처리하고, 추정기 $ \widehat{\sigma}_n^2 $ 를 계산한다.
• 유한요소 이산화에 계층적 하이브리드 메쉬를 적용하여, 반암시적 행렬 조립과 다중 격자 유사한 메쉬 다듬힘을 통해 메모리 사용량을 줄이고 계산 속도를 향상시킨다.
• 알고리즘은 실시간 메쉬 다듬힘을 허용한다: 더 높은 정밀도가 필요할 경우 이전에 계산된 결과를 재사용하여 전체 재계산을 피한다.
• 이 방법은 줄리아로 구현되었으며 무료로 공개되어 있으며, 조각별로 일정한 계수와 다양한 탄성 대비 비율을 갖는 2차원 및 3차원 문제에서 수치적 검증이 이루어졌다.

실험 결과

연구 질문

• RQ1 짧은 척도 기여가 계산을 지배하는 다스케일 표현 방식을 통해 동질화 계수를 구성할 수 있는가? 이때 장기 척도 보정은 효율적으로 근사될 수 있는가?
• RQ2 계층적 하이브리드 메쉬는 랜덤 매질의 동질화 계수 유한요소 계산에서 메모리 사용량과 실행 시간을 어떻게 줄일 수 있는가?
• RQ3 추정기 $ \widehat{\sigma}_n^2 $ 의 수렴 속도는 척도 수 $ n $ 에 따라 어떻게 변화하는가? 이는 이론적 예측과 일치하는가?
• RQ4 해결책에 강한 특이성을 유도하는 모서리 불연속성이 존재하는 문제에서 이 방법은 어떻게 성능을 발휘하는가?
• RQ5 표준 랩탑에서 2차원 문제에서는 수초 내, 3차원 문제에서는 수분 내에 몇 퍼센트의 정밀도로 동질화 계수 추정을 달성할 수 있는가?

주요 결과

• 중간 정도의 탄성 대비 비율을 갖는 2차원 문제에서, 이 방법은 랩탑에서 10초 이내에 몇 퍼센트 이내의 오차로 동질화 계수를 계산한다.
• 3차원 문제에서는 동일한 하드웨어에서 5분 이내에 유사한 정밀도를 달성하며, 경험적 분산은 $ 2^{-3n} $ 으로 감소하여 이론적 예측인 $ \gamma = d = 3 $ 과 일치한다.
• 3차원에서 $ \alpha = 1, \beta = 9 $ 인 경우, $ \widehat{\sigma}^2(n=2, N_{\text{ref}}=4) $ 의 단일 샘플은 동질화 행렬의 95% 신뢰구간 $[3.62, 3.80]$ 내에 포함되며, 약 38분의 계산 시간이 소요된다.
• 이 방법은 2차원에서 약 $ h^{0.1337} $ 의 사전점근 수렴 속도를 보이며, 중간 정도의 메쉬 다듬힘으로도 오차를 급격히 감소시킨다.
• 추정기 $ \widehat{\sigma}_n^2 $ 는 탄성 대비 비율 증가에 대해 뛰어난 내구성을 보이며, 분산은 주로 메쉬 다듬음 수준 $ N_{\text{ref}} $ 와 독립적이고, 감쇠 속도는 $ C3^{-\gamma n} $ 으로 $ \gamma \approx 3.2 $ 로 나타나 이론적 값 $ \gamma = d $ 와 매우 유사하다.
• 이 방법은 실시간 메쉬 다듬힘 기능을 제공한다: 더 높은 정밀도가 필요할 경우 이전에 계산된 결과를 재사용하여 전체 재계산을 피할 수 있다.