CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    这网格会有什么问题

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

      怎么算一点就发散了
      Screenshot from 2022-02-22 09-30-15.png QQ图片20220223135124.png

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

        什么求解器?

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

        S 1 条回复 最后回复 回复 引用
        • S
          shrine @李东岳 最后由 编辑

          @李东岳 interFoam

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

            @shrine 把你log贴全看一下

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

            S 1 条回复 最后回复 回复 引用
            • S
              shrine @李东岳 最后由 编辑

              @李东岳 Screenshot from 2022-02-24 09-40-19.png

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

                @shrine 用代码高亮的方式把log贴全。你这个目前判断不出来。用下面这种方式,不要截图

                ```加一行空格 asdfasdf
                asdfasdf
                asdfasdf
                asdfasdf
                asdfasdf 加一行空格 ```
                

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

                S 1 条回复 最后回复 回复 引用
                • S
                  shrine @李东岳 最后由 shrine 编辑

                  @李东岳 这是最后发散的

                  
                                                  
                  
                  Courant Number mean: 0.00552061 max: 83.2287
                  Interface Courant Number mean: 0.000324483 max: 81.0174
                  deltaT = 1.18982e-53
                  Time = 3.9265
                  
                  MULES: Solving for alpha.water
                  Phase-1 volume fraction = 0.0137961  Min(alpha.water) = -5.68761e-22  Max(alpha.water) = 1
                  MULES: Solving for alpha.water
                  Phase-1 volume fraction = 0.0137961  Min(alpha.water) = -5.6877e-22  Max(alpha.water) = 1
                  DICPCG:  Solving for p_rgh, Initial residual = 0.939399, Final residual = 0.0886856, No Iterations 34
                  time step continuity errors : sum local = 4.07115e-07, global = -7.89166e-08, cumulative = 8.14176e+13
                  DICPCG:  Solving for p_rgh, Initial residual = 0.176142, Final residual = 0.0126346, No Iterations 5
                  time step continuity errors : sum local = 7.05788e-07, global = -7.10764e-08, cumulative = 8.14176e+13
                  DICPCG:  Solving for p_rgh, Initial residual = 0.0465706, Final residual = 8.69174e-07, No Iterations 141
                  time step continuity errors : sum local = 5.76874e-11, global = -1.50599e-12, cumulative = 8.14176e+13
                  smoothSolver:  Solving for epsilon, Initial residual = 4.50714e-08, Final residual = 4.50714e-08, No Iterations 0
                  smoothSolver:  Solving for k, Initial residual = 0.999977, Final residual = 2.51429e-12, No Iterations 1
                  bounding k, min: 1.72065e-45 max: 1.33449e+84 average: 4.66546e+79
                  ExecutionTime = 108.46 s  ClockTime = 109 s
                  
                  Courant Number mean: 1.25464e-05 max: 0.294671
                  Interface Courant Number mean: 5.47167e-07 max: 0.098979
                  deltaT = 1.39171e-53
                  Time = 3.9265
                  
                  MULES: Solving for alpha.water
                  Phase-1 volume fraction = 0.0137961  Min(alpha.water) = -1.12909e-21  Max(alpha.water) = 1
                  MULES: Solving for alpha.water
                  Phase-1 volume fraction = 0.0137961  Min(alpha.water) = -1.41136e-22  Max(alpha.water) = 1
                  DICPCG:  Solving for p_rgh, Initial residual = 0.185686, Final residual = 0.0180753, No Iterations 7
                  time step continuity errors : sum local = 2.01538e-06, global = 8.36086e-09, cumulative = 8.14176e+13
                  DICPCG:  Solving for p_rgh, Initial residual = 1.72997e-14, Final residual = 1.72997e-14, No Iterations 0
                  time step continuity errors : sum local = 2.60495e+41, global = 1.67368e+25, cumulative = 1.67368e+25
                  DICPCG:  Solving for p_rgh, Initial residual = 1.73034e-14, Final residual = 1.73034e-14, No Iterations 0
                  time step continuity errors : sum local = 2.60551e+41, global = 5.1378e+24, cumulative = 2.18746e+25
                  
                  S 1 条回复 最后回复 回复 引用
                  • S
                    shrine @shrine 最后由 编辑

                    @shrine 最开始,中间太多,贴不全

                    /*---------------------------------------------------------------------------*\
                      =========                 |
                      \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
                       \\    /   O peration     | Website:  https://openfoam.org
                        \\  /    A nd           | Version:  8
                         \\/     M anipulation  |
                    \*---------------------------------------------------------------------------*/
                    Build  : 8
                    Exec   : interFoam
                    Date   : Feb 24 2022
                    Time   : 10:28:19
                    Host   : "dyfluid"
                    PID    : 3889
                    I/O    : uncollated
                    Case   : /home/dyfluid/OpenFOAM/wjx/simashan10/simashan10-10
                    nProcs : 1
                    sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
                    fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
                    allowSystemOperations : Allowing user-supplied system call operations
                    
                    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
                    Create time
                    
                    Create mesh for time = 3
                    
                    
                    PIMPLE: No convergence criteria found
                    
                    
                    PIMPLE: Operating solver in transient mode with 1 outer corrector
                    PIMPLE: Operating solver in PISO mode
                    
                    
                    Reading field p_rgh
                    
                    Reading field U
                    
                    Reading/calculating face flux field phi
                    
                    Reading transportProperties
                    
                    Selecting incompressible transport model Newtonian
                    Selecting incompressible transport model Newtonian
                    Selecting turbulence model type RAS
                    Selecting RAS turbulence model kEpsilon
                    RAS
                    {
                        RASModel        kEpsilon;
                        turbulence      on;
                        printCoeffs     on;
                        Cmu             0.09;
                        C1              1.44;
                        C2              1.92;
                        C3              0;
                        sigmak          1;
                        sigmaEps        1.3;
                    }
                    
                    
                    Reading g
                    
                    Reading hRef
                    Calculating field g.h
                    
                    No MRF models present
                    
                    No finite volume options present
                    Restarting alpha
                    DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.43948e-11, No Iterations 234
                    time step continuity errors : sum local = 4.25911e-18, global = -6.53738e-21, cumulative = -6.53738e-21
                    Courant Number mean: 0.00577809 max: 0.443976
                    
                    Starting time loop
                    
                    Courant Number mean: 0.00577809 max: 0.443976
                    Interface Courant Number mean: 0.00156001 max: 0.443976
                    deltaT = 0.00487805
                    Time = 3.00488
                    
                    MULES: Solving for alpha.water
                    Phase-1 volume fraction = 0.0116292  Min(alpha.water) = -5.92388e-18  Max(alpha.water) = 1
                    MULES: Solving for alpha.water
                    Phase-1 volume fraction = 0.011635  Min(alpha.water) = -1.68355e-18  Max(alpha.water) = 1
                    DICPCG:  Solving for p_rgh, Initial residual = 0.0133921, Final residual = 0.000600949, No Iterations 2
                    time step continuity errors : sum local = 3.14837e-05, global = -1.42553e-09, cumulative = -1.42553e-09
                    DICPCG:  Solving for p_rgh, Initial residual = 0.0042361, Final residual = 0.000293054, No Iterations 2
                    time step continuity errors : sum local = 1.5383e-05, global = -1.52887e-09, cumulative = -2.95441e-09
                    DICPCG:  Solving for p_rgh, Initial residual = 0.00109826, Final residual = 9.11591e-07, No Iterations 46
                    time step continuity errors : sum local = 4.78222e-08, global = 2.49579e-09, cumulative = -4.58612e-10
                    smoothSolver:  Solving for epsilon, Initial residual = 0.00269233, Final residual = 1.17508e-07, No Iterations 2
                    smoothSolver:  Solving for k, Initial residual = 0.00185872, Final residual = 1.4954e-07, No Iterations 2
                    ExecutionTime = 0.68 s  ClockTime = 1 s
                    
                    Courant Number mean: 0.00641595 max: 0.484557
                    Interface Courant Number mean: 0.00171789 max: 0.484557
                    deltaT = 0.00502587
                    Time = 3.0099
                    
                    MULES: Solving for alpha.water
                    Phase-1 volume fraction = 0.0116409  Min(alpha.water) = -4.24191e-18  Max(alpha.water) = 1
                    MULES: Solving for alpha.water
                    Phase-1 volume fraction = 0.0116468  Min(alpha.water) = -2.19039e-18  Max(alpha.water) = 1
                    DICPCG:  Solving for p_rgh, Initial residual = 0.0141558, Final residual = 0.000625894, No Iterations 2
                    time step continuity errors : sum local = 3.48527e-05, global = 4.988e-09, cumulative = 4.52939e-09
                    DICPCG:  Solving for p_rgh, Initial residual = 0.00456781, Final residual = 0.000295155, No Iterations 2
                    time step continuity errors : sum local = 1.64706e-05, global = 5.11498e-09, cumulative = 9.64437e-09
                    DICPCG:  Solving for p_rgh, Initial residual = 0.00118887, Final residual = 9.82617e-07, No Iterations 86
                    time step continuity errors : sum local = 5.47957e-08, global = 5.40295e-11, cumulative = 9.6984e-09
                    smoothSolver:  Solving for epsilon, Initial residual = 0.00272111, Final residual = 1.30398e-07, No Iterations 2
                    smoothSolver:  Solving for k, Initial residual = 0.0018956, Final residual = 1.65103e-07, No Iterations 2
                    ExecutionTime = 0.89 s  ClockTime = 1 s
                    
                    Courant Number mean: 0.00661659 max: 0.49063
                    Interface Courant Number mean: 0.0017706 max: 0.49063
                    deltaT = 0.00510359
                    Time = 3.01501
                               
                    
                    1 条回复 最后回复 回复 引用
                    • 李东岳
                      李东岳 管理员 最后由 编辑

                      目前只能看出发散了。库朗数太大。直接发散了。你可以换结构网格试一下。应该会好很多。

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

                      S 1 条回复 最后回复 回复 引用
                      • S
                        shrine @李东岳 最后由 编辑

                        @李东岳

                        背面有圆弧,很难全部用结构网格
                        QQ图片20220226141414.jpg

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

                          你目前画的圆弧的地方就是结构网格啊 :135:
                          目前只能看出发散了

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

                          S 1 条回复 最后回复 回复 引用
                          • S
                            shrine @李东岳 最后由 shrine 编辑

                            @李东岳 把圆弧去掉,全部改成六面体网格,还是发散。
                            从细管流入水池就发散了。
                            附了log文件
                            Screenshot from 2022-03-01 09-14-34.png

                            Screenshot from 2022-03-01 09-12-07.png

                            Screenshot from 2022-03-01 09-12-19.png

                            log

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