Например, Бобцов

ОБЕСПЕЧЕНИЕ ИНФОРМАЦИОННОЙ ДОСТОВЕРНОСТИ И СОГЛАСОВАННОСТИ ТЕРМОДИНАМИЧЕСКИХ ФУНКЦИЙ И БАЗ ДАННЫХ

А.А. Слободов, А.А. Сочагин, Е.Б. Качер, Д.В. Кремнев, М.А. Радин
5 ТЕПЛОФИЗИКА И ТЕОРЕТИЧЕСКАЯ ТЕПЛОТЕХНИКА
УДК 536+541.1
ОБЕСПЕЧЕНИЕ ИНФОРМАЦИОННОЙ ДОСТОВЕРНОСТИ И СОГЛАСОВАННОСТИ ТЕРМОДИНАМИЧЕСКИХ ФУНКЦИЙ
И БАЗ ДАННЫХ
А.А. Слободов, А.А. Сочагин, Е.Б. Качер, Д.В. Кремнев, М.А. Радин
На основании предложенного подхода к ортогональному представлению термодинамических функций веществ разработан метод и построено соответствующее термодинамическое описание температурной зависимости изобарной теплоемкости. Показана его большая эффективность – как с формальных позиций математического и термодинамического описания, так и с точки зрения достоверности получаемых результатов – по сравнению с традиционными подходами. Предложен также метод согласованного описания низко- и высокотемпературных термодинамических данных (теплоемкости, энтальпии и др.). Сочетание подходов позволяет получить наиболее достоверное описание термодинамических свойств веществ. Построена термодинамическая база данных для широкого круга веществ, реализованная в соответствующем программно-информационном комплексе. Ключевые слова: термодинамические функции, информация, моделирование, расчет, база данных, теплоемкость, ортогональное представление, согласование данных, программный комплекс.
Введение
Современный уровень развития термодинамического метода в физикохимических исследованиях дает принципиальную и во многих случаях реальную возможность чисто расчетного изучения самых разнообразных по своей природе систем, исключив необходимость их экспериментального исследования. Применение при этом современного аппарата методов термодинамического моделирования и расчета позволяет рассчитывать на получение количественно достоверных результатов при изучении даже достаточно сложных по составу, свойствам и поведению систем. Этот подход становится особо ценным, когда эксперимент оказывается малоэффективным, труднореализуемым или невозможным – в случаях высокой компонентности исследуемой системы, экстремальности значений параметров состояния, одновременного протекания ряда сложных взаимосвязанных физико-химических превращений и др.
Однако эффективность применения даже абсолютно строгих термодинамических методов расчета во многом ограничена наличием и качеством соответствующих термодинамических данных, определяющих значения параметров применяемых методов. Иными словами, достоверность (как качественная, так и количественная) результатов, получаемых расчетным путем, может быть достигнута только при наличии достаточно адекватной и полноценной информации. Несмотря на кажущееся обилие информации, данные нередко противоречивы, неполны, отсутствуют. Притом даже выявление этих неадекватностей, причин их возникновения (тем более, их устранение) представляет собой самостоятельную, многослойную и отнюдь не простую проблему.
Критерии качества термодинамической информации
Первое требование, которому должны удовлетворять базы и банки термодинамических данных, – достаточная их полнота по спектру содержащихся веществ. Другое условие – полнота набора заданных термодинамических характеристик вещества базы, необходимых для корректного построения его термодинамического описания. И, нако-

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 1(65)

75

ОБЕСПЕЧЕНИЕ ИНФОРМАЦИОННОЙ ДОСТОВЕРНОСТИ И СОГЛАСОВАННОСТИ …

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

Структура термодинамической информации

В целях наиболее всестороннего анализа проблемы и поиска эффективных путей ее решения целесообразно сначала очертить круг термодинамических характеристик, необходимых для построения корректного и полного термодинамического описания произвольной физико-химической системы в допустимой области возможных ее состояний. Для этого может быть использована любая термодинамическая функция вещества – например, его энергия Гиббса G(P, T), которая, используя математический аппарат термодинамики, может быть преобразована к виду:
G0  P,T    f G0  P0,T0   S 0  P0,T0 T T0  

k
          

 СTj 0 T j

P0 ,t

j1 Tj1

t

Tj

dt 

С0 j

T j 1

P 0

,

t

dt



 

T Tj





1

H

0 j

P

P 0

,

T

j

 V0  P0

p,T

dp,

(1)

76 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 1(65)

А.А. Слободов, А.А. Сочагин, Е.Б. Качер, Д.В. Кремнев, М.А. Радин

где

  f G0

P0 ,T0





f

G0 298

,

 S0

P0 ,T0



S0 298

− соответственно энергия Гиббса образова-

ния и энтропия вещества при стандартных (P0 = 1 бар, T0 = 298,15 K) условиях;

Tj



Ttr,

∆H

0 j

(P0,Tj)



∆H

0 tr

– соответственно температуры и стандартные тепловые эф-

фекты фазовых (полиморфных, агрегатных и др.) переходов вещества (если таковые

имеются);

C

0 j

(P0,

t)



C0 Pj

(t)

=

 H0/  t



температурные

зависимости

стандартных

изо-

барных теплоемкостей для каждой из k фаз вещества, образующихся на интервалах

[Tj-1, Tj] (j1:k; T0=To, Tk=T); V0(p, T) – барическая зависимость стандартного мольного

объема на интервале [P0, P] при температуре T.

Анализ даже наиболее авторитетных систем термодинамической информации

[1–5 и др.] вынуждает констатировать невыполнение многих из вышеописанных требований и необходимых характеристик термодинамического описания (1). Описанные

проблемы диктуют необходимость проведения экспертизы и согласования имеющихся

термодинамических данных с целью обеспечения требований, предъявляемых к созда-

ваемым базам данных. Кроме того, в силу нехватки и противоречивости термодинами-

ческой информации остро встает также вопрос о возможностях использования и разра-

ботки методов расчета и оценки недостающих данных.

В связи с этим ниже предлагаются два подхода к описанию термодинамических

данных (функций) и согласованию результатов, полученных при различных условиях

(например, низко- и высокотемпературных), которые открывают новые возможности в решении стоящих проблем качества информации.

Ортогональное представление термодинамических функций

Нехватка и ограниченность возможностей экспериментальных исследований тер-

модинамических свойств многих веществ вынуждает исследователей разрабатывать и

пользоваться различными приближенными, чаще эмпирическими, методами расчета и

оценки. Эта задача особенно актуальна в области повышенных значений давлений и,

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

воречива, ненадежна или вовсе отсутствует.

Полное термодинамическое описание любого вещества или системы (многоком-

понентной) включает, в соответствии со структурой выражения (1), две группы данных

– точечные

значения

(типа

термических

констант

∆fH

0 2

9

8

,

∆fG

0 2

98

,

S

0 29

8

,



(i) tr

H

0

и

др.)

и

функциональные зависимости (по температуре, давлению, концентрации и др.) – изо-

барные

теплоемкости

C

0 P

(T),

изменения

энтальпии

HT–H

0 T

и др. Определение послед-

них – несравненно более сложная задача.

Особо следует остановиться на проблеме приближенного расчета теплоемкости

CP(T) веществ, ее температурной зависимости. В отличие, например, от стандартных энтропии и энтальпии здесь стоит задача оценки функциональных зависимостей (от

температуры), а не отдельных значений, получаемых для стандартных условий (в частности, при фиксированной температуре 298,15 K). Характер температурных зависимостей теплоемкостей очень сильно зависит от природы вещества и далеко не всегда мо-

жет быть с достаточной точностью описан известным трехпараметрическим уравнени-

ем Майера–Келли [6]

CP(T)=a+bT+c/T2

(2)

или ему подобными. Поэтому нередко используются те или иные полиномиальные ва-

рианты его многопараметрического «обобщения», по сути чисто эмпирического. Однако платой за усложнение формы уравнения и повышение точности описания является,

как правило, снижение его универсальности и экстраполяционных характеристик. На

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 1(65)

77

ОБЕСПЕЧЕНИЕ ИНФОРМАЦИОННОЙ ДОСТОВЕРНОСТИ И СОГЛАСОВАННОСТИ …

эти факты неоднократно обращали внимание исследователи [7, 8 и др.], но не более то-

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

Причина отмеченных недостатков лежит в факте взаимозависимости (даже для

простейшего случая (2)) коэффициентов используемого температурного ряда. Поэтому,

добавление (как, впрочем, и исключение) даже одного, например, квадратичного, члена

в уравнение Майера − Келли,

CP(T)=a+bT+c/T2+dT2 ... ,

(3)

влечет за собой не просто его корректировку за счет появления нового (dT2) члена в (3),

но и пересчет вообще всех коэффициентов a, b, c, полученных для исходного уравнения

(2). Таким образом, значения всех коэффициентов варьируются в зависимости от выбо-

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

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

ратурных зависимостей теплоемкостей имеет (причем достаточно ограниченно) прак-

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

структура их одинакова. Физического смысла параметры этих уравнений не имеют во-

все. Имеющаяся же к настоящему времени достаточно обширная, но с различной сте-

пенью точности термохимическая информация [3, 7, 8 и др.] описывается различными

уравнениями как по типу, так и по количеству температурных слагаемых. Поэтому со-

поставление, анализ этих данных и прогнозирование на его основе неизученных тепло-

емкостей ряда соединений, как правило, бесперспективны.

Думается, что выходом из создавшегося положения может служить представле-

ние теплоемкости (и далеко не только теплоемкости!) через систему так называемых

ортогональных (на некотором температурном интервале) функций [9, 10], обладающих

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

CP(T) не зависят от количества используемых функций – при добавлении нового слагаемого «старые» коэффициенты не изменяются. Корректное и эффективное примене-

ние данного подхода – например, для изобарной теплоемкости CP(T) – требует единообразного (по крайней мере, для определенного класса веществ) задания температур-

ных интервалов, на которых определяется CP(T) веществ (и на которых это представление должно быть ортогональным).

Целесообразным и с чисто расчетных позиций, и с точки зрения физического

смысла представляется задание интервала между двумя физически обоснованными, «ре-

перными» температурами. Это, в зависимости от области, для которой определяется

уравнение теплоемкости, либо интервал 298,15 K – точка фазового перехода (плавления

Tm или полиморфного превращения Ttr), либо интервал между двумя соседними точками фазового (полиморфного или агрегатного) перехода. Отметим, что на этой основе целе-

сообразна не только обработка исходных экспериментальных данных, но также может

быть осуществлен «ортогональный пересчет» (т.е. построена искомая система функций),

опираясь на уже имеющиеся зависимости – типа Майера–Келли, например. В основе

расчета и «пересчета» лежит так называемая процедура ортогонализации [9, 10].

Полученное в результате пересчета известных, ранее полученных зависимостей

представление будет иметь то же число слагаемых, что и исходное, однако коэффици-

енты полученных «новых» представлений не будут зависеть от сложности исходной

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

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

турного коэффициента и др.). И тогда сопоставление (с целью, например, выявления

каких-либо закономерностей) соответствующих коэффициентов для соединений опре-

деленного класса будет физически корректным и, как оказалось, очень конструктив-

ным, существенно повышая возможности экспертизы, согласования, сравнительного

анализа, оценивания данных.

78 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 1(65)

А.А. Слободов, А.А. Сочагин, Е.Б. Качер, Д.В. Кремнев, М.А. Радин

Согласованное описание низко- и высокотемпературных данных

Другой достаточно общей задачей построения зависимостей – например, опять для CP(T) – является проблема согласования низко- и высокотемпературных данных. Типична ситуация, когда имеются в целом достаточно надежные данные по теплоемко-

сти

CP(T)

в

высокотемпературной

области

[T1,T2]

и

значение

C

0 P

при фиксированной

температуре T0 (обычно это 298,15 K), но наблюдается рассогласование этих данных,

т.е.

CP(T0)



C

0 P

.

Задачу тогда целесообразно сформулировать так. Необходимо построить новую

функцию

С~P

(T),

проходящую

через

точку

C

0 P

и минимально отклоняющуюся от ис-

ходной зависимости CP(T). При этом в качестве критерия отклонения целесообразно рассматривать изменение энтальпии (которая обычно и измеряется!) на интервале

[T1, T2]. Тогда при использовании квадратичного критерия проблема сводится к решению задачи условной минимизации:

 T2 T1

CP (T

)



CP (T )2

dT



min,

CP (T0)  CP0 .

(4)

Использование указанных выше параметрических представлений изобарных теп-

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

ровать и решить данную задачу. При этом, как можно видеть, использование ортого-

нальных на промежутке [T1, T2] представлений существенно упрощает задачу и позволяет получить наиболее корректные температурные зависимости теплоемкости, учиты-

вающие требуемое согласование низко- и высокотемпературных данных.

Представленная методология была применена для широкого круга соединений.

Анализ результатов показал высокую ее эффективность: ортогональные представления

оказались несравненно более устойчивыми, чем обычные «вариации на тему» полиноми-

альных уравнений типа (2), (3), в плане как описания и экстраполяции данных, так и кор-

реляционных зависимостей внутри классов тех или иных веществ. А совместное рас-

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

веществ построить (особенно на основе ортогональных представлений) согласованное, в

соответствии с критерием (4), их описание, точность и экстраполяционные возможности

которого существенно выше, чем при общепринятых [3, 4, 7, 8 и др.] подходах.

Полученные на основе представленных методов результаты по расчету, оценке,

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

эффективности разработанных методов в таблице представлены результаты расчета ко-

эффициентов A, B, C ортогонального представления CP(T) для оксидов CaO, P2O5 и фосфатов кальция на их основе – Ca(PO3)2, Ca2P2O7, Ca3(PO4)2:

CP (T )  A  Bu2 (x)  Cu3 (x) ,

(5)

где u1(x)1, u2(x), u3(x) – ортогональные функции разложения, получаемые процедурой ортогонализации исходных уравнений типа (2), x – нормированная температурная пе-

ременная (с предварительным указанием температурного интервала). Для метафосфата

Ca(PO3)2 был, кроме того, реализован указанный метод согласованного описания низко- и высокотемпературных данных, что позволило (как подтвердил последующий ана-

лиз) построить в итоге наиболее обоснованную температурную зависимость CP(T). Данные таблицы не только дают указанные важные физические характеристики

соединения (средние на соответствующем температурном интервале значения теплоем-

кости, ее температурного коэффициента и др.). Как нетрудно видеть, для каждого из

соединений ортогональное представление практически удовлетворяют линейно-

аддитивному условию расчета коэффициентов A, B, C (т.е. и CP(T)) фосфатов по исход-

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 1(65)

79

ОБЕСПЕЧЕНИЕ ИНФОРМАЦИОННОЙ ДОСТОВЕРНОСТИ И СОГЛАСОВАННОСТИ …

ным оксидам (аддитивный расчет дает отклонения, практически линейные по составу – отклонения от линейности существенно ниже погрешности эксперимента). Последний факт позволил построить экстраполяционные ортогональные описания CP(T) для ряда экспериментально неизученных фосфатов кальция – CaP4O11, Ca3P4O13, Ca4P2O9, Ca7P10O32 – и убедиться в эффективности такого подхода.

CaO

P2O5

Ca(PO3)2 Ca2P2O7 Ca3(PO4)2

A 54,705 123,084 221,039 259,383 268,491

B 0,444

5,758

–6,775

0,722

24,677

C 1,141

1,220

6,487

5,987

1,492

Таблица. Коэффициенты ортогонального представления CP(T) некоторых фосфатов кальция и составляющих их оксидов

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

Заключение

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

Литература

1. Smith R.M. Critical stability constants / R.M. Smith, A.E. Martell. – N.Y.: Plenum Press, 1976. – V. 4. – 256 p.
2. Wagman D.D. The NBS table of chemical ther-modynamic properties. Selected values for
inorganic and C1 and C2 organic substances in SI units / D.D. Wagman, W.H. Evans,
V.B. Parker, V.B. Schumm, I. Halow, S.M. Bailey, K.L. Churney, R.L. Nuttall // J. Phys. Chem. Ref. Data. – 1982. – V. 11. – Suppl. 2. – 394 p.

80 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2010, № 1(65)

А.Р. Кямяря, Н.В. Серебрянникова, А.А. Залесская

3. Chase M.W. (Ed.) NIST-JANAF Thermochemical Tables: Part I, II. – 4th Ed. / Ed. M.W. Chase, Jr. / J. Phys. Chem. Ref. Data: Monograph № 9. – N.Y.: NSRDS–AChS–AIPh, 1998. – 1151 p.
4. Глушко В.П. (Ред.) Термические константы веществ: Справочник в 10-ти вып. / В.П. Глушко (Ред.) – М.: АН СССР, 1965–1981.
5. Mohr P.J. CODATA Recommended Values of the Fundamental Physical Constants: 2006 / P.J. Mohr, B.N. Taylor, D.B. Newell. – Gaithersburg: NIST, 2007. – 105 p.
6. Maier C.G., Kelley K.K. An equation for the representation of high-temperature heat content data // J. Amer. Chem. Soc. – 1932. – V. 52. – № 8. – P. 3243–3246.
7. Гурвич Л.В., Вейц И.В., Медведев В.А. и дp. Термодинамические свойства индивидуальных веществ. Справочное издание в 4-х тт. – М.: Наука, 1978–1982.
8. Robie R.A., Hemingway B.S. Thermodynamic properties of minerals and related substances at 298.15 K and 1 bar (105 Pascals) pressure and at higher temperatures // U.S. Geol. Surv. Bull. № 2131. – Washington: Dept. Interior, 1995. – 492 p.
9. Суэтин П.К. Классические ортогональные полиномы. – М.: Наука, 1976. – 327 с. 10. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального ана-
лиза. – М.: Наука, 1989. – 486 с.

Слободов Александр Арсеньевич Сочагин Александр Андреевич Качер Екатерина Борисовна Кремнев Дмитрий Валерьевич

– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, доктор химических наук, профессор, aslob@ton.net.ru
– Санкт-Петербургский государственный технологический институт (технический университет), аспирант,
striker112@list.ru – Санкт-Петербургский государственный университет ин-
формационных технологий, механики и оптики, аспирант,
kacherkaterina@list.ru – Санкт-Петербургский государственный технологический
институт (технический университет), аспирант,
spbscout@narod.ru

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2010, № 1(65)

81