CFD中文网

利用DPM_SCALAR_UPDATA

@李东岳 对 有

成功了么？

@lsprxd

$h=\frac{q}{T_w - T_b}$

其中：$q$为热流密度，$T_w$为壁面温度，$T_b$为参考温度

这样可以吗？

回复: fluent源项udf一个吸附组分
fluent表面化学反应中，固体组分分成site species和solid species，我想用udf源项一个吸附成分，随着反应进行被消耗。为什么会想用udf呢，是因为多孔介质表面化学反应的内置反应速率不符合自己的需求，想自己写一个反应速率的udf，无奈一直得不到结果，搞不清楚多孔介质中比表面积的使用，所以想直接用udf定义组分。

在自己的笔记本上跑算例的时候，如果打开Chrome，经常由于内存不足导致网页自动崩溃。回国赶紧换工作站。

我觉得可能是win10系统占用的内存更高了，或者是新版本Chrome内存占用过高。

亏！！

啊哦，第一次发帖好像还把格式弄错了非常抱歉，，
图在这里
1.png

另外请教一下如何发文件，想发例子的压缩包的貌似也没成功

我的16.04也有这个问题

欧拉欧拉那面由于已经在统计学方面已经处理过了，没有噪音。欧拉拉格朗日模拟这面目前有什么文章在处理么？找了找，还没看到。

As an alternative, Lagrangian methods ‘‘discretize” the density function into ‘‘parcels” that are simulated using Monte-Carlo methods. While quite accurate, as in any statistical approach, Lagrangian methods require a relatively large number of parcels to control statistical noise, and thus are computationally expensive.

@mohui 那我感觉应该可以用，我只是知道，没用过，期待你的结果

再继续更新一下最近一周的实验结果，结果是在欧拉双方程模型中开始收敛了，但是在mixture中没法收敛，不知道是什么原因，两个模型的设置都是相同的，离散格式都还是保持的默认一阶格式。
前人文献中应该都是先用mixture+SIMPLEC迭代收敛+欧拉双方程，我也不明白到底哪里错了，会导致发散，从几百万到最高6800万网格都试过了，还是发散。
欧拉双方程模型前20次迭代残差

iter continuity u-air u-sand v-air v-sand w-air w-sand energy-p1 energy-p2 k epsilon vf-sand time/iter 1 1.0000e+00 3.6056e+06 8.4608e+05 3.4192e-09 1.6311e+00 3.4671e-09 1.2974e-08 8.9865e-17 9.1745e-17 3.4616e-02 3.2507e-02 2.1808e-04 0:10:08 19 reversed flow in 12 faces on pressure-outlet 6. 2 1.0000e+00 5.3943e-01 1.9935e+00 1.7837e-01 3.8110e-01 1.9575e-01 3.9505e-01 1.0126e-16 1.0320e-16 7.5627e-03 5.1711e-02 1.6125e-04 0:09:40 18 reversed flow in 2 faces on pressure-outlet 6. 3 1.0000e+00 3.6569e-01 1.7316e-01 9.2265e-02 3.8210e-02 9.7483e-02 4.0064e-02 9.9739e-17 1.0429e-16 3.3266e-03 1.8913e-01 1.8913e-04 0:08:43 17 reversed flow in 2 faces on pressure-outlet 6. 4 1.0000e+00 2.3734e-01 2.1136e-02 4.8203e-02 4.3756e-03 5.0243e-02 4.5314e-03 1.0082e-16 1.0521e-16 9.3765e-04 2.7619e-02 2.7003e-04 0:07:51 16 reversed flow in 5 faces on pressure-outlet 6. 5 8.3879e-01 4.5351e-02 7.4916e-02 1.0088e-02 1.2337e-02 1.0339e-02 1.2760e-02 1.0012e-16 1.0450e-16 2.5818e-04 7.2360e-03 1.6935e-04 0:07:11 15 reversed flow in 2 faces on pressure-outlet 6. 6 7.2566e-01 1.5783e-01 6.2558e-02 2.7641e-02 1.0158e-02 2.8694e-02 1.0502e-02 1.0044e-16 1.0505e-16 1.1932e-04 2.0689e-03 8.0813e-05 0:06:32 14 reversed flow in 1 faces on pressure-outlet 6. 7 6.8356e-01 1.2100e-01 2.3435e-02 1.9995e-02 3.6904e-03 2.0887e-02 3.8252e-03 9.9881e-17 1.0475e-16 6.7527e-05 3.6129e-04 5.6476e-05 0:05:59 13 8 6.1018e-01 8.9233e-02 9.0639e-03 1.4191e-02 1.3958e-03 1.4873e-02 1.4486e-03 1.0016e-16 1.0524e-16 5.4561e-05 8.1376e-05 4.5390e-05 0:05:27 12 9 5.3642e-01 7.4696e-02 9.2564e-03 1.1860e-02 1.3618e-03 1.2386e-02 1.4106e-03 9.9675e-17 1.0550e-16 6.9822e-05 7.6546e-05 1.4094e-04 0:04:51 11 10 4.2423e-01 1.3710e-02 2.2543e-02 7.0280e-03 3.2282e-03 7.1669e-03 3.3568e-03 1.0002e-16 1.0568e-16 5.0774e-05 2.5400e-05 1.0562e-04 0:04:21 10 11 3.4934e-01 4.3842e-02 1.4440e-02 8.0368e-03 2.0472e-03 8.1755e-03 2.1072e-03 1.0028e-16 1.0603e-16 4.8257e-05 3.0984e-05 6.5681e-05 0:03:53 9 iter continuity u-air u-sand v-air v-sand w-air w-sand energy-p1 energy-p2 k epsilon vf-sand time/iter 12 3.1015e-01 4.6611e-02 1.2895e-02 7.6343e-03 1.8324e-03 7.9143e-03 1.9135e-03 1.0023e-16 1.0559e-16 2.7174e-05 2.0044e-05 5.6454e-05 0:03:24 8 13 2.5955e-01 1.9970e-02 5.5753e-03 5.2928e-03 8.0857e-04 5.3187e-03 8.3195e-04 9.9933e-17 1.0536e-16 2.5433e-05 1.9280e-05 5.0605e-05 0:02:57 7 14 2.2851e-01 3.3273e-02 9.6663e-03 5.7408e-03 1.3454e-03 5.8779e-03 1.3995e-03 1.0000e-16 1.0573e-16 2.0854e-05 2.5397e-05 5.6156e-05 0:02:30 6 15 1.8449e-01 7.5023e-03 8.5500e-03 4.3649e-03 1.1736e-03 4.4746e-03 1.2231e-03 9.9658e-17 1.0545e-16 1.6605e-05 1.4800e-05 3.6116e-05 0:02:05 5 16 1.6138e-01 2.4224e-02 9.2798e-03 4.5051e-03 1.2740e-03 4.6164e-03 1.3327e-03 9.9948e-17 1.0539e-16 9.6824e-06 8.5184e-06 2.3386e-05 0:01:39 4 17 1.3852e-01 1.1909e-02 1.8647e-03 3.5484e-03 2.8634e-04 3.5991e-03 2.9363e-04 9.9653e-17 1.0552e-16 1.0118e-05 8.6850e-06 1.2556e-05 0:01:14 3 18 1.2608e-01 1.6858e-02 4.4877e-03 3.4879e-03 6.1594e-04 3.5582e-03 6.4185e-04 9.9543e-17 1.0467e-16 6.7772e-06 7.3624e-06 1.1645e-05 0:00:49 2 Interrupting at the end of timestep. 19 1.0874e-01 8.9217e-03 1.4100e-03 2.9888e-03 2.1968e-04 3.0423e-03 2.2545e-04 9.9718e-17 1.0514e-16 8.4741e-06 8.0399e-06 8.6141e-06 0:00:25 1 20 1.0132e-01 1.4358e-02 4.1734e-03 3.0188e-03 5.6793e-04 3.0721e-03 5.9320e-04 9.9335e-17 1.0506e-16 6.2431e-06 5.2288e-06 4.6646e-06 0:00:00 0

Mixture模型的残差图……每一次都是这样的:zoule:
Snipaste_2019-06-12_10-08-39.png

@李东岳 我最后还是把这个打开了。我模拟的是潮流能水轮机，开之前功率系数不过30%出头，开了之后到了45%左右（TSR=6），功率系数-叶尖速比曲线基本上与实验对上了。而且开了之后对时间步长就不敏感了，时间步从0.0001s到0.002s，CoN从40多一直到700多，功率系数变化在5%以内。

而且，OpenFOAM/tutorials/multiphase/interPhaseChangeDyMFoam/propeller这个算例也是开启了动量预测，我考虑是不是固定转速的旋转机械都需要开这个选项。

最近在用ICEM生成非结构化网格，网格数量大、八叉树生成方法慢，且开始时几何总有各种问题，这个过程中ICEM信息窗格一直在报告进程。在发现问题和解决问题的过程中，我才意识到信息窗格的重要性。想更多了解信息窗格的阅读tips，但没有找到相关资料，希望在这里能和大家交流。相关问题大家踊跃跟楼~hh

现在我有个问题是，信息窗格字体的不同颜色分别代表什么？比如蓝色字体代表什么？

Sorry, Chinese input missing....
0_1476947519913_oie_yyxHcXChIqBJ.png

0_1474848459010_捕获.JPG

充值100送100

这个是哪国人的名字哟

大都忽略重力源项，哪位大神给解释下为什么忽略重力源项

重力对流体的作用就是引起水动力学压力，如果是密闭的流体，其不会引起流动变化。因此可以忽略重力项，这样压力边界条件也可以好给出。

一般考虑重力的主要是重力为一种重要的体积力，如密度不均引起的浮力，其中密度不均进一步由温度导致。

还有一种情况要考虑重力是模拟流体的自由落体。

其他情况可能跟具体的应用有关，如旋转系、海洋热风那种普遍都添加重力。

如果不忽略，怎样恰当的处理？

如果考虑重力，通常求解的压力不是先前的压力p，而是压力p减去水力学压力$\rho g h$。从算法层面理解，重力引致下落，压力对这种效果抵消，采用这种压力计算的时候，可以省去这种来回下落-抵消的求解过程。求解出$p-\rho gh$之后，真实的压力即为$p-\rho gh+\rho gh$

@bishuiyouyou 是指允许的最大体积分数，数值上不会再大了。不影响流动。