@东岳 好的，谢谢李东岳老师，在您的建议下我继续尝试一下，再次表示感谢！

OpenFOAM中如何给初始条件添加合适扰动呢？除了延长入口外（太耗费计算资源），有什么别的方法实现充分发展湍流初始场吗？

0_1537523541221_IMG_20180921_173609.jpg

论坛类似的问题好多。不太好回复。主要是这样：类似这种问题，只要能用函数来表示，在OpenFOAM里面植入通通不是问题。重要的是：

如何用函数表示 用函数表示之后，如何写进入

我觉得你们应该知道如何用函数表示，但如何在OpenFOAM中编程进去，这个几句话说不明白...

grad(U)和wallGrad(U)有关系么？？

从你的描述来看：

wallGradU.boundaryField()[patchi] =-U.boundaryField()[patchi].snGrad();

wallGrad(U)是面法向梯度，grad(U)是梯度。这个区别比较大。如果新版本没有了wallGrad(U)，会不会grad(U)在处理壁面的时候已经当成面法向梯度了，你研究研究这个角度对不对？

@random_ran

我给他代码高亮的时候显示成斜体了，现在改好了

@东岳 谢谢李老师: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: