Математическая модель процесса получения иттриевых сверхпроводников

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

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

Математическая модель процесса получения иттриевых сверхпроводников

Полещук И.Н

К числу важнейших неорганических материалов относятся открытые сравнительно недавно [1] сложные оксиды, обладающие свойствами высокотемпературной сверхпроводимости (ВТСП).

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

Одна из групп веществ, устойчиво демонстрирующие свойства ВТСП, это иттриевые сверхпроводники YBa2Cu3O7-x. Актуальна задача поиска оптимального технологического режима синтеза этого соединения с минимальным содержанием примесей в конечном продукте.

Для получения купрата иттрия и бария необходимое соотношение между иттрием, барием и медью равно 1:2:3 соответственно.

Влияние концентрации металлов на совместную сорбцию целесообразно представить в виде математической модели. Для получения такой модели был использован метод математического планирования и метод крутого восхождения по поверхности отклика. В эксперименте реализован план 23. Использовался центральный композиционный план, отвечающий требованию рототабельности, т.е. план, позволяющий получать модель, способную предсказать значение параметра оптимизации с одинаковой точностью независимо от направления на равных расстояниях от центра плана [2 - 5].

Рототабельность центрального композиционного плана достигается выбором величины “звездного “плеча. Величина “звездного” плеча для ядра, содержащего полный факторный эксперимент, определяется из соотношения:

,

где k -число факторов.

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

сверхпроводник сорбция катионит дисперсия

,

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

При получении центрального плана второго порядка для трех факторов, необходимо к полному факторному эксперименту 23 добавить шесть «звездных» точек с координатами (+; 0; 0), (-; 0; 0), (0; +; 0), (0; -; 0), (0; 0; +), (0; 0; -). Для выполнения эксперимента по сорбции ионов иттрия, бария и меди с к=3 при =6 и = 6 величина «звездного» плеча =1,6818.

Условия опытов о совместной сорбции катионов иттрия, бария и меди на катионите КРФ-10п, матрица плана и статистическая обработка полученных результатов представлены в таблицах 1-4.

Таблица 1. Условия опытов и матрица планирования эксперимента

Таблица 2. Оценка воспроизводимости параметра оптимизации

Таблица 3. Вспомогательная таблица для расчетов коэффициентов регрессии

Таблица 4. Коэффициенты регрессии и оценка их значимости для модели процесса сорбции ионов иттрия, бария и меди катионитом КРФ-10п в Н+-форме

Модель процесса совместной сорбции в работе представлена в виде уравнения регрессии:

1+

Опыты, заданные планом, выполнялись по одному разу. Поэтому дисперсию воспроизводимости по определению параметра оптимизации рассчитываем по результатам шести опытов при одинаковых условиях на основном уровне. Оценка воспроизводимости или дисперсия параметра оптимизации производится по уравнению:

=

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

Расчет коэффициентов регрессии проводился по уравнениям:

; ; ,

где i = 1,2,… k,  - номера факторов; j - номер опыта; N - количество опытов в матрице; уj - значение параметра оптимизации в j - опыте; ,  - кодированные значения (1) факторов i и  в j - опыте.

Оценка значимости коэффициентов регрессии проводилась сравнением абсолютной величины коэффициента с его доверительным интервалом:

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

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

,

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

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

bi =  - доверительный интервал;  =  - средняя квадратичная ошибка в определении коэффициентов регрессии.

Процесс сорбции ионов иттрия, бария и меди на катионите КРФ - 10 п описывается уравнением полинома первой степени в кодированном и натуральном выражении.

Уравнение регрессии:

Уравнение регрессии в кодированном выражении со значимыми коэффициентами:

 = 

Уравнением в кодированной форме неудобно пользоваться, необходимо перевести его в натуральную форму:

где натуральные значения факторов на основном уровне;  натуральные значения факторов; натуральные значения варьирования факторов.

; ; 

Уравнения регрессии, описывающие процесс сорбции:

Гипотезу адекватности полученной модели проверяем по критерию Фишера:

Для расчета дисперсии адекватности используется уравнение:

,

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

Условием адекватности модели процесса сорбции ионов металлов является соблюдение соотношения:

Для иттрия ; 

Для бария ; 

Для меди ; 

Во всех случаях выполняется условие:Следовательно, полученные модели адекватны.

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

Для получения ВТСП гранулята на основе Y-Ba-Cu-O необходимо, чтобы в результате сорбции ионов иттрия, бария и меди на катионите выполнялось мольное соотношение Y:Ba:Си = 1:2:3. Из результатов реализованных опытов видно, что удовлетворяет этим условиям опыт № 17. В нем выполняется соотношение Y :Ва:Си=0,125:0,248:0,375. Полученные уравнения регрессии можно использовать для определения оптимальных условий процесса сорбции методом крутого восхождения. Крутое восхождение начинаем от результатов опыта №17.

Результаты опытов по методу крутого восхождения приведены в табл.5-7.

Таблица 5. Результаты опытов по методу крутого восхождения для иттрия

Наилучшее соотношение получено в опыте №2: Y: Ва: Си = 1: 1,883: 3,083. Результаты реализованного опыта подтверждают расчеты Y:Ва:Си = 1 : 1,866 : 2,984.

Таблица 6. Результаты опытов по методу крутого восхождения для бария

Таблица 7. Результаты опытов по методу крутого восхождения для меди

Таким образом, оптимальные концентрации иттрия, бария и меди соответственно равны СY =0,04 моль/дм3, СВа =0,078 моль/дм3, ССи =0,12 моль/дм3.

Список литературы

1. Bednorz J.G. Muller K.A. Possible high Tc superconductivity in the Ba-La-Cu-O system. Phys. B.-Condensed Matter. 1986. V. 64. p. 189-193.

2. Вознесенский В. А. Статистические методы планирования эксперимента в технико-экономических исследованиях. - М.: Финансы и статистика, 1981, 263 с.

3. Адлер Ю.П. Введение в планирование эксперимента. - М.: Металлургия, 1969, 157 с.

4. Лисенков А.Н. Математические методы планирования многофакторных медико-биологических экспериментов. - М.: Медицина, 1979, 343 с.

5. Рузинов Л.П. Статистические методы оптимизации химических процессов. - М.: Химия, 1972, 199 с.

Размещено на Allbest.ru