CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (https://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   Temperqture field not solved in LES - heat transfer - rhopimplefoam (https://www.cfd-online.com/Forums/openfoam-pre-processing/239800-temperqture-field-not-solved-les-heat-transfer-rhopimplefoam.html)

arun1994 November 24, 2021 15:13

Temperqture field not solved in LES - heat transfer - rhopimplefoam
 
Hello!

I am trying to run a film cooling simulation in openfoam using LES. All the tutorials I have seen so far are cold flows without changes in temperature. I would like to know if there is a 'sample test case' for forced convection problems. I have some nice results with RANS but I am struggling with LES. Is there a solver that works best for these types of problems?



In the case below, I have some velocity field in the flow but no changes in the temperature field with time...

I am pretty sure that I am missing something blatantly wrong. Would be nice if someone could point it out please!


Code:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                |
  \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox
  \\    /  O peration    | Website:  https://openfoam.org
    \\  /    A nd          | Version:  7
    \\/    M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField  uniform (10 0 0);
boundaryField
{
    Body_5.Inlet
    {
        type            surfaceNormalFixedValue;
        refValue        uniform -10;
    }
    Body_5.Inlet_coolant
    {
        type            surfaceNormalFixedValue;
        refValue        uniform -10.68323;
    }
    Body_5.Basewall
    {
        type            noSlip;
    }
    Body_5.side2
    {
        type            noSlip;
    }
    Body_5.Top
    {
        type            noSlip;
    }
    Body_5.Side1
    {
        type            noSlip;
    }
    Body_5.Outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value          uniform (0 0 0);
    }
    Body_5.Wall_coolant
    {
        type            noSlip;
    }
}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                |
  \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox
  \\    /  O peration    | Website:  https://openfoam.org
    \\  /    A nd          | Version:  7
    \\/    M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField  uniform 293;
boundaryField
{
    Body_5.Inlet
    {
        type            fixedValue;
        value          uniform 293;
    }
    Body_5.Inlet_coolant
    {
        type            fixedValue;
        value          uniform 313;
    }
    Body_5.Basewall
    {
        type            zeroGradient;
    }
    Body_5.side2
    {
        type            zeroGradient;
    }
    Body_5.Top
    {
        type            zeroGradient;
    }
    Body_5.Side1
    {
        type            zeroGradient;
    }
    Body_5.Outlet
    {
        type            inletOutlet;
        value          $internalField;
        inletValue      $internalField;
    }
    Body_5.Wall_coolant
    {
        type            zeroGradient;
    }
}


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

Code:

/*--------------------------------*- C++  -*----------------------------------*\
  =========                |
  \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox
  \\    /  O peration    | Website:  https://openfoam.org
    \\  /    A nd          | Version:  7
    \\/    M anipulation  |
/*--------------------------------*- C++ -*----------------------------------*\
  =========                |
  \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox
  \\    /  O peration    | Website:  https://openfoam.org
    \\  /    A nd          | Version:  7
    \\/    M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField  uniform 101325;


boundaryField
{
    Body_5.Inlet
    {
        type            zeroGradient;
        // type            mixed;
        refValue        uniform 101325;
        refGradient    uniform 0;
        valueFraction  uniform 0.3;
    }
    Body_5.Inlet_coolant
    {
        type            zeroGradient;
        // type            mixed;
        refValue        uniform 101325;
        refGradient    uniform 0;
        valueFraction  uniform 0.3;
    }
    Body_5.Basewall
    {
        type            zeroGradient;
    }
    Body_5.side2
    {
        type            zeroGradient;
    }
    Body_5.Top
    {
        type            zeroGradient;
    }
    Body_5.Side1
    {
        type            zeroGradient;
    }
    Body_5.Outlet
    {
        type            fixedValue;
        value          uniform 101325;

        // type            mixed;
        // refValue        uniform 110000;
        // refGradient    uniform 0;
        // valueFraction  uniform 1;
        // type            transonicOutletPressure;
        // U              U;
        // phi            phi;
        // gamma          1.4;
        // psi            psi;
        // pInf            uniform 110000;
    }
    Body_5.Wall_coolant
    {
        type            zeroGradient;
    }
}


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

Code:

/*--------------------------------*- C++  -*----------------------------------*\
  =========                |
  \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox
  \\    /  O peration    | Website:  https://openfoam.org
    \\  /    A nd          | Version:  7
    \\/    M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField  uniform 0;

boundaryField
{
    Body_5.Inlet
    {
        type            calculated;
        value          uniform 0;
    }
    Body_5.Inlet_coolant
    {
        type            calculated;
        value          uniform 0;
    }
    Body_5.Basewall
    {
        type            nutkWallFunction;
        Cmu            0.09;
        kappa          0.41;
        E              9.8;
        value          uniform 0;
    }
    Body_5.side2
    {
        type            nutkWallFunction;
        Cmu            0.09;
        kappa          0.41;
        E              9.8;
        value          uniform 0;
    }
    Body_5.Top
    {
        type            nutkWallFunction;
        Cmu            0.09;
        kappa          0.41;
        E              9.8;
        value          uniform 0;
    }
    Body_5.Side1
    {
        type            nutkWallFunction;
        Cmu            0.09;
        kappa          0.41;
        E              9.8;
        value          uniform 0;
    }
    Body_5.Outlet
    {
        type            calculated;
        value          uniform 0;
    }
    Body_5.Wall_coolant
    {
        type            nutkWallFunction;
        Cmu            0.09;
        kappa          0.41;
        E              9.8;
        value          uniform 0;
    }
}


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

Code:

/*--------------------------------*- C++  -*----------------------------------*\
  =========                |
  \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox
  \\    /  O peration    | Website:  https://openfoam.org
    \\  /    A nd          | Version:  7
    \\/    M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    location    "0";
    object      alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -1 0 0 0 0];

internalField  uniform 0;

boundaryField
{
    Body_5.Inlet
    {
        type            calculated;
        value          uniform 0;
    }
    Body_5.Inlet_coolant
    {
        type            calculated;
        value          uniform 0;
    }
    Body_5.Basewall
    {
        type            compressible::alphatWallFunction;
        Prt            0.85;
        value          uniform 0;
    }
    Body_5.side2
    {
        type            compressible::alphatWallFunction;
        Prt            0.85;
        value          uniform 0;
    }
    Body_5.Top
    {
        type            compressible::alphatWallFunction;
        Prt            0.85;
        value          uniform 0;
    }
    Body_5.Side1
    {
        type            compressible::alphatWallFunction;
        Prt            0.85;
        value          uniform 0;
    }
    Body_5.Outlet
    {
        type            calculated;
        value          uniform 0;
    }
    Body_5.Wall_coolant
    {
        type            compressible::alphatWallFunction;
        Prt            0.85;
        value          uniform 0;
    }
}


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

Code:

/*--------------------------------*- 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      turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType  LES;

LES
{
    LESModel        Smagorinsky;

    turbulence      on;

    printCoeffs    on;

    delta          cubeRootVol;

    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }

    PrandtlCoeffs
    {
        delta          cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }

        smoothCoeffs
        {
            delta          cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }

            maxDeltaRatio  1.1;
        }

        Cdelta          0.158;
    }

    vanDriestCoeffs
    {
        delta          cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }

        smoothCoeffs
        {
            delta          cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }

            maxDeltaRatio  1.1;
        }

        Aplus          26;
        Cdelta          0.158;
    }

    smoothCoeffs
    {
        delta          cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }

        maxDeltaRatio  1.1;
    }
}


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

Code:

/*--------------------------------*- 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      const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
    specie
    {
        nMoles      1;
        molWeight  28.9;
    }
    thermodynamics
    {
        Cp          1007;
        Hf          0;
    }
    transport
    {
        mu          1.8e-05; //1.2 kg/m^3 x 1.5e-05 m^2/s
        Pr        1;
    }
}


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

Code:

/*--------------------------------*- 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    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "(p|rho)"
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance      1e-6;
        relTol          1e-4;
    }

    "(p|rho)Final"
    {
        $p;
        relTol          0;
    }

    "(U|k|nuTilda)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance      1e-6;
        relTol          1e-6;
    }

    "(U|k|nuTilda)Final"
    {
        $U;
        relTol          0;
    }
    e
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance      10;
        relTol          1;
    }

    eFinal
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance      10;
        relTol          1;
    }
}

PIMPLE
{
    momentumPredictor yes;
    nOuterCorrectors 2;
    nCorrectors    2;
    nNonOrthogonalCorrectors 1;
    rhoMin          1.2;
    rhoMax          1.2;
}

relaxationFactors
{
    equations
    {
        ".*"  1;
    }
}

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

Code:

/*--------------------------------*- 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    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default        backward;
}

gradSchemes
{
    default        cellLimited Gauss linear 1;
}

divSchemes
{
    default        none;
    div(phi,U)      Gauss filteredLinear;
    div(phi,e)      Gauss linear;
    div(phi,K)      Gauss linear;
    div(phiv,p)    Gauss filteredLinear;
    //div(phi,k)      Gauss linear;
    div(phi,B)      Gauss linear;
    div(phi,muTilda)Gauss linear;
    div(B)          Gauss linear;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default        Gauss linear corrected;
}

interpolationSchemes
{
    default        linear;
}

snGradSchemes
{
    default        corrected;
}


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

Code:

/*--------------------------------*- 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    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application    rhoPimpleFoam;

//startFrom      startTime;
startFrom    latestTime;

startTime      0;

//stopAt          writeNow;
stopAt        endTime;

endTime        5000000;

deltaT          1e-05;

writeControl    runTime;

writeInterval  4e-04;

purgeWrite      10;

writeFormat    binary;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision  6;

graphFormat    raw;

runTimeModifiable true;

libs ("libLEMOS-2.4.x.so");

functions
{

    /*probes
    {
 
        type probes;
        functionObjectLibs ("libsampling.so");
        enabled true;
        outputControl timeStep;
        outputInterval 1;

        probeLocations
        (
            (0.0738 0.0015 0.0)
            (0.0738 -0.0015 0.0)
        );

        fields
        ( 
            U
        );
    } */


    fieldAverage
        {
            type            fieldAverage;
            functionObjectLibs ("libfieldFunctionObjects.so");
            enabled        true;
            outputControl  outputTime;
       
            fields
            (
                    U
            {
                mean on;
                prime2Mean on;
                base time;
            }
                       
            nut
            {
                mean on;
                prime2Mean off;
                base time;
            }

              );
    }
}

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



All times are GMT -4. The time now is 05:37.