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

小葱小虫

| =========                 |                                                 |
| \\      /  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
{
    grad(U);
}
// ************************************************************************* //

小葱小虫

| =========                 |                                                 |
| \\      /  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         steadyState;
}

gradSchemes
{
     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;
    div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}

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

小葱小虫

0_1490344753819_upload-1d58f4d0-683b-446a-884b-a44fff459031

CFD中文网

谢谢。
进出口边界条件？

小葱小虫

@cfd-china

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

CFD中文网

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

小葱小虫

| =========                 |                                                 |
| \\      /  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
    {
              type            zeroGradient;
    }
    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
    {
          type            zeroGradient;
    }
    ss
    {
    	 type            zeroGradient;
    }
    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
    {	
         type            zeroGradient;
       
    }
    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
    {
        type            zeroGradient;
    }
    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;
    }
    
}

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

random_ran

我注意到你的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

Reading field p

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
Reading field U

Reading/calculating face flux field phi

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

CFD中文网

质量严重不守恒，为何不把网格做到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;
    }
)

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

CFD中文网

为何不用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.
Rerun with cyclic debug flag set for more information.

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

FOAM exiting

小葱小虫

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

小葱小虫

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

CFD中文网

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