这网格会有什么问题

[shrine]

怎么算一点就发散了

[李东岳]

什么求解器？

[shrine]

@李东岳 interFoam

[李东岳]

@shrine 把你log贴全看一下

[shrine]

@李东岳

[李东岳]

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

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

[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

[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
           

[李东岳]

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

[shrine]

@李东岳

背面有圆弧，很难全部用结构网格

[李东岳]

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

[shrine]

@李东岳 把圆弧去掉，全部改成六面体网格，还是发散。
从细管流入水池就发散了。
附了log文件

log