TECHNOLOGY EXPLAINER
マルチスケールシミュレーションとデジタルツイン
― 電子から部材までをつなぐ考え方と、「ツイン」の意味と限界
部材の性質は、その中の組織、原子の動き、電子の振る舞いという、はるかに小さなスケールの現象で決まります。ところが、全部を原子から計算することは、どんな計算機でもできません。そこで、電子→原子→組織→部材と、スケールごとの計算をつなぐのが、マルチスケールシミュレーションです。さらに、デジタルツインは、計算モデルを現物とデータでつなぎ続ける考え方です。この記事では、つなぎ方の型、誤差の伝わり方(本稿の計算)、デジタルツインの意味と、材料に当てはめるときの限界を、公式文書と論文で確認しながら整理します。
- 3行でわかる/読む前に/シリーズの中での位置
- 「スケール」とは:電子・原子・組織・部材
- なぜ、全部を一度には計算できないのか
- つなぎ方の型:値を渡す・往復する・同時に解く
- 下の誤差は、上でどう伝わるか
- デジタルツインとは何か
- 材料のデジタルツインと、その限界
- 結果を信じてよいかを確かめる
- 素材技術者から見て
- 限界と、確認できなかったこと/まとめ
- 評価・標準の読み方
- 用語集
- 参考資料(一次情報)
- 主張と出典の対応表(監査)
事実=公表資料に記載のある内容(出典リンクあり)
本稿の計算=明示した前提をもとに本記事が算出した値
未確定/将来=計画・目標・草案・改訂中・策定中・見通しなど、確定していない事項
このほか、素材技術者向けの読み取りは「解説」として区別しています。推測は太字で示します。
1. 3行でわかる/読む前に/シリーズの中での位置
- マルチスケールシミュレーションは、電子・原子・組織・部材という、大きさの違う段階の計算を「つなぐ」考え方です。つなぎ方には、下の結果を上へ渡す連鎖、往復して収束させる反復結合、同時に1つの式として解く密結合があります事実[[cwa17284]]
- つなぐと、下の段階の小さな誤差が、上の段階で膨らむことがあります。本稿の計算では、活性化障壁が 0.05 eV ずれると、300 K で速度が約6.9倍ずれます本稿の計算
- デジタルツインは、シミュレーションそのものではありません。National Academies の報告書は、「デジタルツインは、モデリングとシミュレーションだけではない」と述べ、物理の相手との双方向のやり取り(データによる更新と意思決定)を要素に挙げます事実[[nasem2024]]。材料では、状態(微細組織)を直接見られない難しさがあります(本記事による解説)
前提は、「原子の動きを計算する方法(分子動力学)」と「電子の振る舞いから計算する方法(DFT)」が、それぞれ別の大きさ・時間を扱う、ということです(本シリーズ「分子動力学」編、「密度汎関数理論の基礎」編)。この記事は単独でも読めるよう書いています。
共通の考え方:材料探索は広大な地図の中から宝のありかを探す作業で、計算は「安く速く候補を出す」、実験は「本当にそうか確かめる」、AIは「探す場所を絞る」という分担です。この記事で扱うのは、計算どうしをつなぐことと、計算と現物をつなぐことです。
「分子動力学」編と「機械学習原子間ポテンシャル」編(原子のスケール)の次に読むのがおすすめです。つないだ結果をどこまで信じてよいかは、「不確かさと適用範囲」編と「MIの結果の評価と再現性の読み方」編で扱います。国や機関の計画は「各国のMI戦略とデータ基盤」編、実験との閉ループは「クローズドループ最適化」編と「自律実験室」編が続きます。
2. 「スケール」とは:電子・原子・組織・部材
同じ金属の部品でも、見る大きさを変えると、別のものが見えます。電子の雲、原子の並び、結晶粒(小さな結晶の集まり)の模様、そして部品全体の形です。それぞれの段階で、計算で使う道具も、扱える大きさと時間も違います。
欧州の規格化団体 CEN の文書(CEN Workshop Agreement 17284、2018年4月)は、材料のモデルを、モデルが記述する「実体」で、まず4つに分けています。電子の実体、原子の実体、メソスケールの実体(原子の集まり:分子、ビーズ、結晶粒など)、連続体の体積要素です事実[[cwa17284]]。同じ材料を、どの実体で記述してもよい、とも述べています。
元のSVG・数値・出典を確認する
元SVGの説明:概念図(ベクター作図)。電子→原子→組織→部材の4段階。「実体」の4分類は CEN CWA 17284 による。スケールの目安は概念的なもので、方法ごとの厳密な範囲ではありません。
- 電子:第一原理計算(DFT)で、エネルギーの面と力を出します。下の段階に渡す値は、原子のあいだの力の決め方(ポテンシャル)です
- 原子:分子動力学で、拡散係数や界面エネルギー、弾性定数、熱伝導率などを計算します(本シリーズ「分子動力学」編)
- 組織(メソスケール):結晶粒や分子の集まりなどを扱います。フェーズフィールド法は、微細組織の変化を扱う方法で、Chen(2002)の総説があります事実[[chen2002]]。転位(結晶の中の線状の欠陥)の動きを扱う方法もあります。熱力学の計算(CALPHAD 法)も、この段階で使われます(CALPHAD の歴史は Spencer、2008)[[spencer2008]]
- 部材:有限要素法などで、部品の変形・熱・電気を計算します。CEN の文書の「連続体の体積要素」に当たります
上の段階は、下の段階の細かい様子をそのまま持ちません。たとえば、部材の計算では、結晶粒一つ一つを描かず、「平均すると、この材料の弾性率はこれ」という値を使います。下の段階の細かい様子を、平均した値にまとめて、上へ渡すことを均質化(アップスケーリング)といいます。逆に、部材にかかる力を下の段階に降ろして、局所の応力を見ることを局所化といいます。Kalidindi ら(2022)は、材料のデジタルツインで、スケールごとの代理モデルを、均質化と局所化の枠組みで用いる、と述べています事実[[kalidindi2022]]。
3. なぜ、全部を一度には計算できないのか
「それなら、部品全体を原子から計算すればよい」と思うかもしれません。計算の量を見積もってみましょう。本シリーズ「分子動力学」編で見たとおり、1 fs(10⁻¹⁵ 秒)のステップで1秒を進めるには 10¹⁵ 回、鉄 1 mm³ には約 8.5×10¹⁹ 個の原子があります本稿の計算。
- 最大級の記録:Germann と Kadau(2008)は1兆個の原子のMDを、Jia ら(2020)は、機械学習ポテンシャルと専用コードで、1億個を超える原子を1日に 1 ns 超進めたと報告しています事実[[germann2008,jia2020]]。1兆個の原子は、鉄なら一辺が約2.3 μm の立方体です本稿の計算
- 見積り:1 mm³ を 1 秒、全原子で計算すると、原子×ステップは約 8.5×10³⁴ です。Jia らの処理量(1日に約 10¹⁴ 以上)で割ると、約 8.5×10²⁰ 日、約 2×10¹⁸ 年です(上限の見積り)本稿の計算。最大級の記録でも、桁がまったく合いません
だから、全部を一度に原子で計算するのではなく、段階ごとに「必要な情報だけ」を取り出して、上の段階へ渡す工夫が要ります。2013年のノーベル化学賞は、「複雑な化学系のための多重スケールモデルの開発」に対して Karplus、Levitt、Warshel に授与されました。その業績は、ニュートンの古典物理学と量子物理学を、並べて使えるようにしたことだとされています事実[[nobel2013]]。
4. つなぎ方の型:値を渡す・往復する・同時に解く
つなぎ方には、大きく3つの型があります。CEN の文書は、次のように用語を定めています(3.14〜3.16)事実[[cwa17284]]。
- 連鎖(linking):モデルAの出力を、モデルBの入力にする。一方向で、流れは開いています。もっとも多い形です。たとえば、DFT で出した値から、分子動力学のポテンシャルを作り、分子動力学の拡散係数を、組織の計算に渡す、といった形です(本記事による解説)
- 反復結合(iterative coupling):AとBが、互いの出力を入力にして、収束するまで往復する。流れは閉じています
- 密結合(tight coupling):AとBの式を集めて、1つの連立方程式として同時に解く。結果は1つにまとまります
「同時」の代表例が、空間的に結合する方法です。
- QM/MM:反応が起きる中心だけを量子力学で、周りを古典力学で計算する方法です。Warshel と Levitt(1976)が、酵素(リゾチーム)の反応の理論的な研究で用いました[[warshel1976]]。2013年のノーベル化学賞の業績に含まれます[[nobel2013]]事実
- 原子と連続体の結合:Rudd と Broughton(1999)は、MEMS(微小な機械部品)の共振器で、中央の領域では数百万個の原子を動かし、周辺の大きな領域は有限要素で計算し、2つの方法が同時に動いて情報を往復させる方法を報告しています。原子の計算が必要なのは、力が非線形になる領域、表面積が大きい領域、欠陥による摩擦が予想される領域だと述べます事実[[rudd1999]]。Broughton らの論文(1999)は、有限要素・分子動力学・強結合の電子状態の結合を扱っています[[broughton1999]]
- 階層的に渡して、実験で確かめた例:Bulatov ら(2006)は、転位の3本のからみ合い(多重接合)の存在を、転位動力学と原子シミュレーションで予測し、モリブデンの透過電子顕微鏡観察で確認しました事実[[bulatov2006]]。複数のスケールの計算を組み合わせ、実験で確かめた例です
ほかに、粗視化という方法もあります。複数の原子を1つの玉にまとめ、大きな系を扱います。高分子の絡み合い(Kremer と Grest、1990)や、生体分子の粗視化力場(Marrink ら、2007、MARTINI)が例で、高分子のマルチスケールの総説(Kremer と Müller-Plathe、2002)があります事実[[kremer1990,marrink2007,kremer2002]]。
材料の分野で、計算を設計・製造に組み込む考え方をICME(Integrated Computational Materials Engineering、統合計算材料工学)と呼びます。米国 National Research Council の2008年の報告書は、ICME を「材料開発を加速し、設計と製造を統合できる、新しい分野」と位置づけています事実[[nrc2008]]。産業での適用例として、鋳造アルミ部品への ICME の適用を題材にした論文があります(Allison ら、2006)[[allison2006a,allison2006b]]。個別の効果の数字は、本稿では確認していません。
5. 下の誤差は、上でどう伝わるか
つなぐと、新しい問題が出ます。下の段階の誤差が、上の段階へ伝わることです。伝わり方には、2つの典型があります。
- 指数での増幅:拡散や反応の速さは、障壁 E に対して exp(−E/k_B T) で決まります。障壁が ΔE ずれると、速度の比は exp(ΔE/k_B T) です。ΔE=0.05 eV で、300 K では約6.9倍、600 K では約2.6倍、1000 K では約1.8倍です。ΔE=0.1 eV なら、300 K で約48倍です本稿の計算。下の段階(DFT)の「小さな」エネルギー誤差が、上の段階(速度)では大きな誤差になります。とくに低温で顕著です
- 段階を重ねる累積:各段階が「95%で2倍以内」だとします。誤差が独立なら、4段階で全体は約4倍(2^√4)に広がります。全てが同じ向きにそろうなら、16倍(2⁴)です本稿の計算。実際の連鎖では、同じ DFT の偏りが上の段階まで伝わるなど、相関がありえます(本記事による解説)
逆です。つないだ全体の精度は、最も弱い段階と、誤差の伝わり方で決まります。最下層に精密な DFT を使っても、途中のモデルが粗ければ、その精度は上へ届きません。さらに、図5のように、誤差は増幅・累積します。つないだ結果は、各段階の検証と、全体の実測との比較で確かめます(第8章)(本記事による解説)。
6. デジタルツインとは何か
次に、デジタルツインです。「ツイン(双子)」の名のとおり、現物(物理の相手)と、計算モデル(仮想の相手)が、データでつながり続ける考え方です。定義を、公式文書で確認します。
- NASA(2012):Glaessgen と Stargel は、将来の NASA と米空軍の機体の認証や保守が、材料特性の統計分布、経験的な設計の考え方、物理試験、試験と使用条件の相似の仮定に大きく頼っていて、過酷な要求には応えられなくなると述べ、その転換としてデジタルツインを提案しました。デジタルツインは、超高精度のシミュレーションを、機体の健全性管理システム、整備履歴、集められる過去のデータと組み合わせて、「飛んでいる双子」の寿命を写すものだとしています事実[[nasa2012]]
- National Academies(2024):デジタルツインを、「自然・工学・社会のシステムの構造・文脈・振る舞いを模した、仮想の情報の集まりで、物理の相手のデータで動的に更新され、予測能力をもち、価値を生む意思決定に情報を与えるもの。仮想と物理の双方向のやり取りが中心」と定義しています(AIAA の定義を改めたもの)。さらに、「デジタルツインは、シミュレーションとモデリングだけではない」(Finding 2-1)と述べます事実[[nasem2024]]
- 国際規格:ISO/IEC 30173:2023(Digital twin — Concepts and terminology)は、デジタルツインの用語と概念(システムの文脈、ライフサイクル、種類、機能の見方、関係者)を定め、2023年11月に発行されています。ISO 23247-1:2021 は、製造向けのデジタルツインの枠組みの概要と一般原則(用語、要件)で、2021年10月に発行されています事実[[iso30173,iso23247]]。本文は有料で、要求事項は確認していません
元のSVG・数値・出典を確認する
元SVGの説明:概念図(ベクター作図)。左:デジタルツインの構造(NASEM の整理にもとづく)。右:統合の度合いによる3段階(Kritzinger らの分類。説明は Barbie と Hasselbring による紹介)。
Kritzinger ら(2018)は、製造分野の論文を調べ、デジタルツインという言葉の使われ方に共通の理解がないとして、統合の度合いで、デジタルモデル、デジタルシャドウ、デジタルツインを区別しました。最も進んだ段階の「デジタルツイン」を扱う論文は少なかった、と報告しています事実[[kritzinger2018]]。3段階の中身(自動のデータ交換の有無と向き)は、Barbie と Hasselbring(2024)の紹介によると、デジタルモデル:自動のデータ交換なし/デジタルシャドウ:物理から仮想への一方向が自動/デジタルツイン:双方向とも自動です[[barbie2024]]。同じ論文は、この分類に対する Grieves の批判と、ISO 23247 には双方向の通信への言及がなく、この定義ではデジタルシャドウにあたる、という見方も紹介しています[[barbie2024]]。本稿は、定義には幅があることを、そのまま示します。
National Academies の報告書は、デジタルツインはシミュレーションとモデリングだけではなく、物理の相手との双方向のやり取り(センサーなどによるモデルの更新と、制御・意思決定)が要素だと述べています[[nasem2024]]。一度作った高精度のモデルを、現物のデータで更新せず、設計の検討だけに使うなら、それは「デジタルモデル」にあたります(Kritzinger らの分類)。名前でなく、更新の方向・頻度・使われ方で見てください(本記事による解説)。
7. 材料のデジタルツインと、その限界
機体や工場のデジタルツインは、センサーで測れる量(ひずみ、温度、振動)が多くあります。では、材料そのもののデジタルツインは、どう考えられているでしょうか。
- 考え方:Kalidindi ら(2022)は、材料そのものを、複数のスケールをもつ複雑な物理システムと見ます。その形(階層をもつ組織の詳細)と働き(外からの刺激への応答、つまり物性)を、デジタルツインで写す。デジタルツインは、加工の履歴と使用環境をふまえた、組織・加工・性能の時間変化を表せる、と主張しています。形はn点空間相関で、働きは均質化と局所化の代理モデルで表し、実験データと物理シミュレーションのデータで校正する、という枠組みです事実[[kalidindi2022]]
- 位置づけ:この論文は、基礎となる概念と枠組みを示すもので、実在の材料での実証を示すものではありません(本稿の読み取り)
- 難しい点:①状態(微細組織など)は、使用中に直接見られないことが多く、測れる量から推定する。②部材ごとに、加工の履歴が違い、個体ごとの更新が要る。③スケールをつなぐ各モデルの誤差が、全体に伝わる(図5)。以上は、本稿の整理(解説)です
National Academies の報告書は、デジタルツインの仮想の側が、決めたい判断と計算資源に合わせた「用途に合う(fit for purpose)」忠実度と解像度で作られるべきだとし(Conclusion 3-1)、高忠実度、簡略化、代理モデルを混ぜて使うこともありうる、と述べます事実[[nasem2024]]。材料の分野でも、何を判断したいかで、必要なスケールと精度が決まります(本記事による解説)。
8. 結果を信じてよいかを確かめる
つないだ計算も、デジタルツインも、結果を信じてよいかを確かめる仕事が、本体の半分を占めます。National Academies の報告書は、検証・妥当性確認・不確かさの定量(VVUQ)が、デジタルツインの全ての要素で役割をもつと述べます。仮想の側では、検証と妥当性確認が信頼性を作り、不確かさの定量が予測の質の目安を与える、とされています事実[[nasem2024]]。また、物理の相手の変化、データの収集方法の変化、データの分布の変化、予測や意思決定の課題の変化、モデルの更新のいずれもが、VVUQ に影響する、とも述べています[[nasem2024]]。
- 検証(verification):式を正しく解いているか(コードの誤り、数値の誤差、時間刻みなど)。本シリーズ「分子動力学」編の、エネルギー保存の確認が、その一例です
- 妥当性確認(validation):正しい式・モデルか。実測と、同じ量・同じ条件で比べる
- 不確かさの定量(UQ):入力やモデルのばらつきが、予測にどれだけ効くか。予測を幅つきで出す
- 規格:ASME の V&V 10-2019 は、計算固体力学の検証と妥当性確認の規格で、2025年に再確認(R2025)され、有効とされています事実[[asmevv10]]。つないだ連鎖の、連続体の段階の確認に関係します。本文は有料で、要求事項は確認していません
9. 素材技術者から見て
各段階の計算で、何を入力し、何を出力するかを、表にします。たとえば、「DFT:構造→活性化障壁」「分子動力学:温度・障壁→拡散係数」「組織の計算:拡散係数・界面エネルギー→組織の変化」「部材の計算:有効な弾性率→応力」。CEN の文書は、シミュレーションの記録様式(MODA)で、使い道、モデル、時間刻み、ソフトウェア、前の段階から受け取る入力を、項目として書き残す形を示しています[[cwa17284]]。どの値が、どの段階から来たかを残すと、後で、どこが怪しいかを追えます。
図5のとおり、指数的に効く量(活性化障壁、界面エネルギー)は、小さな誤差が大きな結果の差になります。まず、結果に最も効く入力を特定し、そこを実測で確かめ(たとえば、拡散の活性化エネルギーを、温度を変えた測定で求める)、モデルの値と比べます。実験の設計(対照、再現性、ロット差)を、効く入力から決めると、無駄が減ります(本記事による解説)。
「デジタルツイン」と大きく構えるより、小さな範囲で、実際に更新できる形から始めるのが現実的です。たとえば、自社の1つの工程(熱処理、成形)で、測れる量(温度履歴、ひずみ、硬さ)と、計算で出す量(組織、強度)を結び、測定値でモデルのパラメータを更新します。更新の頻度と方向が、「モデル」か「シャドウ」か「ツイン」かを分けます(図6)。個体ごとの履歴(ロット、加工条件)を記録する仕組みが、基礎になります(本記事による解説)。
モデルを更新するには、データが要ります。測定値だけでなく、試料の製法、熱履歴、組成の分析値、測定の条件と装置、ロットを付けて、表計算の自由記述ではなく、決まった項目で残してください。うまくいかなかった試験や、外れ値の理由も記録します。モデルが外れた条件は、適用範囲の境界を知る手がかりになります。
10. 限界と、確認できなかったこと/まとめ
- マルチスケールは、電子・原子・組織・部材をつなぐ考え方です。つなぎ方は、連鎖(一方向)、反復結合(往復)、密結合(同時)。QM/MM や、原子と連続体の結合が、同時結合の例です事実
- 全部を原子で計算できないから、つなぎます。1 mm³ を 1 秒、全原子で計算すると、最大級の記録の処理量でも約 2×10¹⁸ 年です本稿の計算
- つなぐと、誤差が伝わります。障壁の 0.05 eV の誤差が、300 K で速度の約6.9倍になり、各段階が「2倍以内」でも、4段階で約4倍〜16倍に広がりえます本稿の計算
- デジタルツインは、シミュレーションそのものではなく、物理の相手との双方向のやり取りを要素とします。材料では、微細組織を直接見られず、個体ごとの履歴が違う難しさがあります。検証・妥当性確認・不確かさの定量が、全ての要素に必要です事実
確認できなかったこと
- 図2のスケールの目安(長さ・時間)は、概念的な整理で、個別の文献で数値を確認していません
- QM/MM の Warshel と Levitt(1976)、Broughton ら(1999)、Allison ら(2006)、Kremer と Müller-Plathe(2002)、Chen(2002)、Spencer(2008)は、題名と、公開されている要旨・記述の範囲で扱いました。本文の詳細は確認していません
- CALPHAD 法とフェーズフィールド法の中身の説明は、一般的な理解にもとづく解説で、個別の文献の本文は確認していません
- Kritzinger らの分類の本文は確認できず、3段階の中身は、Barbie と Hasselbring(2024)の紹介によっています
- 図5の計算は、仮定(アレニウス型、独立な対数正規分布)にもとづく例で、実際の連鎖の誤差を見積もったものではありません
- ISO 23247-1、ISO/IEC 30173、ASME V&V 10 の本文(要求事項)は有料で、確認していません。題名・状態・適用範囲は公式ページで確認しました
- 材料向けのデジタルツインの国際標準は、確認できませんでした。材料でのデジタルツインの実証例の数と成熟度も、本稿では確認していません
11. 評価・標準の読み方
マルチスケール・デジタルツインの「ものさし」は、まず各段階の計算の検証(式を正しく解いているか)、次に実測との比較(妥当性確認)、不確かさ(幅つきの予測)、そして用途への適合です。ものさしが測っていないのは、段階をつないだ全体の、未知の条件での振る舞いです。確認日は2026年10月22日です。
| 評価項目 | 代表的なベンチマーク・評価手法(名称・出典) | 標準・指針(ISO/IEC・NIST・CEN など) | 読み方の注意 |
|---|---|---|---|
| 用語と、つなぎ方の記述 | CEN CWA 17284:実体(電子・原子・メソ・連続体)の分類。連鎖・反復結合・密結合の定義。MODA(記録様式)[[cwa17284]]事実 | CWA 17284:2018(CEN Workshop Agreement、2018年4月)。正式な規格ではない。対象は物理に基づくモデルで、データに基づくモデルは含まない[[cwa17284]]事実 | 記録の様式で、結果の正しさを保証するものではない |
| デジタルツインの用語・枠組み | NASEM(2024)の定義と結論[[nasem2024]]。Kritzinger らの3段階[[kritzinger2018]]事実 | ISO/IEC 30173:2023(用語と概念、2023年11月発行)[[iso30173]]。ISO 23247-1:2021(製造向けの枠組み、概要と一般原則、2021年10月発行)[[iso23247]]事実 | 本文は有料で、要求事項は未確認。ISO 23247 は双方向通信に触れていないとの見方がある(Barbie と Hasselbring)。定義には幅がある |
| 検証・妥当性確認・不確かさの定量(VVUQ) | NASEM(2024):デジタルツインの全ての要素に VVUQ が必要[[nasem2024]]事実 | ASME V&V 10-2019(R2025):計算固体力学の V&V の規格[[asmevv10]]事実 | 規格は連続体(固体力学)の段階が対象。電子・原子の段階の V&V の規格は、確認できなかった |
| つないだ連鎖の誤差の伝わり方 | 本稿の計算(図5):障壁の誤差の増幅、段階の累積[本稿]本稿の計算 | 確認できなかった | 仮定にもとづく例。自社の連鎖では、相関と入力の効き方を実際に見積もる |
| 産業での適用・政策的な位置づけ | NRC(2008):ICME の位置づけ[[nrc2008]]事実 | 確認できなかった(材料向けのデジタルツインの国際標準) | 報告書は政策・提言で、標準ではない。効果の数字は本稿では未確認 |
表の整理は本稿の計算(整理)。番号・題名・状態は ISO・ASME・CEN・NASEM の公式ページで確認(確認日 2026年10月22日)。規格の本文は有料のため、題名・適用範囲にとどめ、要求事項は書いていません。ASTM の規格は、本稿の範囲では該当を確認していません。
読み方の注意(素材技術者向け)
- 名前でなく、中身で判断してください:「デジタルツイン」と名乗っていても、更新が一方向、または手動なら、シャドウやモデルです。導入の効果を、名前でなく、更新の頻度・方向で確認します
- 予測と実験検証の違い:つないだ計算の予測は、各段階の検証と、全体の実測との比較を経て初めて、判断に使えます(Bulatov らの例は、計算の予測を実験で確かめた例)[[bulatov2006]]
- 標準・指針は、推奨か要求か:CWA は正式な規格ではなく、ISO 30173・23247 は用語・枠組み、ASME V&V 10 は自主的な規格で、いずれも、使う側が契約や法令で引用して初めて拘束力をもつのが一般的です(本記事による解説)
- 計画と達成:デジタルツインの導入計画や、国・企業の構想は、目標です。実際に更新が回っているかを、確認してください
12. 用語集
- マルチスケールシミュレーション
- 電子・原子・組織・部材など、大きさの違う段階の計算をつなぐ考え方。
- 実体(エンティティ)
- CEN の文書で、モデルが記述する単位。電子・原子・メソスケール・連続体の4種。
- メソスケール
- 原子と部材の中間(結晶粒、分子の集まりなど)の大きさ。
- 連続体
- 材料を、原子に分けず、連続した物質として扱う見方。有限要素法など。
- 連鎖(linking)
- モデルAの出力をBの入力にする、一方向の結合。
- 反復結合(iterative coupling)
- 2つのモデルが、互いの出力を入力にして、収束するまで往復する結合。
- 密結合(tight coupling)
- 2つのモデルの式を、1つの連立方程式として同時に解く結合。
- QM/MM
- 中心を量子力学(QM)、周りを古典力学(MM)で計算する方法。
- 均質化(アップスケーリング)
- 下の段階の細かい様子を、平均した値にまとめて、上の段階へ渡すこと。
- 局所化
- 部材の条件を、下の段階に降ろして、局所の応力などを見ること。
- 粗視化
- 複数の原子を1つの玉にまとめ、大きな系を扱う方法。
- フェーズフィールド法
- 微細組織の変化を、連続な場として扱う方法。
- CALPHAD
- 熱力学データベースを使い、相の安定性などを計算する方法。
- ICME
- 統合計算材料工学。計算を、材料の開発・設計・製造に組み込む考え方。
- デジタルツイン
- 物理の相手と仮想モデルが、データで双方向につながり続ける考え方。
- デジタルモデル/シャドウ/ツイン
- 自動のデータ交換の有無と向きによる3段階(Kritzinger らの分類)。
- データ同化
- 観測データで、モデルの状態やパラメータを更新すること。
- 代理モデル
- 重い計算を、速く近似する簡易なモデル(機械学習など)。
- 検証/妥当性確認(V&V)
- 式を正しく解いているか/正しい式・モデルか。
- 不確かさの定量(UQ)
- 入力やモデルのばらつきが、予測にどう効くかを定量化すること。VVUQ は3つの総称。
- fit for purpose
- 決めたい判断に足りる精度・解像度であること。
- MODA
- materials modelling data。CEN の文書が示す、シミュレーションの記録様式。
13. 参考資料(一次情報)
- CEN (European Committee for Standardization) "CEN Workshop Agreement CWA 17284:2018, Materials modelling - Terminology, classification and metadata" (2018) — www.cencenelec.eu/media/CEN-CENELEC/CWAs/RI/cwa17284_2018.pdf
- The Royal Swedish Academy of Sciences "The Nobel Prize in Chemistry 2013, press release (9 October 2013)" — www.nobelprize.org/prizes/chemistry/2013/press-release/
- Warshel A., Levitt M. "Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme", J. Mol. Biol. 103, 227–249 (1976), DOI 10.1016/0022-2836(76)90311-9 — doi.org/10.1016/0022-2836(76)90311-9
- Rudd R. E., Broughton J. Q. "Coupling of length scales and atomistic simulation of MEMS resonators", SPIE Proceedings 3680, 104 (arXiv:cond-mat/9812232) (1999), DOI 10.1117/12.341254 — doi.org/10.1117/12.341254
- Broughton J. Q., Abraham F. F., Bernstein N., Kaxiras E. "Concurrent coupling of length scales: Methodology and application", Phys. Rev. B 60, 2391–2403 (1999), DOI 10.1103/PhysRevB.60.2391 — doi.org/10.1103/PhysRevB.60.2391
- Bulatov V. V. ほか "Dislocation multi-junctions and strain hardening", Nature 440, 1174–1178 (2006), DOI 10.1038/nature04658 — doi.org/10.1038/nature04658
- Kremer K., Grest G. S. "Dynamics of entangled linear polymer melts: A molecular-dynamics simulation", J. Chem. Phys. 92, 5057–5086 (1990), DOI 10.1063/1.458541 — doi.org/10.1063/1.458541
- Kremer K., Müller-Plathe F. "Multiscale simulation in polymer science", Mol. Simul. 28, 729–750 (2002), DOI 10.1080/0892702021000002458 — doi.org/10.1080/0892702021000002458
- Marrink S. J. ほか "The MARTINI Force Field: Coarse Grained Model for Biomolecular Simulations", J. Phys. Chem. B 111, 7812–7824 (2007), DOI 10.1021/jp071097f — doi.org/10.1021/jp071097f
- Chen L.-Q. "Phase-Field Models for Microstructure Evolution", Annu. Rev. Mater. Res. 32, 113–140 (2002), DOI 10.1146/annurev.matsci.32.112001.132041 — doi.org/10.1146/annurev.matsci.32.112001.132041
- Spencer P. J. "A brief history of CALPHAD", Calphad 32, 1–8 (2008), DOI 10.1016/j.calphad.2007.10.001 — doi.org/10.1016/j.calphad.2007.10.001
- Kalidindi S. R., Buzzy M., Boyce B. L., Dingreville R. "Digital Twins for Materials", Front. Mater. 9, 818535 (2022), DOI 10.3389/fmats.2022.818535 — doi.org/10.3389/fmats.2022.818535
- National Research Council (U.S.) "Integrated Computational Materials Engineering: A Transformational Discipline for Improved Competitiveness and National Security", The National Academies Press (2008), DOI 10.17226/12199 — doi.org/10.17226/12199
- Allison J., Backman D., Christodoulou L. "Integrated computational materials engineering: A new paradigm for the global materials profession", JOM 58 (11), 25–27 (2006), DOI 10.1007/s11837-006-0223-5 — doi.org/10.1007/s11837-006-0223-5
- Allison J. ほか "Virtual aluminum castings: An industrial application of ICME", JOM 58 (11), 28–35 (2006), DOI 10.1007/s11837-006-0224-4 — doi.org/10.1007/s11837-006-0224-4
- Jia W. ほか "Pushing the limit of molecular dynamics with ab initio accuracy to 100 million atoms with machine learning", arXiv:2005.00223 (SC20) (2020) — arxiv.org/abs/2005.00223
- Germann T. C., Kadau K. "Trillion-atom molecular dynamics becomes a reality", Int. J. Mod. Phys. C 19, 1315–1319 (2008), DOI 10.1142/S0129183108012911 — doi.org/10.1142/S0129183108012911
- Glaessgen E., Stargel D. "The Digital Twin Paradigm for Future NASA and U.S. Air Force Vehicles", NASA Technical Reports Server 20120008178 (AIAA 2012-1818) (2012) — ntrs.nasa.gov/citations/20120008178
- National Academies of Sciences, Engineering, and Medicine "Foundational Research Gaps and Future Directions for Digital Twins (Summary chapter)", The National Academies Press (2024) — nap.nationalacademies.org/read/26894/chapter/2
- Kritzinger W. ほか "Digital Twin in manufacturing: A categorical literature review and classification", IFAC-PapersOnLine 51 (11), 1016–1022 (2018), DOI 10.1016/j.ifacol.2018.08.474 — doi.org/10.1016/j.ifacol.2018.08.474
- Barbie A., Hasselbring W. "From Digital Twins to Digital Twin Prototypes: Concepts, Formalization, and Applications", IEEE Access arXiv:2401.07985 (2024), DOI 10.1109/ACCESS.2024.3406510 — doi.org/10.1109/ACCESS.2024.3406510
- ISO ISO 23247-1:2021 Automation systems and integration — Digital twin framework for manufacturing — Part 1: Overview and general principles (page checked 22 October 2026) (2021) — www.iso.org/standard/75066.html
- ISO/IEC ISO/IEC 30173:2023 Digital twin — Concepts and terminology (page checked 22 October 2026) (2023) — www.iso.org/standard/81442.html
- ASME V&V 10-2019 (R2025) Standard for Verification and Validation in Computational Solid Mechanics (page checked 22 October 2026) (2019) — www.asme.org/codes-standards/find-codes-standards/standard-for-verification-and-validation-in-computational-solid-mechanics
14. 主張と出典の対応表(監査)
| 本文の主張 | 根拠 | 区分 |
|---|---|---|
| CEN CWA 17284:2018(2018年4月):モデルを実体(電子・原子・メソスケール・連続体の体積要素)で分類。連鎖(linking)・反復結合・密結合の定義(3.14〜3.16)。MODA は記録様式。対象は物理に基づくモデル。正式な規格ではない | 出典1 https://www.cencenelec.eu/media/CEN-CENELEC/CWAs/RI/cwa17284_2018.pdf | 事実 |
| 本稿の計算:1 fs のステップで 1 s は 10¹⁵ 回。鉄 1 mm³ は約 8.5×10¹⁹ 個、10¹² 個は一辺約2.3 μm(密度 7.87 g/cm³、モル質量 55.85 g/mol からの概算) | 本稿の計算 | 本稿の計算 |
| 最大級の記録:Germann と Kadau(2008)の1兆個(レナード=ジョーンズ粒子)。Jia ら(2020)の1億個超、1日で 1 ns 超(Summit) | 出典17 https://doi.org/10.1142/S0129183108012911 出典16 https://arxiv.org/abs/2005.00223 | 事実 |
| 1 mm³ を 1 秒、全原子で計算する見積り:原子×ステップ約 8.5×10³⁴、Jia らの処理量(1日に約 10¹⁴ 以上)で約 8.5×10²⁰ 日、約 2×10¹⁸ 年(上限の見積り) | 本稿の計算(Jia らの値は下限なので、年数は上限) | 本稿の計算 |
| 2013年のノーベル化学賞:複雑な化学系のための多重スケールモデルの開発(Karplus、Levitt、Warshel)。古典物理と量子物理を並べて使えるようにした | 出典2 https://www.nobelprize.org/prizes/chemistry/2013/press-release/ | 事実 |
| QM/MM:Warshel と Levitt(1976)、酵素(リゾチーム)の反応の理論的研究 | 出典3 https://doi.org/10.1016/0022-2836(76)90311-9 | 事実 |
| Rudd と Broughton:MEMS 共振器で、中央は数百万個の原子、周辺は有限要素。2つの方法が同時に動き、情報を往復させる。Broughton ら:有限要素・分子動力学・強結合の同時結合(題名) | 出典4 https://doi.org/10.1117/12.341254 出典5 https://doi.org/10.1103/PhysRevB.60.2391 | 事実 |
| Bulatov ら:転位の多重接合を、転位動力学と原子シミュレーションで予測し、モリブデンの透過電子顕微鏡観察で確認 | 出典6 https://doi.org/10.1038/nature04658 | 事実 |
| 粗視化の例:Kremer と Grest(1990、絡み合った高分子溶融体)、Marrink ら(2007、MARTINI)、Kremer と Müller-Plathe(2002、高分子のマルチスケール) | 出典7 https://doi.org/10.1063/1.458541 出典9 https://doi.org/10.1021/jp071097f 出典8 https://doi.org/10.1080/0892702021000002458 | 事実 |
| フェーズフィールド法の総説(Chen、2002)。CALPHAD の歴史(Spencer、2008)。両法の中身の説明は一般的な理解にもとづく解説 | 出典10 https://doi.org/10.1146/annurev.matsci.32.112001.132041 出典11 https://doi.org/10.1016/j.calphad.2007.10.001 | 事実 |
| Kalidindi ら(2022):材料を複数スケールの物理システムと見て、形(n点空間相関)と働き(均質化・局所化の代理モデル)を写す材料のデジタルツインの枠組み。実験データと物理シミュレーションで校正 | 出典12 https://doi.org/10.3389/fmats.2022.818535 | 事実 |
| ICME:NRC(2008)は、材料開発を加速し設計と製造を統合できる新しい分野と位置づけ。Allison ら(2006):ICME の新しい枠組みと、鋳造アルミ部品への適用(題名) | 出典13 https://doi.org/10.17226/12199 出典14 https://doi.org/10.1007/s11837-006-0223-5 出典15 https://doi.org/10.1007/s11837-006-0224-4 | 事実 |
| 活性化障壁の誤差による速度の比 exp(ΔE/k_B T):0.05 eV で 300 K 約6.9倍、600 K 約2.6倍、1000 K 約1.8倍、0.1 eV・300 K で約48倍(k_B=8.617×10⁻⁵ eV/K)。各段階が95%で2倍以内のとき、4段階で独立なら約4倍、完全に同じ向きなら16倍(仮定つき) | 本稿の計算(アレニウス型・独立な対数正規分布の仮定) | 本稿の計算 |
| NASA(2012):認証・保守が材料特性の統計分布や物理試験などに頼ること、デジタルツインは超高精度のシミュレーションと機体の健全性管理・整備履歴・過去のデータを統合し、「飛んでいる双子」の寿命を写す | 出典18 https://ntrs.nasa.gov/citations/20120008178 | 事実 |
| National Academies(2024):デジタルツインの定義(双方向のやり取りが中心)、Finding 2-1(シミュレーションとモデリングだけではない)、Conclusion 3-1(fit for purpose)、VVUQ は全ての要素に役割 | 出典19 https://nap.nationalacademies.org/read/26894/chapter/2 | 事実 |
| Kritzinger ら(2018):デジタルモデル・シャドウ・ツインを統合の度合いで区別。デジタルツインを扱う論文は少ない。3段階の中身(自動のデータ交換の有無と向き)、Grieves の批判、ISO 23247 の評価は Barbie と Hasselbring(2024)の紹介 | 出典20 https://doi.org/10.1016/j.ifacol.2018.08.474 出典21 https://doi.org/10.1109/ACCESS.2024.3406510 | 事実 |
| ISO/IEC 30173:2023(2023年11月、デジタルツインの用語と概念)、ISO 23247-1:2021(2021年10月、製造向けの枠組みの概要と一般原則)。本文は有料で未確認 | 出典23 https://www.iso.org/standard/81442.html 出典22 https://www.iso.org/standard/75066.html | 事実 |
| ASME V&V 10-2019(R2025):計算固体力学の検証と妥当性確認の規格。2025年に再確認、有効 | 出典24 https://www.asme.org/codes-standards/find-codes-standards/standard-for-verification-and-validation-in-computational-solid-mechanics | 事実 |
| 素材技術者の視点(入力と出力の表、効く入力の洗い出し、小さく始める、データの整え方)。均質化・局所化・用途への適合の関連づけ、材料のデジタルツインの難しい点の整理は解説 | 解説(本稿の見解) | |
| 材料向けのデジタルツインの国際標準、図2のスケールの数値の文献、実証例の数と成熟度は確認できなかった | 本稿の確認結果(確認できなかった) |
最終更新:2026年10月22日/出典は一次情報(論文、NASA・National Academies・NRC の報告書、CEN・ISO・ASME・Nobel 財団の公式ページ)に限定しています。数値例は本稿の計算で、実際の連鎖の誤差を示すものではありません。スケールの目安の数値、一部の論文の本文の詳細、Kritzinger らの分類の本文、各規格の本文(要求事項)、材料向けのデジタルツインの国際標準と実証例の数については、確認できなかったため、断定していません。図2〜図8は説明用の概念図または本稿の計算を含む作図で、図1と図9はAI生成イメージ(描く場面の指示のみ記載)です。



