Коэффициент сжимаемости газа

[sergeevcheb 0]

Попытался произвести расчет коэффициента сжимаемости газа использовав данные примера расчета, изложенные в приложении Г ГОСТ 30319.2-96 по модифицированному методу NX19.

Плотность при 0,101325 МПа и 293,15 К:0,6799 кг/м3.

Азот 0,8858 мол.%, диоксид углерода 0,0668 мол.%, давление 2,001 МПа, температура 270,00 К

Используя формулы (17) и (18) ГОСТ 30319.2-96 определяем псевдокритические значения давления и температуры

Рпк=2,9585*(1,608-0,05994*0,6799 (кг/м3)+0,0668 (мол.%)-0,392*0,8858 (молю%) =3,80706 МПа

Тпк=88,25*(0,9915+1,759*0,6799 (кг/м3)-0,0668 (мол.%)-1,681*0,8858 (молю%) =55,7387 К

Долее по формулам (15) и (16) определяем параметры

Ра=0,6714*(2,001 (МПа)/3,80706 (МПа))+0,0147 =0,36758 МПа

Та=0,71892*(270 (К)/55,7387 (К))+0,0007 =3,4832 К

Далее необходимо определить корректирующий множитель F в зависимости от Ра и ДТа где ДТа= Та-1,09 для нашего случая ДТа=3,4832-1,09 =2,393

Здесь возник для меня тупик, так как по ГОСТ ДТа может быть:

1) 0<=ДТа<=0,3

2)-0,25<=ДТа<0

А по расчетам ДТа=2,393 - то, есть не попадает под условия.

Подскажите пожалуйста, где я ошибся

[su215 657]

Ошибка (кстати, распространённая) в том, что Вы в формулы вместо молярных долей подставляете значения процентных концентраций, которые обычно даются в паспортах газораспределительных организаций. Пример: молярная доля компонента составляет 5, общее количество долей всех компонентов 100. Тогда молярная процентная концентрация этого компонента будет (5/100)*100%=5%. (Значение, указанное в паспорте). Молярных же долей будет 5/100=0,05. Это значение и указывается в формулах ГОСТа в виде х_а или х_у.

Для расчёта К_сж я использую программку, которую написал для себя в Exel. Попробуйте, может пригодится. На странице "исходные данные" вводятся плотность и компоненты. На странице "расчёт" - давление, температура и собственно расчёт. Красным цветом выделено то, что вводится с клавиатуры. Остальное, думаю, понятно.

__________.rar

[sergeevcheb 0]

Спасибо за подсказку.

Продолжил расчет. Расчет сравнивал с предложенной программкой.

Возник очередной вопрос.

При расчете Z по формуле (6) ГОСТ 30319.2-96 числитель в квадратных скобках не возводитсься в квадрат, а в формуле расчета по программке возводиться. У меня нет офицального издания вышеуказанного ГОСТа, но я просмотрел несколько электронных версий из разных источников, и там не где вышеуказанный числитель в квадрат не возводиться. Подскажите, что я не прав и я буду искать другой источник, может быть даже официальный.

[su215 657]

Проверил сейчас в двух разных изданиях ГОСТа, одно официальное (печатный вариант), другое скачанное с Интернета. В обоих числитель в квадрате. На всякий случай:

____30319_96.rar

[murik9999 0]

при расчетах не приводил к нормальным условиям только по температуре и давлению не учитывая коэфициент сжимаемости это большая ошибка?