Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Расчет термодинамических свойств газов



Введение

В настоящее время существует довольно разнообразный набор программного обеспечения для расчёта физико-химических свойств газов. Однако, большинство из программных продуктов, либо являются узкоспециализированными приложениями, рассчитанными на определённые вещества или смеси (чаще всего это различные комбинации углеводородов), либо являются коммерческими разработками, стоимость которых весьма значительна. Кроме того, такие коммерческие продукты имеют тот недостаток, что их достаточно сложно встроить в состав разрабатываемого программного обеспечения для моделирования химико-технологических процессов, поскольку они имеют довольно большой объём и требования к оборудованию.

Цель работы:

Создание программного комплекса по расчету ряда термодинамических свойств газов при повышенных давлениях.

Задачи:

1) Расчет сжимаемости газов.

2) Расчет вязкости газов.

3) Сравнение точности используемых методов с литературными данными.

В связи с этим, нами на языке С++ был создан достаточно простой в работе программный комплекс для расчёта сжимаемости и вязкости газов.



1. Литературный обзор

1.1 Уравнение состояния идеального газа и уравнения реальных газов

Механическое состояние вещества в отличие от термического или термодинамического можно описать при наличии известных величин давления, температуры и объема. Поскольку эти параметры связаны между собой посредством, так называемого уравнения состояния f (P, V, Т) = 0, только два из них независимы. Основные уравнения термодинамики также представляют собой соотношения определенных трех параметров, на основе которых можно найти все прочие термодинамические свойства, включая соотношение Р-V-Т. И напротив, данные о соотношении Р-V-Т в сочетании с известной теплоемкостью идеального газа как функцией температуры также составляют одно из основных уравнений термодинамики. Само же уравнение состояния соответствующего вида может применяться для оценки многих важных свойств чистых веществ и их смесей, например, для оценки:

1) плотности жидкой и паровой фаз,

2) давления пара,

3) критических свойств смесей,

4) равновесных отношений между паром и жидкостью,

5) отклонения энтальпии от идеального состояния,

6) отклонения энтропии от идеального состояния.

В настоящее время не существует такого уравнения состояния, которое было бы в равной степени применимо для оценки всех этих свойств какой-либо обширной группы веществ. В то же время работа в этом направлении дала множество положительных результатов, имеющих ограниченную сферу применимости, о чем будет рассказано позднее. В данной же главе будет представлено описание ряда уравнений состояния. Некоторые из этих уравнений представляют только исторический интерес, однако большинство из них имеет как теоретическую, так и практическую ценность.

В 1834 г. французский физик Б. Клапейрон, работавший длительное время в Петербурге, вывел уравнение состояния идеального газа для постоянной массы газа. В 1874 г. Д.И. Менделеев вывел уравнение для произвольного числа молекул.

PV = nRT, (1.1)

Где n - число молей газа,

P - давление газа,

V - объем газа,

T - температура газа, K,

R - газовая постоянная.

Для определения свойств газовых смесей используют уравнения состояния, которые устанавливают связь между температурой, объемом и давлением системы. Термодинамические свойства природных и нефтяных газов и их компонентов значительно отличаются от свойств идеальных газов, особенно при низких температурах и высоких давлениях поэтому уравнение состояния идеальных газов не может быть использовано для определения этих свойств. Для описания поведения реальных газов разработан ряд уравнений состояния. Наибольшее применение для углеводородных систем получили уравнения Бенедикта-Вебба-Рубина и Редлиха-Квонга и их модификации.?



1.2 Уравнение Бенедикта-Вебба-Рубина

Уравнение Бенедикта-Вебба-Рубина имеет вид [1]:

(1.2)

где

P - давление газа, Па,

T - абсолютная температура, K,

R - газовая постоянная, 0,0821 л·атм/моль·K,

с - мольная плотность, моль/мі.

B₀, A₀, C₀, b, a, c, б, г - коэффициенты уравнения, не зависящие от давления, температуры и плотности и постоянные для определенного вещества. Числовые значения коэффициентов являются справочными величинами.

Это уравнение часто используется при расчете таких физико-химических свойств углеводородных смесей, как плотность, вязкость, теплоемкость и др.

Уравнение состояния Бенедикта-Вебба-Рубина является наиболее точным из известных способов определения термодинамических свойств легких углеводородов. Это уравнение слишком сложно, чтобы применять его для непосредственных расчетов равновесия пар - жидкость; для каждой данной серии условий необходимо решать методом последовательных приближений систему нелинейных уравнений. На этих графиках для ряда постоянных давлений приводятся два коэффициента в виде функции температуры и среднемольной температуры кипения паровой и жидкой фаз.

В первоначальном виде уравнение Бенедикта-Вебба-Рубина было выведено для расчета термодинамических свойств и фазового равновесия 12 углеводородов - от метана до н-гептана. В последующем область применения уравнения была расширена до 33 углеводородов и 9 неуглеводородных веществ [2].

Уравнение состояния Бенедикта-Вебба-Рубина применимо и для расчета термодинамических свойств легких углеводородов. В этом случае коэффициенты уравнения БВР для смеси находят по известным коэффициентам чистых компонентов, используя правила смешения [2]:

(1.3)

Популярность этого уравнения привела к появлению его различных модификаций.

Температурные модификации

В ходе обработки экспериментальных данных, некоторые авторы [3, 4] отметили, что при температурах ниже нормальной точки кипения коэффициент С уравнения Бенедикта-Вебба-Рубина лучше заменить на функцию температуры, для того чтобы уравнение более точно описывало давление паров.

Модификация Кауфмана

Кауфман предложил аппроксимацию вида:

(1.4)



где A_i - константы, зависящие от свойств вещества.

Также Кауфман опубликовал константы уравнения состояния Бенедикта - Вебба-Рубина для - 1-бутена, цис-2-бутена, 1-пентена, метилацетилена и 1,3 - бутадиена.

Модификация Орая

Наиболее тщательный количественный анализ проблемы зависимости С(Т) был выполнен Ораем [5]. Он предложил следующую температурную зависимость для С:

C^1 / 2(T) = C^1 / 2(T₀) ? ДC^1 / 2(T), (1.5)

где C(T₀) - значение константы С, а величина ДC^1 / 2(T) - является полиномом 5-ой степени.

ДC^1 / 2(T) = Q₁И² + Q₂И³ + Q₃И⁴ + Q₄И⁵, (1.6)

где И = (T₀ ? T) / T₀ - безразмерный комплекс температуры, Т₀ - реперная температура.

Константу А₀ из уравнения Бенедикта-Вебба-Рубина Орай выражает как:

(1.7)

По простому правилу, впервые предложенному Бенедиктом

(1.8)



Уравнение (1.2.6) будет идентично (1.2.7), если M_ij = 2 для всех ij пар (i ? j). Орай ослабил это ограничение, сделав M_ij настраиваемым параметром, зависящим от температуры, но не зависящим от плотности и состава.

Модификация Старлинга

Старлинг предложил [6, 7] модифицировать уравнение Бенедикта-Вебба-Рубина таким образом, чтобы от температуры зависел не только коэффициент C, но и коэффициент a, получив тем самым уравнение состояния Бенедикта-Вебба-Рубина-Старлинга с одиннадцатью параметрами

(1.9)

Уравнения Бенедикта-Вебба-Рубина и Бенедикта-Рубина-Старлинга нашли широкое применение в практической деятельности. Они более эффективны, чем кубические, например, для любых термодинамических расчетов, особенно для криогенных систем [8]. Соперничать с ними по степени точности расчетов могут разве что уравнения Ли-Кеслера и подобные им. В то же время усложненность формы затрудняет их использование при расчетах дистилляционных процессов в тех случаях, когда необходимо выполнить многократную оценку величин фугитивности, объема и энтальпии, поэтому в перечисленных случаях часто применяют кубические уравнения. Однако ряд ученых не придает большого значения этому ограничению.

Обобщенные модификации

Успешное использование оригинала уравнения Бенедикта-Вебба-Рубина при расчетах волюметрических и термодинамических свойств чистых газов и жидкостей обусловило появление ряда работ, в которых это уравнение или его модификация приводится к обобщенной форме, применимой ко многим типам соединений [9, 10].

Модификация Ли-Кесслера

Ли и Кеслер разработали [11] модифицированное уравнение состояния Бенедикта-Вебба-Рубина, используя трёхпараметрическую корреляцию Питцера [12]. По их методу коэффициент сжимаемости реального вещества связывается со свойствами простого вещества, для которого щ = 0, и н-октана, выбранного в качестве эталона. Для того чтобы рассчитать коэффициент сжимаемости вещества при некоторых значениях температуры и давления, используя критические свойства этого вещества, сначала следует определить приведённые параметры T_r и P_r. Затем рассчитывается идеальный приведённый объём простого вещества по уравнению:

(1.10)

где V_r⁽⁰⁾ = P_kV⁽⁰⁾ / RT_k - идеальный приведенный объем простого вещества;

V⁽⁰⁾ - мольный объем простого вещества, м³/моль;

P_r = P / P_k - приведённое давление;

P_k - критическое давление, Па;

R= 82,057… - м³*атм/(моль*К).универсальная газовая постоянная



 - коэффициенты [13].

После определения V⁽⁰⁾, рассчитывается коэффициент сжимаемости простого вещества:

(1.11)

Далее, используя те же приведенные параметры, определенные ранее, снова решается первое уравнение, относительно V_r⁽⁰⁾, но уже с константами для эталонного вещества. После этого находят коэффициент сжимаемости эталонного вещества:

(1.12)

где Z^(R) - коэффициент сжимаемости эталонного вещества;

V_r^(R) - приведенный объем эталонного вещества.

Коэффициент сжимаемости Z интересующего нас вещества определяется из уравнения:

(1.13)



где щ - фактор ацентричности Питцера исследуемого вещества,

щ^(R) = 0,3978 - фактор ацентричности Питцера эталонного вещества (октана).

Уравнение применяется, в основном, для углеводородов в интервалах значений 0,3 ? T_r ? 4,0 и 0 ? P_r ? 10 для паровой и жидкой фазы, где средняя погрешность составляет менее 2%. [13]

Модификация Нишиуми

По утверждению Хопке [14], как уравнение Бенедикта-Вебба-Рубина, так и уравнение Бенедикта-Вебба-Рубина-Старлинга не дают возможность получить достаточно точные параметры для большинства полярных жидкостей и воды, в частности.

Чтобы устранить этот недостаток Нишиуми разработал [15, 16] обобщающую модификацию уравнения Бенедикта-Вебба-Рубина и привёл данные для 92-х веществ, в том числе воды.

Уравнение Нишиуми для коэффициента сжимаемости Z имеет вид:

(1.14)

где с_r = с / с_k - приведенная плотность, с_k - кинетическая плотность. Все пятнадцать коэффициентов, отмеченные «звездочками», являются функциями коэффициента ацентричности щ, величины и выражают воздействие полярности на свойства паров и жидкостей соответственно.



1.3 Уравнение Редлиха-Квонга

Уравнение Редлиха-Квонга [17, 18] имеет вид:

(1.15)

где P - давление, Па;

T - абсолютная температура, К;

V - мольный объём, мі/моль;

 - универсальная газовая постоянная, Дж/(моль·К);

a и b - некоторые константы, зависящие от конкретного вещества.

Уравнение состояния Редлиха-Квонга значительно проще большинства известных уравнений. Точность данного уравнения в целом соответствует точности определения состава и свойств остатка пластовых смесей [19].

Первоначально уравнение состояния Редлиха-Квонга использовали для определения термодинамических и физических свойств газовых смесей. Применение этого уравнения для расчета свойств жидких углеводородных смесей приводит к большим погрешностям [20].

В практике расчета равновесий широкую известность получили уравнения состояния Редлиха-Квонга, Соаве, Пенга-Робинсона, Бенедикта-Вебба-Рубина, позволяющие рассчитывать равновесие жидкость-пар в однокомпонентных, бинарных и многокомпонентных системах в весьма широком интервале внешних условий. Смеси неполярных веществ обычно с удовлетворительной точностью описывают только на основе данных о чистых компонентах. Параметры уравнений для смесей при этом находят с помощью определенных комбинационных правил на основе констант, характеризующих индивидуальные вещества. В более сложных случаях необходима оценка некоторых бинарных параметров по экспериментальным данным для смесей [21].

Летучесть компонента смеси в газовой фазе рассчитывают по уравнению состояния Редлиха-Квонга, а коэффициент активности по выражению из теории регулярных растворов Гильденбранта. В это выражение входят молярный объем компонента в жидком состоянии и параметр растворимости, который вычисляют по энтальпии испарения компонента смеси [22].

Метод Чао и Сидера рекомендуется использовать в расчетах процессов стабилизации и дебутанизации насыщенного конденсата, а также сепарации природных газоконденсатных смесей при давлениях примерно до 10 МПа и температуре не ниже 263 К - чтобы расширить пределы применимости этого метода, предпринимались попытки его модификации, основанные на использовании различных модификаций уравнения состояния Редлиха-Квонга и теории жидких растворов.

Робинсон и Чао коэффициент летучести компонента в газовой фазе вычисляют по уравнению состояния Редлиха-Квонга с использованием обобщенных коэффициентов Чу и Прауснитца. Для определения летучести чистого вещества подобраны эмпирические выражения, отличающиеся от выражений Чао, Сидера и Шелтона. Коэффициент активности рассчитывают с помощью уравнения Гильдебранта и Скотта, учитывающего отклонение поведения жидких растворов от поведения регулярных растворов.

Особенно эффективной оказалось модификация, предложенная Соаве, которая часто дает хорошие результаты при расчетах равновесия пар-жидкость в углеводородных смесях. Это объясняется, по крайней мере частично, использованием для определения констант уравнения данных по давлениям паров чистых компонентов. В результате уравнение Соаве, надежное почти всегда при расчете констант равновесия К, обычно дает неправильные значения плотности жидкости [23].

В заключение следует отметить, что, несмотря на простоту, осуществленная Соаве модификация уравнения Редлиха-Квонга явилась очень эффективной и позволила улучшить моделирование свойств не только паровой (газовой) фазы, но и газированной жидкой фазы. И все же, свойства паровой фазы моделируются модификациями уравнения состояния Редлиха-Квонга значительно более точно, чем свойства жидкой фазы. Для точного описания фазового равновесия необходимо правильно определять свойства обеих сосуществующих фаз.

1.4 Уравнение Барнера-Адлера

Барнер и Адлер модифицировали уравнение состояния Иоффе [24, 25], чтобы получить P - V - T соотношение, которое является наиболее точным для области насыщенного или слегка перегретого пара. Диапазон применимости уравнения соответствует значениям V_r > 0,6 и T_r < 1,5. [13]

(1.16)

 

 

 

 

 

 

 

 

Чтобы использовать его для смесей, необходимо связать Рсm, Тсm и щm с составом. Для этого рекомендуются следующие уравнения:

(1.17)

Величины Тс_m и Рс_m определяются прямым расчетом, если известны значения k^*_ij для всех возможных пар. Эти параметры k^*_ij характеризуют бинарное взаимодействие компонентов i-j, а не сами компоненты i или j по отдельности.

Для учета взаимодействия пар простых компонентов k^*_ij может быть приближенно скоррелировано с отношением критических объемов чистых i и j. Величина k^*_ij предполагается независимой от давления, температуры и состава.

1.5 Уравнение Суги-Лю

В 1971 г. Суги и Лю [26] предложили уравнение состояния, подобное уравнению Барнера-Адлера и применимое в пределах того же диапазона, т.е. для области насыщенного или слегка перегретого пара. Точность этого соотношения превосходит точность метода Барнера-Адлера. Кроме того, оно дает удивительно хорошие результаты при расчетах объемов насыщенной жидкости.

(1.18)

где

 - давление, Па;

 - абсолютная температура, К;

 - молярный объём, мі/моль;

 - универсальная газовая постоянная, Дж/(моль·К);



 - критическая температура, К;

 - критическое давление, Па;

 - фактор ацентричности Питцера;

 - критический коэффициент сжимаемости;

 - критический объём, мі/моль.

Для этого уравнения состояния предложены следующие правила для определения Тс_m, сс_m и Zс_m:

(1.19)

Параметром бинарного взаимодействия в этом случае является k'_ij, значение которого обычно близко к нулю. Связь его с параметром k^*_ij в уравнении Барнера-Адлера далеко не проста. Значений k'_ij практически не имеется, и эту величину следует либо принимать равной нулю, либо определять по экспериментальным данным. Если k'_ij предполагается равным нулю, то правила для псевдокритических температуры и давления сводятся к соответствующим правилам для уравнения Редлиха-Квонга в оригинальном виде [13].

1.6 Уравнение Ли-Эрбара-Эдмистера

Ли, Эрбар и Эдмистер в 1972-73 годах вывели трехпараметрическое уравнение состояния углеводородных систем [27, 28]. Уравнение разработано таким образом, чтобы определять энтальпии и константы фазового равноавесия углеводородных смесей, а не волюметрических свойств. Данное уравнение является доработанным уравнением состояния, которое было выведено в 1971 году Ли и Эдмистером [29, 30].

Уравнение Ли-Эрбара-Эдмистера имеет вид:

(1.20)

где a = [(0,246105 + 0,02869щ) - (0,037472 + 0,149687щ) T_r + (0,16406 + 0,023727щ) T_r-¹ + (0,04937 + 0,132433щ) T_r-²] R²T_k²/P_k;

b = (0,086313 + 0,002щ) RT_k/P_k;

c = [(0,451169 + 0,00948щ) T_r-^1/2 + (0,387082 + 0,078842щ) T_r-²] R²T_k²/P_k;

 - давление, Па;

 - абсолютная температура, К;

 - молярный объём, мі/моль;

 - универсальная газовая постоянная, Дж/(моль·К);

 - критическая температура, К;

 - критическое давление, Па;

 - фактор ацентричности Питцера;

T_r = T/T_k



Метод Ли-Эрбара-Эдмистера рекомендуется для области температур от 158 до 262°С и давлений, не превышающих 0,8 критического давления смесей. Метод применим для смесей углеводородов (парафиновых, ароматических, нафтеновых) с примесью азота, диоксида углерода и сероводорода. Этот метод не рекомендуется для систем, содержащих 50% или более неуглеводородных газов в жидкой фазе и компоненты, критическая температура которых значительно ниже температуры системы. Среднее отклонение расчетных значений констант фазового равновесия от экспериментальных составляет 6%. [2]

1.7 Фактор ацентричности

Американский физико-химик Питцер установил, что различия между веществами можно установить из анализа температурной зависимости давления насыщенного пара, если ее представить в безразмерной форме, т.е. откладывая по оси абсцисс приведенную температуру, а по оси ординат - приведенное давление.

Оказалось, что для фиксированной приведенной температуры величина приведенного давления связана с геометрической формой молекулы. Для инертных газов, у которых форма молекулы близка к сферической:

- lg P_пр = 1 при Т_пр = 0.7

Поэтому для таких веществ параметр, определяемый соотношением

(1.21)

близок к нулю. Для других веществ расчеты по этому уравнению дают величину щ, отличную от нуля, причем, чем больше форма молекулы отклоняется от сферической, тем больше величина щ. В связи с эти параметр щ, характеризующий не сферичность, асимметрию молекулы, получил наименование фактора ацентричности.

Чтобы получить значение щ, необходимо знать давление паров при T_r = T / T_c = 0,7, а также критическое давление. Сначала эта константа была предложена для описания ацентричности или несферичности молекулы. Для одноатомных газов фактор ацентричности, следовательно, близок к нулю. Для метана он все еще мал. Однако для углеводородов с высокой молекулярной массой значение щ возрастает и резко увеличивается с полярностью. В настоящее время фактор ацентричности широко используется в качестве параметра, который в известной степени характеризует сложность молекулы как в отношении ее геометрии, так и полярности. Применимость корреляций, включающих фактор ацентричности, должна ограничиваться нормальными газами и жидкостями. Не следует использовать такие корреляции для H₂, He, Ne, а также для сильно полярных веществ и соединений с водородными связями. [13]

Методы расчета фактора ацентричности

Метод Эдмистера:

По методу Эдмистера расчет фактора ацентричности [31] производится по формуле:

(1.22)

Метод Ли-Кеслера:

По методу Ли-Кеслера расчет фактора ацентричности производится по формуле:



(1.23)

По критическому коэффициенту сжимаемости:

Между фактором ацентричности щ и критическим коэффициентом Риделя a_c существует следующая зависимость:

(1.24)

По критическому коэффициенту Риделя:

Между фактором ацентричности и критическим коэффициентом Риделя существует следующее соотношение:

(1.25)

1.8 Безразмерные температуры и давления, критические давления и температуры и методы их расчетов

Газы могут быть превращены в жидкость сжатием при условии, что температура не превышает определенного для каждого газа значения. 
Температура, при которой и выше которой данный газ не может быть сжижен никаким повышением давления, называется критической Т_кр. 
Давление, при котором и выше которого повышением температуры нельзя испарить жидкость, называется критическим Р_кр. 
Объем газа, соответствующий критической температуре, называется критическим V_кр, а состояние газа, отвечающее критической температуре, критическому давлению и критическому объему, - критическим состоянием газа. [32, 33]

При изменении величины давления или температуры все газы могут быть превращены в жидкость или из жидкости в пар. Для каждого газа существует определенная температура, выше которой, несмотря на применение любого высокого давления, он не может быть переведен в жидкое состояние. Эта температура называется критической, а давление, необходимое для сжижения газа при этой температуре, критическим.

Иными словами, критическая температура чистых веществ - это та максимальная температура, при которой жидкая и паровая фазы еще могут существовать в равновесии. Давление паров при этой температуре называется критическим давлением, а объем единицы плотности вещества критическим объемом.

Понятие о критической температуре впервые было предложено в 1861 г. Д.И. Менделеевым, который понимал под ней температуру, при которой исчезают силы сцепления между молекулами жидкости, и она, независимо от давления, превращается в газ.

Приведенные выше определения критических параметров справедливы для индивидуальных газов. Технические газы в большинстве случаев представляют собой сложные смеси, критическая температура которых всегда выше критической температуры самого низкокипящего компонента и ниже критической температуры высококипящих компонентов.

Температуру, при превышении которой данный газ не может быть сжижен никаким повышением давления, называют критической температурой газа.

Объем газа, соответствующий критической температуре, называют критическим объемом (F p), а состояние газа, определяемое критическими температурой, давлением и объемом, - критическим состоянием газа. Плотность пара над жидкостью при критическом состоянии становится равной плотности жидкости.

Приведенные критические параметры и закон соответственных состояний.

Отношения параметров газа р, Т и V, характеризующих его состояние, к их критическим значениям, называются приведенными.

Под приведенной температурой понимается отношение абсолютной температуры газа к его критической температуре T_c:

T_r = T / T_c. (1.26)

Приведенным давлением называется отношение абсолютного давления газа к его абсолютному критическому давлению P_c:

P_r = P / P_c. (1.27)

Если параметры газа выразить в безразмерных приведенных величинах, можно установить, что для газов существует так называемое общее уравнение состояния, лишенное величин, характеризующих данный газ:

f (P_r, Т_r, V_r) = 0 (1.28)

Состояние, отвечающее этому условию, называется законом соответственных состояний, по которому два газа будут иметь один и тот же приведенный объем, если у них равны приведенные давление и температура. Значение закона соответственных состояний заключается в том, что он позволяет с достаточной точностью определять коэффициенты сжимаемости газов и другие параметры по графикам для других газов, построенным в приведенных параметрах.



1.9 Вязкость

Вязкость, или внутреннее трение, - это свойство всех текучих тел (жидкостей или газов) оказывать сопротивление перемещению одной их части относительно другой. Молекулярно-кинетическая теория объясняет вязкость тепловым движением и взаимодействием молекул.

Два соприкасающихся элемента жидкости, двигающиеся в одном и том же направлении, но с разными скоростями, воздействуют друг на друга. Сила взаимодействия ускоряет медленно движущийся слой жидкости и замедляет более быстрый. Ньютон предположил, что величина этой силы, называемой силой внутреннего трения, пропорциональна изменению скорости течения жидкости v в направлении, перпендикулярном ее движению, и зависит от площади S соприкосновения слоев жидкости [34]

(1.29)

Коэффициент з пропорциональности называется коэффициентом динамической вязкости.

Величина dv/dy называется градиентом скорости течения жидкости, или быстротой изменения скорости, на единицу длины в направлении, перпендикулярном скорости упорядоченного движения ее слоев.

Теория вязкости газов основывается на кинетической теории газов, которая, в отличие от теории вязкости жидкости, изучена и разработана достаточно подробно. Однако, при повышенных давлениях методы, основанные на теории вязкости, не дают удовлетворительной точности значений коэффициента динамической вязкости. Поэтому в инженерных расчетах эти методы практически не используются. Широкое распространение в математических расчетах получили графические и полуэмпирические методы, предполагающие незначительное влияние давления на значения коэффициентов вязкости. В большинстве случаев при выполнении гидравлических расчетов влиянием давления на величину значения динамической вязкости пренебрегают. Для нахождения вязкости углеводородных газов, а также природных газовых смесей применяют принцип соответственных состояний. Теоретической основой данного принципа является то, что зависимости свойств различных веществ от безразмерных величин, называемых «приведенными параметрами», одинаковы. Отсюда следует, что когда приведенные параметры сравниваемых веществ равны, то и их свойства тоже равны. При определении приведенных параметров соответствующие размерные комплексы данного вещества (такие как, давление, температура и т.д.) относят к их значению в каком-либо характерном соответственном состоянии этого вещества (состоянии сравнения).

Расчет вязкости при умеренных давлениях

Для определения зависимости динамического коэффициента вязкости газа от температуры в области умеренных давлений чаще всего применяется формула Сатерленда:

(1.30)

где м_T⁰ - вязкость газа при температуре T K;

м₀⁰ - вязкость газа при температуре 0 С;

С - постоянная Сатерленда, характерная для данного газа.

Уравнение Сатерленда обычно дает точные результаты, как это видно из сравнения расчетных и экспериментальных значений вязкости этилена:



Значения вязкости этилена по уравнению Сатерленда



2. Теоретическая часть

2.1 Расчет фактора сжимаемости

Фактор сжимаемости учитывает отклонения реального газа от уравнения состояния идеального газа (уравнение Менделеева - Клапейрона) и представляет собой функцию температуры и давления, или функцию приведенных температуры и давления, а для смесей - функцию псевдоприведенной температуры и давления. Для идеальных газов, в которых отсутствует взаимодействие между молекулами, справедлив закон Клайперона, который обобщает законы Бойля-Мариотта и Гей-Люссака:

(P · V) / T = (P0 · V0) / T0 = R, (2.1)

где R ? универсальная газовая постоянная.

Для реальных газов произведение объема на давление (P · V) при постоянной температуре не является постоянной величиной. С ростом давления, это произведение вначале уменьшается до определенного предела, а затем с увеличением давления увеличивается. При практических расчетах для оценки степени указанных отклонений в уравнение Клайперона вводится поправочный коэффициент Z. Уравнение Клайперона - Менделеева с учетом коэффициента сжимаемости имеет вид

PV = ZNRT, (2.2)

где Р - абсолютное давление,

V - объем, занимаемый газом (газоконденсатной смесью) при определенных давлении и температуре,

R - универсальная газовая постоянная,

N - число молей газа.

Для расчета коэффициента сжимаемости, интересующего нас вещества, используем уравнение (по модификации Ли-Кеслера):

(2.3)

Где Z⁽⁰⁾ - коэффициент простого вещества и рассчитывается по формуле:

(2.4)

А Z^(R) - коэффициент эталонного вещества:

(2.5)

2.2 Расчет вязкостей компонентов при низких давлениях

Метод Голубева

Метод Голубева, основанный на преобразовании Трауца:



Размерность искомой величины м - мкПа•с.

В этом выражении м_kr^* определяется по выражению:

где M - молекулярная масса; P_kr - критическое давление, МПа; T_kr - критическая температура, К.

Расчет SO₂ (для T = 317.16 К (40 ^oC))

Рассчитанная величина - 138,241 мкП

Литературные данные - 135 мкП

Погрешность - 2,40%

Метод Тодоса

Соотношение Тодоса для неполярных газов:

Для полярных газов с водородной связью (T_r < 2,0):

Для полярных газов, не образующих водородную связь (T_r < 2,5):



Расчет SO₂ (для T = 317.16 К (40 ^oC))

Рассчитанная величина - 135,196 мкП

Литературные данные - 135 мкП

Погрешность - 0,14%

Метод Райхенберга

Райхенберг предложил другое соотношение, основанное на использовании принципа соответственных состояний:

Также, для органических соединений Райхенберг рекомендует использовать

Расчет ацетилена (для T = 374.16 К (101 ^oC))

Рассчитанная величина - 126,927 мкП

Литературные данные - 126 мкП

Погрешность - 0,74%

Метод Чэпмена

Существующие методики расчета вязкости газа, при давлениях близких к атмосферным, в большинстве случаев, базируются или на теории Чэпмена-Энскога, или на принципе соответственных состояний.

Уравнение Чэпмена-Энскога для вязкости имеет вид:



где M - молекулярная масса; T - температура, K; у - параметр потенциала Леннарда-Джонса, Е; Щ_х - интеграл столкновений.

Чтобы использовать его для расчета вязкостей, необходимо найти параметр потенциала Леннарда-Джонса и интеграл столкновений.

Параметр потенциала Леннарда - Джонса определяется из уравнения (для неполярных газов):

где щ - фактор ацентричности; P_kr - критическое давление, МПа; T_kr - критическая температура, K;

Интеграл столкновений для неполярных газов определяется по выражению Нойфельда:

Для нахождения Щ_х сначала вычисляют величину:



где соотношение е / k определяется по выражению Ти, Готоха и Стьюарта:

где щ - фактор ацентричности; е - параметр потенциала Леннарда-Джонса, эрг; k = 1,3805•10-16, эрг/K - постоянная Больцмана.

2.3 Описание программы

Комплекс создан по модульному принципу и позволяет в будущем расширить его функционал за счёт присоединения дополнительных модулей для расчёта других свойств (теплоёмкость, теплопроводность и пр.), а также модифицировать его для расчёта свойств газовых смесей.

Следует отметить, что модули могут работать, как индивидуально, так и группами, в зависимости от поставленной задачи. Например, для расчёта вязкости компонентов при атмосферном давлении не требуется расчёт факторов сжимаемости, и данный модуль подключается в случае необходимости.

Для расчёта факторов сжимаемости в данном программном комплексе используется уравнение состояния Бенедикта-Вебба-Рубина в модификации Ли-Кесслера. Данное уравнение состояния демонстрирует хорошие результаты, как для неполярных, так и для полярных веществ. Также его можно применять для смесей газов. Выбор именно данного уравнения состояния обусловлен, в частности, и спецификой дальнейшего применения рассчитываемого с его помощью фактора сжимаемости и его составляющих Z⁽⁰⁾ и Z⁽¹⁾ для расчёта прочих термодинамических параметров реальных газов. Также учитывались возможности по дальнейшему развитию данного метода.



Рис. 2.1. Принципиальная схема организации программного комплекса

Фактор сжимаемости Z интересующего нас i-го компонента рассчитывается по уравнению (2.3.1):

Z_i = Z⁽⁰⁾ + (/^R) (Z^(R) - Z⁽⁰⁾) = Z⁽⁰⁾ + Z⁽¹⁾, (2.3.1)

где фактор ацентричности эталонного вещества (н-октана) ^R = 0,3978.

Фактор ацентричности простого вещества равен нулю. Помимо исходными данными для расчёта являются температура, давление и критические параметры (P_c, T_c, V_c). Необходимые параметры уравнения состояния для простого и эталонного веществ брались из [13, табл. 3.9].

В ходе выполнения расчётов требовалось найти решения нелинейного уравнения, имеющего пять корней. Из всех корней выбирается либо минимальный, либо максимальный корень в зависимости от агрегатного состояния вещества. Для паровой фазы применяется наибольший корень.

Для расчёта динамической вязкости при повышенных давлениях использовались метод Райхенберга (МР) и метод Джосси-Стила-Тодоса (ДСТ). Первый из них требует наличие критических параметров T_c, P_c и дипольного момента м _p для каждого вещества_. Второй же в свою очередь - V_c и Z. Что позволяет применять данный модуль отдельно, либо в связке с модулем расчёта сжимаемости. Кроме того, оба метода требуют знания величины вязкости при нормальном давлении (м ⁰), которая рассчитывалась по одному из четырёх методов: Голубева, Тодоса, Райхенберга, Чепмена. Выбор конкретного метода расчёта производится пользователем с учётом особенностей нужного компонента.

Точность расчёта факторов сжимаемости проверялась на основании данных по плотности или удельному объёму газов. Эти величины легко получить из обобщённого уравнения Менделеева-Клапейрона для 1 моль газа (PV_m=ZRT).

Программа состоит из двух расчетных модулей и интерфейса взаимодействия с пользователем.

Первый модуль - модуль расчета фактора сжимаемости. Главная расчетная функция называется Z_Calc.

Входные данные:

Nt, Np-Точек по температуре и давлению

Tr0, Pr0-Массивы температур, давлений (безразм.)

Выходные данные:

Z0, Z1-Коэффициенты фактора сжимаемости Z=Z0+Omega*Z1

Fx(2) - Отклонение от 0 при решении нелинейного уравнения (критерий правильности решения)

Err-Ошибка в расчетах

Второй модуль - модуль расчета вязкости компонентов. В нем содержатся 2 основных функции расчета:

1. Subst_Viscosity0

Ее задача - вычисление вязкостей компонентов при низких давлениях

Входные данные:

Nm-Номер метода 1-4

Варианты:

1-Метод Голубева ур 9.4.14-16

2-Метод Тодоса ур. 9.4.17-20

3-Метод Райхенберга ур 9.4.21-22

4-Метод Чепмена

M - молярная масса г/моль

Tc-критическая температура, К

Pc-критическое давление, атм

Zc-Критический фактор сжимаемости (Тодос)

Т - Температура, К

FlH-Образует ли полярный компонент водородную связь (Только для метода Тодоса)

Dm - Дипольный момент, Дб (Тодос)

SumNC-Сумма групповых составляющих Ci для уравнения Райхенберга

EpsDm-критерий полярности молекулы Если Delta<EpsDm, то молекула неполярна Delta определяется по формуле 9.4.6

Om-фактор ацентричности (Чепман)

Sig, e_k - параметры потенциала Леннарда Джонса (Чепман)

FlLD - Вычислять ли параметры потенциала Леннарда Джонса или принять табличные

Выходные данные:

Mu0 - Искомая вязкость, мкП

Polar - Считается ли молекула полярной

Err - Ошибка

2. Subst_Viscosity

Ее задача - Расчёт вязкости чистых компонентов по формулам Джосси, Стила и Тодоса или Райхенберга

Входные данные:

Nm-Номер метода 1-Джосси, Стила и Тодоса, 2-Райхенберга

Tc, Pc, Vc-критические параметры

Z-Фактор сжимаемости

Mu0-Вязкость рпи нормальном давлении и той же температуре

T, P-Температура, К, Давление атм

FlC-флаг продолжения расчётов даже если метод неприменим

M-Молекулярная масса, кг/кмоль

Polar-Полярный ли компонент

Mp-дипольный момент, Дб

Cf-коэффициент (5%) допуска при применимости формулы 1 метода

Выходные данные

Mu-Искомая вязкость

Err-Сообщения об ошибке



3. Расчетно-практическая часть

В рамках тестирования программного комплекса были проведены расчеты для нескольких веществ:

1) Вода, H₂O

Таблица 3.1. Сравнение расчётных и экспериментальных значений для мольных объёмов воды

Эксперимент, м3/кг [35]	Расчётные значения, м3/кг	|D|, %
P, бар	50,00	100,00	210,00	300,00	600,00	50,00	100,00	210,00	300,00	600,00	50,00	100,00	210,00	300,00	600,00
P, атм	50,66	101,33	212,78	303,98	607,95	50,66	101,33	212,78	303,98	607,95	50,66	101,33	212,78	303,98	607,95
T, К	400	0,057810	0,026450	0,009100	0,002830	0,001638	0,055972	0,025452	0,008610	0,003127	0,002039	3,2	3,8	5,4	10,5	24,5
	600	0,078700	0,038370	0,017190	0,011440	0,004828	0,076527	0,037280	0,016745	0,011179	0,004986	2,8	2,8	2,6	2,3	3,3
	800	0,098030	0,048560	0,022660	0,015620	0,007464	0,095707	0,047490	0,022280	0,015442	0,007603	2,4	2,2	1,7	1,1	1,9
	1000	0,117000	0,058290	0,027560	0,019190	0,009522	0,114381	0,057171	0,027244	0,019110	0,009707	2,2	1,9	1,1	0,4	1,9



Таблица 3.2. Сравнение расчётных и экспериментальных значений для вязкости воды при низких давлениях

	Эксперимент, мкП	Расчётные значения, мкП	|D|, %
T, К	[35]	1	2	3	4	1	2	3	4
400	243±7	190	241	-	241	21,8	3,7	-	3,6
500	284±8	217	279	-	280	23,8	3,4	-	3,2
600	325±10	241	318	-	318	25,8	3,5	-	3,4
700	365±11	264	356	-	356	27,6	3,5	-	3,7

1. Метод Голубева

2. Метод Тодоса

3. Метод Райхенберга (неприменим для воды)

4. Метод Чепмена

Таблица 3.3. Сравнение расчётных и экспериментальных значений для вязкости воды при различных давлениях

Эксперимент [35], мкП	Расчётные значения*, мкП	|D|, %
P, бар	50,00	100,00	210,00	300,00	600,00	50,00	100,00	210,00	300,00	600,00	50,00	100,00	210,00	300,00	600,00
P, атм	50,66	101,33	212,78	303,98	607,95	50,66	101,33	212,78	303,98	607,95	50,66	101,33	212,78	303,98	607,95
T, К	400	250	258	286	458	736	246	253	282	457	741	1,8	2,1	1,4	0,2	0,7
	600	289	295	311	334	485	283	288	303	328	476	2,0	2,2	2,4	1,9	1,9
	800	329	334	346	361	428	321	325	336	351	420	2,4	2,6	2,8	2,8	1,8
	1000	369	374	384	395	458	359	363	371	382	427	2,7	3,0	3,3	3,2	6,7

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

- расчёт вязкости при нормальном давлении производился по методу Тодоса (Метод Чепмена в данном случае даёт практически идентичные результаты). Расчёт при повышенных давлениях - по методу Райхенберга.

Таблица 3.4. Сравнение расчётных и экспериментальных значений для вязкости воды при различных давлениях

Эксперимент [35], мкП	Расчётные значения*, мкП	|D|, %
P, бар	50,00	100,00	210,00	300,00	600,00	50,00	100,00	210,00	300,00	600,00	50,00	100,00	210,00	300,00	600,00
P, атм	50,66	101,33	212,78	303,98	607,95	50,66	101,33	212,78	303,98	607,95	50,66	101,33	212,78	303,98	607,95
T, К	400	250	258	286	458	736	244	249	270	406	660	2,4	3,5	5,6	11,4	10,3
	600	289	295	311	334	485	282	286	296	314	428	2,4	3,2	4,9	6,0	11,7
	800	329	334	346	361	428	320	323	330	341	395	2,7	3,3	4,6	5,6	7,8
	1000	369	374	384	395	458	358	361	366	374	409	3,0	3,6	4,6	5,3	10,8

* - расчёт вязкости при нормальном давлении производился по методу Тодоса. Расчёт при повышенных давлениях - по методу Джосси, Стила и Тодоса (ДСТ).

2) Азот, N₂

Таблица 3.5. Сравнение расчётных и экспериментальных значений для мольных объёмов азота

Эксперимент [36], дм3/кг	Расчётные значения, дм3/кг	|D|, %
P, бар	3	5	10	20	40	3	5	10	20	40	3	5	10	20	40
P, атм	3,04	5,07	10,13	20,27	40,53	3,04	5,07	10,13	20,27	40,53	3,04	5,07	10,13	20,27	40,53
T, К	330	326,5	195,9	97,97	49,02	24,58	318,0	190,8	95,4	47,7	23,9	2,6	2,6	2,6	2,6	2,7
	425	420,8	252,7	126,5	63,49	31,98	409,9	246,1	123,2	61,8	31,1	2,6	2,6	2,6	2,6	2,7
	600	594,3	356,9	178,8	89,79	45,29	578,8	347,6	174,1	87,4	44,1	2,6	2,6	2,6	2,6	2,6
	1000	990,2	594,5	297,8	149,4	75,21	964,5	579,1	290,1	145,5	73,3	2,6	2,6	2,6	2,6	2,6

Таблица 3.6. Сравнение расчётных и экспериментальных значений для сжимаемости азота

Эксперимент [38]	Расчётные значения	|D|, %
P, МПа	3	5	10	20	40	3	5	10	20	40	3	5	10	20	40
P, атм	29,61	49,35	98,69	197,38	394,77	29,61	49,35	98,69	197,38	394,77	29,61	49,35	98,69	197,38	394,77
T, К	350	1,0048	1,0098	1,0273	1,0816	1,2450	1,0042	1,0090	1,0270	1,0850	1,2512	0,06	0,08	0,03	0,31	0,50
	450	1,0114	1,0198	1,0431	1,0973	1,2298	1,0104	1,0185	1,0419	1,0997	1,2401	0,09	0,13	0,12	0,22	0,83
	600	1,0130	1,0220	1,0453	1,0946	1,2014	1,0122	1,0208	1,0441	1,0961	1,2123	0,08	0,12	0,11	0,14	0,91
	1000	1,0105	1,0175	1,0353	1,0711	1,1433	1,0102	1,0171	1,0350	1,0726	1,1518	0,03	0,04	0,03	0,14	0,74

Таблица 3.7. Сравнение расчётных и экспериментальных значений для вязкости азота при низких давлениях

	Эксперимент [37], мкП	Расчётные значения, мкП	|D|, %
T, К	[37]	1	2	3	4	1	2	3	4
323,15	188,0	186,4	189,0	-	186,5	0,86	0,53	-	0,78
348,15	198,5	196,6	199,9	-	196,5	0,98	0,70	-	0,99
423,15	228,0	225,3	230,4	-	224,6	1,19	1,07	-	1,48
523,15	264,0	260,6	267,3	-	258,8	1,30	1,24	-	1,97



1. Метод Голубева

2. Метод Тодоса

3. Метод Райхенберга (неприменим)

4. Метод Чепмена

Таблица 3.8. Сравнение расчётных и экспериментальных значений для вязкости азота при различных давлениях

Эксперимент, мкП [37]	Расчётные значения*, мкП	|D|, %
P, атм	20	50	100	150	300	20	50	100	150	300	20	50	100	150	300
T, К	323,15	190,5	195,5	206,0	218,0	264,5	190,4	194,7	205,3	218,8	269,1	0,0	0,4	0,4	0,4	1,7
	348,15	201,0	205,0	214,5	224,5	265,5	201,1	204,8	213,8	225,3	267,8	0,1	0,1	0,3	0,4	0,9
	423,15	230,0	234,0	240,0	247,5	276,0	231,3	233,9	240,0	247,7	276,2	0,6	0,1	0,0	0,1	0,1
	523,15	265,5	268,0	273,0	278,0	298,0	267,9	269,6	273,9	279,2	298,6	0,9	0,6	0,3	0,4	0,2

* - расчёт вязкости при нормальном давлении производился по методу Тодоса. Расчёт при повышенных давлениях - по методу Райхенберга.

Таблица 3.9. Сравнение расчётных и экспериментальных значений для вязкости азота при различных давлениях

Эксперимент, мкП [37]	Расчётные значения*, мкП	|D|, %
P, атм	20	50	100	150	300	20	50	100	150	300	20	50	100	150	300
T, К	323,15	190,5	195,5	206,0	218,0	264,5	193,3	197,5	207,6	221,2	274,9	1,48	1,01	0,75	1,46	3,94
	348,15	201,0	205,0	214,5	224,5	265,5	204,0	207,8	216,5	228,2	274,6	1,51	1,34	0,94	1,64	3,41
	423,15	230,0	234,0	240,0	247,5	276,0	234,2	237,0	243,3	251,4	283,4	1,84	1,30	1,37	1,56	2,68
	523,15	265,5	268,0	273,0	278,0	298,0	270,7	272,9	277,4	283,0	305,0	1,98	1,81	1,60	1,79	2,36

* - расчёт вязкости при нормальном давлении производился по методу Тодоса. Расчёт при повышенных давлениях - по методу ДСТ.

3) Метанол, CH₃OH

Таблица 3.10. Сравнение расчётных и экспериментальных значений для мольных объёмов метанола

Эксперимент*, м3/кг [39]	Расчётные значения, м3/кг	|D|, %
P, бар	25,00	74,00	90,00	100,00	25,00	74,00	90,00	100,00	25,00	74,00	90,00	100,00
P, атм	25,33	74,98	91,19	101,33	25,33	74,98	91,19	101,33	25,33	74,98	91,19	101,33
T, К	463,15	-	-	0,001674	0,001668	-	-	0,00192	0,00191	-	-	14,4	14,6
	493,15	0,016160	-	0,001911	0,001892	0,01583	-	0,00219	0,00217	2,1	-	14,5	14,9
	513,15	0,019680	0,008899	0,002384	0,002275	0,01856	0,00823	0,00270	0,00260	5,7	7,5	13,1	14,4
	573,15	0,025410	0,015740	0,012080	0,010990	0,02440	0,01484	0,01127	0,01024	4,0	5,7	6,7	6,9



Таблица 3.11. Сравнение расчётных и экспериментальных значений для вязкости метанола при низких давлениях

	Эксперимент [39], мкП	Расчётные значения, мкП	|D|, %
T, К		1	2	3	4	1	2	3	4
463,15	129,3	100,57	126,23	132,38	129,81	22,2	2,4	2,4	0,4
493,15	152,8	117,80	151,07	155,05	153,57	22,9	1,1	1,5	0,5
513,15	162,7	125,16	161,71	164,57	163,58	23,1	0,6	1,2	0,5
573,15	169,2	130,06	168,80	170,85	170,17	23,1	0,2	1,0	0,6

1. Метод Голубева

2. Метод Тодоса

3. Метод Райхенберга (неприменим)

4. Метод Чепмена

Таблица 3.12. Сравнение расчётных и экспериментальных значений для вязкости метанола при различных давлениях

Эксперимент, мкП [39]	Расчётные значения, мкП	|D|, %
P, бар	50,00	75,00	100,00	200,00	400,00	50,00	75,00	100,00	200,00	400,00	50,00	75,00	100,00	200,00	400,00
P, атм	50,66	75,99	101,33	202,65	405,30	50,66	75,99	101,33	202,65	405,30	50,66	75,99	101,33	202,65	405,30
T, К	423,15	ж	ж	1562	1700	1931	-	-	1033	1075	1151	-	-	33,9	36,8	40,4
	473,15	ж	ж	958	1108	1333	-	-	740	802	898	-	-	22,8	27,6	32,6
	523,15	179	197	368	710	943	188	209	366	573	707	4,9	6,1	0,6	19,3	25,1
	543,15	185	191	220	580	827	193	207	238	488	644	4,1	8,2	8,0	15,9	22,2



Таблица 3.13. Сравнение расчётных и экспериментальных значений для вязкости метанола при различных давлениях

Эксперимент, мкП [39]	Расчётные значения*, мкП	|D|, %
P, бар	50,00	75,00	100,00	200,00	400,00	50,00	75,00	100,00	200,00	400,00	50,00	75,00	100,00	200,00	400,00
P, атм	50,66	75,99	101,33	202,65	405,30	50,66	75,99	101,33	202,65	405,30	50,66	75,99	101,33	202,65	405,30
T, К	423,15	ж	ж	1562	1700	1931	-	-	927	963	1031	-	-	40,6	43,3	46,6
	473,15	ж	ж	958	1108	1333	-	-	666	723	810	-	-	30,5	34,7	39,2
	523,15	179	197	368	710	943	183	199	323	513	637	2,3	1,0	12,2	27,8	32,5
	543,15	185	191	220	580	827	189	199	223	433	579	1,9	4,0	1,3	25,3	30,0

* - расчёт вязкости при нормальном давлении производился по методу Тодоса. Расчёт при повышенных давлениях - по методу ДСТ.

4) Оксид углерода, CO

Таблица 3.14. Сравнение расчётных и экспериментальных значений для плотности оксида углерода при различных давлениях

Эксперимент [36], кг/м3	Расчётные значения, кг/м3	|D|, %
P, бар	10,00	40,00	70,00	100,00	10,00	40,00	70,00	100,00	10,00	40,00	70,00	100,00
P, атм	10,13	40,53	70,93	101,33	10,13	40,53	70,93	101,33	10,13	40,53	70,93	101,33
T, К	463,15	11,2640	45,34	79,43	113,00	11,563	46,496	81,359	115,584	2,7	2,6	2,4	2,3
	493,15	8,4020	33,31	57,67	81,33	8,628

Страница:

•  1 
•  2 

дипломная работа "Расчет термодинамических свойств газов" скачать

Подобные документы

• Термодинамические свойства реальных газов

  Уравнение Менделеева–Клапейрона - самое простое, надежное и известное уравнение состояния идеального газа. Межмолекулярное взаимодействие в реальных газах, приводящее к конденсации (образование жидкости). Среднее значение его потенциальной энергии.
  
  презентация [1,2 M], добавлен 13.02.2016
  

• Реальные газы, жидкости и твердые тела

  Описание реальных газов в модели идеального газа. Особенности расположения молекул в газах. Описание идеального газа уравнением Клапейрона-Менделеева. Анализ уравнения Ван-дер-Ваальса. Строение твердых тел. Фазовые превращения. Диаграмма состояния.
  
  реферат [1,1 M], добавлен 21.03.2014
  

• Теоретические методы описания эффекта Джоуля-Томсона в рамках уравнения Ван-дер-Ваальса и сопоставление их с экспериментальными данными для простых газов

  Использование уравнения состояния для описания свойств реальных газов в термодинамике. Уравнение Ван-Дер-Ваальса, связывающее давление, молярный объем и температуру. Физическая природа эффекта Джоуля-Томсона. График инверсии по теоретическим данным.
  
  курсовая работа [1014,0 K], добавлен 27.09.2013
  

• Разработка рациональных методов сжигания газа

  Основы теории диффузионного и кинетического горения. Анализ инновационных разработок в области горения. Расчет температуры горения газов. Пределы воспламенения и давления при взрыве газов. Проблемы устойчивости горения газов и методы их решения.
  
  курсовая работа [794,4 K], добавлен 08.12.2014
  

• Термодинамические процессы в реальных газах и парах. Влажный воздух

  Отклонение свойств реального газа от идеального. Расчет свойств реальных газов. Процесс перехода твердого вещества непосредственно в пар. Испарение жидкости в ограниченном пространстве. Определение массы сухого пара во влажном и массы влажного пара.
  
  реферат [246,1 K], добавлен 24.01.2012
  

• Идеальный газ в молекулярно-кинетической теории

  Определение и модель идеального газа. Микроскопические и макроскопические параметры газа и формулы для их расчета. Уравнение состояния идеального газа (уравнение Менделеева-Клайперона). Законы Бойля Мариотта, Гей-Люссака и Шарля для постоянных величин.
  
  презентация [1008,0 K], добавлен 19.12.2013
  

• Законы идеальных газов

  Уравнение состояния газа Ван-дер-Ваальса, его сущность и краткая характеристика. Влияние сил молекулярного притяжения на стенки сосуда. Уравнение Ван-дер-Ваальса для произвольного числа молей газа. Изотермы реального газа и правило фаз Максвелла.
  
  реферат [47,0 K], добавлен 13.12.2011
  

• Расчет параметров газа

  Взаимоотношение объема и давления, оценка влияния изменения объема на значение давления. Уравнение давления при постоянном значении массы газа. Соотношение массы и температуры по уравнению Менделеева-Клапейрона. Скорость при постоянной массе газа.
  
  контрольная работа [544,5 K], добавлен 04.04.2014
  

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

  Изучение сущности, вероятностных характеристик идеального газа, выведение его уравнения. Рассмотрение понятий теплообмена и температуры. Ознакомление с плотностью равновесного распределения молекул в потенциальном силовом поле и распределением Максвелла.
  
  курс лекций [86,0 K], добавлен 29.03.2010
  

• Молекулярная физика

  Содержание основных газовых законов. Свойства классического идеального газа, реальных газов и жидкостей. Понятие и принципы создания тепловой машины. Распределение Максвелла и распределение Больцмана. Сущность вероятности состояния. Перенос в газах.
  
  учебное пособие [569,9 K], добавлен 20.01.2011
  

Другие документы, подобные "Расчет термодинамических свойств газов"

• главная
• рубрики
• по алфавиту
• вернуться в начало страницы
• вернуться к началу текста
• вернуться к подобным работам