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

Fluid properties as function of Temperature with polynomial using chtMultiRegionFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 24, 2018, 05:40
Default Fluid properties as function of Temperature with polynomial using chtMultiRegionFoam
  #1
New Member
 
Quentin
Join Date: Aug 2018
Posts: 1
Rep Power: 0
qferon is on a distinguished road
Hello everyone,

I'm trying to calculate natural convection in a small slot surrounded by a solid with temperature from 110 K and 173 K. I'm using chtMultiRegionFaom to solve it.

When I'm using constant properties for the fluid everything is doing good, but as soon as try to implement polynomial transport properties, my case run few iterations and then crash with the following message (see picture attached below).

I try to run the simulation with deltaT equals to 5e-04 but it still crash. I'm trying with delta T equals to 5e-08 but it takes a very long time. My checkMesh is OK. There is my fvScheme and fvSolution for the fluid attached below.

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

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;

div(phi,U) Gauss linearUpwind grad(U);
div(phi,K) Gauss linearUpwind grad(U);
div(phi,h) Gauss linearUpwind grad(U);
div(phi,k) Gauss linearUpwind grad(U);
div(phi,epsilon) Gauss linearUpwind grad(U);
div(phi,omega) Gauss linearUpwind grad(U);
div(phi,R) Gauss linearUpwind grad(U);
div(R) Gauss linearUpwind grad(U);
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}
// ************************************************** *********************** //


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

solvers
{
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-6;
relTol 0.1;
}

rhoFinal
{
$rho;
tolerance 1e-5;
relTol 0;
}

p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-3;
relTol 0.01;
maxIter 2000;
}

p_rghFinal
{
$p_rgh;
tolerance 1e-3;
relTol 0;
}

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

"(U|h|k|epsilon|omega|R)Final"
{
$U;
tolerance 1e-7;
relTol 0;
}
}

PIMPLE
{
momentumPredictor yes;
nCorrectors 10;
nOuterCorrectors 10;
nNonOrthogonalCorrectors 1;
}

relaxationFactors
{
equations
{
"h.*" 0.95;
"p_rgh.*" 0.7;
"rho.*" 0.7;
"U.*" 1;
"omega.*" 0.7;
"epsilon.*" 0.7;
}
}

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



Thanks for everyone that can help !
Quentin
Attached Images
File Type: jpg screenshot_divergence.jpg (130.2 KB, 6 views)
qferon is offline   Reply With Quote

Old   November 3, 2018, 07:14
Default Loved It
  #2
New Member
 
James Carrol
Join Date: Nov 2018
Posts: 1
Rep Power: 0
Dembir is on a distinguished road
I am a coder and dont love the codding while at job for work but I found few error in your codding though which is minnor and could be resolve soon as currently we are essay writers for hire to get the job done for the work and surely you will enjoy the same essay and assignment for college student
Dembir is offline   Reply With Quote

Old   November 3, 2018, 14:41
Default
  #3
Senior Member
 
Join Date: Sep 2013
Posts: 166
Rep Power: 9
Bloerb is on a distinguished road
Your settings are no good. A final pressure residual of 1e-3 does not suffice for heat transfer. And you are using linearUpwind on every scheme which is no good for complex heat transfer simulations. In addition you are relaxing the final iteration even though this is an unsteady simulation.
Bloerb is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, polynomial properties

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 07:20
Setting the height of the stream in the free channel kevinmccartin CFX 10 July 9, 2015 22:36
Temperature Dependent Fluid Properties eliwide OpenFOAM Running, Solving & CFD 4 December 15, 2012 20:21
error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Native Meshers: blockMesh 2 March 14, 2012 10:56
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11


All times are GMT -4. The time now is 22:46.