IAPWS95/viscosity.h マニュアル

(The documentation of IAPWS95/viscosity.h)
Last Update: 2023/9/22


IAPWS95/viscosity.hでは 粘性係数を計算する関数が定義されている。
Functions to compute the viscosity are defined in IAPWS95/viscosity.h.
粘性係数の計算式はWagner and Pruss (2002)には含まれておらず、 Huber et al. (2009)において別途定義されている。 2023年9月21日現在、 NIST Chemistry WebBook において同論文の式が用いられている。
The equation of the viscosity is not included in Wagner and Pruss (2002) but is defined separately in Huber et al. (2009). The equation in that paper is used in NIST Chemistry WebBook as of 21 September 2023.
Huber et al. (2009)において粘性係数\(\mu\)は \[\begin{equation} \frac{\mu}{\mu_{ref}} =\bar{\mu}_0(\tau) \bar{\mu}_1(\delta,\tau) \bar{\mu}_2(\delta,\tau) \label{eq.mu} \end{equation}\] の形で与えられる(同論文(2)式)。 ここで\(\mu_{ref}\equiv 1\times 10^{-6}\) [Pa s]は基準となる粘性係数を表す。 また\(\delta\equiv\rho/\rho_c\)は無次元化した密度、 \(\tau\equiv T_c/T\)は無次元化した温度の逆数であり、 これらはIAPWS-95全体の表記に合わせてある (そのためHuber et al. (2009)の表記とは異なる)。
Huber et al. (2009) gives the viscosity \(\mu\) in the form of Eq. (\ref{eq.mu}), which is from Eq. (2) of the paper, where \(\mu_{ref}\equiv 1\times 10^{-6}\) [Pa s] is a reference viscosity; \(\delta\equiv\rho/\rho_c\) is a nondimensional density and \(\tau\equiv T_c/T\) is the inverse of a nondimensional temperature, defined consistently with the notation of IAPWS-95. Note that the notation is different from that of Huber et al. (2009).
\(\bar{\mu}_0\)は低密度極限での粘性係数を表し、 \[\begin{equation} \bar{\mu}_0(\tau)= \frac{100}{\sqrt{\tau}\sum_{i=0}^3 H_i \tau^i} \label{eq.mu0} \end{equation}\] で与えられる(論文(11)式)。ここで\(H_i\)は表1のように与えられる。
The factor \(\bar{\mu}_0\) represents the viscosity in the limit of zero density and is given as Eq. (\ref{eq.mu0}), which is from Eq. (11) of the paper. The values of \(H_i\) are given in Table 1.
表1. (\ref{eq.mu0})式の係数\(H_i\)。 Huber et al. (2009)の表2に基づく。
Table 1. The coefficients \(H_i\) of Eq. (\ref{eq.mu0}) from Table 2 of Huber et al. (2009).
\(i\) \(H_i\)
0 1.67752
1 2.20462
2 0.6366564
3 -0.241605

\(\bar{\mu}_1\)は密度増加の粘性係数への寄与を表し、 \[\begin{equation} \bar{\mu}_1(\delta,\tau)= \exp\left[ \delta \sum_{i=0}^5 (\tau-1)^i \sum_{j=0}^6 H_{ij} (\delta-1)^j \right] \label{eq.mu1} \end{equation}\] で与えられる(同論文(12)式)。 ここで\(H_{ij}\)は表2のように与えられる。
The factor \(\bar{\mu}_1\) represents the contribution to the viscosity due to increasing density and is given as Eq. (\ref{eq.mu1}) which is from Eq. (12) of the paper. The values of \(H_{ij}\) are given in Table 2.
表2. (\ref{eq.mu1})式の係数\(H_{ij}\)。 Huber et al. (2009)の表3に基づく。
Table 2. The coefficients \(H_{ij}\) of Eq. (\ref{eq.mu1}) from Table 3 of Huber et al. (2009).
\(i=0\) \(i=1\) \(i=2\) \(i=3\) \(i=4\) \(i=5\)
\(j=0\) 5.20094e-01 8.50895e-02 -1.08374 -2.89555e-01 0 0
\(j=1\) 2.22531e-01 9.99115e-01 1.88797 1.26613 0 1.20573e-01
\(j=2\) -2.81378e-01 -9.06851e-01 -7.72479e-01 -4.89837e-01 -2.57040e-01 0
\(j=3\) 1.61913e-01 2.57399e-01 0 0 0 0
\(j=4\) -3.25372e-02 0 0 6.98452e-02 0 0
\(j=5\) 0 0 0 0 8.72102e-03 0
\(j=6\) 0 0 0 -4.35673e-03 0 -5.93264e-04

\(\bar{\mu}_2\)は臨界点近傍での振る舞いを表し、 論文(33)式, (20)-(23)式, (26)-(28)式, および(34)式の直後の注釈により \[\begin{equation} \bar{\mu}_2(\delta,\tau)=\exp\left[\chi_{\mu} Y(\delta,\tau)\right] \label{eq.mu2} \end{equation}\] \[\begin{equation} Y(\delta,\tau)= \begin{cases} Y_1(\delta,\tau) & (\xi\leq 0.3817016416\times 10^{-9}\mbox{ [m]}) \\ Y_2(\delta,\tau) & (\xi>0.3817016416\times 10^{-9}\mbox{ [m]}) \end{cases} \label{eq.Y} \end{equation}\] \[\begin{equation} Y_1(\delta,\tau)= \frac{1}{5}q_Cq_D^5\xi(\delta,\tau)^6 \left[1-q_C\xi(\delta,\tau)+q_C^2\xi(\delta,\tau)^2 -\frac{765}{504}q_D^2\xi(\delta,\tau)^2\right] \label{eq.Y1} \end{equation}\] \[\begin{eqnarray} Y_2(\delta,\tau) &=& \frac{1}{12}\sin\left[3\psi_D(\delta,\tau)\right] -\frac{1}{4q_C\xi(\delta,\tau)}\sin\left[2\psi_D(\delta,\tau)\right] \nonumber \\ & & +\frac{1}{q_C^2\xi(\delta,\tau)^2} \left[1-\frac{5}{4}q_C^2\xi(\delta,\tau)^2\right] \sin\left[\psi_D(\delta,\tau)\right] \nonumber \\ & & -\frac{1}{q_C^3\xi(\delta,\tau)^3}\left\{ \left[1-\frac{3}{2}q_C^2\xi(\delta,\tau)^2\right]\psi_D(\delta,\tau) -|q_C^2\xi(\delta,\tau)^2-1|^{3/2}l(\delta,\tau) \right\} \label{eq.Y2} \end{eqnarray}\] \[\begin{equation} l(\delta,\tau)= \begin{cases} \ln\frac{1+w(\delta,\tau)}{1-w(\delta,\tau)} & (q_C\xi(\delta,\tau)>1) \\ 2\tan^{-1}|w(\delta,\tau)| & (q_C\xi(\delta,\tau)\leq 1) \end{cases} \label{eq.l} \end{equation}\] \[\begin{equation} w(\delta,\tau)= \sqrt{\left|\frac{q_C\xi(\delta,\tau)-1}{q_C\xi(\delta,\tau)+1}\right|} \tan\frac{\psi_D(\delta,\tau)}{2} \label{eq.w} \end{equation}\] \[\begin{equation} \psi_D(\delta,\tau)= \cos^{-1}\frac{1}{\sqrt{1+q_D^2\xi(\delta,\tau)^2}} \label{eq.psi_D} \end{equation}\] \[\begin{equation} \xi(\delta,\tau)= \xi_0\left[\frac{\Delta\bar{\chi}(\delta,\tau)}{\Gamma_0}\right]^{\nu/\gamma} \label{eq.xi} \end{equation}\] \[\begin{equation} \Delta\bar{\chi}(\delta,\tau)= \max\left\{0, \delta\frac{P_c}{\rho_c}\left[ \frac{1}{\left(\PartialDiff{P}{\rho}\right)_T} -\frac{1}{\left.\left(\PartialDiff{P}{\rho}\right)_T\right|_{\tau=\tau_R}} \frac{\tau}{\tau_R} \right] \right\} \label{eq.Delta_chi} \end{equation}\] と書ける。ここで\(P\)は圧力、\(P_c\)は臨界点での圧力であり、 他の定数は表3のように与えられる。
The factor \(\bar{\mu}_2\) represents the behaviour of the viscosity near the critical point, given as Eqs. (\ref{eq.mu2})-(\ref{eq.Delta_chi}) according to Eqs. (33), (20)-(23), (26)-(28), and the note immediately after Eq. (34) of the paper. In these equations, \(P\) is a pressure and \(P_c\) is the pressure at the critical point. The other constants in these equations are given in Table 3.
表3. (\ref{eq.mu2})-(\ref{eq.Delta_chi})式の定数。 Huber et al. (2009)の表5に基づく。
Table 5. The constants in Eqs. (\ref{eq.mu2})-(\ref{eq.Delta_chi}) from Table 5 of Huber et al. (2009).
定数
Constant 		値
Value
\(\chi_{\mu}\) 0.068
\(q_C\) \(\frac{1}{1.9\times 10^{-9}}\) [m\(^{-1}\)]
\(q_D\) \(\frac{1}{1.1\times 10^{-9}}\) [m\(^{-1}\)]
\(\nu\) 0.630
\(\gamma\) 1.239
\(\xi_0\) \(0.13\times 10^{-9}\) [m]
\(\Gamma_0\) 0.06
\(\tau_R\) \(\frac{1}{1.5}\)

このヘッダファイル内で定義されている関数を以下に示す。 各関数の詳細は関数名をクリックしてリンク先を参照のこと。
Functions defined in this header file are listed below. For details of individual functions, click the links.
関数名
Function name 		機能・用途
Purpose
IAPWS95_calculate_viscosity_mu0 \(\bar{\mu}_0(\tau)\)を計算する。
Compute \(\bar{\mu}_0(\tau)\).
IAPWS95_calculate_viscosity_mu1 \(\bar{\mu}_1(\delta,\tau)\)を計算する。
Compute \(\bar{\mu}_1(\delta,\tau)\).
IAPWS95_calculate_viscosity_mu2 \(\bar{\mu}_2(\delta,\tau)\)を計算する。
Compute \(\bar{\mu}_2(\delta,\tau)\).
IAPWS95_calculate_viscosity Huber et al. (2009)の計算式を用いて 指定された密度・温度における粘性係数を計算する。
Compute the viscosity at given density and temperature using the equations of Huber et al. (2009).


◆検証 (Validation)

関数IAPWS95_calculate_viscosity を用いて粘性係数を計算し、 Huber et al. (2009)の表6, 表7と厳密に一致する結果が得られることを確認した。
The viscosities computed by function IAPWS95_calculate_viscosity were exactly identical to the numerical values shown in Tables 6 and 7 of Huber et al. (2009).

◆引用文献 (References)