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

icoPoly8ThermoPhysics limits density below 2 kg/m3

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 16, 2012, 05:28
Default icoPoly8ThermoPhysics limits density below 2 kg/m3
  #1
New Member
 
Sinisa Majer
Join Date: Oct 2009
Location: Zagreb, CROATIA
Posts: 8
Rep Power: 7
smajer is on a distinguished road
Hello all,

I am trying to use chtMultiRegionSimpleFoam with liquid (water). Since it is important to take in account buoyant forces I have specified thermophysical properties as:
thermoType hRhoThermo<pureMixture<icoPoly8ThermoPhysics>>;

Unfortunately, it seems that a density is limited to 2 kg/m3.

To demonstrate this, I have changed the thermophysicalProperties file in multiRegionHeater case (tutorial for chtMultiRegionSimpleFoam).
For density (in this case constant, but it is just to demonstrate the problem) lower then 2 solver is using the value from thermophysicalProperties file. For density above 2 solver is fixing density to 2.
--------------
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant/topAir";
object thermophysicalProperties;
}

thermoType hRhoThermo<pureMixture<icoPoly8ThermoPhysics>>;

mixture
{
specie
{
nMoles 1;
molWeight 28.9;
}
equationOfState
{
// rhoCoeffs<8> ( 1.999 0 0 0 0 0 0 0 ); //If you use this, density is constant but lower then 2
rhoCoeffs<8> ( 2.1 0 0 0 0 0 0 0 ); //If you use this, density is constant and always 2 (it is same even if you put other coefficients for polynom)
}
thermodynamics
{
Hf 0;
Sf 0;
CpCoeffs<8> ( 1000 0 0 0 0 0 0 0 );
}
transport
{
muCoeffs<8> ( 1.8e-05 0 0 0 0 0 0 0 );
kappaCoeffs<8> ( 0.0242 0 0 0 0 0 0 0 );
Pr 0.7;
}
}

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

---------------


Why is that so? In OF documentation (http://www.openfoam.com/features/thermophysical.php) we can find this statement:
"icoPolynomial - Incompressible polynomial equation of state, e.g. for liquids"

Can someone explain this?

Best regards,

Sinisa
smajer is offline   Reply With Quote

Old   March 16, 2012, 06:41
Default
  #2
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 154
Rep Power: 7
Aurelien Thinat is on a distinguished road
Hi,

Check out your fvSolution file. You may have limited the density value range.
Aurelien Thinat is offline   Reply With Quote

Old   March 16, 2012, 07:11
Default
  #3
New Member
 
Sinisa Majer
Join Date: Oct 2009
Location: Zagreb, CROATIA
Posts: 8
Rep Power: 7
smajer is on a distinguished road
Well, that was quick and efficient!!
Thank you very much!!
That was the problem. I searched everything before posting here, but I didn't notice that density is limited in fvSolution.

Once again,

thanks!

Sinisa
smajer is offline   Reply With Quote

Reply

Tags
chtmultiregionsimplefoam, constant rho, icopoly8thermophysics, liquid

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
A problem about density in liquid air definition alloveyou CFX 2 June 14, 2012 14:20
Variable Density Function ryzd FLUENT 1 August 25, 2011 14:16
REAL GAS UDF brian FLUENT 6 September 11, 2006 08:23
Warning 097- AB CD-adapco 6 November 15, 2004 05:41
variable density water Atit Koonsrisuk CFX 2 July 24, 2003 03:07


All times are GMT -4. The time now is 16:23.