研究概要 |
われわれの研究対象は,スピノ-ダル分野と呼ぶ,合金など空間的に均一な組成が時間発展とともに空間的に不均一に変化してゆく,いわゆる異常拡散現象の数値解析である.その現象のモデルを記述する方程式として,Cahn-Hilliard方程式が知られている.この方程式は非線形項をもつ拡散方程式の形をしているが,拡散項の係数が負であることおよび非線形項の非線形性がかなり強いことにより,数値計算によって正確な解を求めることがかなりむずかしい. われわれは数学的考察と物理的考察を同時に行うことによって,Cahn-Hilliard方程式の数値解を求める差分スキームに対する安定性の条件を得た.それらの安定性を,F-安定および最適F-安定と名付けた.そして,F-安定および最適F-安定の条件を満たす差分パラメータの典型的な組合せを計算によって求め,その組合せの幾組かを使って実際に数値解を計算し,十分満足する結果を得た,以上の解析は空間変数2次元以上の場合にもそのまま適用でき,実際2次元問題のいくつかの例題に対しても十分安定な解を得ることができた. 次に,われわれはこの数値解の安定性をさらに詳細に調べるために,考えている物理系のリアプノフ関数を利用することを考えた.その解析過程で,より安定な差分スキームを得るためにはむしろリアプノフ関数を積極的に利用すべきであることに気付いた.そこで,このリアプノフ関数を出発点にとり,それを直接離散化し,それからCahn-Hilliard方程式に対応する差分スキームを合理的に導出した.こうして得られた差分スキームを使えば,数値解はいずれの場合も安定に計算することができることが確かめられた,次年度以降は,これらの結果に対する数学的に裏付けを中心に研究を続けたい.
|