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 Search this Thread 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, 20 views)
qferon is offline   Reply With Quote

Old   November 3, 2018, 14:41
Default
  #2
Senior Member
 
Join Date: Sep 2013
Posts: 353
Rep Power: 20
Bloerb will become famous soon enough
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 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
Setting the height of the stream in the free channel kevinmccartin CFX 12 October 13, 2022 22:43
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 07:20
Temperature Dependent Fluid Properties eliwide OpenFOAM Running, Solving & CFD 4 December 15, 2012 20:21
[blockMesh] error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Meshing & Mesh Conversion 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:02.