我在昨天发现了这个错误，怀疑也是内存出的问题。仔细回想了一下，昨天在程序里面写了几行：

scalar* zeta = new double [2*n];

后面忘了删除。昨晚上计算的时候三个算例全部killed：

diagonal: Solving for M11, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for M12, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for M13, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for M14, Initial residual = 0, Final residual = 0, No Iterations 0 PIMPLE: iteration 1 GAMG: Solving for p, Initial residual = 0.0004130904583986163, Final residual = 2.405751520795531e-06, No Iterations 4 time step continuity errors : sum local = 8.330982304837603e-10, global = 1.267524541545882e-12, cumulative = -2.029919160042075e-06 Killed

无任何消息。然后我在源代码中将这个new出来的数组delete之后，就不会出错了。

另外，我测试了一下，不删除数组计算的时候内存会一直上升。因此，动态数组记得释放！

@东岳

好的，感谢东岳教授的及时回复。

若有新进展，再来更新！

Peng提到在空化区域的壁面应用slip wall boundary，并提到这与超音速流类似
MODELLING INJECTOR FLOW INCLUDING CAVITATION EFFECTS FOR DIESEL APPLICATIONS
为什么超音速流要用slip wall boundary，有什么具体说法吗？

各位大佬，我做了一个 T 形管的稳态模拟，该T形管左端进口速度2m/s，下端进口速度1m/s，右端为outlet，用的是simpleFoam，在我计算了10多步，就出现问题了，发现bounding k和bounding eplison的值特别大，我猜想这是出错的原因，，但是我不知道自己该如何进行弥补，因此我附上了0文件夹下的k、eplison、p、U、nut以及log文件的截图，求大佬们帮助一下，谢谢！0_1537274990785_k.PNG 0_1537274998945_eplison.PNG 0_1537275008082_p.PNG 0_1537275013088_U.PNG 0_1537275020063_nut.PNG 0_1537275025878_log.PNG

想在某些特定的时刻赋予某些物理场，比如速度场压力场等一个初值，请问除了在特定时刻手动setFields之外，有没有办法能实现自动周期性初始化，多次尝试失败备受困扰，希望有大神和老师能帮忙解答一下，跪谢

看一下Versteeg那本书TVD章节，TVD里面一些格式需要制定梯度格式，迎风，中心这种的特别简单不需要梯度，因此不需要指定 :wanan:

那个地址是运行程序临时分配的地址，也不是完整的物理地址，应该是某个段地址的偏移地址。如果要看，得用gdb调试，设置断点再看。
除非是相当底层的调试，一般用不到这种调试方法。

这是个ddtFoam()函数引起的浮点溢出错误，就排查排查程序参数有没有不合理的。

@东岳 参照此方法已解决问题！感谢！

@izumi 谢谢您！^_^

@dzw05 嗯，谢谢你，这种方法的确比较简单，但我更想用上面那种方法来实现，因为之后会加非稳态风场等等，后续会比较方便一点:xinlei:

@dzw05 谢谢！！

静候大作！

用旧版本的paraview可以，我现在用新版的也有这个问题。paraFoam确实慢。
没有解决办法，只能心中自有云图咯 :mihu:

耦合求解和分离求解的结果有什么不同呢？

俩种求解策略。高超那面比较用户耦合求解。不可压缩低速这面SIMPLE PISO近似垄断。这不好说哪个结果好哪个结果不好。你看SCI，发表出来的，结果都很好。算法上的区别，耦合求解构建一个大矩阵，压力速度放在一起求解，一锅粥出来了。分离求解分开构建矩阵，比如压力方程、速度方程。

@东岳 谢谢李老师！

@linhan-ge 改过来也没有发生变化，据说难以保证气泡的形状精准。

附上温度限制：

temperatureLimit { type limitTemperature; active true; limitTemperatureCoeffs { selectionMode all; min <Tmin>; max <Tmax>; } }

div改成upwind有所改善，但不能根本上解决问题。
如果温度使用uniformFixedValue将压力从小到大随时间增大，为何实际算的时候也不起作用？
初始压力极小如绝对压力为1，即使温度也很低，则能正常计算，压力和温度同事增大或减小才可以计算，这是什么原因？

谢谢各位指教

@东岳 李老师好，目前发散问题已解决，但是由于epsilon会偏离正常值范围，所以会导致流场混乱，且层流算出来的流场和湍流k-e模型算出来的基本相似，尝试用低雷诺数模型SA计算，则无法观察到流场中的涡。

对于场，IOobject类型里面如果是AUTO_WRITE就是输出。对于颗粒，目前不确定是不是IOobject ..