MMS（Method of Manufactured Solutions）とは何ですか？

MMS（Method of Manufactured Solutions）は、数値解析コードやソルバーの検証（Verification）に用いられる手法です。あらかじめ解析的に微分・積分可能な「製造解」（例：温度分布 T(x,y)=sin(2πx)・cos(2πy)+5）を設定し、これを支配方程式に代入して対応するソース項（例：熱源Q）を逆算します。この製造解とソース項をソルバーに入力し、得られた数値解が製造解にどれだけ一致するかを評価することで、コードの実装誤差を検証します。

MMSはどのような問題の検証に使えますか？

MMSは、熱伝導、流体力学、構造解析など、支配方程式が定式化できる幅広いCAE（数値シミュレーション）問題のコード検証に適用できます。特に、非線形性や移流項が支配的な現象（例：乱流モデル）の検証にも有効です。使用する製造解は、検証対象の物理的性質（非線形性、特異点、境界層など）を反映している必要があり、これによりソルバーや離散化手法の特定の挙動を系統的にテストできます。

MMSで使用する「製造解」の選び方の基準は？

MMSで用いる製造解の選択には、主に2つの重要な基準があります。第一に、解は解析的に（何度でも）微分および積分可能であること。これは、支配方程式に代入してソース項を正確に計算するために必須です。第二に、解が検証対象の問題の物理的性質（例えば、非線形性、移流項の支配性、特異挙動など）を適切に反映していることです。これにより、コードが実問題に対して正しく動作するかをより現実的に検証できます。

MMSの具体的な実施手順を教えてください。

MMSの実施手順は以下の通りです。1. 検証対象の支配方程式（例：定常熱伝導方程式 ∇・(k∇T)+Q=0）を設定する。2. 解析的に扱える「製造解」（例：T(x,y)=sin(2πx)・cos(2πy)+5）を定義する。3. この製造解を支配方程式に代入し、必要なソース項（この例では熱源 Q = -∇・(k∇T) ）を解析的に計算する。4. 計算されたソース項と製造解（境界条件としても使用）を検証対象の数値ソルバーに入力してシミュレーションを実行する。5. 得られた数値解と元の製造解を比較し、誤差を定量評価してコードの妥当性を検証する。

MMS — CAE用語解説

カテゴリ: 用語集 | 2026-01-15

CAE visualization for mms - technical simulation diagram

MMS

🧑‍🎓

先生、MMS（Method of Manufactured Solutions）って何ですか？

理論と物理

MMSの基本概念

🧑‍🎓

MMSって、具体的に何をする手法なんですか？「製造された解」という言葉がよく出てきますが、どういう意味ですか？

🎓

良い質問だ。MMS（Method of Manufactured Solutions）は、数値解法やコードの検証（Verification）に使う手法だ。まず、自分で「解」を決めてしまう。例えば、2次元熱伝導の温度分布を

$$ T(x, y) = \sin(2\pi x) \cdot \cos(2\pi y) + 5 $$

と「製造」する。この解は、解析的に微分できるので、支配方程式に代入して、対応する「ソース項」を逆算するんだ。

🧑‍🎓

ソース項を逆算？普通はソース項から解を求めるんじゃないですか？

🎓

その通り。通常の順問題と逆の手順を踏むのがMMSの肝だ。定常熱伝導方程式

$$ \nabla \cdot (k \nabla T) + Q = 0 $$

があったとする。ここに先ほどの製造解Tを代入すると、

$$ Q = -\nabla \cdot (k \nabla T) $$

というソース項Qが計算できる。このQと製造解Tのセットを、ソルバーに入力して計算させ、出てきた数値解が製造解Tにどれだけ近いかで、コードの正しさを検証する。

🧑‍🎓

なるほど。でも、そんな都合のいい解を考えるのは難しそうですね。解を選ぶ時の基準はありますか？

🎓

いくつか重要な基準がある。第一に、解は解析的に微分・積分可能であること。第二に、解が問題の物理的性質（例えば、非線形性、移流項の支配性など）を反映していること。例えば、乱流モデルの検証には、速度勾配が大きい解を選ぶ。第三に、境界条件を厳密に満たすこと。例えば、ディリクレ条件

$$ T|_{\Gamma} = g $$

を課すなら、製造解Tが境界Γ上で正確にgになるように設計する。

数値解法と実装

離散化と誤差評価

🧑‍🎓

MMSを使って誤差を評価する時、具体的に何を見るんですか？単に数値解と製造解の差を見るだけ？

🎓

差を見るのは第一歩だが、それだけでは不十分だ。重要なのは「収束率」を確認すること。メッシュサイズhを系統的に細かくしていき（例えば、1mm, 0.5mm, 0.25mm）、各メッシュでの誤差ノルム（L2ノルムやH1ノルム）を計算する。理論的な収束率（例えば、線形要素ならL2誤差が

$$ O(h^2) $$

）と一致するかどうかを対数グラフで確認する。これが一致しなければ、コードの実装にバグがある可能性が高い。

🧑‍🎓

誤差ノルムの計算では、製造解をメッシュの節点で評価するだけでいいんですか？

🎓

節点値の差（ポイントワイズ誤差）も見るが、より厳密には、数値解と製造解の関数としての差を積分で評価する。例えば、L2ノルムは

$$ ||e||_{L^2} = \sqrt{\int_{\Omega} (u_{num} - u_{exact})^2 d\Omega } $$

だ。この積分は、数値求積法（ガウス積分など）を使って各要素ごとに計算する。商用ソフトのAbaqusやANSYSの検証マニュアルでも、同様の手順が記載されている。

🧑‍🎓

複雑な非線形問題、例えば弾塑性解析にもMMSは使えますか？降伏条件や硬化則が入ると、製造解を作るのが非常に難しそうですが。

🎓

可能だが、高度な技巧が必要だ。例えば、等方硬化則を持つ弾塑性材料を考える。まず、変位場

$$ u(x, y) $$

を製造する。そこからひずみ、応力を計算し、それが降伏面

$$ f(\sigma, \bar{\epsilon}^p) = 0 $$

を満たすように、あらかじめ塑性ひずみや硬化パラメータを調整して「製造」する。実は、NIST（米国標準技術研究所）が公開しているMMSベンチマーク問題集には、いくつかの非線形材料モデルの例が含まれている。

実践ガイド

ワークフローとチェックリスト

🧑‍🎓

自分でMMS検証をやってみたいです。最初に何から手を付けるべきですか？

🎓

まずは最も単純な線形定常問題から始めることを勧める。手順はこうだ：1) 1次元Poisson方程式を選ぶ。2) 解として3次多項式

$$ u(x) = x^3 + 2x^2 - x + 1 $$

を製造する。3) 境界条件（x=0でu=1, x=1でu=3）とソース項

$$ f(x) = -6x - 4 $$

を計算する。4) 自作FEMコードまたは商用ソフトのユーザーサブルーチン機能でこれを解かせる。5) 誤差がマシンイプシロン（~1e-15）レベルであることを確認する。

🧑‍🎓

商用ソフト、例えばANSYS MechanicalでMMSを実行するにはどうすればいいですか？ソース項を設定する機能がありますか？

🎓

直接「MMSモード」はないが、APDLやUser Defined Function (UDF) を使って実現できる。熱解析なら、`BFE`コマンドで要素ごとに体積発熱を定義できる。あるいは、より一般的には、`ANSYS UPFs` (User Programmable Features) を使って、方程式レベルでソース項を追加するカスタム要素を開発する。ただし、これは上級者向けだ。検証目的なら、オープンソースのFEMライブラリ（FEniCS, deal.II）で始める方が敷居は低い。

🧑‍🎓

MMSで検証が完了したとみなす基準はありますか？誤差が小さければそれでOK？

🎓

誤差の絶対値が小さいだけでは不十分だ。先述した「収束率」が理論値と統計的に有意に一致することが最も強力な証拠になる。例えば、メッシュを4段階細かくし、各段階の誤差をプロットしたとき、その傾きが期待される収束次数（例：2）と一致するか。さらに、境界条件や保存則（例えば、全熱流量のバランス）がメッシュの粗さに依存せず満たされているかも併せてチェックする。ASME V&V 20規格でも、コード検証にはこのような系統的な収束性試験が推奨されている。

ソフトウェア比較

各ソフトでの実装アプローチ

🧑‍🎓

COMSOL Multiphysicsは「方程式ベースのモデリング」が売りですが、MMSとは相性が良さそうですか？

🎓

非常に良い指摘だ。COMSOLはMMSと本質的に親和性が高い。GUIで「係数型PDE」インターフェースを選び、製造解を「解析関数」として定義する。その後、ソース項の欄に、その製造解を微分して得た式を直接入力できる。例えば、拡散方程式のソース項fを `-d(d(u_manufactured,x),x)-d(d(u_manufactured,y),y)` と定義するだけで済む。この手軽さは、Abaqusのユーザーサブルーチン（Fortran）やANSYSのAPDLと比べて圧倒的にハードルが低い。

🧑‍🎓

Abaqusでは、どのモジュールを使ってMMSを実施するのが現実的ですか？

🎓

標準的な熱伝導や構造解析モジュールでは、任意のソース項を定義するのは難しい。主な方法は二つだ。一つは、`DLOAD`や`DFLUX`などのユーザーサブルーチンを使って、要素ごとにソース項を定義する方法。もう一つは、`Abaqus/Standard`の`USER ELEMENT` (UEL) インターフェースを使って、要素行列と右辺ベクトルをスクラッチから作成する方法だ。後者は柔軟性が極めて高いが、FEMの深い理解が必須で、開発コストが大きい。Dassault Systèmesのカスタマーサポートは、UELを用いた検証事例をいくつか公開している。

🧑‍🎓

オープンソースの環境ではどうでしょうか？

🎓

オープンソースはMMSの本丸と言える。Pythonライブラリの`FEniCS`や`Firedrake`は、弱形式をそのままコードに落とし込めるので、製造解とソース項の定義が数学的表記に近い形で可能だ。例えば、変分形式の中で `f = -div(grad(u_manufactured))` と書ける。また、`MFEM`や`deal.II`のチュートリアルには、MMSを用いた収束率検証のサンプルコードが豊富に含まれている。研究開発の現場では、これらのツールを使って新しい解法アルゴリズムの検証を行うのが一般的だ。

トラブルシューティング

よくあるエラーと対策

🧑‍🎓

MMSをやってみたのですが、メッシュを細かくしても誤差がまったく減りません。考えられる原因は？

🎓

まず疑うべきは「境界条件の設定ミス」だ。製造解が、数値計算で課した境界条件と完全に一致しているか確認してほしい。例えば、ノイマン条件（勾配指定）で数値の符号を間違えていないか。第二に、ソース項の実装エラー。自作コードなら、ソース項の積分（右辺ベクトルの組立）にバグがないか。商用ソフトなら、ユーザーサブルーチン内の座標変換（全体座標→自然座標）や積分点の重み付けが正しいか。第三に、製造解が「多項式の次数が高すぎる」ケース。要素の積分精度（ガウス積分点の数）が低くて、ソース項の積分自体に誤差が出ている可能性がある。

🧑‍🎓

収束率をプロットしたら、傾きが理論値より低くなりました（例えば、理論O(h^2)なのにO(h^1.5）。これはコードが間違っている証拠ですか？

🎓

必ずしもコードの根本的な間違いとは限らないが、何らかの「問題」があるサインだ。まず、誤差ノルムの計算が正しいか確認する。特に、H1ノルム（勾配の誤差）を計算する時、数値解の勾配をどのように計算しているか（節点値から差分か、要素内部の導関数か）。次に、メッシュが十分に細かいか。漸近収束領域に入っていない初期の粗いメッシュデータを含めてしまうと、傾きが低く見える。最後に、問題に特異性がないか。製造解自体は滑らかでも、領域の形状（リエントラントコーナーなど）や材料定数の不連続性が、数値解の正則性を損ない、収束率を低下させることがある。

🧑‍🎓

非線形問題で、ソルバーが収束しないか、別の解に落ち着いてしまいます。MMSで非線形問題を検証するのは無理ですか？

🎓

無理ではないが、初期条件とソルバー設定に細心の注意が必要だ。非線形ソルバー（ニュートン・ラフソン法）は初期推定値が解から遠いと発散する。MMSでは「解が分かっている」ので、初期条件を製造解に極めて近い値に設定できる。にもかかわらず収束しないなら、それはヤコビアン行列（接線剛性行列）の実装に誤りがある強力な証拠になる。また、材料パラメータ（例えば、弾塑性モデルの硬化係数）を極端な値にすると数値的に不安定になる。まずは線形弾性の範囲で検証し、少しずつ非線形性を強くしていく段階的アプローチが確実だ。