[논문 리뷰] Multilevel Monte Carlo methods for highly heterogeneous media

Aretha L. Teckentrup|arXiv (Cornell University)|2012. 06. 07.

Probabilistic and Robust Engineering Design참고 문헌 17인용 수 5

한 줄 요약

이 논문은 다층 몬테카를로(MLMC) 방법을 고도로 이방성이고 비균일하게 강성 조건을 갖는 무작위 계수를 가진 타원형 PDE에 확장한다—예를 들어 지하류 모델링에서 흔한 로그노멀 투과도 필드와 같은 경우이다. 무작위 필드의 정규성에 대한 최소한의 가정 하에, 텐서값 계수와 압력 및 다리 플럭스의 점 평가에 대해 수렴성과 계산 비용의 경계를 증명한다.

ABSTRACT

We discuss the application of multilevel Monte Carlo methods to elliptic partial differential equations with random coefficients. Such problems arise, for example, in uncertainty quantification in subsurface flow modeling. We give a brief review of recent advances in the numerical analysis of the multilevel algorithm under minimal assumptions on the random coefficient, and extend the analysis to cover also tensor--valued coefficients, as well as point evaluations. Our analysis includes as an example log--normal random coefficients, which are frequently used in applications.

연구 동기 및 목표

• 지하류에서 고도로 이방성 투과도를 갖는 경우에 표준 몬테카를로 방법의 높은 계산 비용 문제를 해결한다.
• 표준 몬테카를로 이론 분석을 비균일하게 강성 조건이나 유계성이 없는 경우, 예를 들어 로그노멀 필드와 같은 경우로 확장한다.
• 지하수 흐름에서 이방성 매체를 모델링하는 데 관련된 텐서값 계수를 다룰 수 있도록 MLMC 프레임워크를 일반화한다.
• 압력 또는 다리 플럭스의 점별 평가를 위한 MLMC 추정기의 수렴성과 계산 비용 경계를 확립한다.
• 해의 경로별 정규성과 유한요소 근사의 모멘트 경계를 포함한 최소한의 가정 하에 엄밀한 오차 분석을 제공한다.

제안 방법

• 지하류에서 투과도를 모델링하는 데 사용되는 무작위 텐서값 계수 $\mathbf{A}$ 를 가진 타원형 PDE $-\mathrm{div}(\mathbf{A}\nabla u) = f$ 에 다층 몬테카를로 알고리즘을 적용한다.
• 메쉬 크기가 $h_\ell$ 인 중첩된 메쉬 계층에서 해 $u$ 를 유한요소 방법으로 근사하고, 기능 $M$ 에 대해 추정기 $Q_{h_\ell} = M(u_{h_\ell})$ 를 정의한다.
• 예상값 $\mathbb{E}[Q]$ 를 수준 간 차이 $\mathbb{E}[Q_{h_\ell} - Q_{h_{\ell-1}}]$ 를 사용한 텔레스코프 합으로 분해하여 분산을 감소시킨다.
• 해의 해 $u(\omega,\cdot)$ 에 대한 소볼레프 임베딩과 샤투르 유형 정규성 추정을 사용하여 편향 $\mathbb{E}[|Q - Q_{h_\ell}|]$ 의 수렴 속도를 확립한다.
• 무작위 필드의 호일더 정규성과 적분 가능성에 대한 가정을 사용하여 해와 그 유한요소 오차의 모멘트 경계를 유도한다.
• 최소한의 가정 하에 루트 평균 제곱 오차 $\varepsilon$ 를 달성하기 위한 계산 비용이 $\mathcal{O}(\varepsilon^{-2})$ 이내로 경계지어지며, 이는 표준 몬테카를로의 최적 비율을 유지하지만 상수가 훨씬 작아진다.

실험 결과

연구 질문

• RQ1 비균일하게 강성 조건을 갖거나 유계성이 없는 경우, 다층 몬테카를로 방법이 타원형 PDE에 대해 엄밀히 분석 가능할 수 있는가?
• RQ2 이방성 매체에서 사용되는 텐서값 무작위 투과도 필드에 대해 MLMC 방법이 여전히 최적의 복잡도를 유지하는가?
• RQ3 전역 기능이 아닌 압력 또는 다리 플럭스의 점 평가에 대해 수렴성과 계산 비용 경계를 확립할 수 있는가?
• RQ4 해가 거친 무작위 계수로 인해 낮은 공간 정규성을 갖는 경우, MLMC 추정기의 수렴 속도는 어떻게 되는가?
• RQ5 경로별 정규성 추정(예: 샤투르 이론에서의 것들)은 MLMC 프레임워크에서 편향과 분산 감소에 어떻게 기여하는가?

주요 결과

• 최소한의 무작위 계수 가정 하에 다층 몬테카를로 추정기는 루트 평균 제곱 오차 $\varepsilon$ 를 달성하기 위해 계산 비용이 $\mathcal{O}(\varepsilon^{-2})$ 이내로 경계지어지며, 이는 표준 몬테카를로의 최적 비율을 유지하지만 상수가 훨씬 작다.
• 압력 또는 다리 플럭스의 점 평가에 대해 편향은 $\alpha = 1 + s - d/2$ (여기서 $s < 1$) 또는 샤투르 이론에서의 $\alpha = 1 + t$ 의 속도로 감소하고, 분산은 $\beta = 2\alpha$ 의 속도로 감소한다.
• 분석은 유계성이 없는 로그노멀 무작위 필드에 대해서도 적용 가능하며, 오직 리프시츠 연속성 공분산만 요구한다. 이는 이전 결과에서 요구한 균일한 경계 조건을 필요로 하지 않는다.
• 유한요소 오차의 수렴 속도는 $\|u - u_h\|_{L^p(\Omega, W^{1,\infty}(D))} \lesssim h^{s - d/2}$ (여기서 $s < 1$) 이며, 샤투르 정규성 조건 하에 향상된 경로별 속도 $h^t$ 도 가능하다.
• 이 방법은 균일한 강성 조건이나 유계성이 없더라도 텐서값 계수, 예를 들어 이방성 투과도에 대해 적용 가능하다.
• 이론적 프레임워크는 투과도 필드가 고도로 공간적 변동성을 보이고 정규성이 낮은 실질적인 지하류 시뮬레이션에서 MLMC의 사용을 지원한다.