thumb|物理情報ニューラルネットワークを用いてナビエ–ストークス方程式を扱う概念図

物理情報ニューラルネットワーク(ぶつりじょうほうニューラルネットワーク、、PINN)は、既知の微分方程式や初期条件、境界条件などの物理的・数学的制約をニューラルネットワークの学習過程に組み込む手法である。主として偏微分方程式や常微分方程式で記述される系の数値解析に用いられ、既知のモデルから解を求める順問題だけでなく、観測データから未知の物理パラメータなどを推定する逆問題にも利用される。

典型的なPINNでは、空間座標や時間を入力し、速度、圧力、温度などの未知の物理量を出力するニューラルネットワークを構築する。ネットワークの出力を自動微分することで支配方程式の残差を求め、方程式の残差、初期・境界条件との不一致、必要に応じて観測データとの不一致を組み合わせた損失関数を最小化する。このため、正解となる数値解を大量の教師データとして用意しなくても、既知の方程式自体を学習時の制約として利用できる。

日本語では「物理情報ニューラルネットワーク」の表記が学術文献で用いられている。また、「Physics-Informed Neural Network」の語をそのまま用いる例や、「物理法則を組み込んだニューラルネットワーク」と説明する例もある。

歴史 ニューラルネットワークを微分方程式の数値解法に用いる研究は、PINNという名称が定着する以前から存在した。1994年、M・W・M・G・ディサナヤケとN・ファン=ティエンは、ニューラルネットワークによる普遍近似器と選点法を組み合わせ、偏微分方程式を解く問題を最小化問題として扱う手法を報告した。

1998年には、I・E・ラガリス、A・リカス、D・I・フォティアディスが、初期条件や境界条件を構造的に満たす試行関数の一部としてニューラルネットワークを用い、常微分方程式および偏微分方程式を解く方法を提案した。

2019年、Maziar Raissi、Paris Perdikaris、George Em Karniadakisは、非線形偏微分方程式によって表現される物理法則をニューラルネットワークの学習に組み込み、順問題と逆問題を共通の枠組みで扱う「physics-informed neural networks」を体系化した。以降、PINNは物理モデルと機械学習を組み合わせる科学機械学習(scientific machine learning)の代表的手法の一つとして研究されている。

定式化 PINNが対象とする微分方程式を一般的に

とする。ここで、は求める物理量、は微分演算子、は方程式に含まれる物理パラメータ、は空間領域である。さらに、初期条件または境界条件を

と表す。

PINNでは、未知関数を、パラメータをもつニューラルネットワーク

によって近似する。このとき支配方程式の残差は

となる。の空間・時間微分は一般に自動微分によって計算される。

代表的な損失関数は

のように表される。は領域内部の選点(コロケーション点)における微分方程式の残差、は初期条件・境界条件との不一致、は観測データとの不一致を表し、、、はそれぞれの項の重みである。観測データを用いない順問題では、データ項を含めず、方程式と初期・境界条件のみからネットワークを学習させることもできる。

順問題 順問題では、支配方程式に含まれる物理パラメータや必要な初期・境界条件を既知として、未知の解を求める。従来の有限差分法、有限要素法、有限体積法などとは異なり、典型的なPINNでは要素間の接続関係を持つ計算格子を必須とせず、領域内に配置した選点で方程式残差を評価する。この意味でPINNはしばしばメッシュフリーの手法と説明される。

学習後のネットワークは座標から物理量への連続的な近似関数として利用できるため、学習時の選点とは異なる位置でも値を評価できる。ただし、計算領域、支配方程式、物理パラメータなどを変更した場合に、通常のPINNがそのまま利用できるとは限らず、再学習が必要となる場合がある。

また実運用面では、PINNを用いたソルバーはTensorFlowなどのテンソル演算フレームワーク上で動作するため、GPUの活用が容易であったり一般的なソルバーより高速に動作することが期待されるというメリットが存在する。

逆問題 逆問題では、観測された物理量を利用して、状態だけでなく、支配方程式に含まれる未知のパラメータなどもニューラルネットワークのパラメータと同時に推定する。これにより、疎な観測データから物理場を再構築したり、未知の係数や物理パラメータを推定したりすることができる。

例えば、菊本英紀らは、PINNを用いて航空機客室の縮尺模型から得られた疎な平均流速データと流体の支配方程式を組み合わせ、換気流れ場を再構築する方法を検討している。またプラズマ中に現れるソリトンに対してPINN解析を行い、それがKdV方程式に従うことを定量評価した研究も存在する。

ただし、観測データを利用できる逆問題で境界情報の不足を補える場合があることは、一般の順問題が必要な初期条件や境界条件なしに一意に解けることを意味しない。問題の適切性は支配方程式や与えられる条件に依存する。

また逆問題は一般的に順問題より難しく、支配方程式の係数を精度よく推論するためにはより適切な数と質を持つデータセットが必要となる。

特徴と応用 PINNは、既知の物理モデルと観測データを同一の最適化問題に組み込める点に特徴がある。特に、観測データが疎またはノイズを含む場合の状態推定や逆問題、通常の数値計算から得られたデータを利用する代理モデルの構築などが主要な利用分野となっている。

応用は流体力学、熱輸送、固体力学、材料科学などに広がっている。日本でも換気流れ場の再構築や、プラズマ物理における順問題・逆問題への応用が紹介されている。

PINNには多数の拡張手法が提案されている。例えば、計算領域を複数の部分領域に分割して複数のニューラルネットワークを学習させるcPINN(conservative PINN)やXPINN(extended PINN)、弱形式・変分形式を利用する手法、初期・境界条件をニューラルネットワークの構造によって厳密に満たす手法、不確実性を扱うベイズ型PINNなどが研究されている。

限界 PINNは従来の数値解法を一般的に置き換える万能な方法ではない。ニューラルネットワークの学習そのものが非線形な最適化問題であり、支配方程式、境界条件、観測データなど複数の損失項を同時に小さくする必要があるため、学習が困難になる場合がある。損失項の大きさや収束速度が異なることによって勾配が不均衡となり、特定の制約だけが優先される現象も報告されている。

標準的なPINNは、移流、反応、拡散を含む一部の方程式で、問題が複雑になると適切な解を学習できなくなる場合がある。例えば、KdV方程式の1ソリトン解が一様な背景のもとで移流している様子を逆問題解析させたとする。1ソリトン解は「初期値がソリトン形の移流方程式の解」ともみなすことができるため、PINNによる逆問題解析で支配方程式が移流方程式かKdV方程式かを見極めることは原理的に不可能である。Krishnapriyanらは、こうした失敗がニューラルネットワークの表現能力だけによるものではなく、微分方程式をソフトな制約として加えた結果、損失関数の最適化が困難になることにも起因すると報告している。

また、急峻な勾配、不連続性、複数の空間・時間スケールを含む問題などでは、標準的なPINNの学習の難易度が上がり、必要とされるデータ点が大きくなることが知られている。従来の数値解法と比べて学習に長い計算時間を要する場合もあり、特に高次元問題や三次元問題では計算負荷が問題となる。

PINNの近似誤差や汎化誤差については理論解析も進められているが、収束性は問題の種類や解の正則性、安定性、最適化精度などの条件に依存する。De RyckとMishraによる数値解析のレビューでは、近似誤差・汎化誤差・訓練誤差を分けて解析したうえで、特に訓練誤差がPINNの性能を制限する主要な要因の一つであると指摘されている。

応用面で、気象予想では2026年現在はPINNによる予報は第一原理計算による予報を精度面で上回っていないとする報告もある。

関連項目 ニューラルネットワーク ディープラーニング 偏微分方程式 数値解析 逆問題 自動微分 有限要素法 計算流体力学

脚注

外部リンク PINNs - RaissiらによるPINNの公開実装