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