CFD中文网

@xiaofei6538567 谢谢，确实没有找到办法直接编辑，最后非常窄的地方都单独分出来不加边界层了

再继续更新一下最近一周的实验结果，结果是在欧拉双方程模型中开始收敛了，但是在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

有些站友提出问题之后，其他站友们给了追问或者讨论之后，就没有消息了.....
方案能不能用，倒是回复一下啊
完全不能理解这种人提出问题就再也不管的人，不靠谱

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

@bestucan 先说声谢谢，我去仔细看看

如楼上所说，不指定就是六个自由度，重叠网格也不会出现网格变形的问题

成功了么？

计算域是否对称，可否切出1/3生成网格后旋转

如果能实现顶盖驱动MAC算法，simple也没有问题啊~可以先试试MAC，速度边界条件都是第一类边界条件，其余都为0.只有顶盖有速度。

fluent中的压力出口边界条件，资料上写着

压力出口边界条件在流场出口边界上定义静压，而静压的值仅在流场为亚音速时使用。如果出口边界上流场达到超音速，那么边界上的压力将从流场内部通过插值得到。
其他流场变量均从流场内部插值获得。

这里说的插值是怎样的插值，导数为零？还是二阶导数为零？

半新手，看其他论坛上可以实现这样的操作，即把一个复杂的几何体分成多部分分别划分网格之后再进行拟合，但是我自己试了下建块merge后导入到fluent中出错，请问该如何操作呢？不胜感激

是我应该感谢cfd中文网的各位老师，对我的帮助很大.
:xiexie: :xiexie: :xiexie: :xiexie: :xiexie: :xiexie: :xiexie: :xiexie:
另外，我测试的这个构造函数和析构函数只承担了new和delete一个double[]的任务，如果是功能更复杂的构造函数，就不好说了。

哦，这样啊，先看完这本书，再学新的。

我也是想了解一下

Sep 25, 2018 , Jan 08, 2019, Required Reviews Completed

替代文字

@东岳

just for fun:

They used to say if man could fly, he'd have wings. But he did fly. He discovered he had to. Do you wish that the first Apollo mission hadn't reached the moon, or that we hadn't gone on to Mars and then to the nearest star? That's like saying you wish that you still operated with scalpels and sewed your patients up with catgut like your great-great-great-great-grandfather used to. I'm in command. I could order this. But I'm not because, Doctor McCoy is right in pointing out the enormous danger potential in any contact with life and intelligence as fantastically advanced as this. But I must point out that the possibilities, the potential for knowledge and advancement is equally great. Risk. Risk is our business. That's what the CFD is all about. That's why we're aboard her. You may dissent without prejudice. Do I hear a negative vote? (silence around the table) Engineer, solving that imcompressible N.S. equations using PISO algorithm.

捕获.JPG

重写了一个，以后用这个更新

用的Adobe AI

谢谢李老师的回复！因为雷诺应力不仅出现在动量方程里面，也出现在turbulent kinetic energy equation 与 turbulent frequency equation中，我的想法是如果我分别修改湍流模型的文件（kOmegaSSTBase）与动量方程的文件（linearEddyViscosity）就会要定义两次非线性项，所以我选择就在kOmegaSSTBase中更新turbulence->divDevRhoReff(U),在其中加入非线性项，参考的资料是OpenFoam里面的nonlineareEddyViscosity。由于这样修改需要在kOmegaSSTBase.C文件的constructors中加入eddyViscosity 也就是上图中蓝色后面的部分，我直接在原来的代码Turbulence后面加这一项会报错，想请问这怎么解决？非常感谢李老师耐心看我的问题，不知道我表述清楚没有？

直接相乘即可： $\mathbf S_b \cdot (\phi \mathbf u)_b$