24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jianchaoyv2651

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 131
帖子: 29
在线:
虫号: 657814
注册: 2008-11-19

[交流] 【求助】求一模拟程序（fortran）

本人初次编程解决实际问题，感觉力不从心！！请求各位高手帮忙！
问题如下：用fortran语言编程即时计算和展现十个氢原子在模拟盒中的运动情况。初速度满足麦克斯韦分布，采取硬球完全弹性碰撞

[ Last edited by lei0736 on 2009-11-25 at 11:58 ]

回复此楼

» 猜你喜欢

心脉受损已经有5人回复
博士读完未来一定会好吗已经有15人回复
Springer期刊投稿求助已经有4人回复
读博已经有3人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复

1楼2009-04-12 15:12:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jianchaoyv

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 964.9
帖子: 198
在线: 94.8小时
虫号: 657809
注册: 2008-11-19
专业: 通信理论与系统

请问我是按所给程序运行的，所得结果怎样转化为我在其他文献上看到此程序所得的总能量随步数演化的图形、模拟中动能的分布图以及其他一些物理量的值等，请问应该怎样去操作？
下面是我运行结果的前面的一部分，后面是从100到3000部分我未写出
MOLECULAR DYNAMICS SIMULATION PROGRAM

NUMBER OF PARTICLE IS          256
SIDE LENGTH OF THE BOX IS 6.752840
CUT OFF IS                2.500000
REDUCED TEMPERATURE IS    0.7220000
BASIC TIME STEP IS       6.4000003E-02

TOTAL LINEAR MOMENTUM IN X DIRECTION IS  6.7786448E-02
TOTAL LINEAR MOMENTUM IN Y DIRECTION IS -6.6548891E-02
TOTAL LINEAR MOMENTUM IN Z DIRECTION IS  0.1201641

VELOCITY ADJUSTMENT
TEMPERATURE BEFORE SCALING IS 45.41020
SCALE FACTOR IS  0.1260933
      2    270.799927 -1699.370728 -1428.570801
      4    252.771790 -1681.374146 -1428.602295
      6    218.253845 -1646.222656 -1427.968750
      8    171.139252 -1598.697998 -1427.558716
   10    132.214249 -1559.054932 -1426.840698
   12    116.545746 -1543.839722 -1427.293945
   14    119.052315 -1546.832153 -1427.779785
   16    127.827530 -1555.922119 -1428.094604
   18    134.569000 -1562.949707 -1428.380737
   20    137.285446 -1565.422729 -1428.137329
   22    137.336853 -1565.384644 -1428.047852
   24    136.563614 -1564.932617 -1428.369019
   26    134.083801 -1562.593628 -1428.509766
   28    129.720490 -1558.355103 -1428.634644
   30    126.987198 -1555.609863 -1428.622681
   32    127.910660 -1556.748169 -1428.837524
   34    131.486847 -1560.345459 -1428.858643
   36    135.715729 -1564.671997 -1428.956299
   38    138.538544 -1567.548828 -1429.010254
   40    139.642487 -1568.478882 -1428.836426
   42    140.040283 -1568.885742 -1428.845459
   44    140.190887 -1569.073975 -1428.883057
   46    139.712982 -1568.837280 -1429.124268
   48    138.952988 -1568.239868 -1429.286865
VELOCITY ADJUSTMENT
TEMPERATURE BEFORE SCALING IS  0.3634446
SCALE FACTOR IS 1.409449
   50    276.164978 -1568.357178 -1292.192139
   52    271.862701 -1564.260986 -1292.398315
   54    259.602783 -1552.116089 -1292.513306
   56    242.570480 -1535.268799 -1292.698364
   58    225.266388 -1518.148560 -1292.882202
   60    211.759369 -1504.485107 -1292.725708
   62    205.603088 -1498.357300 -1292.754150
   64    209.019333 -1501.927979 -1292.908691
   66    213.457825 -1506.205566 -1292.747803
   68    215.538666 -1508.412109 -1292.873413
   70    217.871109 -1510.706299 -1292.835205
   72    223.228973 -1516.197266 -1292.968262
   74    229.282425 -1522.213135 -1292.930664
   76    230.661316 -1523.383911 -1292.722656
   78    227.155609 -1519.937866 -1292.782227
   80    221.844513 -1514.623657 -1292.779175
   82    215.889771 -1508.402832 -1292.513062
   84    211.333252 -1503.536133 -1292.202881
   86    211.333344 -1503.780518 -1292.447144
   88    216.985687 -1509.352051 -1292.366333
   90    227.358459 -1519.550293 -1292.191895
   92    235.516296 -1527.578247 -1292.062012
   94    236.259933 -1528.119873 -1291.859985
   96    232.384796 -1524.102417 -1291.717651
   98    227.971985 -1519.898438 -1291.926514
VELOCITY ADJUSTMENT
TEMPERATURE BEFORE SCALING IS  0.5878155
SCALE FACTOR IS 1.108276
   100    276.164978 -1516.690308 -1240.525391