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



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



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



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



  • @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.
        Cell to face addressing 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
    {
        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-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
        {
                  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;
        }
        
    }
    
    // ************************************************************************* //


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


  • 质量严重不守恒,为何不把网格做到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)这么给。当时从这里看见的0_1490441899214_cyclic.jpg

    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-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 这个进口条件的问题。
    我也在做类似的研究,可以交流一下吗


Log in to reply