[silviliril]

I have a solver based on OpenFOAM-2.2.0 (link). that can simulate boiling and condensation. The thermo-physical properties can be written in "transportProperties" file as shown below:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.6                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Numerical simulation of bubble growth in film boiling using a coupled
//level-set and volume-of-fluid method
//Autors: G. Tomar and etc

// for critical water at near critical pressure Pr=0.99



phaseChangeTwoPhaseMixture Hardt;

pSat             pSat       [1 -1 -2 0 0 0 0]   1E5;
       			 	
TSat		TSat	    [0 0 0 1 0]	373.15;//add saturation temprature
R               R   [0 2 -2 -1 0 0 0]     461.0; //add
TSatLocal       false;

Hfg		Hfg	[0 2 -2 0 0] 2257000;//add  J/kg

sigma           sigma [1 0 -2 0 0 0 0] 0;

phase1
{
    transportModel Newtonian;
    nu              nu  [0 2 -1 0 0 0 0]    2.94e-7; //46.7e-6/402.4
    rho             rho [1 -3 0 0 0 0 0]    958.35;
    C		    C   [0 2 -2 -1 0 0 0]   4217;//add
    K		    K   [1 1 -3 -1 0 0 0]   0.6791;//add
}

phase2
{
    transportModel Newtonian;
    nu              nu  [0 2 -1 0 0 0 0]   20.51e-6;//1.33e-07;//32.38e-6/242.7
    rho             rho [1 -3 0 0 0 0 0]   0.598141;
    C		    C   [0 2 -2 -1 0 0 0]  2077;//look
    K		    K   [1 1 -3 -1 0 0 0]  0.0251;//add
}
HardtCoeffs
{
   Cv            Cv   [0 0 0 0 0 0 0]    1; 
   rv            rv   [0 0 0 0 0 0 0]    1;
   rc            rc   [0 0 0 0 0 0 0]    0; 
}

LeeCoeffs
{
  r            r   [0 0 -1 0 0 0 0]    1000;//100; //add
}
FourierCoeffs {}
// ************************************************************************* //

As, you can see, the constant properties can be easily supplied in "transportProperties" folder. But, I want to supply variable properties as function of temperature in polynomial equation. such as

Surface tension as function of temperature:
σ = 0.09805856 - 1.845 * 10^-5*T - 2.3 * 10^-7*T^2

and density of phase 1:
ρ = 859.0083 + 1.252209*T + 0.0026429*T^2

Please help me, how can I write any property as polynomial function of temperature.

[mldrmz]

Hello
could you apply your changes on sigma calculation? would you please share your achievement?