Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

MATHEMATICAL SIMULATION OF THE SODIUM HYPOPHOSPHITE PRODUCING PROCESS FROM THE PHOSPHOROUS SLUDGE

Koltsova E.M. 1 Еgorkin A.S. 2
1 Moscow university of chemical technology of Russia, Moscow
2 Moscow university of chemical technology of Russia, Moscow
1409 KB
Under observation is the mathematical model construction of the phosphorus sludge to sodium hypophosphite processing. Mathematical model constructed as a system of ordinary differential equations based on experimentally substantiated sodium hypophosphite synthesis mechanism. Subsequently the mathematical model solved using a semi-implicit finite difference scheme with n-time decomposition. The stability of the finite difference scheme checked by tests. After reducing the system of equations to dimensionless form, the kinetic constants searched with a random search method. Finally the adequacy of the resulting mathematical model checked using Fisher criterion. As a result there is an adequate mathematical model that can be used for phosphorous sludge processing optimization.
sodium hypophosphite
phosphorous sludge
mathematical simulation
1. Lotosh V.E. Pererabotka othodov prirodopol’zovanija. Ekaterinburg: Poligrafist, 2007. pp. 503.
2. Melihov I.V., Merkulova M.S. Sokristallizacija. Moskva, 1975. pp. 279.
3. Morgunova Je.M., Averbuh T.D. Izuchenie processa sinteza gipofosfita natrija // Zhurnal prikladnoj himii. 1967. Vol. 40, pp. 274–284.
4. Morgunova Je.M., Averbuh T.D. K voprosu o mehanizme sinteza rastvorov gipofosfitov // Zhurnal prikladnoj himii. 1967. Vol. 40, pp. 1660–1667.
5. Latimer V.M. Okislitel’noe sostojanie jelementov i ih potencialy v vodnyh rastvorah. Moskva, Izdatel’stvo IL, 1954. pp. 402.
6. Strezhnev I.V., Morgunova Je.M., Gabova E.L. Poluchenie gipofosfita natrija ih zheltogo fosfora i edkogo natra v prisutstvii gidrookisi barija // Zhurnal prikladnoj himii. 1963. no. 36, pp. 1873–1832.
7. Lisicina V.V. Modelirovanie i optimizacija processa poluchenija gipofosfita natrija: Dis. … kand. teh. nauk. M., 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