好的，谢谢李老师，我试一下

0时间步的U，是用湍流生成方法直接生成的，可以计算得到div(U)。然后就可以计算div(phi)： U -> phi = fvc::flux(U) -> div(phi)。测试结果是，这种方式计算的div(U)和div(phi)结果基本一样的，区别是小数点后五六位。此时div(U)=div(phi)=0

那这个算法很屌。不错。

@李东岳 在 CodedFixedValue边界设置问题 中说：

http://dyfluid.com/code.html 来这查一下

好的，谢谢东岳哥

velocity.,txt
(
(0 (4 0 0))
(0.5 (4 0.0698 0))
(1 (4 0.0698 0))
(1.5 (4 0.0698 0))
(2 (4 0.0698 0))
(2.5 (4 0.0698 0))
(3 (4 0.0698 0))
(3.5 (4 0.0698 0))
(4 (4 0.1397 0))
(4.5 (4 0.2096 0))
(5 (4 0.2797 0))
(5.5 (4 0.3499 0))
(6 (4 0.4204 0))
(6.5 (4 0.4911 0))
(7 (4 0.5622 0))
(7.5 (4 0.6335 0))
(8 (4 0.7053 0))
(8.5 (4 0.7775 0))
(9 (4 0.8502 0))
(9.5 (4 0.9235 0))
(10 (4 0.9973 0))
(10.5 (4 1.0718 0))
(11 (4 0.9973 0))
(11.5 (4 0.9235 0))
(12 (4 0.8502 0))
(12.5 (4 0.7775 0))
(13 (4 0.7053 0))
(13.5 (4 0.6335 0))
(14 (4 0.5622 0))
(14.5 (4 0.4911 0))
(15 (4 0.4204 0))
(15.5 (4 0.3499 0))
(16 (4 0.2797 0))
(16.5 (4 0.2096 0))
(17 (4 0.1397 0))
(17.5 (4 0.0698 0))

)
入口设置
{
type uniformFixedValue;
uniformValue tableFile;
uniformValueCoeffs
{
file "velocity.txt";
outOfBounds repeat;
}
value $internalField;
}

流场内的模型（stl文件）已经创建好，想在模型外面创建一个矩形的外流场，类似snappyHexMesh先利用blockMesh生成矩形网格。cfmesh有这个功能吗？

好的
我试试

@dxl 您好，在你的 https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/47 这篇文章中我也遇到了相同的情况可以讨论一下吗

@学流体的小明 感谢感谢:146: :146:

如题，在一个三维计算域中，均匀提取空间中各个位置的数据，除了使用probesDict文件添加监测点外，有没有其他方法实现？

@李东岳 谢谢李老师的回复，期待:xiexie:

多谢老师！

谢谢啦

多谢老铁，$\frac{1}{\Delta t}$已更正为$\frac{\rho_\rP^{*}}{\Delta t}$

想复现溃坝一物理模型的数据，需要保持上游水头一直不变，请问各位大神这怎么操作

@麦迪文 你好，请问能参考下您的代码吗:xiexie:

@chengan-wang 怎么修正啊

@dzw05 @溯光 @李东岳 您好，我想请问一下这个interface里面的C表征的是流体对粒子一个什么样的力呀？我最近在验证波浪作用下一个粒子的输移运动，发现改变不同的C会使得粒子在水平和垂向位移上与stokes质量输移有较大的偏差

image.png

各位老师好，想请教一个问题，我现在在controlDict文件里写了一个function，里面需要读取一个字典中的值。字典中的值为一个标量，读取之后想改变这个值并重新写入字典中，这样下次读取之后就是新的值。
其中字典文件定义如下
1.png
function里面读取的代码如下
2.png
我尝试用write函数来改写，但是报错是这样
3.png
请问各位老师这个该怎么解决呢？谢谢！

@Wayne 你好，请问一下这个问题怎么解决？

@李东岳 谢谢师兄，因为模型有辐射模型，流固耦合，还有多物种传输，所以一直觉得是自己设置有问题，或者部分模型需要重新编译，忘了最基本的，刚才把所有收敛系数都调小了.可以跑了，应该就是前期不稳定。但是，也只是多跑了3步，我回头再看看