1 定義と基本的考え方
カーン・ヒリアード型は、保存される量の空間分布が時間発展する様子を記述する偏微分方程式の枠組みである。主に濃度や組成の不均一が、拡散と熱力学的な駆動力によってどう変化するかを表し、相分離や界面形成を扱う際の基礎モデルとして用いられる。
1.1 モデルの目的
このモデルの目的は、混合状態の中に生じた小さなゆらぎが、どのように増幅されて領域分離へ至るかを定式化することにある。単なる拡散ではなく、系全体のエネルギーを考慮しながら、構造の発達や粗大化の過程を追跡できる点が特徴である。
1.2 保存量と秩序変数
カーン・ヒリアード型では、対象となる量が系の内部で保存される。これを表すために、秩序変数として濃度場や組成場を導入し、その値が時空間で変化する様子を解析する。
1.2.1 濃度場の扱い
濃度場は、各点における成分の割合を連続的に表す。局所的な増減は起こっても、系全体としての総量は保たれるため、変化は主として流れではなく拡散として現れる。
1.2.2 自由エネルギーとの関係
秩序変数の分布は、自由エネルギー汎関数によって評価される。系はより安定な状態へ向かうため、エネルギーが低くなる方向に進み、その結果として相分離や界面の生成が生じる。
1.3 勾配流としての解釈
この方程式は、自由エネルギーの勾配に沿って状態が変化する勾配流として理解できる。時間発展は、エネルギーを最小化する方向に進むが、保存則の制約があるため、単純な緩和とは異なる特徴を示す。
2 基本方程式
標準的なカーン・ヒリアード方程式は、秩序変数の時間変化を化学ポテンシャルの空間的な不均一と結びつけて表す。通常は四階の空間微分を含み、界面を滑らかに保ちながら相の分離を表現する。
2.1 標準形
典型的な形では、秩序変数の時間微分が、化学ポテンシャルのラプラシアンに比例する。化学ポテンシャル自体は自由エネルギーの変分微分として定義され、保存量の拡散を駆動する役割を果たす。
2.1.1 化学ポテンシャル
化学ポテンシャルは、局所的な状態がどれだけ不安定かを示す量である。自由エネルギーを秩序変数に関して変分すると得られ、系の変化方向を決める主要な要素となる。
2.1.2 拡散項
拡散項は、濃度の急な変化をならしつつ、全体としては相分離を促進する働きを持つ。通常の拡散方程式と異なり、高次の空間微分が現れるため、界面の厚みや波長選択を表現しやすい。
2.2 境界条件
方程式を解くには、領域の端での振る舞いを定める必要がある。境界条件は保存量の出入りや勾配の制約を規定し、解の形や安定性に大きく影響する。
2.2.1 ノイマン型条件
ノイマン型条件では、境界を横切る流束をゼロにすることが多い。これにより、領域外への物質移動がない閉じた系として扱える。
2.2.2 周期境界条件
周期境界条件では、領域の端と端をつなげたように扱う。数値計算で用いられることが多く、有限領域でも無限に繰り返す系の近似として有効である。
2.3 初期条件
初期条件は、計算開始時点での濃度分布を与える。均一状態に微小な乱れを加える設定がよく用いられ、そこから相分離がどのように進むかを観察する。
3 理論的性質
カーン・ヒリアード型は、保存則を保ちながらエネルギーを下げていくという明確な構造をもつ。そのため、理論解析では散逸性、収束性、解の存在といった性質が重要になる。
3.1 エネルギー散逸
このモデルでは、時間とともに自由エネルギーが減少する。散逸構造が組み込まれているため、発展は無秩序な増幅ではなく、一定の方向をもった緩和として現れる。
3.1.1 自由エネルギーの減少
適切な条件の下では、自由エネルギーは単調に減少する。これは系がより安定な配置へ向かうことを意味し、理論的な安定性解析の基盤となる。
3.1.2 平衡解への収束
長時間後には、変化が小さくなり、平衡状態あるいは準安定状態へ近づくことが多い。実際には局所最小に留まる場合もあり、複数の終状態が存在しうる。
3.2 保存則
秩序変数の総量は、外部とのやり取りがなければ保存される。これはこの方程式の中心的性質であり、単純な拡散系と区別される重要な点である。
3.3 解の存在と一意性
数学的には、初期値や境界条件に応じて、解が存在するかどうか、また一意に定まるかが検討される。高次の非線形性を含むため、解析は容易ではないが、適切な関数空間では多くの結果が知られている。
4 応用と拡張
カーン・ヒリアード型は、相分離を伴う現象の記述だけでなく、複合場の相互作用を扱う拡張へも発展している。材料、流体、反応系、確率過程などに広く応用される。
4.1 材料科学への応用
材料科学では、混合物の成分分布や相構造の形成を調べる道具として活用される。合金、ポリマー、薄膜などの微視的な組織変化を記述する際に有用である。
4.1.1 相分離の解析
相分離の解析では、初期の微小ゆらぎが成長し、異なる相が形成される過程を追う。ドメインの拡大速度や界面の挙動を調べることで、実験結果との比較も行われる。
4.1.2 微細構造の予測
微細構造の予測では、時間発展の計算から最終的な模様や粒子配置を見積もる。これにより、材料特性の設計や加工条件の検討に役立つ。
4.2 流体との結合
流体が関与する場合には、濃度場だけでなく速度場も考慮する必要がある。対流と拡散が同時に働くため、構造形成はより複雑になる。
4.2.1 カーン・ヒリアード・ナビエ・ストークス系
この結合系では、カーン・ヒリアード方程式が流体の運動方程式と連立する。相分離する混相流や界面運動を扱う際に使われ、相の移動と流体力学的な効果を同時に記述できる。
4.3 反応拡散系への拡張
反応項を加えると、保存されない生成・消滅過程を取り込める。これにより、生体材料や化学反応を伴う場でのパターン形成を表現しやすくなる。
4.4 確率的モデル
実際の系では、熱ゆらぎや外乱が無視できないことがある。そのため、決定論的な方程式に確率項を加えたモデルが用いられる。
4.4.1 ゆらぎの導入
ゆらぎの導入は、微視的な不規則性や温度効果を反映する。これにより、同じ初期条件でも異なるパターンが現れる可能性を表せる。
4.4.2 数値計算上の工夫
確率項を含む場合、数値計算では安定性や保存則の扱いに注意が必要である。離散化の方法や時間刻みの選択を工夫することで、物理的整合性を保ちながら計算を進める。