物性値の動粘度や密度は、温度の関数として与える方法

OpenFOAMでは物性値（例えば動粘度や密度）を温度の関数として設定することが可能です。これを実現するために、constantディレクトリのtransportPropertiesファイルで関数形式のプロパティを定義します。具体的には、temperatureFunctionやpolynomialといったモデルを使用することができます。

以下に、温度の関数として動粘度と密度を設定する例を示します。

1. 動粘度を温度の関数として設定する例
例えば、動粘度を温度 T の関数として設定する場合、次のようにtransportPropertiesファイルを編集します。

plaintext
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*  See: https://www.openfoam.com/documentation/source-code/documentation/    */
/*---------------------------------------------------------------------------*/

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}

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

transportModel  Newtonian;

nu              nu [0 2 -1 0 0 0 0] polynomial
(
    (1.0e-06)  // Coefficient for T^0
    (2.0e-08)  // Coefficient for T^1
    (-1.0e-10) // Coefficient for T^2
);

rho             rho [1 -3 0 0 0 0 0] 998.2;  // 密度 (kg/m^3)

// ************************************************************************* //
Public code references from 5 repositories
2. 密度を温度の関数として設定する例
密度を温度の関数として設定する場合も同様に、transportPropertiesファイルを編集します。

plaintext
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*  See: https://www.openfoam.com/documentation/source-code/documentation/    */
/*---------------------------------------------------------------------------*/

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}

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

transportModel  Newtonian;

nu              nu [0 2 -1 0 0 0 0] 1.0e-06;  // 動粘度 (m^2/s)

rho             rho [1 -3 0 0 0 0 0] polynomial
(
    (998.2)     // Coefficient for T^0
    (-0.2)      // Coefficient for T^1
    (0.001)     // Coefficient for T^2
);

// ************************************************************************* //
Public code references from 5 repositories
3. 物性値の関数モデル
上記の例では多項式モデルを使用しましたが、他にも linear や exponential などの関数モデルが利用できます。詳細な設定方法については、OpenFOAMの公式ドキュメントやユーザーガイドを参照してください。