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

tutorial>heatTransfer>buouyantPimpleFoam>hotRoomBo ussinesq

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

Like Tree3Likes
  • 1 Post By Yann
  • 1 Post By Yann
  • 1 Post By Yann

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 25, 2022, 02:42
Question tutorial>heatTransfer>buouyantPimpleFoam>hotRoomBo ussinesq
  #1
Member
 
Amal Chummar
Join Date: May 2021
Posts: 31
Rep Power: 3
CHUMMAR is on a distinguished road
Hi .... I am new to openFOAM and i wanted to run a basic turbulent case using the hotroomboussinesq tutorial...... i wanted to implement that same turbulence properties into a blockmesh created by me. Is this possible. The files i changed for this were the '0' folder files, and the blockMeshDict file as :


scale 1;

vertices
(
(-0.2 0 -0.1)
( 1 0 -0.1)
( 1 1 -0.1)
(-0.2 1 -0.1)
(-0.2 0 0.1)
( 1 0 0.1)
( 1 1 0.1)
(-0.2 1 0.1)
);

blocks
(
hex (0 1 2 3 4 5 6 7) (48 40 8) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
top
{
type wall;
faces
(
(3 7 6 2)
);
}
hotwall
{
type patch;
faces
(
(0 4 7 3)
);
}
coldwall
{
type patch;
faces
(
(2 6 5 1)
);
}
bottom
{
type wall;
faces
(
(1 5 4 0)
);
}
back
{
type wall;
faces
(
(0 3 2 1)
);
}
front
{
type wall;
faces
(
(4 5 6 7)
);
}
);

mergePatchPairs
(
);



..............................................
from this :


convertToMeters 1;

vertices
(
(0 0 0)
(10 0 0)
(10 5 0)
(0 5 0)
(0 0 10)
(10 0 10)
(10 5 10)
(0 5 10)
);

blocks
(
hex (0 1 2 3 4 5 6 7) (20 10 20) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
floor
{
type wall;
faces
(
(1 5 4 0)
);
}
ceiling
{
type wall;
faces
(
(3 7 6 2)
);
}
fixedWalls
{
type wall;
faces
(
(0 4 7 3)
(2 6 5 1)
(0 3 2 1)
(4 5 6 7)
);
}
);

mergePatchPairs
(
);




Thank you
CHUMMAR is offline   Reply With Quote

Old   June 28, 2022, 05:46
Smile I used salome instead
  #2
Member
 
Amal Chummar
Join Date: May 2021
Posts: 31
Rep Power: 3
CHUMMAR is on a distinguished road
i made the polymesh file using salome-ideasUnvToFoam-checkMesh and got mesh non-orthonogality as 56.38. i am trying to obtain a turbulent analysis using buoyantPimpleFoam from the hotRoomBousinessq tutorial.


now after running pimpleFoam i am getting such an error message. Pls help
Attached Images
File Type: jpg Error Message.jpg (69.8 KB, 6 views)
CHUMMAR is offline   Reply With Quote

Old   June 29, 2022, 04:26
Default
  #3
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 563
Rep Power: 19
Yann will become famous soon enough
Hi Amal,

You have a fatal error stating the keyword div(((rho*nuEff)*dev2(T(grad(U))))) is missing in system/fvShemes. You need to add something like this in the divSchemes in the fvShemes file:

Code:
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
Cheers,
Yann
CHUMMAR likes this.
Yann is offline   Reply With Quote

Old   June 30, 2022, 01:52
Exclamation Got the case open in paraFoam after creating a log file but not writing cases
  #4
Member
 
Amal Chummar
Join Date: May 2021
Posts: 31
Rep Power: 3
CHUMMAR is on a distinguished road
after using salome and finally giving checkMesh.......i edited the Allrun file and '#' out the blockMesh creation and straightly went for calling the application bouyantPimpleFoam. A log file is now gettting written but the case is not writing any time intervals. only the 0 folder is present and the logfile gets created after i give Allrun. Even if i edit the controlDict file for 1000 iterations with 100 writing interval nothing gets written . What is the issue now???



log file:
/*---------------------------------------------------------------------------*\
========= |
\\ / 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 : Jun 30 2022
Time : 08:28:22
Host : "chummar-GL553VE"
PID : 4565
I/O : uncollated
Case : /home/chummar/OpenFOAM/chummar-7/run/hotRoomBq
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


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Reading thermophysical properties

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo eConst;
equationOfState Boussinesq;
specie specie;
energy sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
RASModel kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}


Reading g

Reading hRef
Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

Radiation model not active: radiationProperties not found
Selecting radiationModel none
No finite volume options present
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Time = 2

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 9.87889e-07, No Iterations 479


controlDict File:
/*--------------------------------*- 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 controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application buoyantPimpleFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 1000;

deltaT 2;

writeControl adjustableRunTime;

writeInterval 100;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep no;

maxCo 0.5;

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


thanks
CHUMMAR is offline   Reply With Quote

Old   June 30, 2022, 04:09
Default
  #5
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 563
Rep Power: 19
Yann will become famous soon enough
Hi,

Can you post your full log file for buoyantPimpleFoam please?

Yann
CHUMMAR likes this.
Yann is offline   Reply With Quote

Old   July 3, 2022, 23:52
Question Not getting Ux, Uy in graphs
  #6
Member
 
Amal Chummar
Join Date: May 2021
Posts: 31
Rep Power: 3
CHUMMAR is on a distinguished road
Hi Yann,
I am sorry for the late reply. I was actually trying to edit the hotRoomBoussinesq case within the tutorials>heatTransfer>BuoyantPimpleFoam directory to view a turbulent heat transfer case within a differentially heated cube. I tried all things but finally when i simply ran the hotRoomBoussinesq tutorial case it was not writing case files at required intervals itself without any of my modifications. I then tried the hotRoom case and it was writing case files at correct intervals properly. So i then decided to edit the hotRoom case instead. Now i got it running and was able to see results in parafoam after running for 20000 iterations.



My problem now is i am able to get the temperature variation but i thought there would be variation for Ux and Uy within the differentially heated cubical enclosure but instead its showing zero. Do u know why?



I also have another very basic doubt.....When i run blockMeshDict and then checkMesh the mesh seems too coarse. how can i change it to a more medium or fine mesh.
I tried with salome>ideasUnvtoFoam and got a very good mesh but dont knw how to implement it into this case neglecting blockMesh.


Thankyou
Attached Images
File Type: jpg ux error.jpg (96.0 KB, 3 views)
File Type: jpg Salome mesh quality.jpg (73.1 KB, 1 views)
CHUMMAR is offline   Reply With Quote

Old   July 4, 2022, 05:13
Default
  #7
Senior Member
 
Yann
Join Date: Apr 2012
Location: France
Posts: 563
Rep Power: 19
Yann will become famous soon enough
Hi Amal,

Quote:
Originally Posted by CHUMMAR View Post
My problem now is i am able to get the temperature variation but i thought there would be variation for Ux and Uy within the differentially heated cubical enclosure but instead its showing zero. Do u know why?
I have no idea about what and how you modified the original case so I cannot really help if you don't share more details about the case you have run. (eg. setup files, log files, etc...)

Quote:
Originally Posted by CHUMMAR View Post
I also have another very basic doubt.....When i run blockMeshDict and then checkMesh the mesh seems too coarse. how can i change it to a more medium or fine mesh.
I tried with salome>ideasUnvtoFoam and got a very good mesh but dont knw how to implement it into this case neglecting blockMesh.
If you only want to refine the mesh, you can just modify the block(s) definition in blockMeshDict to had more cells: https://www.openfoam.com/documentati...ckmesh-utility
CHUMMAR likes this.
Yann is offline   Reply With Quote

Reply

Tags
heat transfer, pimple. openfoam

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



All times are GMT -4. The time now is 19:53.