CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    怎样限制流场中最大速度

    OpenFOAM
    2
    2
    446
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 小
      小考拉 最后由 编辑

      各位前辈好!

      我在compressibleInterFoam中引入了空化源项,现在调用自己写的非线性SGS模型计算的时候,流场初始速度10m/s,计算一段时间后,maxU会突然增大,导致流场发散掉,下面是速度暴增的过程

      MULES: Solving for alpha.water
      Liquid phase volume fraction = 0.997517  Min(alpha.water) = 0.000412402  Max(alpha.water) = 1
        Min(alpha.vapor) = -8.20677e-13  Max(alpha.vapor) = 0.999588
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      smoothSolver:  Solving for T, Initial residual = 4.32259e-05, Final residual = 2.38488e-09, No Iterations 2
      min(T) 290.144
      GAMG:  Solving for p_rgh, Initial residual = 4.54434e-06, Final residual = 7.05483e-10, No Iterations 3
      max(U) 13.9186
      min(p_rgh) 1843.8
      GAMGPCG:  Solving for p_rgh, Initial residual = 3.62757e-07, Final residual = 4.88638e-10, No Iterations 1
      max(U) 13.9184
      min(p_rgh) 1843.67
      ExecutionTime = 275.32 s
      
      
       Integral of alpha.water = 0.00190802  [0 3 0 0 0 0 0]
      fieldAverage fieldAverage1 output:
          Calculating averages
      
      Courant Number mean: 0.00487524 max: 0.135574
      Time = 0.011952
      
      MULES: Solving for alpha.water
      Liquid phase volume fraction = 0.997517  Min(alpha.water) = 0.000413145  Max(alpha.water) = 1
        Min(alpha.vapor) = -4.79172e-13  Max(alpha.vapor) = 0.999587
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      smoothSolver:  Solving for T, Initial residual = 4.324e-05, Final residual = 2.42709e-09, No Iterations 2
      min(T) 290.146
      GAMG:  Solving for p_rgh, Initial residual = 5.42286e-06, Final residual = 7.32723e-10, No Iterations 3
      max(U) 26.57
      min(p_rgh) 1813.11
      GAMGPCG:  Solving for p_rgh, Initial residual = 3.66982e-07, Final residual = 5.21535e-10, No Iterations 1
      max(U) 25.4251
      min(p_rgh) 1817.97
      ExecutionTime = 275.6 s
      
      
      fieldAverage fieldAverage1 output:
          Calculating averages
      
      Courant Number mean: 0.00487537 max: 0.135663
      Time = 0.011953
      
      MULES: Solving for alpha.water
      Liquid phase volume fraction = 0.997517  Min(alpha.water) = 0.000412369  Max(alpha.water) = 1
        Min(alpha.vapor) = -8.14016e-13  Max(alpha.vapor) = 0.999588
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      smoothSolver:  Solving for T, Initial residual = 4.32537e-05, Final residual = 2.41625e-09, No Iterations 2
      min(T) 290.148
      GAMG:  Solving for p_rgh, Initial residual = 5.31025e-06, Final residual = 6.64697e-10, No Iterations 3
      max(U) 20.2828
      min(p_rgh) 1797.39
      GAMGPCG:  Solving for p_rgh, Initial residual = 3.95326e-07, Final residual = 5.53382e-10, No Iterations 1
      max(U) 20.6401
      min(p_rgh) 1803.66
      ExecutionTime = 275.88 s
      
      
      fieldAverage fieldAverage1 output:
          Calculating averages
      
      Courant Number mean: 0.00487547 max: 0.135752
      Time = 0.011954
      
      MULES: Solving for alpha.water
      Liquid phase volume fraction = 0.997516  Min(alpha.water) = 0.000412357  Max(alpha.water) = 1
        Min(alpha.vapor) = -2.22555e-12  Max(alpha.vapor) = 0.999588
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      smoothSolver:  Solving for T, Initial residual = 4.32478e-05, Final residual = 2.41079e-09, No Iterations 2
      min(T) 290.15
      GAMG:  Solving for p_rgh, Initial residual = 4.91279e-06, Final residual = 6.3296e-10, No Iterations 3
      max(U) 30.2344
      min(p_rgh) 1778.9
      GAMGPCG:  Solving for p_rgh, Initial residual = 4.09173e-07, Final residual = 3.01626e-10, No Iterations 1
      max(U) 29.9168
      min(p_rgh) 1779.32
      ExecutionTime = 276.17 s
      
      
      fieldAverage fieldAverage1 output:
          Calculating averages
      
      Courant Number mean: 0.00487556 max: 0.135842
      Time = 0.011955
      
      MULES: Solving for alpha.water
      Liquid phase volume fraction = 0.997516  Min(alpha.water) = 0.000412479  Max(alpha.water) = 1
        Min(alpha.vapor) = -4.44089e-16  Max(alpha.vapor) = 0.999588
      diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
      smoothSolver:  Solving for T, Initial residual = 4.32198e-05, Final residual = 2.70013e-09, No Iterations 2
      min(T) 290.153
      GAMG:  Solving for p_rgh, Initial residual = 1.71158e-05, Final residual = 7.83321e-10, No Iterations 4
      max(U) 140.097
      min(p_rgh) 1138.55
      GAMGPCG:  Solving for p_rgh, Initial residual = 1.088e-06, Final residual = 4.13587e-10, No Iterations 2
      max(U) 147.712
      min(p_rgh) 1226.71
      ExecutionTime = 276.46 s
      

      我仔细检查发现很可能是由于求解U方程时,某个网格上面速度突然变得很大,最终导致流场发散,目前想在U方程中 或者 P方程中限制最大速度不超过某一个值,由于U是矢量,我现在不太清楚怎么限制流场中的U不超过某一个值呢?是分别限制U某一方向的速度 还是 直接限制max(mag(U).value)呢?

      麻烦哪位前辈能够指点一下呢?谢谢

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        @小考拉 maxU会突然增大只是发散的迹象,不是发散的原因,所以这种问题很难debug :135:

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

        CFD高性能服务器 http://dyfluid.com/servers.html

        1 条回复 最后回复 回复 引用
        • First post
          Last post