[論文レビュー] Equation-free optimal switching policies for bistable reacting systems using coarse time-steppers

Antonios Armaou, Ioannis G. Kevrekidis|ArXiv.org|Oct 21, 2004

Numerical methods for differential equations参考文献 36被引用数 6

ひとこと要約

本稿では、確率的力学に基づく二安定反応系において、安定状態間の最適スイッチング方策を同定するための、方程式を必要としない計算フレームワークを提示する。反応系のマクロスケールの運動方程式が不明な状況においても、反応系の期待状態のマクロスケールの時間発展を短時間のマイクロスケール（反応キネティクスモンテカルロ）シミュレーションに基づいて推定する粗大タイムステッパーを用いることで、Hooke-Jeeves法などの導出フリー最適化手法と統合し、明示的なマクロスケール方程式を必要とせずに、共存する定常状態間を最小パラメータでスイッチングする軌道を同定可能である。100×100格子系においては、4,000秒未満で収束を達成した。

ABSTRACT

We present a computer-assisted approach to locating approximate coarse optimal switching policies between stationary states of chemically reacting systems described by microscopic/stochastic evolution rules. The ``coarse time-stepper" constitutes a bridge between the underlying kinetic Monte Carlo simulation and traditional, continuum numerical optimization techniques formulated in discrete time. The approach is illustrated through two simple catalytic surface reaction models, implemented through kinetic Monte Carlo: NO reduction on Pt, and CO oxidation on Pt. The objective sought in both cases is to switch between two coexisting stable stationary states by minimal manipulation of a macroscopic system parameter.

研究の動機と目的

• マクロスケール方程式が得られない化学反応系において、共存する安定定常状態間の最適スイッチング方策を同定する計算手法の開発。
• 反応キネティクスモンテカルロ（KMC）シミュレーションによるマイクロスケールのシミュレーションと、マクロスケール最適化手法とを、粗大タイムステッパーを介して橋渡しする。
• 複雑で非線形な挙動を示し、複数のアトラクタを持つ系に対して、時間変化する制御方策の動的最適化を可能にする。
• 粗大スケールの観測量を用いて、導出フリー最適化アルゴリズム（例：Hooke-Jeeves）を適用する可能性を実証する。
• 粗い時間グリッドからの解を再利用するマルチグリッド的リファインメント戦略を用いて、計算コストを低減する。

提案手法

• 粗大タイムステッパーを用いて、短時間のマイクロスケール（反応キネティクスモンテカルロ）シミュレーションの結果に基づき、系の期待状態のマクロスケール時間発展をシミュレートする。
• 最適制御問題を離散時間で定式化し、時間変化する制御パラメータ（例：温度や表面被覆率依存の反応定数）を最適化すべき意思決定変数として扱う。
• 大域的変化や長時間のスイッチングをペナルティとするコスト関数を最小化する目的で、導出フリー最適化アルゴリズム（Hooke-Jeeves）を適用する。
• マルチグリッド的リファインメントプロセスにおいて、最適パラメータプロファイルを粗い時間グリッドから細かい時間グリッドに三次スプライン補間を用いて移行する。
• 最適化プロセス中、系の応答を候補となるパラメータ軌道に対して評価するために、反応系の運動方程式を明示的に必要とせずに、繰り返し粗大タイムステッパーを呼び出す。
• 本手法は、二酸化窒素還元反応とプラチナ表面における一酸化炭素酸化反応の二つの触媒表面反応モデルを用いて検証された。両者とも二安定性を示す。

実験結果

リサーチクエスチョン

• RQ1 マクロスケールの時間発展方程式が不明な状況において、二安定反応系の共存する安定定常状態間の最適スイッチング方策を計算可能か？
• RQ2 粗大タイムステッパーは、確率的マイクロスケールシミュレーションとマクロスケール最適化アルゴリズムを効果的に橋渡しできるか？
• RQ3 導出フリー最適化と組み合わせたマルチグリッドリファインメントの計算効率は、このようなスイッチング問題に対してどの程度高いか？
• RQ4 反応キネティクスモンテカルロシミュレーションを用いた実際の触媒反応系に適用した場合、収束性と精度の観点で本手法はどの程度の性能を示すか？
• RQ5 本手法は、高次元の粗大スケール状態空間や空間的に分布するダイナミクスを持つ系へ拡張可能か？

主な発見

• 本手法は、温度や被覆率依存反応定数といった単一のマクロスケールパラメータを操作することで、一つの安定定常状態から別の定常状態へ遷移させる最適スイッチング方策を効果的に同定した。
• 200回のKMC実現を用いた100×100格子系のCO酸化反応系において、Hooke-Jeeves法は時間ステップT=0.1で3,983秒で最適方策に収束した。
• 最適軌道は不安定なサドル状態の分離子を通過しており、ポテンシャルランドスケープ内での制御された遷移が、臨界点を通じて実現されていることを確認した。
• 三次スプライン補間を用いたマルチグリッドリファインメント戦略は、繰り返しの時間グリッドのリファインメントにおいて、収束速度と解の精度を顕著に向上させた。
• 粗大タイムステッパーを用いることで、複雑で非線形なダイナミクスを示す系においても、閉形式のマクロスケール方程式を導出することなく、効果的な最適化が可能であった。
• 本手法は、マイクロスケールシミュレーションから抽出可能なマクロスケール観測量を持つ他の系へ一般化可能であり、希少事象のダイナミクスや空間的に分布する系に対しても適用可能である。