[논문 리뷰] Uncertainty quantification for hyperbolic conservation laws with flux coefficients given by spatiotemporal random fields

Andrea Barth, Franz Fuchs|arXiv (Cornell University)|2015. 06. 25.

Probabilistic and Robust Engineering Design참고 문헌 30인용 수 3

한 줄 요약

이 논문은 유변법의 초순수 보존법칙에서 유속 계수를 시공간 랜덤 필드로 모델링한 경우의 불확실성 정량화를 위한 몬테카를로 유한체적 방법을 제안한다. 공간적 성분은 상관된 가우시안 랜덤 필드이며, 시간적 성분은 옴스타인-울렌벡 확률미분방정식(Ornstein–Uhlenbeck SDEs)으로 표현된다. 이 방법은 해의 통계적 모멘트(평균 및 분산)를 안정적으로 계산할 수 있게 하며, 자기유도와 선형 음향학에서 배경 속도 필드가 랜덤인 경우에 적용하여 수렴성과 복잡한 해 구조의 포착을 입증한다.

ABSTRACT

In this paper hyperbolic partial differential equations with random coefficients are discussed. We consider the challenging problem of flux functions with coefficients modeled by spatiotemporal random fields. Those fields are given by correlated Gaussian random fields in space and Ornstein-Uhlenbeck processes in time. The resulting system of equations consists of a stochastic differential equation for each random parameter coupled to the hyperbolic conservation law. We define an appropriate solution concept in his setting and analyze errors and convergence of discretization methods. A novel discretization framework, based on Monte Carlo Finite Volume methods, is presented for the robust computation of moments of solutions to those random hyperbolic partial differential equations. We showcase the approach on two examples which appear in applications: The magnetic induction equation and linear acoustics, both with a spatiotemporal random background velocity field.

연구 동기 및 목표

• 유속 계수가 시간 및 공간에 따라 변화하는 랜덤 필드로 모델링된 초순수 보존법칙에서의 불확실성 정량화 문제를 다루기.
• 랜덤 매개변수를 기술하는 초순수 PDE와 확률미분방정식(SDE)의 결합 시스템에 대한 해 개념을 개발하기.
• 고도로 정규화된 확률 변수가 필요하지 않은 비침습적이고 안정적인 수치 프레임워크를 제공하여 해의 통계적 모멘트(평균 및 분산)를 계산하기.
• 물리적으로 관련된 문제들인 자기유도 방정식과 랜덤 배경 속도 필드를 갖는 선형 음향학에 대해 방법을 검증하기.
• 이러한 확률적 초순수 시스템의 맥락에서 제안된 이산화 체계의 수렴성과 오차 행동을 확립하기.

제안 방법

• 유속 계수를 시공간 랜덤 필드로 모델링하며, 공간 성분은 상관된 가우시안 랜덤 필드로, 시간 성분은 옴스타인-울렌벡 SDE의 해로 표현한다.
• 결합 시스템을 랜덤 유속 함수를 갖는 초순수 보존법으로 기술하며, 각 계수는 시공간 노이즈에 의해 구동되는 SDE를 통해 진화한다.
• 몬테카를로 유한체적(MC FV) 프레임워크를 구현: 랜덤 필드의 실현값을 샘플링하고, 각 실현값에 대해 일阶 안정적인 업윈드 스킴을 사용하여 결정론적 초순수 PDE를 푼다.
• M개의 몬테카를로 샘플에 대한 앙상블 평균을 사용하여 해의 통계적 모멘트(평균 및 분산)를 근사한다.
• 속도 필드를 지배하는 옴스타인-울렌벡 과정의 시간 이산화에 밀스타인 스킴을 사용한다.
• 자기유도 방정식에 대해 주기적 경계 조건과 수렴이 없는 초기 조건을 적용하여 시뮬레이션의 물리적 일관성을 확보한다.

실험 결과

연구 질문

• RQ1 시공간 랜덤 필드에 의해 구동되는 초순수 보존법칙에 대해 일관된 해 개념을 어떻게 정의할 수 있는가?
• RQ2 이러한 확률적 초순수 시스템에 몬테카를로 유한체적 방법을 적용했을 때의 수렴 행동은 어떠한가?
• RQ3 물리적 응용에서 랜덤 배경 속도 필드 하에서 해의 통계적 모멘트(평균 및 분산)는 어떻게 변화하는가?
• RQ4 제안된 방법이 확률적 매개변수에 대해 해의 정규성 부족을 효과적으로 다룰 수 있는가? 이는 gPC 및 콘테이션 메서드의 제약을 피하는가?
• RQ5 자기유도 및 음향학과 같은 시스템에서 스트로스틱 운반에 의해 해의 분산에 어떤 구조적 특징이 나타나는가?

주요 결과

• 제안된 몬테카를로 유한체적 방법은 첫 번째 및 두 번째 통계적 모멘트 모두에서 기대되는 수렴 속도를 달성하며, 두 번째 모멘트는 더 큰 오차 상수를 보인다.
• 해의 평균 구조는 단일 결정론적 시뮬레이션과 매우 유사하여, 앙상블 평균이 주요 파동 전파 행동을 잘 포착하고 있음을 시사한다.
• 해의 분산은 복잡하고 비자명한 공간적 구조를 보이며, 이는 랜덤 속도 필드가 파동 왜곡 및 운반에 미치는 영향을 반영한다.
• t=0.75에서 256×256 격자에서 자기유도 방정식의 수치 결과는 랜덤 속도 필드로 인한 초기 자기장의 명확한 운반 및 왜곡을 보여준다.
• 이 방법은 강건하고 비침습적이며, 일반 다항식 혼합(gPC)과 같은 침투적 방법의 높은 복잡성에서 벗어나며, 확률적 매개변수의 정규성 저하가 있는 문제에 적합하다.
• 이 프레임워크는 3차원 문제 및 더 복잡한 SDE 또는 SPDE를 갖는 랜덤 필드로의 확장이 가능하며, 향후 효율성을 높이기 위해 GPU 가속 다중격자 솔버를 구현할 계획이다.