diff --git a/examples/case04/Readme.html b/examples/case04/Readme.html index 22d9a9b..fded8aa 100644 --- a/examples/case04/Readme.html +++ b/examples/case04/Readme.html @@ -211,14 +211,15 @@ MathJax = {

2.1 质量

- - - - + + + +
天体质量 / kg与地球质量之比
太阳$1.989 \times 10^{30}$$333\,000$
地球$5.972 \times 10^{24}$$1$
月球$7.35 \times 10^{22}$$0.0123 \;(1/81.3)$
天体质量 / kg模拟质量(月球=1)
太阳$1.989 \times 10^{30}$$27\,000$
地球$5.972 \times 10^{24}$$81$
月球$7.35 \times 10^{22}$$1$
- 质量比:$M_\odot : M_\oplus : M_\text{moon} \approx 27\,100\,000 : 81.3 : 1$ + 质量比:$M_\odot : M_\oplus : M_\text{moon} \approx 27\,100\,000 : 81.3 : 1$
+ 模拟中采用缩放比例 $27\,000 : 81 : 1$(太阳质量缩至 $1/1000$ 以保持数值稳定)。

2.2 轨道参数

@@ -310,8 +311,8 @@ MathJax = { - - + + @@ -321,12 +322,12 @@ MathJax = {
参数说明
gravity_strength100.0万有引力强度(已缩放)
box_a20.0盒子半边长
gravity_strength100.0万有引力强度
box_a30.0盒子半边长
methodleapfrog蛙跳法(辛积分器,能量守恒)
T_total10.0 s总模拟时间
NSTEP2抽帧间隔(密采帧)
- - - + + +
天体质量位置 $(x,y,z)$速度 $(v_x,v_y,v_z)$固定约束
太阳1$(0,0,0)$$(0,0,0)$全部固定
地球1$(4,0,0)$$(0,0,4)$
月球0.1$(5,0,0)$$(0,0,6)$
太阳$27\,000$$(0,0,0)$$(0,0,0)$全部固定
地球$81$$(10,0,0)$$(0,0,520)$
月球$1$$(11.5,0,0)$$(0,0,526)$
-

地球在 $z$ 方向获得初速度,产生绕太阳的轨道运动;月球在地球基础上获得额外的 $z$ 方向速度,形成绕地球的轨道。

+

地球在 $z$ 方向获得初速度 $v=520$,产生绕太阳的轨道运动;月球在地球基础上附加 $v=6$ 的绕地速度,形成绕地球的微小轨道。速度由圆形轨道公式 $v = \sqrt{G_{\text{eff}} M / r}$ 计算,其中 $G_{\text{eff}} = \text{gravity\_strength} = 100$。

五、使用方法

@@ -357,8 +358,8 @@ MathJax = {

七、已知局限

diff --git a/examples/case04/input/coord.txt b/examples/case04/input/coord.txt index 7b5d01a..036ade1 100644 --- a/examples/case04/input/coord.txt +++ b/examples/case04/input/coord.txt @@ -1,4 +1,4 @@ n mass radius x y z vx vy vz fix_x fix_y fix_z -1 1 0.28 0 0 0 0 0 0 1 1 1 -2 1 0.28 4 0 0 0 0 4 0 0 0 -3 0.1 0.18 5 0 0 0 0 6 0 0 0 \ No newline at end of file +1 27000 2.0 0 0 0 0 0 0 1 1 1 +2 81 0.5 10 0 0 0 0 520 0 0 0 +3 1 0.3 11.5 0 0 0 0 526 0 0 0 \ No newline at end of file diff --git a/examples/case04/input/input.txt b/examples/case04/input/input.txt index 6d70e56..1bf6fbd 100644 --- a/examples/case04/input/input.txt +++ b/examples/case04/input/input.txt @@ -23,7 +23,7 @@ force_calc: 0 # 强制重新计算:1=跳过缓存强算,0=自动使用 engine: python # 默认使用 Python 引擎 # ── 盒子 ────────────────────────────────────── -box_a: 20.0 # 立方体半边长,粒子被限制在 [-box_a, box_a]³ 内 +box_a: 30.0 # 立方体半边长,粒子被限制在 [-box_a, box_a]³ 内 # ── 初始构型 ────────────────────────────────── # 坐标文件格式: