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



  • 我需要给定进口总温总压,进口气流速度角度,以及出口静压请问,这么给是否有问题,总是挂。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中文网 的连接断开,我们正在尝试重连,请耐心等待