[논문 리뷰] Approximation of probability density functions for SPDEs using truncated series expansions

Giacomo Capodaglio, Max Gunzburger외 1명|arXiv (Cornell University)|2018. 10. 02.

Probabilistic and Robust Engineering Design인용 수 3

한 줄 요약

이 논문은 스토하스틱 편미분방정식(SPDEs)의 함수적 요소들의 확률밀도함수(PDF)를 근사하기 위해 잘라낸 그람-샤를리에르 및 에지워스 급수 전개를 제안한다. 몬테카를로 및 스토하스틱 갈레르킨 방법을 통해 계산된 모멘트로부터 유도된 세 moments를 활용하여, 다항식 전개를 통해 비정규 분포를 효율적으로 모델링하는 방법으로, 샘플링 기반 밀도 추정에 비해 계산 비용이 저렴한 대안을 제공한다.

ABSTRACT

The probability density function (PDF) of a random variable associated with the solution of a stochastic partial differential equation (SPDE) is approximated using a truncated series expansion. The SPDE is solved using two stochastic finite element (SFEM) methods, Monte Carlo sampling and the stochastic Galerkin method with global polynomials. The random variable is a functional of the solution of the SPDE, such as the average over the physical domain. The truncated series are obtained considering a finite number of terms in the Gram-Charlier or Edgeworth series expansions. These expansions approximate the PDF of a random variable in terms of another PDF, and involve coefficients that are functions of the known cumulants of the random variable. To the best of our knowledge, their use in the framework of SPDEs has not yet been explored.

연구 동기 및 목표

• 스토하스틱 편미분방정식(SPDEs)의 해에서 유도된 기능적 요소들의 확률밀도함수(PDF)를 정확하게 근사하는 데 도전하는 것. 이러한 PDF는 종종 비정규성을 나타낸다.
• 고전적인 급수 전개인 그람-샤를리에르 및 에지워스 급수의 SPDEs에 대한 적용을 탐색하는 것. 이 분야에서는 이론적 잠재력은 있지만 아직 탐색이 부족한 분야이다.
• 샘플링 기반 방법과 비교하여, 잘라낸 급수 전개가 SPDE 해의 PDF의 형태와 꼬리 행동을 얼마나 잘 포착하는지 평가하는 것.
• 몬테카를로 샘플링 및 전역 다항식을 사용한 스토하스틱 갈레르킨 방법을 포함한 두 가지 스토하스틱 유한요소 방법에서 유도된 결과를 기반으로 한 모멘트 기반 근사의 효과를 평가하는 것.
• 특히 샘플링이 매우 비용이 많이 드는 경우에, 몬테카를로 샘플링에 비해 계산 비용이 저렴한 SPDEs의 PDF 추정을 위한 대안을 제공하는 것.

제안 방법

• 이 방법은 SPDE 해의 기능적 요소(예: 도메인 전체에 대한 공간 평균)의 확률밀도함수(PDF)를 잘라낸 그람-샤를리에르 및 에지워스 급수 전개를 통해 근사한다.
• 확률변수의 모멘트는 몬테카를로 샘플링 및 전역 다항식 카오스 전개를 사용한 스토하스틱 갈레르킨 방법을 통해 SPDE 해로부터 계산된다.
• 급수 전개는 알려진 기준 PDF(일반적으로 정규분포)를 기반으로 하며, 확률변수의 모멘트를 기반으로 한 수정 항을 포함한다.
• 급수의 항 수를 유한한 수로 잘라내어 정확도와 계산 비용의 균형을 맞추며, 고차항은 비대칭성과 첨도 효과를 포착한다.
• 결과로 도출된 근사 PDF는 몬테카를로 샘플링을 통해 확보된 기준 PDF와 비교하여 정확도를 평가한다.
• 진짜 해의 PDF가 알려져 있거나 잘 근사되는 테스트 케이스에서 이 방법을 검증함으로써 근사 오차의 정량적 평가가 가능해진다.

실험 결과

연구 질문

• RQ1 잘라낸 그람-샤를리에르 및 에지워스 급수는 비정규적 특성을 띠는 SPDE 해의 기능적 요소의 PDF를 효과적으로 근사할 수 있는가?
• RQ2 이 급수 근사의 정확도와 수렴성은 몬테카를로 샘플링에 비해 PDF의 형태와 꼬리 행동을 얼마나 잘 포착하는가?
• RQ3 몬테카를로 샘플링과 스토하스틱 갈레르킨 방법에서 유도된 모멘트를 사용할 경우, PDF 근사의 품질에 어떤 영향을 미치는가?
• RQ4 급수 전개의 고차항은 어떤 상황에서 근사 품질을 크게 향상시키며, 언제는 저차항 근사로도 충분한가?
• RQ5 이 급수 전개는 SPDEs의 PDF 추정에서 전체 샘플링에 비해 얼마나 효율적인 대안이 될 수 있는가?

주요 결과

• 잘라낸 그람-샤를리에르 및 에지워스 급수 전개는 광범위한 샘플링 없이도 비정규 PDF를 효율적으로 근사할 수 있는 실용적이고 계산 비용이 저렴한 방법을 제공한다.
• 사용된 모멘트의 수가 많을수록 근사 품질이 향상되며, 특히 비정규 분포에서 중요한 비대칭성과 첨도를 잘 포착한다.
• 동일한 계산 비용에서 스토하스틱 갈레르킨 방법은 몬테카를로 샘플링보다 더 정확한 모멘트 추정치를 제공하여 더 나은 PDF 근사 결과를 이끈다.
• 이 방법은 높은 밀도 영역의 PDF 형태를 효과적으로 포착하며, 꼬리 영역에서도 합리적인 성능을 보이지만, 매우 비대칭적이거나 꼬리가 두꺼운 분포에서는 정확도가 떨어진다.
• 특히 고해상도 PDF 추정이 필요한 경우, 몬테카를로 샘플링에 비해 뚜렷한 계산적 절감 효과를 보인다.
• 스토하스틱 갈레르킨 프레임워크 내에서의 전역 다항식 카오스의 사용은 부드럽고 안정적인 모멘트 계산을 가능하게 하며, 이는 신뢰할 수 있는 급수 전개 계수를 확보하는 데 필수적이다.