私は、それぞれのソース項を介して2つのadv-diff結合ドメインがあるプロジェクトに取り組んでいます(1つのドメインは質量を加算し、他のドメインは質量を減算します)。簡潔にするために、定常状態でモデル化しています。方程式は、次のようなソース項を持つ標準の移流拡散輸送方程式です。
ここで、は種iの拡散および移流フラックスであり、Q iは種iのソース項です。
Newton-Raphson法を使用して問題のソルバーを記述でき、ブロック質量行列を使用して2つのドメインを完全に結合しました。
項は、ヤコビ行列を決定し、c 1とc 2の両方を更新するために使用されます。
または
速度を上げるために、反復ごとにヤコビアンを計算することはありません。現在、5回ごとに遊んでいます。これは十分に機能し、解を安定させているようです。
問題は、両方のドメインが2D / 2.5Dであるより大きなシステムに移動することであり、ヤコビ行列を計算すると、使用可能なコンピューターリソースがすぐに使い果たされます。後で最適化設定で使用するためにこのモデルを構築しているので、減衰係数の調整などの反復ごとにハンドルを握ることはできません。
私の問題に対してより堅牢でアルゴリズムを求めて他の場所を探すのは正しいのですか、それともこれで問題ないのですか?私は準線形化について少し調べましたが、それが自分のシステムにどの程度適用できるかはわかりません。
ヤコビアンを頻繁に再計算せずに非線形方程式系を解くことができる、見逃した他の巧妙なアルゴリズムはありますか?