Fluent风沙两相流求解发散



  • 模型长度100m,高10m,宽1.8,对数速度入口、压力出口,中间有个挡板,边界层网格20层,第一层网格高0.0001m,网格数400wimage_2.png
    挡板的平面结构:
    2da8c93c-20e9-4b0a-adb2-5ac5d1529c8f-image.png
    单相流收敛之后换成双欧拉,发散。d97078dd-5c98-4ceb-9430-e6d2fd529551-image.png
    Mixture和双欧拉模型都没法收敛,时间步长0.001,迭代50次
    1489d110-9a20-44b3-9517-410bbc1fa9c8-image.png
    请问我该怎么调整算例和网格呢?



  • 你这个直接结构网格搞起啊,为何用四面体呢



  • 那个沙障周围结构化太难了吧,得分好多块啊



  • @东岳 之前用过ICEM尝试过,最后感觉分块太复杂了,然后就转成非结构了,混合网格的话我没试过,我估计还是沙障周围和后边一段网格达不到计算要求,沙障后边会出现很长的低速区,再细化网格电脑就没法算了……
    单相流收敛情况现在还行,基本上满足计算要求了,两相流特别容易发散。
    最近上了超算,那个超算也只给我提供32G内存的虚拟机,最多2000W网格左右。



  • @veenFluent风沙两相流求解发散 中说:

    那个沙障周围结构化太难了吧,得分好多块啊

    把这块附近的网格截图看看



  • @东岳 当前算例中网格是这样的:
    e3227c49-ad6b-4255-aae6-33bb337a8551-image.png
    里边截面和这张类似:
    8bcf6432-1b64-4f87-86e9-1ac3dcbbe539-image.png
    还尝试过这种局部加密,但是从单相流的情况来看这个网格还不太合理
    546d1393-57ab-4c85-90e6-8d259a5e5653-image.png
    网格的Element Quality平均大概在0.33左右,是不是太低了?扭曲度平均0.22,从扭曲度来看还可以接受。



  • 又看了看网格质量,质量差的网格竟然都是边界层的网格,那我可能考虑不要边界层了:zoule:地面和沙障周围网格质量都特别差 359d1044-68f3-4afb-853c-66694ea16377-image.png



  • 取消了边界层之后加密了一下网格(加密区域如下,网格尺寸0.1m),两相流还是没有成功收敛,不过现在单相流计算收敛性很好,湍粘度过大和回流都没有出现,看来还得继续调整调整,记得之前算过一个不准确的两相流,调整了一周:zoule: image_6.png



  • :huahua: 网格关系很大 继续搞 正在往收敛的步骤走 如果要发SCI 画一星期搞个结构网格也是值得的!



  • @东岳 好的,谢谢东岳老师!刚才又想了想,结构化网格也是能画的:xinxin2:



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


Log in to reply