by 管理人 » 2026年2月21日(土) 11:15
1粒子系では自由落下のような単純な運動しか検証できません。したがって、まず2粒子系(粒子―粒子、粒子―壁)の運動を対象とし、解析解と数値解が一致するかどうかを確認します。減衰振動など、力学の教科書に掲載されている典型例で確認します。
3粒子以上になると一般に解析解は存在しないため、可視化によって不自然な挙動が生じていないかを確認します。さらにエネルギー収支を逐次チェックすることで、時間積分が正しく実装されているかどうかを検証できます。
モデルが独自のものでない場合には、信頼できるオープンソースのプログラムと同じ系で計算結果を比較し、丸め誤差の範囲内で一致することを確認するのが望ましいと考えます。
1粒子系では自由落下のような単純な運動しか検証できません。したがって、まず2粒子系(粒子―粒子、粒子―壁)の運動を対象とし、解析解と数値解が一致するかどうかを確認します。減衰振動など、力学の教科書に掲載されている典型例で確認します。
3粒子以上になると一般に解析解は存在しないため、可視化によって不自然な挙動が生じていないかを確認します。さらにエネルギー収支を逐次チェックすることで、時間積分が正しく実装されているかどうかを検証できます。
モデルが独自のものでない場合には、信頼できるオープンソースのプログラムと同じ系で計算結果を比較し、丸め誤差の範囲内で一致することを確認するのが望ましいと考えます。