#discrete element method

[fruitsoftheweb]

EDEM Simulation of Potato Harvester: Assess Grading, Segregation Performance & Crop Damage

[fruitsoftheweb]

“Ship model cruising over an ice layer. DEM results obtained with DEMPACK“

[fruitsoftheweb]

“Particles settling under gravity. Yellow arrows represent forces exerted by fluid on particles.“

[fruitsoftheweb]

“Adding wall movement. Particle assembly is compressed by walls moving with constant velocity.“

[tamachanishere]

t=0.145付近でも時間刻み幅の切り直しが観測できました。

この場合、最大粒子速度が上昇を続け対応するvmaxに対する比が0.01を超えたところで時間刻み幅を切り直し。その後、最大の粒子速度は漸減して切り直し前程度に戻ったところでまたほぼ一定の速度になりました（計算している系が重力落下なのでそのぶん緩やかに加速しています）。

t=0.09付近でも速度が急に落ちる場所がありますが、こちらは時間刻みの調節とは関係ないようです。

[tamachanishere]

時間刻みの調節アルゴリズムをもう少し意味のあるものに直したいのですが、その前に現行のものの動作を確認。

横軸が時刻tで、それに対応する時間刻みΔtが紫でplotしてあります。初期値が0.00001で、刻み幅を大きくしながら進行していってt=0.0002を過ぎたあたりで後戻りしています。ここで時間刻み幅が大きすぎたのでrollbackしてから細かな幅で進め直しています。このあとはしばらくこの刻み幅で進んでいます。

各ステップで最も速い粒子の速度をplotしたものが一番上の緑で、t=0.0002あたりまで速度が増え続け、時間刻み幅の調節後はほぼ同じ値を維持しています。

青のラインは許容できる最大速度に対する上の最速粒子速度の比で、今回の計算ではこれが0.01を越えると時間刻み幅を小さく切り直すことにしてあります。これにもrollbackしている様子が見えています。この許容できる最大速度はシステムの特徴長さとその時点での時間刻み幅で決まっており、1ステップの最大移動距離の特徴長さに対する比と捉えることもできます。

[tamachanishere]

デモ計算：混合精度版DEM計算5（2倍スケール）

88200particles,  184820steps,  2871sec by nVidia Tesla K80 x 4

[tamachanishere]

デモ計算：混合精度版DEM計算5

10800particles,  130775steps,  414sec by nVidia Tesla K80 x 4