进口马赫数0.7左右,用rhosimpleFoam不收敛的问题



  • 与上一个帖子一样,http://www.cfd-china.com/topic/744/请教关于进口条件的设置的问题/17
    还是算的平面叶栅,只是改成可压缩的,但是总是挂掉,我贴出所有边界条件以及网格文件,希望大家给些建议。![替代文字]0_1494768632866_计算区域.png

    0_1494769166888_湍流.png
    弦长是40mm,不知道这是不是有问题Turbulence length scale,要不要成个0.07。

    | =========                 |                                                 |
    | \\      /  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;
    
    boundaryField
    {
      hub
        {
         type symmetryPlane;
        }
        inlet
        {
            type            turbulentIntensityKineticEnergyInlet;
            intensity       0.05;
            value           uniform 194;
        }
        outlet
        {
            type            zeroGradient;
        }
        per1
        {
               type            cyclicAMI;
       	 value           uniform 190;
        }
        per2
        {
           type            cyclicAMI;
    	  value           uniform 190; 
       	
        }
        ps
        {
            type            kqRWallFunction;
            value           uniform 190;
        }
        ss
        {
        	type            kqRWallFunction;
            value           uniform 190;
        }
        sym
        { 
        type symmetryPlane;
        }
        
    }
    
    
    // ************************************************************************* //
    ```/*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  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      epsilon;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 2 -3 0 0 0 0];
    
    internalField   uniform 6043;
    
    boundaryField
    {   
    
       hub
        {
         	type            epsilonWallFunction;
            value           uniform 6043;
        }
        inlet
        {
            
            type           turbulentMixingLengthDissipationRateInlet;
            mixingLength    0.05;
            value           uniform 6043;
        }
        outlet
        {
           type            zeroGradient;
        }
        per1
        {
            type            cyclicAMI;
       	
        }
        per2
        {
        	   type            cyclicAMI;
       	
        }
        ps
        {
           type            epsilonWallFunction;
            value           uniform 6043;
        }
        ss
        {
         	type            epsilonWallFunction;
            value           uniform 6043;
    ;
        }
        sym
        { 
        	type symmetryPlane;
        }
    	
    }
    
    
    // ************************************************************************* //
    ```/*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  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      500;
    
    dimensions      [0 0 -1 0 0 0 0];
    
    internalField   uniform $omegaInlet;
    
    boundaryField
    {
    
         hub
        {
            type            omegaWallFunction;
            value           $internalField;
        }
          inlet
        {
           type            turbulentMixingLengthFrequencyInlet;
            mixingLength    0.05;
            k               k;
            value           $internalField;
        }
        outlet
        {
            type            zeroGradient;
        }
         per1
        {
            type            cyclicAMI;
              value           $internalField; 
       	
        }
        per2
        {
        	   type            cyclicAMI;
               value           $internalField; 
               
       	
        }
    
        ps
        {
            type            omegaWallFunction;
            value           $internalField;
        }
        ss
        {
            type            omegaWallFunction;
            value           $internalField;
        }
        sym
        { 
           type symmetryPlane;
        }
        
    }
    
    // ************************************************************************* //
    ```/*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  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      [1 -1 -2 0 0 0 0];
    
    internalField   uniform 103325;
    
    boundaryField
    {
        hub
        {
             type            zeroGradient;
        }
        
        inlet
        {
               type            totalPressure;
               p0              uniform 119000;
               gamma           1.4; 
               value           uniform 119000;
        }
        outlet
        {
        	
            type            fixedValue;
            value           uniform 101325;
       
        }
        
         per1
        {
           type            cyclicAMI;
          
        }
        per2
        {
         type            cyclicAMI;
           
        }
        ps
        {
              type            zeroGradient;
        }
        ss
        {
        	 type            zeroGradient;
        }
        sym
        { 
       	type symmetryPlane;
           
        }
           
    }
    
    
    // ************************************************************************* //
    ```/*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  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      T;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 0 0 1 0 0 0];
    
    internalField   uniform 298;
    
    boundaryField
    {
      
        hub
        {
             type            zeroGradient;
        }
        inlet
        {
           type	totalTemperature;
            T0	uniform 300;
    	gamma 	1.4;
        }
        outlet
        {
        	
            type            zeroGradient;
        }
        per1
        {
          type            cyclicAMI;
          
        }
        per2
        {
        	type            cyclicAMI;
           
        }
        ps
        {
              type            zeroGradient;
        }
        ss
        {
        	 type            zeroGradient;
        }
        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      T;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 0 0 1 0 0 0];
    
    internalField   uniform 298;
    
    boundaryField
    {
      
        hub
        {
             type            zeroGradient;
        }
        inlet
        {
           type	totalTemperature;
            T0	uniform 300;
    	gamma 	1.4;
        }
        outlet
        {
        	
            type            zeroGradient;
        }
        per1
        {
          type            cyclicAMI;
          
        }
        per2
        {
        	type            cyclicAMI;
           
        }
        ps
        {
              type            zeroGradient;
        }
        ss
        {
        	 type            zeroGradient;
        }
        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       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 (126 113 0);
        }
        outlet
        {	
                  type            inletOutlet;
                  inletValue      $internalField;
                  value           uniform (63 55.5 0);
           
        }
        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;
        }
    }
    
    
    // ************************************************************************* //


  • 介质就是空气

    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  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;
        location    "constant";
        object      thermophysicalProperties;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    thermoType
    	{
    	    type            hePsiThermo;
        mixture         pureMixture;
        transport       sutherland;
        thermo          hConst;
        equationOfState perfectGas;
        specie          specie;
        energy          sensibleInternalEnergy;
    	}
    mixture
    {
     
        specie
        {
            nMoles      1;
            molWeight   28.9;
        }
        thermodynamics
        {
            Cp          1007;
            Hf          0;
        }
        transport
        {
            As          1.4792e-06;
            Ts          116;
        }
    }
    
    
    // ************************************************************************* //
    


  • | =========                 |                                                 |
    | \\      /  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;
        location    "constant";
        object      turbulenceProperties;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    simulationType RAS;
    
    RAS
    {
        // Tested with kEpsilon, realizableKE, kOmega, kOmegaSST, v2f,
        // ShihQuadraticKE, LienCubicKE.
        RASModel        kOmegaSST;
    
        turbulence      on;
    
        printCoeffs     on;
    }
    
    
    // ************************************************************************* //


  • 网格checkmesh OK

    | =========                 |                                                 |
    | \\      /  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          16560;
            startFace       2399630;
        }
        inlet
        {
             type patch;
            nFaces          3430;
            startFace       2416190;
        }
        outlet
        {
            type patch;
            nFaces          2842;
            startFace       2419620;
        }
        per1
        {
            type            cyclicAMI;
            nFaces          9163;
            startFace       2422462;
            neighbourPatch   per2;
            matchTolerance    0.01;
    	transform         translational;
    	separationVector (0 0.022 0);
        }
        per2
        {  
    
           type            cyclicAMI;
           nFaces          9163;
           startFace       2431625;
           neighbourPatch   per1;
           matchTolerance    0.01;
           transform         translational;
           separationVector (0 -0.022 0);
        }
        ps
        {
            type            wall;
            nFaces          5831;
            startFace       2440788;
        }
        ss
        {
             type            wall;
            nFaces          5831;
            startFace       2446619;
        }
        sym
        {
            type symmetryPlane;
            nFaces          16560;
            startFace       2452450;
        }
    )
    
    // ************************************************************************* //


  • | =========                 |                                                 |
    | \\      /  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      controlDict;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    application     rhoSimpleFoam;
    
    startFrom        latestTime;
    
    startTime       0;
    
    stopAt          endTime;
    
    endTime         100;
    
    deltaT          1e-4
    
    writeControl    timeStep;
    
    writeInterval   1000;
    
    purgeWrite      0;
    
    writeFormat     ascii;
    
    writePrecision  6;
    
    writeCompression off;
    
    timeFormat      general;
    
    timePrecision   6;
    
    graphFormat     raw;
    
    runTimeModifiable true;
    
    adjustTimeStep  on;
    
    // ************************************************************************* //


  • | =========                 |                                                 |
    | \\      /  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
    {
        default          none;
        
        div(phi,U)       bounded Gauss upwind;
        div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
        div(phi,e)      bounded Gauss upwind;
        div(phi,epsilon) bounded Gauss upwind;
        div(phi,k)       bounded Gauss upwind;
        div(phi,Ekp)     bounded Gauss upwind;
        div(phi,K)       bounded Gauss upwind;
        div(phi,h)       bounded Gauss upwind;
        div(phi,omega)   bounded Gauss upwind;
        div(phid,p)      bounded Gauss upwind;
    }
    
    laplacianSchemes
    {
        default         Gauss linear corrected;
    }
    
    interpolationSchemes
    {
        default         linear;
    }
    
    snGradSchemes
    {
        default         corrected;
    }
    wallDist
    {
        method meshWave;
    }
    fluxRequired
    {
      default         no;
        p               ;
    }
    
    // ************************************************************************* //


  • | =========                 |                                                 |
    | \\      /  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;
            tolerance       1e-05;
            relTol          0.1;
            smoother        GaussSeidel;
            nPreSweeps      0;
            nPostSweeps     2;
            nFinestSweeps   2;
            cacheAgglomeration true;
            nCellsInCoarsestLevel 20;
            agglomerator    faceAreaPair;
            mergeLevels     1;
        }
    
        "(U|omega|e|k|epsilon)"
        {
           solver           smoothSolver;
            smoother         GaussSeidel;
            tolerance        1e-8;
            relTol           0.1;
            nSweeps          1;
        }
    }
    
    SIMPLE
    {
        nNonOrthogonalCorrectors 1;
        rhoMin          0.1;
        rhoMax          2;
    //   consistent      yes;
      transonic       yes;
    
        residualControl
        {
            p               1e-5;
            U               1e-5;
            e               1e-5;
    
            // possibly check turbulence fields
            "(k|epsilon|omega)" 1e-5;
        }
    }
    
    relaxationFactors
    {
      fields
        {
            p               0.3;
            rho             0.7;
        }
        equations
        {
            U               0.7;
            "(k|epsilon)"   0.7;
            e               0.7;
        }
    }
    // ************************************************************************* //

  • 管理员

    你的速度不是正方向?并且和进口压力有关?试过简单的进出口边界条件不?直觉你这个边界条件不稳定哈。





  • 我给的进口总温总压,出口背压,速度让它自己计算,速度是有方向的,不知道如何给稳定些



  • omega飙升

    0_1494809875821_upload-4d2ab3d3-7660-4854-96eb-7690a9048342

    Time = 0.0275
    Duration of pickling 0.0523500442505 too long. Extending frequency from 1.0 to 2.61750221252
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00152996, Final residual = 0.000123741, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.00222516, Final residual = 0.000169915, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00285202, Final residual = 0.000244214, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00798671, Final residual = 0.000494355, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0160753, Final residual = 0.0013806, No Iterations 4
    GAMG:  Solving for p, Initial residual = 0.00143423, Final residual = 0.000143128, No Iterations 94
    time step continuity errors : sum local = 0.32827, global = 0.0097095, cumulative = 1.35785
    rho max/min : 1.34304 1.20262
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 1.08444e-49, No Iterations 1
    bounding omega, min: -5.71564e+59 max: 7.78487e+77 average: 1.01456e+72
    smoothSolver:  Solving for k, Initial residual = 8.8899e-07, Final residual = 1.0766e-40, No Iterations 1
    bounding k, min: -8.90538e-08 max: 713.841 average: 2.87127
    ExecutionTime = 3087.11 s  ClockTime = 3100 s
    
    Time = 0.0276
    Duration of pickling 0.0545029640198 too long. Extending frequency from 1.0 to 2.72514820099
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0015129, Final residual = 0.000121723, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.00220558, Final residual = 0.000167395, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00284466, Final residual = 0.000244908, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00802377, Final residual = 0.000497091, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0172986, Final residual = 0.00170545, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.00174044, Final residual = 0.000172569, No Iterations 68
    time step continuity errors : sum local = 0.326715, global = 0.0096409, cumulative = 1.36749
    rho max/min : 1.34304 1.20279
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 2.9831e-72, No Iterations 1
    bounding omega, min: -1.96379e+59 max: 1.91134e+73 average: 2.35575e+67
    smoothSolver:  Solving for k, Initial residual = 1.89765e-06, Final residual = 1.06009e-29, No Iterations 1
    bounding k, min: -8.15996e-06 max: 718.357 average: 2.79046
    ExecutionTime = 3096.76 s  ClockTime = 3109 s
    
    Time = 0.0277
    Duration of pickling 0.0812048912048 too long. Extending frequency from 1.0 to 4.06024456024
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0014993, Final residual = 0.000119976, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.00218975, Final residual = 0.000165293, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00283678, Final residual = 0.000245435, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.0080669, Final residual = 0.00049988, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0190055, Final residual = 0.00162944, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.00165444, Final residual = 0.000165306, No Iterations 58
    time step continuity errors : sum local = 0.325202, global = 0.0095477, cumulative = 1.37704
    rho max/min : 1.34304 1.20296
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 1.65172e-63, No Iterations 1
    bounding omega, min: -3.89243e+67 max: 4.18007e+77 average: 6.1041e+71
    smoothSolver:  Solving for k, Initial residual = 4.19584e-08, Final residual = 1.05679e-23, No Iterations 1
    bounding k, min: -2.76439e-08 max: 720.659 average: 2.71377
    ExecutionTime = 3105.59 s  ClockTime = 3118 s
    
    Time = 0.0278
    Duration of pickling 0.0585401058197 too long. Extending frequency from 1.0 to 2.92700529099
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00149028, Final residual = 0.000118581, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.00218168, Final residual = 0.000163814, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00282467, Final residual = 0.000245514, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00812588, Final residual = 0.000503407, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0211776, Final residual = 0.00167735, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.00169416, Final residual = 0.000167178, No Iterations 52
    time step continuity errors : sum local = 0.323742, global = 0.0094258, cumulative = 1.38646
    rho max/min : 1.34304 1.18013
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 5.93357e-72, No Iterations 1
    bounding omega, min: -1.81139e+71 max: 4.20153e+75 average: 5.17821e+69
    smoothSolver:  Solving for k, Initial residual = 1.6527e-05, Final residual = 5.61923e-22, No Iterations 1
    bounding k, min: -6.39732e-08 max: 711.177 average: 2.6393
    ExecutionTime = 3113.97 s  ClockTime = 3127 s
    
    Time = 0.0279
    Duration of pickling 0.0739150047302 too long. Extending frequency from 1.0 to 3.69575023651
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00149312, Final residual = 0.000117778, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.0021903, Final residual = 0.000163573, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00283009, Final residual = 0.00024665, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00820162, Final residual = 0.000508201, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0246477, Final residual = 0.00198208, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.00200497, Final residual = 0.000198396, No Iterations 55
    time step continuity errors : sum local = 0.322345, global = 0.00924313, cumulative = 1.39571
    rho max/min : 1.34304 1.17683
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 1.13229e-69, No Iterations 1
    bounding omega, min: -4.25636e+69 max: 8.12744e+67 average: -5.13238e+63
    smoothSolver:  Solving for k, Initial residual = 2.60154e-07, Final residual = 7.34386e-26, No Iterations 1
    bounding k, min: -5.22454e-09 max: 713.695 average: 2.56462
    ExecutionTime = 3122.56 s  ClockTime = 3135 s
    
    Time = 0.028
    Duration of pickling 0.067948102951 too long. Extending frequency from 1.0 to 3.39740514755
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00150484, Final residual = 0.000117816, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.00220887, Final residual = 0.000164596, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00285433, Final residual = 0.000249085, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00832458, Final residual = 0.00051524, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0284494, Final residual = 0.00250348, No Iterations 4
    GAMG:  Solving for p, Initial residual = 0.00262877, Final residual = 0.000259279, No Iterations 92
    time step continuity errors : sum local = 0.321011, global = 0.00897452, cumulative = 1.40468
    rho max/min : 1.34305 1.09793
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 3.60599e-62, No Iterations 1
    bounding omega, min: -1.30817e+62 max: 5.18859e+60 average: -1.70377e+56
    smoothSolver:  Solving for k, Initial residual = 0.650092, Final residual = 5.24055e-24, No Iterations 1
    bounding k, min: -5.18956e-11 max: 710.995 average: 2.49681
    ExecutionTime = 3134.49 s  ClockTime = 3147 s
    
    Time = 0.0281
    Duration of pickling 0.0619690418243 too long. Extending frequency from 1.0 to 3.09845209122
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00152832, Final residual = 0.000117872, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.00223761, Final residual = 0.000166456, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00292184, Final residual = 0.000253548, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00846883, Final residual = 0.000523429, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0354473, Final residual = 0.00331063, No Iterations 6
    GAMG:  Solving for p, Initial residual = 0.00364502, Final residual = 0.000362752, No Iterations 125
    time step continuity errors : sum local = 0.319631, global = 0.00861118, cumulative = 1.41329
    rho max/min : 1.34305 1.06718
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 9.78348e-48, No Iterations 1
    bounding omega, min: -4.71803e+50 max: 1.27516e+67 average: 1.57147e+61
    smoothSolver:  Solving for k, Initial residual = 0.993077, Final residual = 6.53872e-31, No Iterations 1
    bounding k, min: -1.47294e-10 max: 799.109 average: 2.43261
    ExecutionTime = 3149.37 s  ClockTime = 3162 s
    
    Time = 0.0282
    Duration of pickling 0.0754199028015 too long. Extending frequency from 1.0 to 3.77099514008
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00160754, Final residual = 0.000122563, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.00234167, Final residual = 0.000174383, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.00305302, Final residual = 0.000264702, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00890675, Final residual = 0.000544697, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0456857, Final residual = 0.00450887, No Iterations 11
    GAMG:  Solving for p, Initial residual = 0.00532448, Final residual = 0.000530223, No Iterations 148
    time step continuity errors : sum local = 0.318517, global = 0.00804085, cumulative = 1.42133
    rho max/min : 1.34304 0.944974
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 5.03085e-61, No Iterations 1
    bounding omega, min: -6.12926e+56 max: 1.89257e+63 average: 2.58001e+57
    smoothSolver:  Solving for k, Initial residual = 0.0986294, Final residual = 2.34926e-17, No Iterations 1
    bounding k, min: -3.64697e-08 max: 1561.98 average: 2.37097
    ExecutionTime = 3166.71 s  ClockTime = 3180 s
    
    Time = 0.0283
    Duration of pickling 0.0689640045166 too long. Extending frequency from 1.0 to 3.44820022583
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00188211, Final residual = 0.00013498, No Iterations 5
    smoothSolver:  Solving for Uy, Initial residual = 0.0026032, Final residual = 0.000188799, No Iterations 5
    smoothSolver:  Solving for Uz, Initial residual = 0.0033437, Final residual = 0.000287587, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.00997634, Final residual = 0.000591898, No Iterations 3
    GAMG:  Solving for p, Initial residual = 0.0732098, Final residual = 0.00729724, No Iterations 26
    GAMG:  Solving for p, Initial residual = 0.0102857, Final residual = 0.00102317, No Iterations 160
    time step continuity errors : sum local = 0.316687, global = 0.00719629, cumulative = 1.42853
    rho max/min : 1.368 0.265827
    smoothSolver:  Solving for omega, Initial residual = 1, Final residual = 2.76757e-55, No Iterations 1
    bounding omega, min: -2.79393e+51 max: 1.89919e+65 average: 2.72917e+59
    smoothSolver:  Solving for k, Initial residual = 0.00338243, Final residual = 8.55618e-21, No Iterations 1
    bounding k, min: -2.59751e-09 max: 2425.35 average: 2.31181
    ExecutionTime = 3186.64 s  ClockTime = 3200 s
    
    Time = 0.0284
    Duration of pickling 0.0688588619232 too long. Extending frequency from 1.0 to 3.44294309616
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00844988, Final residual = 0.000751562, No Iterations 3
    smoothSolver:  Solving for Uy, Initial residual = 0.00963742, Final residual = 0.000907189, No Iterations 3
    smoothSolver:  Solving for Uz, Initial residual = 0.00423285, Final residual = 0.000343988, No Iterations 5
    smoothSolver:  Solving for e, Initial residual = 0.24105, Final residual = 0.013761, No Iterations 2
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::sigFpe::sigHandler(int) at ??:?
    #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
    #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
    #5  ? at ??:?
    #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #7  ? at ??:?
    Floating point exception (core dumped)
    cc@ubuntu:~/OpenFOAM/cc-3.0.1/run/tutorials/compressible/rhoSimpleFoam$


  • 最近看你一直在算这个算例,好像卡在这了?你这个就是发散了,发散的原因非常多,不能简单用一种药治疗。
    把你的算例发给我试试?网格数量最好少于1万,或者5万以下。
    li.dy@dyfluid.com



  • @小葱小虫

    要不算个初始场?

    在fvSolution中加

    solvers
    {
    ...
         Phi
         {
              solver PCG;
              preconditioner PIC;
              tolerance 1e-12;
              relTol 0;
         }
    }
    
    potentialFlow
    {
        nCorrectors 10;
    }
    

    然后运行potentialFlow搞个速度初始场试试?

    或者在controlDict中加一些debugSwitches

    DebugSwitches
    {
    lduMatrix 2;
    ...// 去$FOAM_ETC/controlDict下找和湍流有关的。
    }
    

    还有就是加functionObject : fieldMinMax看看是哪个位置的场出了问题

        fieldMinMax1
        {
            type        fieldMinMax;
            libs ("libfieldFunctionObjects.so");
            write       yes;
            log         yes;
            location    yes; // 关键是这个!
            mode        magnitude;
            fields
            (
                U
                p
                omega
            );
        }
    


  • @程迪 感谢!正在尝试


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待