如果问题确实不稳定，则舍入误差可能足以引发现象，尽管可能需要一段时间。如果有一个小的稳定区域，使您陷入稳定状态，您可以尝试从稳定的解决方案重新启动流动，并在温度中产生小的随机或正弦扰动。我更喜欢这个而不是搞乱速度和压力，因为通常需要担心的连续性方程没有一致性。可能有类似的论点不干扰物种浓度，因为您不想担心质量的一致性。

编辑补充：如果扰动你的温度场让你担心与你的状态方程的一致性，你总是可以在你的能量方程中添加一个右手边的函数，并将扰动函数放在那里几个时间步，而不是依赖于扰动的初始条件。

添加高斯函数很烦人，因为它没有得到紧凑的支持。凹凸函数解决了这个问题，但具有需要分辨率的令人讨厌的高阶导数。我花了一点时间研究如何为您描述的目的获得更好的凹凸功能。写在 http://www.scribd.com/mobile/doc/61560918/Behavior-of-the-bump-function-s-derivatives