[velikaDjoka]

Hello everyone,

Does anybody already encoded in openfoam the thermophysical properties of the Ti6Al4V?

I have done it but it does not seem correct. As I launched a simulation with the default material in openfaom, the simulation worked . While my new material ( Ti6Al4V), the simulation crashes immediately ( the diffusion number is negative ).

Here the thermophysical properties that I have written :

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            heSolidThermo; // properties for solid
    mixture         pureMixture;
    transport       polynomial;// constIso
    thermo          hPolynomial; //hConst
    equationOfState icoPolynomial; // rhoConst
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    // ############ Properties for Ti6Al4V ##############
    // ATTENTION: values are mainly correct from 275K -> 1600K
    
    specie
    {
        molWeight   446.0697; // [g/mol] => ex: carbon = 12, copper = 63.546 
    }

    transport
    {
        muCoeffs<8>     ( 0 0 0 0 0 0 0 0 ); //[Pa s]
        kappaCoeffs<8>  ( -4.561 0.09348 -0.0002819 3.932e-07 -2.226e-10 3.478e-14 4.816e-18 0 ); // [W/mK]
    }

    thermodynamics
    {
        Hf              0; // no reactions => 0
        Sf              0; // no reactions => 0
        CpCoeffs<8>     ( -800.1667 12.7983 -0.0474 -0.0474 -8.8041e-08 4.7729e-11 -1.3240e-14 1.4709e-18 ); // [J/kgK]
    }

    equationOfState
    {
        rhoCoeffs<8>    ( 4.5243e+03 -0.4901 3.9646e-04 6.2064e-07 -1.5296e-09 1.0303e-12 -2.2831e-16 0 ); // [kg/m3]

    }
}

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