C语言化学反应模拟程序while循环提前终止问题排查

• 未初始化变量引发未定义行为
  代码中reactionRate、nitrogenCurrent、waterCurrent、waterProduced四个浮点变量未做初始化，第一次执行while条件判断时，reactionRate读取的是内存中的随机垃圾值，同时未初始化的浓度变量会让后续浓度计算完全偏离真实反应轨迹，氨浓度下降速度被异常放大，很快就低于初始值的一半（1.0mol/L）。
• 循环逻辑运算符使用错误
  你的需求是循环需要在「氨浓度高于初始值一半」且「反应速率高于0.04mol/(L*s)」两个条件同时成立时持续运行，而代码中while条件用了||（逻辑或），只要两个条件任意一个不满足，后续只要另一个条件也失效就会直接终止循环，提前了退出时机。
• 反应量计算逻辑错误
  计算各物质浓度变化时，你乘的系数是(time-0.38)（从0.38s到当前的累计时间），而非单步迭代的时间间隔0.01s，相当于把历史总反应量每轮都重复叠加一次，直接让反应物消耗速度比真实值快几十倍，这是循环仅运行到0.66s就终止的核心原因。
• 反应速率判断时机错误
  代码先判断reactionRate是否符合循环条件，再在循环体内计算当前步的新反应速率，导致判断的永远是上一轮的反应速率，最后一轮符合终止条件的速率没有被正确计入，进一步提前了循环终止时间。

内容的提问来源于stack exchange，提问作者Rish Tamhane