[論文レビュー] Gaussian process modeling for stochastic multi-fidelity simulators, with application to fire safety

Rémi Stroh, Julien Bect+3 more|arXiv (Cornell University)|May 9, 2016

Probabilistic and Robust Engineering Design被引用数 4

ひとこと要約

本論文は、メッシュ分解能の異なる火炎動態シミュレーションを対象として、ベイジアン多忠実度ガウス過程モデリングを確率的シミュレータに拡張する。平均と分散の両方をガウス過程としてモデル化し、非定常相関構造を組み込むことで、低忠実度シミュレーションを用いて高忠実度の火炎挙動を正確に予測可能となり、計算コストを著しく削減しつつも高忠実度モデルと同等の予測性能を達成する。

ABSTRACT

To assess the possibility of evacuating a building in case of a fire, a standard method consists in simulating the propagation of fire, using finite difference methods and takes into account the random behavior of the fire, so that the result of a simulation is non-deterministic. The mesh fineness tunes the quality of the numerical model, and its computational cost. Depending on the mesh fineness, one simulation can last anywhere from a few minutes to several weeks. In this article, we focus on predicting the behavior of the fire simulator at fine meshes, using cheaper results, at coarser meshes. In the literature of the design and analysis of computer experiments, such a problem is referred to as multi-fidelity prediction. Our contribution is to extend to the case of stochastic simulators the Bayesian multi-fidelity model proposed by Picheny and Ginsbourger (2013) and Tuo et al. (2014).

研究の動機と目的

• 高忠実度火炎シミュレーションの高い計算コストを、安価な低忠実度シミュレーションを活用することで低減すること。
• 元々決定論的な多忠実度ガウス過程モデルを、固有のばらつきを有する確率的シミュレータに拡張すること。
• 粗いメッシュからのデータを用いて、細かいメッシュ解像度における火炎挙動の正確な予測を可能にすること。
• 火災安全評価における臨界温度閾値を超える確率を推定すること。
• 計算負荷を軽減しつつも、予測精度を維持する計算効率の良いサーモンモデルの開発

提案手法

• 出力のシミュレータを、平均関数と分散関数を持つ確率過程としてモデル化し、出力がガウス分布に従うと仮定する。
• 理想シミュレータと、メッシュサイズが0に近づくにつれて減少する決定論的誤差項に平均関数を分解する。
• 入力とメッシュサイズに依存する分離可能な共分散構造を持つ非定常ガウス過程事前分布を平均関数に実装する。
• 分散をメッシュサイズのみの関数としてモデル化し、メッシュが細かくなるにつれて減少すると仮定する。
• 平均と分散に共同ガウス過程事前分布を適用し、階層的事前分布を用いた完全ベイジアン推論を可能にする。
• 低忠実度データが事前分布を決定し、高忠実度データが事後予測分布を精緻化する多忠実度フレームワークを採用する。

実験結果

リサーチクエスチョン

• RQ1 固有のばらつきを有する確率的シミュレータを扱えるように、ベイジアン多忠実度ガウス過程モデルを拡張できるか？
• RQ2 低忠実度モデルは、平均と分散の両面で高忠実度火炎シミュレーションの結果をどの程度正確に予測できるか？
• RQ3 提案手法の予測精度は、臨界温度閾値の推定において高忠実度モデルと比較してどの程度か？
• RQ4 故障確率推定におけるモデルの不確実性は、高忠実度基準モデルと比較してどの程度か？
• RQ5 限定的な計算予算のもとで、提案モデルは信頼できる安全評価を達成できるか？

主な発見

• 提案された多忠実度非定常モデルは、20 cmメッシュ解像度における温度予測において、高忠実度モデルと同等の予測性能を達成した。
• 多忠実度定常モデルも良好なフィットを示したが、わずかに高温出力を低く見積もっていた。
• 正規化残差の事後標準偏差は1に近く、予測分散が適切にキャリブレーションされていることを示している。
• 高忠実度モデルは60°Cを超える確率の事後分布が最も狭いが、多忠実度モデルの結果も基準と整合的であった。
• 多忠実度定常モデルは小さな事後不確実性を達成したが、そのサポートは高忠実度モデルからわずかに逸脱していた。
• 提案モデルは、正確性と計算効率のバランスを図り、シミュレーションコストを削減しつつ信頼できる安全評価を可能にした。