@东岳 谢谢李老师:huahua:

@d-benjamin 这个是正则表达式，是用来匹配字符串的。"."代表任意字符，"*"代表出现多次（包括0次），合起来".*"代表匹配任意字符。

"epsilon.*"可以匹配以epsilon开头的任意字符（包括本身），如epsilon，epsilonFinal等。其中epsilonFinal是用于PIMPLE算法最后一个迭代循环的控制参数。

以上规则同样适用于"U.*"和"k.*"。

不想patch被分割，可以用这个：？

preservePatches (xxx xxx)

@D-Benjamin

这是代码高亮，大片大片的代码用esc下面那个，前面点击三次，后面点击三次，把代码包围起来，你可以试试下面这个（分为三行）

``` 代码```

你在/opt/下面编译，应该没有权限。这面需要sudo，建议把你的东西拷贝到别的地方，修改Make文件夹中files：

LIB = $(FOAM_LIBBIN)/libxxx

LIB = $(FOAM_USER_LIBBIN)/libxxx

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

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 改过来也没有发生变化，据说难以保证气泡的形状精准。