CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Temperqture field not solved in LES - heat transfer - rhopimplefoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 24, 2021, 15:13
Default Temperqture field not solved in LES - heat transfer - rhopimplefoam
  #1
Member
 
Arun subramanian
Join Date: Jun 2016
Location: Florence,Italy
Posts: 48
Rep Power: 8
arun1994 is on a distinguished road
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;
            }

               );
    }
}

// ************************************************************************* //
arun1994 is offline   Reply With Quote

Reply

Tags
heat transfer, les, openfoam 7, rhopimplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Solvers in OpenFOAM for LES + heat transfer arun1994 Main CFD Forum 1 November 26, 2021 07:57
Heat transfer coefficient in regenerator gartz89 FLUENT 0 May 26, 2018 01:23
Problem with total heat transfer rate aswathy_raghu FLUENT 0 July 26, 2016 07:39
Solver for Incompressible flow with heat transfer using LES Tarak OpenFOAM 2 December 20, 2013 08:43
Heat transfer from a heated plate using fins pathakamit FLUENT 1 April 30, 2013 04:07


All times are GMT -4. The time now is 11:47.