Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ПОЛУЧЕНИЯ ГИПОФОСФИТА НАТРИЯ ИЗ ФОСФОРНОГО ШЛАМА

Кольцова Э.М. 1 Егоркин А.С. 2
1 ОУ ВПО «Российский химико-технологический университет имени Д.И. Менделеева»,Москва
2 ГОУ ВПО «Российский химико-технологический университет имени Д.И. Менделеева», Москва
Рассматривается построение математической модели процесса переработки фосфорного шлама с получением гипофосфита натрия. На основании экспериментальных данных строится схема механизма синтеза гипофосфита натрия, используемая далее для построения математической модели в виде системы обыкновенных дифференциальных уравнений. Затем приводится решение разработанной математической модели с помощью полу-явной разностной схемы с n-разбиением по времени. Устойчивость разностной схемы проверяется методом тестов. После приведения системы уравнений к безразмерному виду кинетические константы находятся методом случайного поиска. Адекватность полученной математической модели проверяется с помощью критерия Фишера. Делается вывод о том, что полученная модель может быть использована для оптимизации рассматриваемого процесса.
гипофосфит натрия
фосфорный шлам
математическое моделирование
1. Лотош В.Е. Переработка отходов природопользования. – Екатеринбург: Полиграфист, 2007. – 503 с.
2. Мелихов И.В., Меркулова М.С. Сокристаллизация. – М., 1975. – 279 с.
3. Моргунова Э.М., Авербух Т.Д. Изучение процесса синтеза гипофосфита натрия // Журнал прикладной химии. – 1967. – №. 40. – С. 274–284.
4. Моргунова Э.М., Авербух Т.Д. К вопросу о механизме синтеза растворов гипофосфитов // Журнал прикладной химии. – 1967. – №. 40. – С. 1660–1667.
5. Латимер В.М. Окислительное состояние элементов и их потенциалы в водных растворах. – М.: Изд-во ИЛ, 1954. – 402 с.
6. Стрежнев И.В., Моргунова Э.М., Габова Е.Л. Получение гипофосфита натрия их желтого фосфора и едкого натра в присутствии гидроокиси бария // Журнал прикладной химии. – 1963. – №. 36. – С. 1873–1832.
7. Лисицина В.В. Моделирование и оптимизация процесса получения гипофосфита натрия: дис. … канд. техн. наук. – М., 1994.

Одним из способов утилизации отходов фосфорной промышленности, основу которых представляет фосфорный шлам, является его обработка щелочным раствором с последующим получением гипофосфита натрия [1].

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

На разных этапах исследования взаимодействия фосфора с гидроксильными группами высказывались различные мнения о протекании данных реакций. На основании работ [2-6] и экспериментальных исследований [7] схему механизма синтеза гипофосфита натрия можно записать с помощью уравнений реакций следующим образом:

(1)

 (2)

 (3)

 (4)

 (5)

 (6)

 (7)

 (8)

(9)

Наличие реакций (1) и (2) проверялось экспериментально [7]. В связи с тем, что вещество NaOH играет роль катализатора в реакции (1), скорость реакции W1 была представлена в виде:

где C1, С2, С3 - объемные концентрации в растворе, кг/м3; С1 - фосфора; С2 - гидроксида натрия; С3 - гидроксида кальция.

Протекание (3), (4), (7) и (8) реакций подтверждается экспериментальными результатами [7]. Показано наличие в осадках синтеза фосфита кальция и двойной соли NaCa(H2PO2)3. Реакции (5) и (6) экспериментально исследованы и подтверж- дены [3, 4].

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

Синтез гипофосфита натрия протекает при наличии двух фаз: твердой - обусловленной присутствием в фосфорном шламе неорганической составляющей, в которой фосфор находится в расплавленном виде, и жидкой. Математическая модель представляет собой систему обыкновенных дифференциальных уравнений, составленных по каждому веществу, участвующему в процессе синтеза. Обозначим следующим образом объемные концентрации участвующих веществ [кг/м3]:

С1 = [P4], C2 = [NaOH],

C3 = Ca(OH)2, C4 = [NaH2PO2],

C5 = [Na2HPO3], C6 = [Ca(H2PO2)2],

C7 = [CaHPO3], C8 = [NaCa(H2PO2)3].

Запишем выражение скоростей протекающих реакций для принятой схемы:

?

 

 

 

 

где W1-W9 - скорости реакций (1)-(9).

В выражении W8, как и в W1, учитывается влияние гидроксида натрия как катализатора, поэтому концентрация гидроксида натрия в выражении W8 имеет вторую степень.

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

(10)

где - движущая сила процесса извлечения фосфора, кг/кг; m1 - масса фосфора в частице шлама, кг; k0 - константа массоотдачи, кг/(м2 с); S - площадь поверхности частицы шлама (м2), которая рассчитывается следующим образом:

;

Vш - объем минеральной и водной части шлама, м3; - плотность фосфора, кг/м3; - массовая концентрация фосфора в частице фосфорного шлама, кг/кг;
mш - масса минеральной и водной части частицы фосфорного шлама, кг; - массовая концентрация фосфора в растворе, кг/кг.

Объемная концентрация фосфора в растворе С1 связана  следующим образом:

где r - плотность раствора, кг/м3.

Изменение концентрации фосфора в растворе С1 будет выглядеть таким образом:

 (11)

где N - число частиц шлама в единице объе­ма реактора.

Первый член уравнения описывает количество фосфора, приходящее из фосфорного шлама в раствор, а последующие четыре члена описывают расход фосфора в ходе химических реакций (1), (3), (5) и (6). Расходование щелочной суспензии, состоящей из гидроксидов натрия и кальция, представлено следующими двумя уравнениями:

Изменение концентрации гидроксида натрия:

(12)

Уравнение описывает расход гидроксида натрия по (1), (2), (5), (6), (8) реакциям, и образование гидроксида натрия по реакции (7) и поступление щелочи при приливах суспензий, где υq - расход суспензии, м3/с; V - объем реакционной массы, м3; - концентрация гидроксида натрия в подаваемой суспензии, кг/м3.

Изменение концентрации гидроксида кальция:

(13)

где - концентрация гидроксида кальция в подаваемой суспензии, кг/м3.

Изменение остальных компонентов реакций описывается следующими уравнениями:

(14)

(15)

(16)

(17)

(18)

Принятая схема протекания реакций и математическая модель содержит пять обратных связей по гидроксидам натрия (реакция (7)) и кальция (реакция (2)), гипофосфиту натрия (реакции (4) и (8)), фосфиту натрия (реакции (7) и (9)) и гипофосфиту кальция (реакции (2), (4), (9)).

В математической модели коэффициенты γij (где i = 1, 9 - номер, соответствующий индексу вещества; j = 1, 9 - соответствует номеру реакции) находили следующим образом. Записывали материальный баланс образования и расходования всех компонентов, участвующих в реакции. Коэффициенты γij рассчитывали пропорционально массовой доле каждого компонента в каждой рассматриваемой реакции.

Для реакции (1):

Σ прихода массы = Σ расхода массы = Σ массы

 (19)

?

 

 

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

Материальный баланс системы уравнений (1) и (9) с использованием рассчитанных коэффициентов γij полностью сохраняется.

(20)

В окончательном виде математическая модель процесса синтеза гипофосфита натрия выглядит следующим образом:

(21)

(22)

(23)

(24)

(25)

(26)

(27)

(28)

(29)

где

.

Начальные условия системы

 i = 2, 3, ..., 8.

Для решения разработанной математической модели процесса синтеза гипофосфита натрия используем полу-явную разностную схему с n-разбиением по времени. Производную по времени от Сi концентрации в момент времени tn аппроксимировали «правой» конечной разностью:

(30)

Аппроксимацию правых частей системы уравнений (43) производили следующим образом: если рассматривается изменение во времени концентрации Сi, и в правой части уравнения имеются члены, наименьшие по Ci-концентрации (например член , то этот член в правой части уравнения представляется в виде произведения: (то есть используем метод «замораживания»). Таким образом, схема становится линейной, разрешимой относительно  и более устойчивой по сравнению со схемой, где член  аппроксимируется соотношением ,

Разностная схема для системы уравнений (30) имеет вид:

?

Так как полу-неявная разностная схема, аппроксимирующая систему уравнений (30) с первым порядком по времени, является, в целом, нелинейной, то для проверки ее устойчивости аналитического метода (типа метода гармоник) не существует.

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

не превышала 5% для каждого n, где - решение разностного уравнения; - истинное решение дифференциального уравнения; i = 1, 8 - индекс, соответствующий компоненту.

Такому требованию отвечали значения шагов Δt: 0,001; 0,004; 0,008. Так как, в дальнейшем разностная схема должна использоваться для поиска кинетических констант процесса синтеза, для целей оптимизации режимов синтеза, и мы заинтересованы в более быстром счете, то выбрали из полученного набора значений шага наиболее крупный шаг Δt - 0,008.

Для дальнейшей работы система дифференциальных уравнений была приведена к безразмерному виду для того, чтобы при поиске кинетических констант можно было работать с выбранным нормированным шагом Δt = 0,008. Константы математической модели искались стандартным методом случайного поиска.

При проведении расчетов использовали экспериментальный материал 22 синтезов [7]. Температура синтезов поддерживалась в интервале T = 85-95°C, использовался шлам, содержащий 68% фосфора. Количество молей воды на 1 г/атом фосфора вводили оптимально-необходимое для протекания реакций, а именно 18-20 молей.

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

(31)

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

Результатом расчетов являлся следующий набор констант:

k′0 = 0,5, k1 = 3,8, k2 = 120,0,

k′3 = 14,45, k4 = 1,88 k5 = 9,6,

k′6 = 15,0, k7 = 3,0, k8 = 0,012,

k′9 = 1,5

Зная характерные значения констант, найденные ранее, определяем истинные значения констант.

k0 = 0,5 10-6 г/(см2•мин);

k1 = 64000,0 см9/(г3•мин);

k2 = 20000,0 см6/(г2•мин);

k3 = 2400,0 см6/(г2•мин);

k4 = 3,0 см3/(г•мин);

k5 = 1600,0 см6/(г2•мин);

k6 = 2500,0 см6/(г2•мин);

k7 = 5,0 см3/(г•мин);

k8 = 2,0 см6/(г2•мин);

k9 = 2,5 см3/(г•мин).

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

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

где N - число опытов; l - число вычисляемых параметров (в данном случае l = 10); fl = N - l - число степеней свободы дисперсии адекватности.

Дисперсия воспроизводимости рассчитывалась по следующей зависимости:

где m - число параллельных опытов (брали m = 3);

f2 = m - 1 - число степеней свободы дисперсии воспроизводимости.

Используем критерий Фишера:

Для уровня значимости p = 0,05 и f1 = 12, f2 = 2 табличное значение критерия Фишера равно 19,4; расчетное - 11,2. Fcalc < Ftabl, следовательно математическая модель процесса синтеза гипофосфита натрия адекватна описываемому процессу.

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

Работа выполнена в рамках государственного контракта № 11.519.11.5005 с Министерством образования и науки РФ.

Рецензенты:

  • Софиев А.Э., д.т.н., профессор, заведующий кафедрой «Техническая киберенетика и автоматика», ФГБ ОУ ВПО «Московский государственный университет инженерной экологии», г. Москва;
  • Урманчеев С.Ф., д.физ.-мат.н., профессор, директор Учреждения Российской академии наук Института механики Уфимского научного центра РАН, г. Уфа.

Работа поступила в редакцию 19.03.2012


Библиографическая ссылка

Кольцова Э.М., Егоркин А.С. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ПОЛУЧЕНИЯ ГИПОФОСФИТА НАТРИЯ ИЗ ФОСФОРНОГО ШЛАМА // Фундаментальные исследования. – 2012. – № 6-1. – С. 117-122;
URL: http://fundamental-research.ru/ru/article/view?id=29907 (дата обращения: 22.04.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074