スタガード連成手法
理論：交互に解くことで生じる時間の遅れ
概要
熱解析と構造解析を交互に回して連成させていますが、それで正しい答えになりますか？
結合が弱ければ、ほぼ正しい。交互に解く方法では、1つの時間ステップの中で、先に解く側はもう一方の前のステップの値を使う。この時間の遅れが誤差になり、結合が弱いうちは小さいが、結合が強いと大きくなる。さらに重要なのは安定性で、それぞれのソルバーが無条件に安定な陰解法でも、交互に値を渡す部分が陽解法のようにふるまうので、時間刻みが大きいと連成全体が発散することがある。2つの変数が結合した簡単な方程式で、その様子を確かめてみよう。
結合した2変数のモデル
u と v は2つの物理を代表する量。b と c が結合の強さ。弱い結合（b = c = 0.5）と、振動的な強い結合（b = 5、c = −5）を比べた。初期値 u = 1、v = 0、時刻2までの解を厳密解と比べる。
交互に解く方法（ガウス・ザイデル型）
それぞれの式は後退オイラー法（陰解法）で解く。ヤコビ型は v の式でも u^n（前の値）を使う。G は1ステップでの増幅行列で、そのスペクトル半径 ρ が1より小さければ安定。
ソルバーをつなぐ時代
流体、構造、熱、電磁界のソルバーは、それぞれの分野で何十年もかけて別々に発達してきた。1990年代以降、それらを組み合わせて連成問題を解く需要が高まり、既存のソルバーをそのまま使えるスタガード（分離型）の連成が広く使われるようになった。一方、ファーハットらは航空機の空力弾性の研究で、交互に解く方法の精度と安定性を詳しく解析し、時間の遅れを補う予測子や、ステップ内で繰り返す強い連成の手法を提案した。
計算例
例1：弱い結合（b = c = 0.5）での誤差（時刻2）
| 方法 | Δt = 0.2 | Δt = 0.1 | Δt = 0.05 | Δt = 0.025 | 観測次数 |
|---|---|---|---|---|---|
| 一体型 | 2.34×10⁻² | 1.19×10⁻² | 6.00×10⁻³ | 3.01×10⁻³ | 0.99 |
| ガウス・ザイデル型 | 2.34×10⁻² | 1.21×10⁻² | 6.16×10⁻³ | 3.11×10⁻³ | 0.99 |
| ヤコビ型 | 3.47×10⁻² | 1.79×10⁻² | 9.12×10⁻³ | 4.60×10⁻³ | 0.99 |
例2：振動的な強い結合（b = 5、c = −5）での安定性
| 時間刻み | 一体型 | ガウス・ザイデル型 | ヤコビ型 |
|---|---|---|---|
| 0.5 | 0.329 | 0.577 | 1.555（発散） |
| 1.0 | 0.180 | 3.283（発散） | 2.082（発散） |
| 2.0 | 0.093 | 6.122（発散） | 2.595（発散） |
表の値は増幅行列のスペクトル半径。安定の限界の時間刻みは、ガウス・ザイデル型で0.567、ヤコビ型で0.130。系の固有値は −1.5 ± 4.98i で、物理的には減衰する振動。
それぞれを陰解法で解いているのに、交互に解くと発散するんですね。
それぞれのソルバーの中では陰解法でも、相手の値は前のステップから持ってくるので、結合の部分は陽解法と同じ扱いになっている。陽解法には時間刻みの上限があり、結合が強いほどその上限は小さい。だから結合が弱い問題では交互に解いてもまったく問題ないが、結合が強い問題では時間刻みを小さくするか、1つのステップの中で交互の計算を収束するまで繰り返す（強い連成）必要がある。繰り返して収束させれば、結果は一体型と同じになり、安定性も戻る。どれくらい結合が強いかは、それぞれの物理の時間の速さと、互いにどれだけ影響し合うかで決まるので、まず両方を一度に解いた小さな問題で安定な時間刻みを確かめておくとよい。
連成計算の進め方
- それぞれの物理の時間の速さと、互いへの影響の大きさを見積もる。
- 結合が弱ければ一方向や交互の連成、強ければ強い連成や一体型を選ぶ。
- 時間刻みを変えて結果が変わらないこと、発散しないことを確かめる。
- 強い連成では、ステップ内の繰り返しの収束判定と緩和を設定する。
- 受け渡す量（温度、熱流束、変位、力）が保存されているか確認する。
「時間刻みを大きくしたら振動し始めた」
ある電磁界と熱の連成解析で、計算時間を短くしようと時間刻みを2倍にしたところ、温度が時間ステップごとに上下に振れ始め、やがて発散した。それぞれのソルバーは陰解法で、単独なら大きな刻みでも安定だった。原因は、温度で電気抵抗が大きく変わる材料で、温度と発熱の結合が強かったことだった。1ステップの中で電磁界と熱を3回ほど繰り返して収束させるようにすると、2倍の刻みでも安定に解け、全体の計算時間も短くなった。
よくある間違い
間違いと対策
| 間違い | 起きること | 対策 |
|---|---|---|
| 陰解法なら安定と考える | 連成で発散 | 連成全体の安定性を確認 |
| ヤコビ型で受け渡す | 安定範囲が狭い | ガウス・ザイデル型 |
| 強い結合で1回だけ受け渡す | 誤差と発散 | ステップ内で繰り返す |
| 時間刻みの確認をしない | 遅れの誤差を見逃す | 刻みを変えて比較 |
| 受け渡しで保存が崩れる | エネルギーが増減 | 保存型の写像 |
関連する内容も知りたいです。
関連するページは連成計算の収束、コシミュレーション、異なるメッシュ間の写像、モデル次数低減、マルチフィジックスのデジタルツインが参考になる。
なった
詳しく
報告