Определение динамической вязкости воды и водяного пара

Пример применения функции calc_ws_dvisc

Автор: Сергей Медведев, 2020

import numpy as np
import matplotlib.pyplot as plt
#Импортируем фцнкцию calc_ws_dvisc
from calcwsdvisc import calc_ws_dvisc

Протестируем функцию calc_ws_dvisc на корректность выдаваемых результатов расчёта динамической вязкости для нескольких пар входных параметров, которые приведены для проверки в http://www.iapws.org/relguide/visc.pdf

#Температура, К
Ts = np.array([298.15, 298.15, 373.15, 433.15, 433.15, 873.15, 873.15, 873.15, 
               1173.15, 1173.15, 1173.15])
#Температура, C
ts = Ts - 273.15
#Плотность, кг/м3
denss = np.array([998, 1200, 1000, 1, 1000, 1, 100, 600, 1, 100, 400])
#Динамическая вязкость, Па*с, правильные ответы
dviscs0 = np.array([889.7351, 1437.649467, 307.883622, 14.538324, 217.685358, 32.619287, 35.802262,
             77.430195, 44.217245, 47.640433, 64.154608]) * 1e-6

#Находим с помощью функции calc_ws_dvisc значения динамической вязкости
#для каждой тестовой пары входных значений
dviscs = []
for i in range(Ts.shape[0]):
    dviscs.append(calc_ws_dvisc(ts[i], denss[i]))
dviscs = np.array(dviscs)
dviscs

array([8.89735100e-04, 1.43764947e-03, 3.07883622e-04, 1.45383245e-05,
       2.17685358e-04, 3.26192870e-05, 3.58022617e-05, 7.74301952e-05,
       4.42172445e-05, 4.76404331e-05, 6.41546078e-05])

#Рассчитываем относительные ошибки расчётных значений, %
err = (dviscs - dviscs0) / dviscs0 * 100
err

array([ 1.68374208e-08, -2.16756201e-08,  1.10926244e-07,  3.34136474e-06,
        1.21783182e-07, -7.97673192e-08, -7.76836891e-07,  2.93532069e-07,
       -1.09749859e-06,  1.70164286e-07, -2.36364193e-07])

#Максимальная по модулю относительная ошибка, %
np.abs(err).max()

3.3413647385231315e-06

Область применимости и погрешность расчётной методики

http://www.iapws.org/relguide/visc.pdf