请教关于进口条件的设置的问题

• 我需要给定进口总温总压，进口气流速度角度，以及出口静压请问，这么给是否有问题，总是挂。1_1490324845490_U 0_1490324845490_p

• 你要很小心设置压力进出口边界条件，从理论来讲这并不是最稳定的边界条件，当然是可以算的并且可以算正确。你用的哪个求解器？网格类型是什么？网格质量如何？把边界条件贴一下？

• @cfd-china 您好，我算的马赫数大约0.6，是一个平面的压气机叶栅，目前用的是simpleFoam，网格用pointwise画的，检查质量ok，但总是挂掉，K，Omega总是超级大。

• @cfd-china 用的不可压，只是想尝试算个结果，曾经在fluent中算过，用的komega-sst，湍流模型，看曲线k和omega特别大，达到了10e35量级，2_1490340702077_omega 1_1490340702077_k 0_1490340702077_epsilon 错误信息0_1490340438293_PyFoamRunner.simpleFoam.logfile

1. 你这个网格哪里调用了交界面边界？抱歉看到了。
2. 结果从第一步就没有保证连续性
3. 将网格贴图（而非显示面）
4. chechMesh结果
5. fvSolution
6. fvScheme

注意是贴哦不是上传附近如：

Time = 2001

GAMG:  Solving for Ux, Initial residual = 0.00638252, Final residual = 0.000400813, No Iterations 1
GAMG:  Solving for Uy, Initial residual = 0.00114759, Final residual = 5.63775e-05, No Iterations 1
GAMG:  Solving for Uz, Initial residual = 0.0283682, Final residual = 0.000925109, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.867854, Final residual = 0.0670769, No Iterations 4
time step continuity errors : sum local = 43.0388, global = -0.914839, cumulative = -0.914839
GAMG:  Solving for omega, Initial residual = 0.0563256, Final residual = 0.00245155, No Iterations 1
GAMG:  Solving for k, Initial residual = 1, Final residual = 0.0431817, No Iterations 1
ExecutionTime = 32.24 s  ClockTime = 33 s

Time = 2002

GAMG:  Solving for Ux, Initial residual = 0.00830621, Final residual = 0.000239452, No Iterations 1
GAMG:  Solving for Uy, Initial residual = 0.00165236, Final residual = 4.55896e-05, No Iterations 1
GAMG:  Solving for Uz, Initial residual = 0.0188184, Final residual = 0.000526957, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.837434, Final residual = 0.0711894, No Iterations 2
time step continuity errors : sum local = 148.659, global = -0.294784, cumulative = -1.20962
GAMG:  Solving for omega, Initial residual = 0.0143057, Final residual = 8.19713e-06, No Iterations 1
bounding omega, min: -2.58031e+07 max: 3.30734e+11 average: 4.70801e+09
GAMG:  Solving for k, Initial residual = 0.11034, Final residual = 0.000207143, No Iterations 1
ExecutionTime = 39.79 s  ClockTime = 40 s

Time = 2003

GAMG:  Solving for Ux, Initial residual = 0.00392685, Final residual = 7.6632e-05, No Iterations 1
GAMG:  Solving for Uy, Initial residual = 0.00110286, Final residual = 1.82308e-05, No Iterations 1
GAMG:  Solving for Uz, Initial residual = 0.0143816, Final residual = 0.000210414, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.791161, Final residual = 0.0685184, No Iterations 2
time step continuity errors : sum local = 111.167, global = -0.185086, cumulative = -1.39471
GAMG:  Solving for omega, Initial residual = 1, Final residual = 0.0313884, No Iterations 3
bounding omega, min: -1.11841e+23 max: 5.89039e+26 average: 3.21879e+22
GAMG:  Solving for k, Initial residual = 0.999599, Final residual = 4.49893e-17, No Iterations 1
ExecutionTime = 47.54 s  ClockTime = 48 s

Time = 2004

GAMG:  Solving for Ux, Initial residual = 0.0028955, Final residual = 7.92784e-05, No Iterations 1
GAMG:  Solving for Uy, Initial residual = 0.0009558, Final residual = 2.10947e-05, No Iterations 1
GAMG:  Solving for Uz, Initial residual = 0.0218884, Final residual = 0.000373211, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.777876, Final residual = 0.065162, No Iterations 3
time step continuity errors : sum local = 78.0342, global = 0.432795, cumulative = -0.961914
GAMG:  Solving for omega, Initial residual = 0.538686, Final residual = 3.00162e-09, No Iterations 1
bounding omega, min: -9.64586e+19 max: 1.32289e+42 average: 7.89521e+35
GAMG:  Solving for k, Initial residual = 0.178897, Final residual = 1.16811e-20, No Iterations 1
ExecutionTime = 54.91 s  ClockTime = 56 s

Time = 2005


另外，这是一个极为简单的几何，建议采用全六面体。

• | =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-119cac7e8750
Exec   : checkMesh
Date   : Mar 23 2017
Time   : 23:51:13
Host   : "ubuntu"
PID    : 41427
Case   : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/nasa65
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points:           1876500
faces:            5504237
internal faces:   5380555
cells:            1814132
faces per cell:   6
boundary patches: 8
point zones:      0
face zones:       0
cell zones:       0

Overall number of cells of each type:
hexahedra:     1814132
prisms:        0
wedges:        0
pyramids:      0
tet wedges:    0
tetrahedra:    0
polyhedra:     0

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch               Faces    Points   Surface topology
hub                 30748    31275    ok (non-closed singly connected)
inlet               5546     5700     ok (non-closed singly connected)
outlet              4543     4680     ok (non-closed singly connected)
per1                15340    15660    ok (non-closed singly connected)
per2                15340    15660    ok (non-closed singly connected)
ps                  10266    10500    ok (non-closed singly connected)
ss                  11151    11400    ok (non-closed singly connected)
sym                 30748    31275    ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-0.0925 -0.0130406 -1.48492e-15) (0.13 0.0278965 0.02)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (-1.04226e-16 4.59207e-16 4.83911e-16) OK.
***High aspect ratio cells found, Max aspect ratio: 10084.9, number of cells 119355
<<Writing 119355 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 9.85732e-14. Maximum face area = 3.49678e-06.  Face area magnitudes OK.
Min volume = 9.85915e-19. Max volume = 2.44324e-09.  Total volume = 9.66559e-05.  Cell volumes OK.
Mesh non-orthogonality Max: 54.3818 average: 6.13782
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.752069 OK.
Coupled point location match (average 0) OK.

Failed 1 mesh checks.

• | =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       dictionary;
object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p
{
solver          GAMG;
smoother        GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator    faceAreaPair;
mergeLevels     1;

tolerance       1e-06;
relTol          0.05;
}

pFinal
{
$p; tolerance 1e-06; relTol 0; } "(U|k|epsilon)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-05; relTol 0.1; } "(U|k|epsilon)Final" {$U;
tolerance       1e-05;
relTol          0;
}
}

PIMPLE
{
nOuterCorrectors 4;
nCorrectors     1;
nNonOrthogonalCorrectors 0;
pRefCell        0;
pRefValue       0;
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
residualControl
{
p               1e-2;
U               1e-3;
"(k|epsilon)"   1e-3;
}
}

relaxationFactors
{
fields
{
p               0.3;
}
equations
{
U               0.7;
k               0.7;
"epsilon.*"     0.7;
}
}

cache
{
}
// ************************************************************************* //

• | =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       dictionary;
object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
}

{
default         Gauss linear;
}

divSchemes
{
efault         none;
div(phi,U)      bounded Gauss limitedLinearV 1;
div(phi,k)      bounded Gauss limitedLinear 1;
div(phi,epsilon) bounded Gauss limitedLinear 1;
div(phi,R)      bounded Gauss limitedLinear 1;
div(R)          Gauss linear;
div(phi,nuTilda) bounded Gauss limitedLinear 1;
}

laplacianSchemes
{
default         Gauss linear corrected;
}

interpolationSchemes
{
default         linear;
}

{
default         corrected;
}

wallDist
{
method meshWave;
}

// ************************************************************************* //

• 谢谢。
进出口边界条件？

• @cfd-china

Failed 1 mesh checks.
这个是不是网格不好的意思

• 这不是个很重要的报错，可以容忍。

• | =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       volScalarField;
object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 122825;

boundaryField
{
hub
{
}
inlet
{
type		totalPressure;
p0 		uniform 132700;
value 		uniform 132700;
gamma 		1.4;
}
outlet
{

type            fixedValue;
value           uniform 111500;

}
per1
{
type            cyclicAMI;

}
per2
{
type             cyclicAMI;

}
ps
{
}
ss
{
}
sym
{
type symmetryPlane;

}

}

• @cfd-china 您看有什么问题吗

| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       volVectorField;
object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (126 113 0);

boundaryField
{

hub
{
type            fixedValue;
value           uniform (0 0 0);

}
inlet
{
type            pressureDirectedInletVelocity;
inletDirection 	uniform	(0.743 0.669 0);
value           uniform (0 0 0);
}
outlet
{

}
per1
{
type            cyclicAMI;

}
per2
{
type            cyclicAMI;

}
ps
{
type            fixedValue;
value           uniform (0 0 0);
}
ss
{
type            fixedValue;
value           uniform (0 0 0);
}
sym
{
type symmetryPlane;
}
}

// ************************************************************************* //

• | =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       volScalarField;
location    "0";
object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 194.6;

boundaryField
{
hub
{
type            kqRWallFunction;
value           uniform 194.6;
}
inlet
{
type            turbulentIntensityKineticEnergyInlet;
intensity       0.05;
U               U;
phi             phi;
value           uniform 194.6;
}
outlet
{
}
per1
{
type            cyclicAMI;

}
per2
{
type            cyclicAMI;

}
ps
{
type            kqRWallFunction;
value           uniform 194.6;
}
ss
{
type            kqRWallFunction;
value           uniform 194.6;
}
sym
{
type symmetryPlane;
}

}

// ************************************************************************* //

• | =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       volScalarField;
object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

omegaInlet      630.0;

dimensions      [0 0 -1 0 0 0 0];

internalField   uniform $omegaInlet; boundaryField { hub { type omegaWallFunction; value$internalField;
}
inlet
{
type            turbulentMixingLengthFrequencyInlet;
mixingLength    0.025;
k               k;
value           $internalField; } outlet { type zeroGradient; } per1 { type cyclicAMI; } per2 { type cyclicAMI; } ps { type omegaWallFunction; value$internalField;
}
ss
{
type            omegaWallFunction;
value           $internalField; } sym { type symmetryPlane; } } // ************************************************************************* // • 我注意到你的log 文件里面从第二个时间步长开始， specific dissipation rate (omega)就已经发散了(10^9)。 omega 文件里 FoamFile { version 2.0; format ascii; class volScalarField; object omega; } 你的 internalField uniform$omegaInlet = 630 这个神奇的数字是怎么算的？
从你的log file 看出，你是从2000 s 开始算的，贴出2000s的初始条件或许能帮助解决问题？

对于你的的进，出口边界条件不是特别理解。用进口速度边界条件，出口用压力边界条件不能满足你的需求吗？气动问题大家通常都是用这个方法。

• @random_ran 我开始没有开湍流模型，想算个初场，2000部之前都是没有湍流模型的，omega用这个公式算的
omega = k^0.5/(C^0.25 * l) c=0.09
我的是进口条件是进口已知总压，和气流方向出口知道静压，算出来速度大约是0.67Ma

• 这是最新算的几部，麻烦老师们给看一下

| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-119cac7e8750
Exec   : simpleFoam
Date   : Mar 25 2017
Time   : 02:33:53
Host   : "ubuntu"
PID    : 44090
Case   : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/nasa65
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

SIMPLE: convergence criteria
field U	 tolerance 1e-05
field p	 tolerance 1e-05
field "(k|epsilon|omega)"	 tolerance 1e-05

AMI: Creating addressing and weights between 15340 source faces and 15340 target faces
AMI: Patch source sum(weights) min/max/average = 0.999911, 1.00059, 1.00001
AMI: Patch target sum(weights) min/max/average = 0.999839, 1.00646, 1.00006

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
kOmegaSSTCoeffs
{
alphaK1         0.85;
alphaK2         1;
alphaOmega1     0.5;
alphaOmega2     0.856;
gamma1          0.555556;
gamma2          0.44;
beta1           0.075;
beta2           0.0828;
betaStar        0.09;
a1              0.31;
b1              1;
c1              10;
F3              false;
}

No MRF models present

No finite volume options present

Starting time loop

Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0813711, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0714855, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0771402, No Iterations 4
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0876464, No Iterations 5
time step continuity errors : sum local = 10.8537, global = 4.01707, cumulative = 4.01707
smoothSolver:  Solving for omega, Initial residual = 0.259035, Final residual = 0.0221475, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.0890898, No Iterations 5
ExecutionTime = 34.31 s  ClockTime = 36 s

Time = 2

smoothSolver:  Solving for Ux, Initial residual = 0.886786, Final residual = 0.058124, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.979345, Final residual = 0.0631154, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.810242, Final residual = 0.0670691, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.986904, Final residual = 0.0565546, No Iterations 4
time step continuity errors : sum local = 1345.89, global = -10.0278, cumulative = -6.0107
smoothSolver:  Solving for omega, Initial residual = 0.0850778, Final residual = 0.00220061, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.814193, Final residual = 0.0238861, No Iterations 2
ExecutionTime = 45.36 s  ClockTime = 47 s

Time = 3

smoothSolver:  Solving for Ux, Initial residual = 0.0193303, Final residual = 0.00152468, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.011814, Final residual = 0.00102405, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.422504, Final residual = 0.0419746, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.791662, Final residual = 0.0642592, No Iterations 4
time step continuity errors : sum local = 957.279, global = -143.421, cumulative = -149.432
smoothSolver:  Solving for omega, Initial residual = 0.00274129, Final residual = 4.02132e-05, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.169936, Final residual = 0.0157838, No Iterations 1
ExecutionTime = 56.05 s  ClockTime = 58 s

Time = 4

smoothSolver:  Solving for Ux, Initial residual = 0.0983965, Final residual = 0.00679746, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0657851, Final residual = 0.00457164, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.206983, Final residual = 0.0189299, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.980449, Final residual = 0.0634521, No Iterations 4
time step continuity errors : sum local = 1261.26, global = -92.2404, cumulative = -241.672
smoothSolver:  Solving for omega, Initial residual = 0.000609973, Final residual = 3.54731e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.10451, Final residual = 0.0098634, No Iterations 1
ExecutionTime = 66.79 s  ClockTime = 69 s

Time = 5

smoothSolver:  Solving for Ux, Initial residual = 0.0261376, Final residual = 0.00214716, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0187723, Final residual = 0.00145964, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.0626224, Final residual = 0.00595206, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.921146, Final residual = 0.0735434, No Iterations 4
time step continuity errors : sum local = 1305.03, global = -174.17, cumulative = -415.842
smoothSolver:  Solving for omega, Initial residual = 0.000273909, Final residual = 9.03935e-06, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0349887, Final residual = 0.00218001, No Iterations 3
ExecutionTime = 78.4 s  ClockTime = 81 s

Time = 6

smoothSolver:  Solving for Ux, Initial residual = 0.0566235, Final residual = 0.00458913, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0400099, Final residual = 0.00302616, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.0384566, Final residual = 0.00324406, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.916166, Final residual = 0.0603851, No Iterations 4
time step continuity errors : sum local = 668.506, global = -94.6009, cumulative = -510.443
smoothSolver:  Solving for omega, Initial residual = 0.000292328, Final residual = 1.62378e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.046623, Final residual = 0.00262596, No Iterations 2
ExecutionTime = 89.53 s  ClockTime = 92 s

Time = 7

smoothSolver:  Solving for Ux, Initial residual = 0.233958, Final residual = 0.0176323, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.128911, Final residual = 0.00936957, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.0442786, Final residual = 0.00376373, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.952315, Final residual = 0.0728294, No Iterations 4
time step continuity errors : sum local = 1263.9, global = -195.294, cumulative = -705.737
smoothSolver:  Solving for omega, Initial residual = 0.000206228, Final residual = 1.22213e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0167436, Final residual = 0.00109875, No Iterations 4
ExecutionTime = 100.52 s  ClockTime = 103 s

Time = 8

smoothSolver:  Solving for Ux, Initial residual = 0.064942, Final residual = 0.00538914, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0316677, Final residual = 0.00284742, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.0769102, Final residual = 0.00692843, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.791177, Final residual = 0.0625584, No Iterations 6
time step continuity errors : sum local = 645.649, global = -319.458, cumulative = -1025.19
smoothSolver:  Solving for omega, Initial residual = 0.000258477, Final residual = 1.59141e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0181236, Final residual = 0.00171421, No Iterations 3
ExecutionTime = 113.31 s  ClockTime = 116 s

Time = 9

smoothSolver:  Solving for Ux, Initial residual = 0.219725, Final residual = 0.0183762, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.158121, Final residual = 0.0111089, No Iterations 7
smoothSolver:  Solving for Uz, Initial residual = 0.303528, Final residual = 0.0257045, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.94797, Final residual = 0.064793, No Iterations 4
time step continuity errors : sum local = 1113.34, global = -309.45, cumulative = -1334.64
smoothSolver:  Solving for omega, Initial residual = 0.000122954, Final residual = 1.00576e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0197785, Final residual = 0.00176744, No Iterations 3
ExecutionTime = 126.36 s  ClockTime = 129 s

Time = 10

smoothSolver:  Solving for Ux, Initial residual = 0.0821976, Final residual = 0.00721542, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.03459, Final residual = 0.00247223, No Iterations 7
smoothSolver:  Solving for Uz, Initial residual = 0.1429, Final residual = 0.0115741, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.861997, Final residual = 0.084008, No Iterations 5
time step continuity errors : sum local = 1083.47, global = 389.776, cumulative = -944.868
smoothSolver:  Solving for omega, Initial residual = 0.000136441, Final residual = 2.43071e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0235868, Final residual = 0.00172782, No Iterations 3
ExecutionTime = 138.69 s  ClockTime = 141 s

Time = 11

smoothSolver:  Solving for Ux, Initial residual = 0.163168, Final residual = 0.0136899, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0793402, Final residual = 0.00752288, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.386931, Final residual = 0.0294783, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.870129, Final residual = 0.0711586, No Iterations 6
time step continuity errors : sum local = 919.618, global = 516.016, cumulative = -428.851
smoothSolver:  Solving for omega, Initial residual = 0.000106145, Final residual = 1.01946e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0196961, Final residual = 0.0018465, No Iterations 3
ExecutionTime = 151.28 s  ClockTime = 154 s

Time = 12

smoothSolver:  Solving for Ux, Initial residual = 0.274151, Final residual = 0.0245354, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.186447, Final residual = 0.0165338, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.156144, Final residual = 0.0136401, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.887478, Final residual = 0.070432, No Iterations 6
time step continuity errors : sum local = 1119.32, global = 661.637, cumulative = 232.786
smoothSolver:  Solving for omega, Initial residual = 0.000100293, Final residual = 2.24457e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0248731, Final residual = 0.00212404, No Iterations 3
ExecutionTime = 164.01 s  ClockTime = 167 s

Time = 13

smoothSolver:  Solving for Ux, Initial residual = 0.205525, Final residual = 0.0187076, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0923904, Final residual = 0.00894338, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.196898, Final residual = 0.0175523, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.817651, Final residual = 0.0790382, No Iterations 8
time step continuity errors : sum local = 1128.57, global = 728.044, cumulative = 960.83
smoothSolver:  Solving for omega, Initial residual = 0.00010952, Final residual = 2.80076e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0218195, Final residual = 0.00158256, No Iterations 4
ExecutionTime = 178.57 s  ClockTime = 181 s

Time = 14

smoothSolver:  Solving for Ux, Initial residual = 0.209248, Final residual = 0.0174402, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.234617, Final residual = 0.0220479, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.330718, Final residual = 0.0281098, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.871195, Final residual = 0.0759581, No Iterations 10
time step continuity errors : sum local = 1172.79, global = 858.035, cumulative = 1818.86
smoothSolver:  Solving for omega, Initial residual = 0.000111802, Final residual = 3.26755e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0256925, Final residual = 0.00171392, No Iterations 4
ExecutionTime = 193.33 s  ClockTime = 196 s

Time = 15

smoothSolver:  Solving for Ux, Initial residual = 0.173376, Final residual = 0.0142151, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.141625, Final residual = 0.0112524, No Iterations 7
smoothSolver:  Solving for Uz, Initial residual = 0.203415, Final residual = 0.0168199, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.850352, Final residual = 0.0820699, No Iterations 8
time step continuity errors : sum local = 1261.82, global = 852.091, cumulative = 2670.96
smoothSolver:  Solving for omega, Initial residual = 0.000119946, Final residual = 3.04546e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0523867, Final residual = 0.00451596, No Iterations 2
ExecutionTime = 206.27 s  ClockTime = 209 s

Time = 16

smoothSolver:  Solving for Ux, Initial residual = 0.189402, Final residual = 0.0147132, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.172538, Final residual = 0.0148313, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.216101, Final residual = 0.0175322, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.831111, Final residual = 0.0800892, No Iterations 13
time step continuity errors : sum local = 1124.03, global = -888.775, cumulative = 1782.18
smoothSolver:  Solving for omega, Initial residual = 0.000132821, Final residual = 2.71341e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0429148, Final residual = 0.00279213, No Iterations 3
ExecutionTime = 220.4 s  ClockTime = 223 s

Time = 17

smoothSolver:  Solving for Ux, Initial residual = 0.240906, Final residual = 0.0185245, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.481738, Final residual = 0.0406032, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.306464, Final residual = 0.0227411, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.829137, Final residual = 0.0797512, No Iterations 5
time step continuity errors : sum local = 1176.95, global = 495.085, cumulative = 2277.27
smoothSolver:  Solving for omega, Initial residual = 0.000148683, Final residual = 1.27859e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0542324, Final residual = 0.0033722, No Iterations 3
ExecutionTime = 231.68 s  ClockTime = 235 s

Time = 18

smoothSolver:  Solving for Ux, Initial residual = 0.206106, Final residual = 0.0161844, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.255402, Final residual = 0.0223374, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.341787, Final residual = 0.0275786, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.827099, Final residual = 0.082525, No Iterations 7
time step continuity errors : sum local = 1134.39, global = 819.278, cumulative = 3096.54
smoothSolver:  Solving for omega, Initial residual = 9.46476e-05, Final residual = 2.25188e-06, No Iterations 2
smoothSolver:  Solving for k, Initial residual = 0.0505738, Final residual = 0.00348548, No Iterations 3
ExecutionTime = 242.31 s  ClockTime = 246 s

Time = 19

smoothSolver:  Solving for Ux, Initial residual = 0.199931, Final residual = 0.0167285, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.361021, Final residual = 0.0275666, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.234612, Final residual = 0.0186392, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.829252, Final residual = 0.0781392, No Iterations 18
time step continuity errors : sum local = 1182.87, global = 1048.32, cumulative = 4144.86
smoothSolver:  Solving for omega, Initial residual = 0.000103443, Final residual = 9.28303e-06, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.160832, Final residual = 0.0122526, No Iterations 1
ExecutionTime = 255.21 s  ClockTime = 259 s

Time = 20

smoothSolver:  Solving for Ux, Initial residual = 0.203207, Final residual = 0.0177404, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.328249, Final residual = 0.0269088, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.224519, Final residual = 0.0164899, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.838448, Final residual = 0.0774968, No Iterations 16
time step continuity errors : sum local = 1282.19, global = 1084.48, cumulative = 5229.34
smoothSolver:  Solving for omega, Initial residual = 0.000137159, Final residual = 1.18779e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.105212, Final residual = 0.00733032, No Iterations 2
ExecutionTime = 268.59 s  ClockTime = 273 s

Time = 21

smoothSolver:  Solving for Ux, Initial residual = 0.195891, Final residual = 0.0166656, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.295498, Final residual = 0.0264212, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.256556, Final residual = 0.019235, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.820159, Final residual = 0.0809549, No Iterations 27
time step continuity errors : sum local = 1210.13, global = -1110.62, cumulative = 4118.72
smoothSolver:  Solving for omega, Initial residual = 0.000209472, Final residual = 1.33081e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.0963909, Final residual = 0.00779254, No Iterations 2
ExecutionTime = 283.95 s  ClockTime = 288 s

Time = 22

smoothSolver:  Solving for Ux, Initial residual = 0.210533, Final residual = 0.0180219, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.451871, Final residual = 0.0368704, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.267514, Final residual = 0.0242343, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.828406, Final residual = 0.0758993, No Iterations 30
time step continuity errors : sum local = 1183.8, global = 1080.35, cumulative = 5199.06
smoothSolver:  Solving for omega, Initial residual = 0.00040392, Final residual = 1.58342e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.116199, Final residual = 0.00839545, No Iterations 2
ExecutionTime = 299.9 s  ClockTime = 304 s

Time = 23

smoothSolver:  Solving for Ux, Initial residual = 0.223249, Final residual = 0.0207695, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.479522, Final residual = 0.0438371, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.287509, Final residual = 0.0233311, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.819938, Final residual = 0.0816971, No Iterations 21
time step continuity errors : sum local = 1364, global = -1206.8, cumulative = 3992.26
smoothSolver:  Solving for omega, Initial residual = 0.000308337, Final residual = 2.31198e-05, No Iterations 1
smoothSolver:  Solving for k, Initial residual = 0.183739, Final residual = 0.0103415, No Iterations 2
ExecutionTime = 314.31 s  ClockTime = 319 s

• 质量严重不守恒，为何不把网格做到10万以下，轻量级跑跑试试？

• @cfd-china 好的，我这么设置周期边界没有问题吧，两个面的距离是22mm

| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
version     2.0;
format      ascii;
class       polyBoundaryMesh;
location    "constant/polyMesh";
object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

8
(
hub
{
type wall;
nFaces          30748;
startFace       5380555;
}
inlet
{
type patch;
nFaces          5546;
startFace       5411303;
}
outlet
{
type patch;
nFaces          4543;
startFace       5416849;
}
per1
{
type            cyclicAMI;
nFaces          15340;
startFace       5421392;
neighbourPatch   per2;
matchTolerance    0.01;
transform         translational;
separationVector (0 0.022 0);
}
per2
{
type            cyclicAMI;
nFaces          15340;
startFace       5436732;
neighbourPatch   per1;
matchTolerance    0.01;
transform         translational;
separationVector (0 -0.022 0);
}
ps
{
type            wall;
nFaces          10266;
startFace       5452072;
}
ss
{
type            wall;
nFaces          11151;
startFace       5462338;
}
sym
{
type symmetryPlane;
nFaces          30748;
startFace       5473489;
}
)

// ************************************************************************* //

• 为何不用cyclic？

per1
{
type            cyclic;
nFaces          15340;
startFace       5421392;
neighbourPatch   per2;
}
per2
{
type            cyclic;
nFaces          15340;
startFace       5436732;
neighbourPatch   per1;
}



• @cfd-china 我感觉这个我理解的可能有问题，是应该给平移距离的大小（0.022m），还是（0,1,0）这么给。当时从这里看见的

separationVector (0 -0.022 0);

• @cfd-china 使用cyclic会有这个问题，因为两侧网格并不完全匹配

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-119cac7e8750
Exec   : simpleFoam
Date   : Mar 25 2017
Time   : 05:33:37
Host   : "ubuntu"
PID    : 48984
Case   : /home/cc/OpenFOAM/cc-3.0.1/run/tutorials/incompressible/simpleFoam/cynasa
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

--> FOAM FATAL ERROR:
face 47 area does not match neighbour by 10.4554% -- possible face ordering problem.
patch:per1 my area:1.63789e-10 neighbour area:1.47515e-10 matching tolerance:0.1
Mesh face:5421439 fc:(0.0411928 0.00446869 5e-07)
Neighbour fc:(0.040347 0.0264296 5e-07)
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.

From function cyclicPolyPatch::calcTransforms()
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 221.

FOAM exiting

• @cfd-china 我的算例可以算了，但是不知道为什么，壁面的分离算的特别小

• @小葱小虫 我只是换了一套网格，其他没动

• @cfd-china 老师您好，由于我的算例进口马赫数为0.7，所以要考虑可压缩性，但是算着算着就发散了，用的rhosimplefoam，请问，能否帮我调试一下，由于我的上面是对称面，下面是壁面，所以网格10万很难满足壁面要求，能否帮忙调试，谢谢

• @小葱小虫 你好！我在计算压缩机时，也出现了time step continuity errors过大导致发散的情况，参考苏军伟的博客，我调整了非正交修正次数，目前没发散，得到了结果。你可以参考一下。
主要就是修改fvsolution文件。

SIMPLE
{
//nNonOrthogonalCorrectors 0;
nNonOrthogonalCorrectors 1;
residualControl
{
p               1e-2;
U               1e-3;
"(k|epsilon)"   1e-3;
}
}

• @Izumi 是增加nNonOrthogonalCorrector这个数吗，谢谢我试试

• @小葱小虫
@Izumi 应该说的就是nNonOrthogonalCorrectors

另外你两个网格之间差别是什么，我也是用pointwise画网格，出现aspect ratio那么大可能是手误也有可能是中间有breakpoint，你用examine aspect ratio看看那些畸形的网格，用connector->distribution->delete breakpoint去掉。因为如果一个方向的网格意外的小，虽然你的求解器是simpleFoam没有CFL数，但也对稳定性有影响吧？

• 我也跟你遇到了相同的问题，不知道是不是计算设置的问题，总压入口是不是还有别的需要注意的地方

• 你好，你现在解决这个问题了吗？我现在也遇到了这个问题，不过是因为pressureDirectedInletVelocity 这个进口条件的问题。
我也在做类似的研究，可以交流一下吗