这个codeStream写的扰动，可以保证$\nabla\cdot\bfU=0$么

我只看到codeStream跟perturbU的代码是一致的，具体的公式在EugeneDeVilliers博士论文5.1.2章也能找到

谢谢。@学流体的小明

我的计算域没有入口出口之分，沿流动方向的两个边界都是周期性的边界条件，展向也是周期性的边界条件。所以无法回答你的问题。

我的算例也是周期边界的槽道流，其实应该是internalField初始场的影响。
前面我采用的是uniform的初始速度场，计算湍流雷诺数远达不到目标值。最近按你算例的codeStream代码加入了初始脉动速度，计算雷诺数比较接近目标值了。
以Re_tau=1000为例，我现在计算得到的雷诺数数值是910

是的，用的就是EugeneDeVilliers在他博士论文中提到的初始化槽道流流场的方法。

EugeneDeVilliers大佬的论文提到了这种湍流需要加入初始扰动才能发展，这个初始化方法也有一个 perturbU的开源代码。不过codeStream的实现更方便

@学流体的小明 感谢此贴，感谢几位的研究，解答了我的很多疑问。
并请教：槽道流LES计算的入口inflow速度场影响大吗？
因为我近期也开展了多个雷诺数Re180~5200的槽道流计算，但计算得到的壁面摩擦速度偏小（大约只有目标值的一半）。
对比了一下你的channel1000案例文件，发现一个主要的差别是初始场0/U，看起来你是采用了一种生成脉动速度的入口条件吗？

@coolhhh 楼主有可分享的结论么？我理解 x+和z+如果太大，那近壁面的网格长宽比就太大了，对LES计算还是有影响的

@李东岳 谢谢李老师。1. solver.pcorr改成PCG，虽然还是出现iteration 1000，但确实快了很多。2.尝试将pcorr的tolerance从 1e-06改成1e-05，还是会出现 iteration 1000。这个有没有可能跟网格有关呢？整体计算速度慢的原因还有P_rgh也挺费事的

@韬智tz ESI版本的OFv2012之后有一个ccmToFoam工具但默认没编译,在这个网站搜一下能找到一个编译的教程

@李东岳 李老师，借楼问一下：滑移网格计算，每三个时间步会出现pcorr iteration 1000，可能出现的原因在哪方面？网格checkMesh无问题，可能会是GAMGPCG这些格式吗？

@微凉风 我认为普通动网格应该就是会变形，运动幅值很小都发散的话可能先检查一下壁面的边界条件以及dynamic设置吧

@五好青年 抱歉我记错了，OF没有这个选项。好像这个库有个loadBalance的研究，但我测试感觉也不是很稳定，也会报错。
https://github.com/ElsevierSoftwareX/SOFTX_2018_143

印象中把各核unbalance的参数提高一点，可以起到稍微改善作用