Математическое моделирование систем и явлений / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1986.

( I) где - плотистьатомовсортаЯ; Е- полнаяэнергияединицыобъемасмеси; V - среднем£.оссЕаяскоростьмолекулгаза; Нд. - вектордиффузииатомовсортаЯ; ? - векторпотокагнергии; f - плотностьгаза, Равновесныйсоставсмесиие температуравэтихусловияхопределяется поформуле: гденоі - концентрациямолекулсортаснаі-м энергетическомуровне, кото ромусоответствуетвнутренняяэнергия £оі истатистическийвес шс- масса молекулы сортас; кАс" чиолоатомовсорта Авмолекулесортас; П - постояннаяПланка; к ~ постояннаяБольцмана. Параметрыj^(v=l,..,r) иравновеснаятемператураТопределяютсяизуравне нийбалансаатомовиэнерги Этиравенствасовместно с (2) представляютсобой системутрансцендентных уравненийдляопределенияпараметровравновесногораспределения. Система былазаписанавбезразмерномвидеирешаласьметодомНьютона. Предваритель нодоказанаположительностьЯкобианаидругиеусловияприменимостиметода Ньютона. Предложеннымметодомрассчитаны: 1) равновеснаядиссоциациядвухатомныхгомоядерныхмолеул(Н2 , Ьі2 ; С2, к2. Fg, Na2, s2, сі2> к2, Br2 , l2) вширокомдиапазонеатомарныхкон центрацийианергий; 2) равновесие смесей, состоящихизтрехатомарныхкомпонент: паровга логенаидвух сортов щелочныхметаллов, двухгалогеновиодногощелочного металла; предполагалось, чтоможетобразовыватьсявосемьвидовразличныхве ществ; варьировались как атомарнаяплотность, такиэнергетическиеусловия; 3) равновесие сложной смеси с тремяатомарнымисоставляющими, вкоторой полагалось наличие12 веществв газовойфазе. Прификсированнойатомарной плотности(NcstfH sM0:NL = 5:7:8:10 , где»L - числоЛошмицта) изменялась полнаяэнергия системы. Вычисленияпроводилиськаксиспользованиемпрос тейшеймолекулярной модели - гармоническогоосциллятора, такисиспользо ваниемуниверсальной модели, учитывающейвращениямолекулы, ангармонизмко лебаний, электроные степенисвободы привысокихтемпературах /2/. Нари- оункеI приведенырезультатырасчетовравновесногосоставацляуниверсаль ноймодели /2/ (вверху) идлягармоническогоосциллятора(внизу). Всекон центраци даны в ML , равновеснаятемпературавабсолютнойшкале. Номерасо- 1? - тензорнапряжений. Ыоі“ •^(2tfmck'r)j/2 exp - йоіАТ}» (2 ) ci’ (3 ) 63