個別要素法への陰解法の導入による数値移動床の改良

抄録

　数値移動床の基幹ツールである個別要素法のコードには一般に陽解法が用いられるので，複雑な境界条件下で流動抵抗が増加すると，粒子の重なりが過大となって異常反発が発生することがある．この対応策として，粒子間接触力の計算ルーチンに陰解法を導入し，抜本的なコードの改変を試みる．新しいコードはシリンダー堆積層崩壊の物理実験と比較して，従来の陽解法と同レベルの良好な再現性を示し，さらに，陽解法では安定した解が得られない急勾配斜面上でも安定に動作した．