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 ..

@东岳 谢谢东岳老师，本来我的solver是基于OpenFOAM4.0版本，根据您的建议，我换用了OpenFOAM-extend4.0来编写solver，成功通过编译。
如您所说，错误确实很有可能是因为粘弹性的库不匹配OF的库函数。
谢谢东岳老师的提示。

Zalesak disk可以参考isoAdvector里面的算例配置：
https://github.com/isoAdvector/isoAdvector/tree/master/OpenFOAM/run/prescribedU/notchedDiscInSolidBodyRotation