View Issue Details

Jump to NotesJump to History
IDProjectCategoryView StatusDate SubmittedLast Update
0004096OpenFOAMFeaturepublic2024-06-12 08:252024-06-12 09:20
ReporterRaghav Assigned Tohenry  
PriorityimmediateSeveritymajorReproducibilityunable to reproduce
Status closedResolutionno change required 
PlatformGNU/LinuxOSUbuntuOS Version22.04
Product Version11 
Summary0004096: How to simulate conduction and forced convection in a heater heatsink model?
DescriptionI have added my files below...
I tried to run a simulation using a heater heatsink model with an fluid enclosure.
But when the simulation was done I created a .foam file and opened it in paraview expecting to view a contour temperature distribution over the surface of the heatsink, I got it but according to it the temperature is uniformly at 297.25 K everywhere, which is not possible if conduction has occurred. Am I missing something in the boundary conditions or is it the wrong solver or anything other issue.. I don't know please help...
TagsNo tags attached.

Activities

Raghav

2024-06-12 08:25

reporter  

Allclean (169 bytes)    
#!/bin/sh
cd "${0%/*}" || exit 1

. $WM_PROJECT_DIR/bin/tools/CleanFunctions

cleanCase

#------------------------------------------------------------------------------
Allclean (169 bytes)   
Allrun (429 bytes)    
#!/bin/sh

cd ${0%/*} || exit 1

. $WM_PROJECT_DIR/bin/tools/RunFunctions

rm -rf constant/polyMesh/sets

#runApplication blockMesh
#runApplication snappyHexMesh -overwrite

ideasUnvToFoam ../heatsinkmesh2.unv
transformPoints "scale=(2 2 2)"
runApplication splitMeshRegions -cellZones -overwrite

runApplication decomposePar -allRegions
runParallel $(getApplication)
runApplication reconstructPar -allRegions
#paraFoam -touchAll
Allrun (429 bytes)   
snappyHexMeshDict (2,940 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    object      snappyHexMeshDict;
}

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

#includeEtc "caseDicts/mesh/generation/snappyHexMeshDict.cfg"

castellatedMesh true;
snap            true;
addLayers       false;

geometry
{
    heatedDuct
    {
        type triSurfaceMesh;
        file "heatedDuct.stl";

        regions
        {
            metalInlet     { name metalInlet;     }
            heaterInlet    { name heaterInlet;    }
            fluidInlet     { name fluidInlet;     }
            metalOutlet    { name metalOutlet;    }
            heaterOutlet   { name heaterOutlet;   }
            fluidOutlet    { name fluidOutlet;    }
            metalExternal  { name metalExternal;  }
            heaterExternal { name heaterExternal; }
        }
    }

    metalToHeater
    {
        type triSurfaceMesh;
        file "metalToHeater.stl";
    }

    fluidToMetal
    {
        type triSurfaceMesh;
        file "fluidToMetal.stl";
    }
};

castellatedMeshControls
{
    refinementSurfaces
    {
        heatedDuct
        {
            level (0 0);
            regions
            {
                metalInlet     { level (0 0); patchInfo { type patch; } }
                heaterInlet    { level (0 0); patchInfo { type patch; } }
                fluidInlet     { level (0 0); patchInfo { type patch; } }
                metalOutlet    { level (0 0); patchInfo { type patch; } }
                heaterOutlet   { level (0 0); patchInfo { type patch; } }
                fluidOutlet    { level (0 0); patchInfo { type patch; } }
                metalExternal  { level (1 1); patchInfo { type wall;  } }
                heaterExternal { level (1 1); patchInfo { type wall;  } }
            }
        }

        fluidToMetal
        {
            level (1 1);
            faceZone fluidToMetal;
            cellZone metal;
            mode        insidePoint;
            insidePoint (0.025 0.0025 0.075);
        }

        metalToHeater
        {
            level (1 1);
            faceZone metalToHeater;
            cellZone heater;
            mode        insidePoint;
            insidePoint (0.025 -0.005 0.075);
        }
    }

    nCellsBetweenLevels 1;

    insidePoint (0.025 0.025 0.075);
}

addLayersControls
{
    relativeSizes       true;
    minThickness        1;
    finalLayerThickness 1;
    expansionRatio      1;
    layers
    {}
}

// ************************************************************************* //
snappyHexMeshDict (2,940 bytes)   
fvSolution (683 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

PIMPLE
{
}

// ************************************************************************* //
fvSolution (683 bytes)   
decomposeParDict (1,034 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    note        "mesh decomposition control dictionary";
    location    "system";
    object      decomposeParDict;
}

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

numberOfSubdomains 8;

method          simple;

simpleCoeffs
{
    n           (3 1 2);
    delta	0.001
    order       xyz;
}
hierarchicalCoeffs
{
    n           (1 1 1);
    delta       0.001;
    order       xyz;
}
manualCoeffs
{
    dataFile    "";
}
distributed     no;

roots           ( );

// ************************************************************************* //
decomposeParDict (1,034 bytes)   
cellMaxHeaterTemp (722 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
-------------------------------------------------------------------------------
Description
    Writes out the maximum cell value for one or more fields.

\*---------------------------------------------------------------------------*/

#includeEtc "caseDicts/postProcessing/volFieldValue/cellMax.cfg"
region heater;
fields  (T);

// ************************************************************************* //
cellMaxHeaterTemp (722 bytes)   
controlDict (1,711 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     foamMultiRun;

regionSolvers
{
    fluid           fluid;
    heatsink        solid;
    heater          solid;
}

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         50;

deltaT          1;

writeControl    adjustableRunTime;

writeInterval   50;

purgeWrite      0;

writeFormat     binary;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

//maxCo           1.0;

//maxDi           10.0;

//adjustTimeStep  yes;
function
{
	#includeFunc cellMaxHeaterTemp
	#includeFunc cellMaxHeatsinkTemp
	#includeFunc volAvgHeatsinkTemp
	#includeFunc volAvgHeaterTemp
	
	averageTempHeaterHeatsinkInterface
	{
		type		surfaceFieldValue;
		surfaceFormat	off;
		libs		("libfieldFunctionObjects.so");
		enabled		true;
		writeControl	timeStep;
		writeInterval	1;
		log		true;
		writeFields	false;
		region		heater;
		regionType	patch;
		name		heater_to_heatsink;
		operation	areaAverage;
		mode		magnitude;
		fields
		(
			T
		);	
	}
}

// ************************************************************************* //
controlDict (1,711 bytes)   
volAvgHeaterTemp (729 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
-------------------------------------------------------------------------------
Description
    Writes out the volume-weighted average of one or more fields.

\*---------------------------------------------------------------------------*/

#includeEtc "caseDicts/postProcessing/volFieldValue/volAverage.cfg"
region heater;
fields  (T);

// ************************************************************************* //
volAvgHeaterTemp (729 bytes)   
cellMaxHeatsinkTemp (724 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
-------------------------------------------------------------------------------
Description
    Writes out the maximum cell value for one or more fields.

\*---------------------------------------------------------------------------*/

#includeEtc "caseDicts/postProcessing/volFieldValue/cellMax.cfg"
region heatsink;
fields  (T);

// ************************************************************************* //
cellMaxHeatsinkTemp (724 bytes)   
volAvgHeatsinkTemp (731 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
-------------------------------------------------------------------------------
Description
    Writes out the volume-weighted average of one or more fields.

\*---------------------------------------------------------------------------*/

#includeEtc "caseDicts/postProcessing/volFieldValue/volAverage.cfg"
region heatsink;
fields  (T);

// ************************************************************************* //
volAvgHeatsinkTemp (731 bytes)   
meshQualityDict (737 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system";
    object      meshQualityDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#includeEtc "caseDicts/mesh/generation/meshQualityDict.cfg"

// ************************************************************************* //
meshQualityDict (737 bytes)   
fvSolution-2 (950 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system/heater";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    e
    {
        solver           GAMG;
        smoother         symGaussSeidel;
        tolerance        1e-6;
        relTol           0.1;
    }

    eFinal
    {
        $e;
        relTol           0;
    }
}

PIMPLE
{
    nNonOrthogonalCorrectors 2;
}

// ************************************************************************* //
fvSolution-2 (950 bytes)   
decomposeParDict-2 (19 bytes)    
../decomposeParDict
decomposeParDict-2 (19 bytes)   
fvSchemes (989 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system/heater";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

// ************************************************************************* //
fvSchemes (989 bytes)   
fvSolution-3 (1,432 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system/fluid";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    rho
    {
        solver          diagonal;
    }

    rhoFinal
    {
        $rho;
    }

    p_rgh
    {
        solver           GAMG;
        smoother         symGaussSeidel;
        tolerance        1e-7;
        relTol           0.01;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol           0;
    }

    "(U|h)"
    {
        solver           PBiCGStab;
        preconditioner   DILU;
        tolerance        1e-7;
        relTol           0.1;
    }

    "(U|h)Final"
    {
        $U;
        relTol           0;
    }
}

PIMPLE
{
    momentumPredictor   yes;
    nNonOrthogonalCorrectors 2;
}

relaxationFactors
{
    equations
    {
        h               0.7;
        U               0.3;
        T		0.9;
    }
}

// ************************************************************************* //
fvSolution-3 (1,432 bytes)   
decomposeParDict-3 (19 bytes)    
../decomposeParDict
decomposeParDict-3 (19 bytes)   
fvSchemes-2 (1,145 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system/fluid";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;

    div(phi,U)      Gauss upwind;
    div(phi,K)      Gauss linear;
    div(phi,h)      Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

// ************************************************************************* //
fvSchemes-2 (1,145 bytes)   
fvConstraints (916 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    object      fvConstraints;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

/*energySource
{
    type            heatSource;

    select          all;

    q               200;
}*/
limitTemperature
{
    type	limitTemperature;
    active	true;
    limitTemperatureCoeffs
    {
        select		all;
        min		200;
        max		400;
    }
}


// ************************************************************************* //
fvConstraints (916 bytes)   
fvSolution-4 (982 bytes)    
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  11
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system/heater";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    e
    {
        solver           GAMG;
        smoother         symGaussSeidel;
        tolerance        1e-6;
        relTol           0.1;
    }

    eFinal
    {
        $h;
        relTol           0;
    }
}

PIMPLE
{
    nNonOrthogonalCorrectors 0;
    nNonOrthogonalCorrectors 2;
}

// ************************************************************************* //
fvSolution-4 (982 bytes)   
decomposeParDict-4 (19 bytes)    
../decomposeParDict
decomposeParDict-4 (19 bytes)   

Raghav

2024-06-12 08:28

reporter   ~0013258

This is the mesh file I am using for this simulation

henry

2024-06-12 08:47

manager   ~0013259

User support request

Issue History

Date Modified Username Field Change
2024-06-12 08:25 Raghav New Issue
2024-06-12 08:25 Raghav File Added: Allclean
2024-06-12 08:25 Raghav File Added: Allrun
2024-06-12 08:25 Raghav File Added: snappyHexMeshDict
2024-06-12 08:25 Raghav File Added: fvSolution
2024-06-12 08:25 Raghav File Added: decomposeParDict
2024-06-12 08:25 Raghav File Added: cellMaxHeaterTemp
2024-06-12 08:25 Raghav File Added: controlDict
2024-06-12 08:25 Raghav File Added: volAvgHeaterTemp
2024-06-12 08:25 Raghav File Added: cellMaxHeatsinkTemp
2024-06-12 08:25 Raghav File Added: volAvgHeatsinkTemp
2024-06-12 08:25 Raghav File Added: meshQualityDict
2024-06-12 08:25 Raghav File Added: fvSolution-2
2024-06-12 08:25 Raghav File Added: decomposeParDict-2
2024-06-12 08:25 Raghav File Added: fvSchemes
2024-06-12 08:25 Raghav File Added: fvSolution-3
2024-06-12 08:25 Raghav File Added: decomposeParDict-3
2024-06-12 08:25 Raghav File Added: fvSchemes-2
2024-06-12 08:25 Raghav File Added: fvConstraints
2024-06-12 08:25 Raghav File Added: fvSolution-4
2024-06-12 08:25 Raghav File Added: decomposeParDict-4
2024-06-12 08:28 Raghav Note Added: 0013258
2024-06-12 08:47 henry Assigned To => henry
2024-06-12 08:47 henry Status new => closed
2024-06-12 08:47 henry Resolution open => no change required
2024-06-12 08:47 henry Note Added: 0013259