Теплофизические свойства двуокиси углерода - Алтунин В.В.
Скачать (прямая ссылка):
7 31 1,18966+28 —5,4564 Ы8 9,4:] -4 Ь4 -51+20 —69+52
8 5 1,19006+5 -5,535: -4 12,4- = 1 —ЗІ- -2 265+И — 1292+43
9 5 1,19007+5 —5,539- :5 12,6- -2 -35: :4 299+31 -1474+159
10 5 1,19006+5 -5,534: :7 12,3- Ез -27= -7 209+76 —857+498
11 5 1,19006-1-5 -5,531- -8 12,1= :5 —21: 43 116±170 —60+1397
15
6 6 5
31 точка; максимальная плотность ^А =
6 43 1,18934+39 -5,409+22 7,9 + 4
7 13 1,189934-12 -5,504 + 9 11,1+2
8 5 1,19007 + 5 -5,540 + 5 12,6 + 2
9 5 1,19006 + 5 -5,533+6 12,3+3
10 5 1,19006+5 -5,533+8 12,2+4
28 точек максимальная плотность
6 14 1,18988 + 13 -5,489+9 10,4+2
7 5 1,19005 + 5 —5,530±4 12,1+2
8 5 1,19006+5 —5,535 + 6 12,4 + 3
9 5 Г, 19006+5 —5,535 + 8 12,4 + 5
492 (0)^2,08, /?<519 атм)
31+3 —157+13
—9+3 80+16
—35+3 296 + 20
—27 + 5 207+51
—26+11 204+137
»,=387 (со<1,64,
3 + 2 -24+2 —30+5 —30 + 12
: 181 атм)
— 18 + 12 191 + 16 252 + 51 247 + 154
1,18986+13 1,19006+5 1,19006+6 1,19006+5
12 точек; максимальная плотность ^Л = 92 ((0^0,38, р^7\ атм) —5,491+6 10,5±1
-5,523+5 -5,525+10 -5,5б6±22
И,5±1 11,6+5 15,3+19
-7+1 —9+9 -123+58
10+50 1458+730
250+23 —474+47 —1436+85 —897 + 298 -875+ 1036
—7,2+27 —875+60 -1224+282 -1183+1152
-6435+3239
Примечание. В таблице приводятся также вероятные допуски к коэффициентам, причем последняя цифра допуска соответствует последней цифре коэффициента.
сматриваемого интервала плотностей. Однако из-за большой Неопределенности в коэффициентах В{ при степенях р выше Четвертой, которая обусловлена неизбежной неточностью опытных данных, определение вириальных коэффициентов выше четвертого только по измерениям сжимаемости лишено смысла *. Таким образом, в нашем случае может идти речь об уравнении состояния, в котором вириальными будут лишь два-три младших коэффициента, а остальные являются простыми температурными функциями. По этой причине использованный в [1.88] критерий вириальности полинома становится излишним, так как гораздо важнее иметь подходящий критерий вириальности каждого коэффициента полиномиального ряда.
По определению, вириальный коэффициент В{ = 1(Т). Поэтому коэффициент Вг можно отождествить с вириальным только в том случае, когда значение искомого коэффициента в полиномах (1.10) будет сохраняться неизменным (разумеется, в пределах ошибки его определения) при вариации степени полинома г и интервала плотности Ар. Это позволяет сформулировать по крайней мере два критерия вириальности коэффициентов полиномиального ряда: 1) при фиксированной степени полинома г с изменением интервала плотностей на изотерме коэффициент полинома не должен зависеть от плотности; 2) при фиксированном интервале плотностей Д$ на изотерме с изменением степени полинома его коэффициенты не должны зависеть от плотности. Последующие расчеты показали, что предпочтительнее первый критерий вириальности В{, в соответствии с которым в [1.26] разработана машинная методика последовательного выделения коэффициентов ви-риального разложения.
Блок-схема программы представлена на рис. 3. Для заданной начальной степени полинома (в нашем случае г0=2) выбирается начальный интервал плотностей, который включает 3—4 опытные точки, и по методу наименьших квадратов определяются значения коэффициентов полиномам Из найденных коэффициентов запоминается коэффициент В1 при р' (аналог второго вириального коэффициента). Далее интервал плотностей расширяется на величину Ар (1—2 опытные точки) и процедура вычислений повторяется. В результате серми расчетов для фиксированной степени полинома отыскивается зависимость коэффициента Вх от плотности. Затем аналогичные расчеты выполняются для полиномов более высоких степеней и вновь запоминается только В^
__Из рис. 4 видно, что зависимость В1=/(р) распадается на
* Выполненное позднее [1.73] на том же объекте исследование ортогональных полиномов дало практически такие же результаты, как н в обсуждаемой здесь работе [1.88].
29
три характерных участка: начальный участок, на котором обнаруживается зависимость В{ от р (по-видимому, эта зависимость обусловлена влиянием погрешности исходных данных при малом интервале плотности); средний участок — участок стабильности, на котором коэффициент полинома практически не зависит от плотности; конечный участок — на котором обнаруживается резкая зависимость коэффициента полинома от плотности (эта зависимость обусловлена низкой степенью по-
Тнач и г
Начальный \ А интервал ) Др,, плотности )
•
ОпоШлщиг, коэффициентов полинома г = 1+?Н1 Р1
о. зонд °1 -к накопления с! V/
*
щиен ?к = Ртах V/ к.
Расширение интервала плотности Дрк — Дрк+*
Л\
*
Т_т
*
^3 3^ Выделение участка стабилизации
*
Усреднение В1 = {В1,к)
¦
Аппроксимация В[=В[(Т)
Фиксация В[=В[(Т)
¦
X 1? 1 ^ гпах
ОстаноВ
Рис. 3. Блок-схема программы метода последовательного выделения коэффициентов вириального разложения
линома для рассматриваемого интервала плотности). Таким образом, в качестве истинного вириального коэффициента должно быть принято среднее значение коэффициента полинома на участке стабильности. Такая же картина получается и для следующих вириальных коэффициентов.
