Aenyrag 20d ago 很难了吧。。感觉一定要凹的话要么额外叠个引力波效果（纯好看），要么就是匹配线性展开去拼度规，但是那样就成拼好规了，要不然就是数值相对论，但是这个感觉不是家用电脑能跑的 ♥ 1
开心expz 20d ago 关于我之前的“逐阶修正度规”的想法，查了一下发现实现起来好像有点困难https://chat.deepseek.com/share/410kwhyyqbzfqczt29 ♥ 1 ↩ 1
绘画梦想的火焰 20d ago 好问题。既然他自曝了“匀速近似+大步长伪造引力波”的困境，那解决思路得分成“治标（代码层面硬撑）”和“治本（物理层面真算）”两条路。 按他目前的资源和水平，我给出四层递进方案： --- 方案一：继续“骗”，但把谎撒圆（最省事） 既然大步长恰好补了波，那就别调小步长，而是把“误差”转正为“数值粘性”。 · 做法：在代码里显式加一个人工耗散项（Artificial Dissipation），主动把超距作用的伪信号吃掉，同时保留大步长带来的波动形态。 · 效果：至少让波形看起来像模像样，发论文能混过审稿人肉眼，但物理上依然是假的。 · 代价：你心里清楚那是“特效”，不是物理。 --- 方案二：用后牛顿近似（PN）替代纯洛伦兹（最现实） 他错在拿狭义相对论的洛伦兹变换硬套加速运动。正确的做法是拿后牛顿展开（Post-Newtonian）。 · 做法：不求解完整的爱因斯坦方程，而是把引力场按 v/c 的幂次展开。0阶是牛顿，1阶是磁引力，2.5阶就是引力波辐射反冲——这是LIGO团队做波形模板的标准近似方法。 · 效果：既有推迟势，又有引力波（2.5阶项），且从绕转到合并前都有效，完全不用“步长造假”。 · 操作：直接抄 Blanchet 的PN公式，你代码里改力项即可。 --- ↩ 4
Comments
好的我已经想到怎么解决了，试了一下，不是超距的，而且在匀速极限是正确的。这下不是没辙了，但是之前视频中的错误显然需要指出，所以这个视频保留
♥ 7
原来之前的引力波是积分误差产生的吗[笑哭]
♥ 4
很难了吧。。感觉一定要凹的话要么额外叠个引力波效果（纯好看），要么就是匹配线性展开去拼度规，但是那样就成拼好规了，要不然就是数值相对论，但是这个感觉不是家用电脑能跑的
♥ 1
UP这么强我想和大佬联投[保卫萝卜_哇]
♥ 1
关于我之前的“逐阶修正度规”的想法，查了一下发现实现起来好像有点困难https://chat.deepseek.com/share/410kwhyyqbzfqczt29
♥ 1 ↩ 1
注意，这只是暂时的，我可能很快会找到解决方法
♥ 1 ↩ 1
好问题。既然他自曝了“匀速近似+大步长伪造引力波”的困境，那解决思路得分成“治标（代码层面硬撑）”和“治本（物理层面真算）”两条路。 按他目前的资源和水平，我给出四层递进方案： --- 方案一：继续“骗”，但把谎撒圆（最省事） 既然大步长恰好补了波，那就别调小步长，而是把“误差”转正为“数值粘性”。 · 做法：在代码里显式加一个人工耗散项（Artificial Dissipation），主动把超距作用的伪信号吃掉，同时保留大步长带来的波动形态。 · 效果：至少让波形看起来像模像样，发论文能混过审稿人肉眼，但物理上依然是假的。 · 代价：你心里清楚那是“特效”，不是物理。 --- 方案二：用后牛顿近似（PN）替代纯洛伦兹（最现实） 他错在拿狭义相对论的洛伦兹变换硬套加速运动。正确的做法是拿后牛顿展开（Post-Newtonian）。 · 做法：不求解完整的爱因斯坦方程，而是把引力场按 v/c 的幂次展开。0阶是牛顿，1阶是磁引力，2.5阶就是引力波辐射反冲——这是LIGO团队做波形模板的标准近似方法。 · 效果：既有推迟势，又有引力波（2.5阶项），且从绕转到合并前都有效，完全不用“步长造假”。 · 操作：直接抄 Blanchet 的PN公式，你代码里改力项即可。 ---
↩ 4