verlet-integration

    2

    1答えて

    周期境界条件または反射境界条件の2次元ボックスで、lennard-jones流体の分子動力学シミュレーションを作成しようとしています。シミュレーションは反射境界でうまくいくように見えますが、何らかの理由で周期的なボックス内の粒子が数百の積分ステップを経て余分な速度を得るようになり、最終的に粒子が箱に残っていません。 私は暴走速度がタイムステップを小さすぎるように設定することによって引き起こされる可

    1

    1答えて

    2つの物体の運動方程式を解くための私の休暇中のコードがありましたが、その結果は粒子が走っていてエラーがどこで見つけられなかったのですか import numpy as np import matplotlib.pyplot as plt plt.style.use('ggplot') DIM = 2 N = 2 ITER = 1000 def acc(r, v, a):