射流算例速度发散了,有谁知道怎么解决吗



  • FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        location    "system";
        object      controlDict;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    application     myFoam;
    
    startFrom       latestTime;
    
    startTime       0;
    
    stopAt          endTime;
    
    endTime         1.8e-04;
    
    deltaT          1e-09;
    
    writeControl    adjustableRunTime;
    
    writeInterval   1e-07;
    
    purgeWrite      0;
    
    writeFormat     binary;
    
    writePrecision  8;
    
    writeCompression uncompressed;
    
    timeFormat      general;
    
    timePrecision   6;
    
    runTimeModifiable yes;
    
    adjustTimeStep  yes;
    
    maxCo           0.5;
    maxAlphaCo      0.5;
    
    maxDeltaT       1e-09;
    
    
    // ************************************************************************* //
    


  • *--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.4.0                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        location    "system";
        object      fvSchemes;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    ddtSchemes
    {
        default         Euler;
    }
    
    gradSchemes
    {
        default         Gauss linear;
    }
    
    divSchemes
    {   
        default          Gauss upwind;
        div(phi,alpha)  Gauss vanLeer;
        div(phirb,alpha) Gauss linear;
        div(rhoPhi,U)  Gauss upwind;
        div(phi,thermo:rho.fuel) Gauss upwind;
        div(phi,thermo:rho.air) Gauss upwind;
        div(rhoPhi,T)  Gauss upwind;
        div(rhoPhi,K)  Gauss upwind;
        div(phi,p)      Gauss upwind;
        div(phi,k)      Gauss upwind;
        div((muEff*dev2(T(grad(U))))) Gauss linear;
    }
    
    laplacianSchemes
    {
        default         Gauss linear uncorrected;
    }
    
    interpolationSchemes
    {
        default         linear;
    }
    
    snGradSchemes
    {
        default         uncorrected;
    }
    
    fluxRequired
    {
        default         no;
        p_rgh;
        pcorr;
    }
    
    
    // ************************************************************************* //


  • /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.4.0                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        location    "system";
        object      fvSolution;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    solvers
    {
        alpha.fuel
        {
            nAlphaCorr      1;
            nAlphaSubCycles 1;
            cAlpha          1;
        }
    
        ".*(rho|rhoFinal)"
        {
            solver          diagonal;
        }
    
        pcorr
        {
            solver          GAMG;
            preconditioner
            {
                preconditioner  GAMG;
                tolerance       1e-05;
                relTol          0;
                smoother        GaussSeidel;
                nPreSweeps      0;
                nPostSweeps     2;
                nFinestSweeps   2;
                cacheAgglomeration true;
                nCellsInCoarsestLevel 10;
                agglomerator    faceAreaPair;
                mergeLevels     1;
            }
            tolerance       1e-05;
            relTol          0;
            maxIter         100;
        }
    
        p_rgh
        {
            solver          GAMG;
            tolerance       1e-07;
            relTol          0.01;
            smoother        GaussSeidel;
            nPreSweeps      0;
            nPostSweeps     2;
            nFinestSweeps   2;
            cacheAgglomeration true;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels     1;
        }
    
        p_rghFinal
        {
            solver          GAMG;
            preconditioner
            {
                preconditioner  GAMG;
                tolerance       1e-07;
                relTol          0;
                nVcycles        2;
                smoother        GaussSeidel;
                nPreSweeps      2;
                nPostSweeps     2;
                nFinestSweeps   2;
                cacheAgglomeration true;
                nCellsInCoarsestLevel 10;
                agglomerator    faceAreaPair;
                mergeLevels     1;
            }
            tolerance       1e-07;
            relTol          0;
            maxIter         20;
        }
    
        "(U|T|k|B|nuTilda).*"
        {
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1e-08;
            relTol          0;
        }
    "e.*"
        {
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1e-8;
            relTol          0;
            minIter         1;
        }
    }
    
    PIMPLE
    {
        momentumPredictor yes;
        transonic       yes;
        nOuterCorrectors 5;
        nCorrectors     1;
        nNonOrthogonalCorrectors 2;
        residualControl
        {
            U
            {
              tolerance 1e-07;
              relTol    0;
            }
            p_rgh
            {
              tolerance 1e-07;
              relTol    0;
            }
            T
            {
              tolerance 1e-07;
              relTol    0;
            }
            k
            {
              tolerance 1e-07;
              relTol    0;
            }
    }
    
    relaxtionFactors
    {
        fields
        {
           p_rgh 0.4;
           p_rghFinal 1;
        }
        equations
        {
           "(U|T|k|B|nuTilda)" 0.4;
           "(U|T|k|B|nuTilda)Final" 1;
        }
    }
    // ************************************************************************* //
    


  • Much appreciated if you provide more information and enforce the operation depicted in this thread.



  • 0_1464929999007_QQ截图20160603130508.png
    大口进,小口出



  • 速度

    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.4.0                                 |
    |   \\  /    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 (0 0 0);
    
    boundaryField
    {
        INLET
        {
            type            zeroGradient;
        }
    
        WALL
        {
            type            zeroGradient;
        }
    
        OUTLET
        {
            type            zeroGradient;
        }
    }
    
    // ************************************************************************* //
    


  • /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.4.0                                 |
    |   \\  /    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 3e+06;
    
    boundaryField
    {
        INLET
        {
            type            fixedValue;
            value           uniform 100e+6;
        }
    
        WALL
        {
            type            zeroGradient;
        }
    
        OUTLET
        {
            type            zeroGradient;
        }
    }
    
    // ************************************************************************* //
    


  • /*------------none--------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.4.0                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volScalarField;
        object      p_rgh;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [1 -1 -2 0 0 0 0];
    
    internalField   uniform 3e+06;
    
    boundaryField
    {
        INLET
        {
            type            fixedValue;
            value           uniform 100e+06;
        }
    
        WALL
        {
            type            zeroGradient;
        }
    
        OUTLET
        {
            type            zeroGradient;
        }
    }
    
    // ************************************************************************* //
    


  • /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2.4.0                                 |
    |   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volScalarField;
        object      alpha.fuel;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 0 0 0 0 0 0];
    
    internalField   uniform 0;
    
    boundaryField
    {
        INLET
        {
            type            fixedValue;
            value           uniform 1;
        }
    
        WALL
        {
            type            zeroGradient;
        }
    
        OUTLET
        {
            type            zeroGradient;
        }
    }
    
    // ************************************************************************* //
    


  • 速度和压力设置不对?



  • 请把log文件的信息贴出来。


Log in to reply