関数IAPWS95_forward_calculation マニュアル

(The documentation of function IAPWS95_forward_calculation)

Last Update: 2023/11/29


◆機能・用途(Purpose)

Wagner and Pruss (2002)の表6.3の計算式を用いて 指定された密度・温度における熱力学変数を計算する。
Compute thermodynamic properties under given density and temperature using formulas in Table 6.3 of Wagner and Pruss (2002).


◆形式(Format)

#include <IAPWS95/forward.h>
inline struct IAPWS95_state IAPWS95_forward_calculation
(const double rho,const double T)


◆引数(Arguments)

rho 密度ρ [kg m3]。
The density ρ [kg m3].
T 温度T [K]。
The temperature T [K].


◆戻り値(Return value)

引数で指定した密度・温度における熱力学変数の値。 戻り値のメンバの値は以下のようになる。
The values of thermodynamic variables at the density and temperature specified by the arguments. The values of the members of the return value are as follows.

メンバ
Member

Quantity
計算式(Wagner and Pruss (2002)の表6.3に基づく)
Formula (from Table 6.3 of Wagner and Pruss (2002))
p 圧力 [Pa]
Pressure [Pa]
ρRT(1+δϕδr)
s 単位質量あたりのエントロピー [J kg1 K1]
Entropy per unit mass [J kg1 K1]
R[τ(ϕτo+ϕτr)ϕoϕr]
u 単位質量あたりの内部エネルギー [J kg1]
Internal energy per unit mass [J kg1]
RTτ(ϕτo+ϕτr)
h 単位質量あたりのエンタルピー [J kg1]
Enthalpy per unit mass [J kg1]
RT[1+τ(ϕτo+ϕτr)+δϕδr]
f 単位質量あたりのヘルムホルツの自由エネルギー [J kg1]
Helmholtz free energy per unit mass [J kg1]
RT(ϕo+ϕr)
g 単位質量あたりのギブスの自由エネルギー [J kg1]
Gibbs free energy per unit mass [J kg1]
RT(1+ϕo+ϕr+δϕδr)
cv 定積比熱 [J kg1 K1]
Isochoric heat capacity [J kg1 K1]
Rτ2(ϕττo+ϕττr)
cp 定圧比熱 [J kg1 K1]
Isobaric heat capacity [J kg1 K1]
Rτ2(ϕττo+ϕττr)+R(1+δϕδrδτϕδτr)21+2δϕδr+δ2ϕδδr
w 音速 [m s1]
Sound speed [m s1]
RT[1+2δϕδr+δ2ϕδδr(1+δϕδrδτϕδτr)2τ2(ϕττo+ϕττr)]
dp_drho_fixT (pρ)T RT(1+2δϕδr+δ2ϕδδr)

ここでδ=ρ/ρcは無次元化した密度、 τ=Tc/Tは無次元化した温度の逆数、 ρcは臨界点での密度、Tcは臨界点での温度、 Rは単位質量あたり気体定数である。 ϕoは理想気体項、ϕrは残差項であり、 これらの関数やその導関数(下付き添字で微分する変数を表す)は関数 IAPWS95_calculate_idealIAPWS95_calculate_residual を内部で用いて計算している。
Here, δ=ρ/ρc is a nondimensional density, τ=Tc/T is the inverse of a nondimensional temperature, ρc and Tc are the density and temperature at the critical point, respectively, and R is the gas constant per unit mass; ϕo and ϕr are the ideal and residual parts, respectively, which and whose derivatives (where the subscripts show the variables to differentiate) are computed by internally calling functions IAPWS95_calculate_ideal and IAPWS95_calculate_residual.


◆使用例(Example)

const double rho=20.0;
const double T=600.0;
struct IAPWS95_state values=IAPWS95_forward_calculation(rho,T);
printf("Pressure = %e [Pa]\n",values.p);
printf("Entropy = %e [J/(kg K)]\n",values.s);