Теория фазового состояния водно-спиртовой смеси: температура, давление, концентрация

Тут у Стабникова смотри какие модельки есть:
может пытаться приближать состояние жидкости не какой то абстрактной экспонентой,
а построить формулу на основании какой то известной модели из термодинамики
а потом уже в ней твоей чудо прогой коефициенты пытаться подбирать

Rudy, давай я тебе все 32 графика выгружу каждый отдельным файликом?
Это чем то поможет?
Я думаю не с того конца мы её интерполируем.
Формула T(P,C) слишком сложная. Надо сначала её как то разбить на более простые элементы.
Что говорит термодинамика о фазовом равновесии жидкость-пар?

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

Я просто грубо прикинул в чем ошибка и, опять же грубо, ее подправил. А точность стала лучше - посмотри, там приведена абсолютная погрешность, а не относительная. Сотые градуса уже вполне прилично, точнее не нужно.

Я потихоньку сварганю нужные массивы и скорее всего, сделаю интерполяцию по ним. А аппроксимация нужна только для того, чтобы считать на пролете, это я так, для интереса.

Rudy, позырь в каких интересных координатах они графики строят

https://ru.m.wikipedia.org/wiki/Законы_Рауля

http://chem21.info/...69089002182203/

Там спирта нет. Да и точность от балды.
Честно говоря, я себе верю больше. Я попробовал некоторые функции, но они дают худшие результаты. Может каких-то более правильных и не нашел, то точность в сотые градуса вполне устраивает.

Я попробовал некоторые функции, но они дают худшие результаты.Rudy, 25 Нояб. 16, 19:10

А теперь смотри сюды!

Получи свою простую формулу. Ошибка под градус.

Круто!
А почему +273,16 а не +273,15,
Ты этим параметром тоже варьировал?

Нет, но это даст постоянную ошибку в тысячные градуса и рояля не играет.

Кстати, ты видел эти посты?

Добавлено через 9дн. 12ч. 39мин.:

Блин, добавил к калькулятору расчет вакуумной перегонки по новым данным и обалдел. Вижу, у меня температура в дефлегматоре низковата -77.81. Залез в калькулятор выяснить при каком давлении должна быть такая температура - получил 745.52. Полез в инет выяснить какое же сейчас давление - а оно рухнуло и оказалось 745.8. Теперь можно по температуре в дефлегматоре давление измерять

добавил к калькулятору расчет вакуумной перегонки по новым данным и обалделRudy, 26 Нояб. 16, 09:08

Руди, с этого места пожалуйста поподробнее...
Где лежит вымя за которое нужно щупать?  ::)
Может быть там и кнопки местами поменялись?  ;D

Котище - Это еще не проверенный вариант, но попробуй. Файл в приложении

Для Чапаев1945 - вот картинка для Krect=A+B*exp(C*x) для объемных процентов. Коэффициенты - на рисунке.
Погрешности, на мой взгляд, достаточно малы. Но важная зона - от 0 до 30% не прописана

Rudy, спасибо большое.

Но важная зона - от 0 до 30% не прописанаRudy, 09 Дек. 16, 13:13

так для этой хоны есть всего 2 значения. Для 1% и для 25%. А дальше только 30%. Как тут апроксимировать? Вот это вообще непонятный диапазон... Может попробовать продлить существующую апроксимацию  прогнозом. Или еще поколдовать.
Посмотрю вечером. Еще раз спасибо.

есть всего 2 значения. Для 1% и для 25%. А дальше только 30%.Чапаев1945, 09 Дек. 16, 13:29

Но это всё равно будет интерполяция хоть и плохонькая.
А экстраполяция без опорных точек в любом случае хуже.

Но это всё равно будет интерполяция хоть и плохонькая.Kotische, 09 Дек. 16, 13:50

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

так для этой хоны есть всего 2 значения. Для 1% и для 25%.Чапаев1945, 09 Дек. 16, 13:29

Так дай их. Это позволит заметно уточнить параметры.

Rudy,

Так дай их.Rudy, 09 Дек. 16, 22:24

Вот
1   35,00
25   5,55
26   
27   
28   
29   
30   3,00

Если найду точнее и недостающее, то выложу. Хочется на весь диапазон общее уравнение...

Что-то ты совсем меня запутал. У Цыганкова по твоей ссылке (стр. 329, приложение 3) только формулы.
А откуда твои данные?

Rudy, У Цыганкова там формулы по расчету К испарения. А зная К исп уже считаем К рект.
А данные что привел (табличка) из Стабникова (Перегонка и ректификация этилового спирта стр 234), а значение для 1 % спирта (35) выведено обратным перечсетом из К рект. В таблице данные по К испарения. Извини что запутал. Из них уже получал К рект. Просто данных по К рект гораздо меньше чем по К испарения. Вот и делал так. Так что интерполяции и уравнения все по К испарения, а не по Крект. К рект потом отдельно считается. Крект примеси = Кисп примеси/Кисп спирта...

Добавлено через 1дн. 1ч. 52мин.:

Rudy, Вот данные.

Для изоамилового спирта достаточно простой аппроксимации