Not obtaining temperature contours
 

Hello Guys,


I am quite new to openfoam and was trying a turbulent heat transfer case. I edited the blockmesh file for the tutorials>heatTransfer>buoyantPimpleFoam>hotRoomBo ussinesq. Introduced a cubical block with a hot wall (673K) on the left side and a cold wall (193K) on the right side wall. when i ran the case through openFoam for 10000 iterations it is not writing a folder at 100 interval mark which i gace in the controlDict and also no heat transfer is seen between the hot and cold wall.


WHAT AM I DOING WRONG :confused: ..... Pls Help

Hi Amal,

Does the solver run properly? Can you post the log files here please?

Thanks,
Yann

log files
 

these are the log files...... I blocked the log.setFields file from being written.


the other two log files are here


log.blockMesh
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 7-1ff648926f77
Exec : blockMesh
Date : Apr 07 2022
Time : 11:13:08
Host : "chummar-GL553VE"
PID : 4186
I/O : uncollated
Case : /home/chummar/OpenFOAM/chummar-7/run/realone/turbtrial
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Creating block mesh from
"/home/chummar/OpenFOAM/chummar-7/run/realone/turbtrial/system/blockMeshDict"
Creating block edges
No non-planar block faces defined
Creating topology blocks
Creating topology patches

Creating block mesh topology

Check topology

Basic statistics
Number of internal faces : 0
Number of boundary faces : 6
Number of defined boundary faces : 6
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1
Block 0 cell size :
i : 0.5
j : 1
k : 0.5

Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (0 0 0) (10 10 10)
nPoints: 4851
nCells: 4000
nFaces: 12800
nInternalFaces: 11200
----------------
Patches
----------------
patch 0 (start: 11200 size: 400) name: floor
patch 1 (start: 11600 size: 400) name: ceiling
patch 2 (start: 12000 size: 200) name: hotwall
patch 3 (start: 12200 size: 200) name: coldwall
patch 4 (start: 12400 size: 400) name: fixedWalls

End



log.bouyantPimpleFoam


/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 7-1ff648926f77
Exec : buoyantPimpleFoam
Date : Apr 07 2022
Time : 11:13:09
Host : "chummar-GL553VE"
PID : 4187
I/O : uncollated
Case : /home/chummar/OpenFOAM/chummar-7/run/realone/turbtrial
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0



--> FOAM FATAL IO ERROR:
keyword PIMPLE is undefined in dictionary "/home/chummar/OpenFOAM/chummar-7/run/realone/turbtrial/system/fvSolution"

file: /home/chummar/OpenFOAM/chummar-7/run/realone/turbtrial/system/fvSolution from line 22 to line 65.

From function const Foam::dictionary& Foam::dictionary::subDict(const Foam::word&) const
in file db/dictionary/dictionary.C at line 708.

FOAM exiting





THANKS:)

Thank you Amal,

Have a look at your buoyantPimpleFoam log file. The solver did not run because it misses the PIMPLE dictionary in the fvSolution file.

You have to fix this error.

Yann

fvSolution and fvSchemes
 

in the controlDict bouyantPimpleFoam is called but


fvSolution


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

solvers
{
p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-08;
relTol 0.01;
}

"(U|e|k|epsilon)"
{
solver PBiCGStab;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;

pRefCell 0;
pRefValue 0;

residualControl
{
p_rgh 1e-2;
U 1e-4;
e 1e-2;

// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}

relaxationFactors
{
fields
{
p_rgh 0.7;
}
equations
{
U 0.3;
e 0.5;
"(k|epsilon|R)" 0.7;
}
}


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






fvSchemes


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

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;

div(phi,U) bounded Gauss upwind;
div(phi,e) bounded Gauss upwind;

div(phi,k) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;

div(phi,Ekp) bounded Gauss linear;

div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}


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


Is this cause of mentioning simple in fvSolution instead of pimple :confused:

Yes, you are probably using a fvSolution file from a steady case (like buoyantSimpleFoam).
If your case is based on the hotRoomBoussinesq tutorial (tutorials>heatTransfer>buoyantPimpleFoam>hotRoomB oussinesq), try to replace your fvSolution file by the one from the tutorial.

Yann