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