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

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の公式ドキュメントやユーザーガイドを参照してください。