Математическое моделирование процесса горения капли водоэмульсионного топлива Текст научной статьи по специальности «Математика»

Бюллетень науки и практики /Bulletin of Science and Practice Т. 5. №7. 2019

https://www.bulletennauki.com DOI: 10.33619/2414-2948/44

ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ / PHYSICAL & MATHEMATICAL SCIENCES

УДК 51.73+62.642 https://doi.org/10.33619/2414-2948/44/01

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССА ГОРЕНИЯ КАПЛИ ВОДОЭМУЛЬСИОННОГО ТОПЛИВА

©Абдалиев У. К., ORCID: 0000-0002-8994-722X, канд. техн. наук, Институт природных ресурсов Южного отделения Национальной академии наук Киргизской Республики,

г. Ош, Кыргызстан, abdaliev.u@mail.ru ©Ташполотов Ы., SPIN-код: 2425-6716, д-р физ.-мат. наук, Институт природных ресурсов Южного отделения Национальной академии наук Киргизской Республики; Ошский государственный университет; Ошский государственный социальный университет,

г. Ош, Кыргызстан, itashpolotov@mail.ru ©Садыков Э., канд. техн. наук, Институт природных ресурсов Южного отделения Национальной академии наук Киргизской Республики; Ошский государственный университет; Ошский государственный социальный университет, г. Ош, Кыргызстан, sadykov.erkinbai@mail.ru

MATHEMATICAL PROCESS MODELING OF BURNING A DROP OF AQUEOUS EMULSION FUEL

©Abdaliev U., ORCID: 0000-0002-8994-722X, Ph.D., Institute of Natural Resources in the Southern Branch of the National Academy of Sciences of the Kyrgyz Republic, Osh, Kyrgyzstan, abdaliev.u@mail.ru ©Tashpolotov Y., SPIN-code: 2425-6716, Dr. habil., Institute of Natural Resources in the Southern Branch of the National Academy of Sciences of the Kyrgyz Republic, Osh State University,

Osh State Social University, Osh, Kyrgyzstan, itashpolotov@mail.ru ©Sadykov E., Ph.D., Institute of Natural Resources in the Southern Branch of the National Academy of Sciences of the Kyrgyz Republic, Osh State University, Osh State Social University,

Osh, Kyrgyzstan, sadykov.erkinbai@mail.ru

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

Abstract. In this article considered mathematical design of process of burning of drop of water-emulsion fuel. Drop of emulsion of type water-oil is the difficult system, consisting of fuel, in that evenly as drops microdrops are up-diffused waters being into the drop of emulsion. At burning isolated microdrops of water with emulsion, there is an irreversible exothermic reaction. For such reactions considered nonlinear task, with the use of theory of spatially-homogeneous thermal

Бюллетень науки и практики /Bulletin of Science and Practice https://www.bulletennauki.com

Т. 5. №7. 2019 DOI: 10.33619/2414-2948/44

explosions in the process of burning of fuel. This theory shows a few features, not considered to it. These features include appearance more than of one external area, more than one (internal) layer and nonlinear transformation of compression-expansion. As a result, an asymptotic decision for an incompressible liquid allows to get the hydrodynamic fields of explosion in a vicinity steam bubble drop, and also in the vicinity of temperature front.

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

Keywords: mathematical modeling, process, burning, drop of water emulsion fuel, explosion, exothermic reaction, asymptotic decision.

Капля эмульсии типа вода-масло представляет собой сложную систему, состоящую из топлива, в котором равномерно в виде микрокапель распределены капельки воды. Благодаря этому — микрокапли воды, находящиеся внутри капли эмульсии, в процессе ее прогрева быстрее превращаются в парообразные состояние и образуют паровые пузырьки, чем пленка топлива, которая обволакивает эти пузырьки пара [1-2]. При этом пленка топлива вследствие испарения с поверхности капли непрерывно уменьшается по толщине. В момент, когда давление водяных паров внутри частицы превысит силы поверхностного натяжения пленки, произойдет разрушение поверхности капли, т. е. взрыв, или микровзрыв.

При взрыве частиц эмульсионного топлива непосредственно в топочном объеме происходит дополнительное перемешивание паров топлива с кислородом воздуха вследствие того, что они разлетаются в различном направлении. Это ускоряет процесс горения и возможно само горение эмульсии протекает более бурно и за меньший промежуток времени, чем горение безводного топлива [3-4].

Рассмотрим нелинейную задачу, с использованием теории пространственно-однородных тепловых взрывов в процессе горения топлива [5-6]. Данная теория показывает несколько особенностей, не рассмотренные ранее в исследованных работах. Эти особенности включают появление более чем одной внешней области, более одного (внутреннего) слоя и нелинейного преобразования сжатия-расширения.

С таких позиций рассмотрим изолированную микрокаплю воды с эмульсией, в которой происходит необратимая экзотермическая реакция типа А^Б. Математическую модель задачи представим в следующем виде:

— = - AYe<E/RT ), dt

= QAYe~(E/RT),

dt c

(1)

(2)

с начальным условием Y(0)=Yo, ^0)=^.

Здесь время (^ — независимая переменная, Yи T — соответственно массовая доля реагента A и температуры, c — удельная теплоемкость, Q — теплота реакции, A — предэкспоненциальный фактор, E — энергия активации и R — универсальная газовая постоянная. Все эти величины положительные.

Из (1) и (2) имеем:

Бюллетень науки и практики / Bulletin of Science and Practice Т. 5. №7. 2019

https://www.bulletennauki.com DOI: 10.33619/2414-2948/44

— = -Q,^ dT = -QdY ^ T -T = -Q(Y - Y ),T ,Y - const dY c c c

Q

Qi

Отсюда

T + QY = T + QY0 ^ QY = T + QY -T

c

c c

c

QY = T0 + -T

Учитывая c c , уравнению (2) можно записать в виде

dT = J* Q

dt

T0 + QY0 - Tie -(E/RT),

T (0)=T0

(3)

Используя Т0 в характеристической температуре получим:

t = cT0R e (E/RT0)

0 AQY0E

Уравнение (3) представим в виде

— = -(1 + — - T )e^J (0) = 1 dt — У

здесь p=QY0 /(cT0), s=RT0/E.

Исследуем задачу Коши

^ = -(1 + — - T )e dt -У ^

T-1 —

(4)

T(0)= 1.

(5)

Для начала попробуем найти точное решение. Для этого уравнению (1) запишем в виде

dt =

Р

—(1+Р- T)

1 1

—I—

e s —TdT,

или

, Р -1 1

dt = — e s — — T

1+ Р

1+ Р- T

Л ±

-1

e—TdT,

, Р -1

dt = — e s

1+ Р

\ 1_

л =—

VJ

1 1

t(1+р-t) t

e—TdT,

J

1 -1- 1 + - -1+1- ^

-e sesTdT + —,-^e s sTdT

T T(1 + Р- t)

отсюда имеем:

t В I I 1 -1— 1 T \ + R 1, 1

idt = —\ f|--e sessds 1+ f , +Р чe s+ssds\

J s|JI T I J s(1 + — -s)

c

\

s

s

Бюллетень науки и практики / Bulletin of Science and Practice Т. 5. №7. 2019

https://www.bulletennauki.com DOI: 10.33619/2414-2948/44

t (T )=£<

s

1 T

1 1

—e ssd I — 1 V ss

V ss

— e

>7(1+3)

I т

d

1___

ss s(1 + /)

/

ss

s(1 + /).

л 1

Если ввести обозначение Et(л) = PV I—eydy, то получим

' y

—ад-7

' (t )=

3

s

a'/ts !

s I — eydy — e

1/ y

i i

1 — sT s(l+/) s s(l+/) "

1

=3

s

— f — Ts — s ,

H I - J ) —/dy

V —ад —ад У ✓

e^ dy — e

J — eydy > =

/ y

s(l+/)

f^__^ —__

1 1 sT s(l+/) s s(l+/)

J — J

s(l+3)

y

e y dy

Или

t (T ) = 7^e ^ Et ^ j — e ^+s(l+3) E

/f^ f—s+3)rf1

—sie s E lsJ—e

1 1 j sT s(l + /)y

1

(6)

ss

(i+/),

Мы получили решение задачи (4)-(5) в неявном виде.

Попробуем теперь построить явное асимптотическое решение задачи (4) и (5) по малому параметру в при в^-0.

Сначала построим внешнее асимптотическое решение (5) и будем искать в виде:

где

т = i+st +s2T2 + ...+ sT +...,

т=т (t\т = T (t).

(7)

T—1

sT

Подставляя (7) в выражение e s имеем:

— 1 1

еsT = eT + seT (T2 — T2) + sV1 (T3 — 2TJ2 + T3 + ~T22 — TT2 + 4)+

(8)

+sV1 (T — 2TT — T + 3T2T — T + T3T2 — 2TT + 3T2T — T3T — T +

+-T3 — - T,2T2 + - TT4 — - T6 У+•••

6 2 2 2 1 2 2 1 6

Подставляя (7) и (8) в (4) получим:

T-+ ••• = - (/ — sT- — sT — s3T —•••)( eTl +seTl (T — T-2) +

1

1

+в V1 (Т - 2ТТ + Т3 + -Т2 - ТТ2 + -Т4)+ •••)

Отсюда, приравнивая коэффициенты при одинаковых степенях малого параметра в и учитывая начальное условие (5), получим:

11

1

T

1

1

e

e

s

1

e

1

—ад

—ад

Бюллетень науки и практики / Bulletin of Science and Practice https://www.bulletennauki.com

Т. 5. №7. 2019 DOI: 10.33619/2414-2948/44

^ = TJ 0) = 0 dt

dTL = T -T2)-eTT, T(0) = 0

dT

dt

3 = eT

T -2TT + T3 +1T2 -TT2 +1T41 -eT T(T2 -T2)-^e'1, T(0) = 0

'Г2 1

L2-M

2 1

T -

(9) (10) (11)

Решение задачи (9) представим в виде:

1

1

T = ln-,eT =-, 0 < t < 1.

1 -1

1 -1

Подставляя это выражение в (10) имеем:

T =-T -

-ln

ln-

1 -1 2 1 -1 1 -1 Р(1 -1) 1 -1

интегрируя последнее выражение получим:

rr 1 2 1 о 7 1 о 1 7 1 1 С -,1

T = ln--+ 2 ln--+ 2 + — ln--1---1--, c = -2--.

1 -1 1 -1 - 1 -1 - 1 -1 Заметим, что справедливы асимптотические оценки:

T = ln-^— ,t — 1 1 1 -1

Р

T = O

1-t

,t — 1

Учитывая эти асимптотические оценки, получим следующие соотношения:

( ! Л

T = о

f 1 (1-t)2

,t — 1, Tk = о

1

(1-t)k

,t — 1,к > 3.

Следовательно, справедливо разложение

T = \ + elл — + e2—\тl+—t+...-\-t 1 -t 1-t -

1 -1

при s—>0.

(12)

Ряд (12) является асимптотическим только при 8< 1 — t, т. е. I <1-8. В малой окрестности точки Т=1 теряется свойство асимптотичности. Поэтому в окрестности точки Т=1 введем растянутую переменную а.

Пусть 1 — t = е-а/8, о < а < 1 тогда ^ = 1 е-°'8^ст.

8

Подставляя эти выражения в задачу (4) имеем:

s di e° / s=—(1 + — -ф)

da -V '

ф-1

(13)

или

^ =1 (1+ — -ф) e' da -V '

Чтобы получить ограниченное решение требуем выполнения соотношения

1

1

1

1

2

Бюллетень науки и практики / Bulletin of Science and Practice https://www.bulletennauki.com

Т. 5. №7. 2019 DOI: 10.33619/2414-2948/44

1---с = O(s),8 ^ 0

Ф

Поэтому, решение (13) ищем в виде:

1

Ф:

1 — с

+ 8ф + s Ф + .

(14)

Подставляя (14) в (13), имеем:

(1 — с)

2 +8Ф;+82Ф2 + •••=-

/

1

л 11—с—"^Т^з

1 + / —"--8фl —Гф2 — • • •

1 — с у

--+8ф1+8 ф2+-

V 1—с

так как

1—с—

1 2

-+sфl+s Ф2+^ •

1—с

= е

(1—с) Ф1

+ е(1—с)2ф1(1 — с)2 (ф—ф2 + ф12с)8 +

1 2

+1 ф (—12ффс2 + 4ффс3 + 20фф2с3 — 10с4ф2фЗ + 2с5ф2фЗ +

+10сф2ф2 — 4фф + 12сфф — 20с ф^ — 4сф3 — 8с фз + 2с ф3 + +2с4ф3 + 6с2ф2 + 15с4ф4 — 4с3ф2 — 6с5ф4 + с4ф2 + с6ф4 — 8сф3 + 12с 2ф3 + +ф2 — 2фф2 — 6сф4 + 15с 2ф4 — 4сф2 — 20с3ф4 + 2ф + 2ф3 + ф4)2 + O(s3)

поэтому

■+8ф;+82ф;2 +

•••=1 и+/-L.

/V 3 1—с

8Фl — 8 ф2

(1 — с)

х( е(1-°)2 ф1 + е(1-°)2 ф1 (1 — с)2 (ф — фЗ + фЗс)8 + • • • ), и приравнивая коэффициенты при одинаковых степенях 8, получим:

х

1 1

= 7711 + / —■

1

,(1—с) ф

(1 — с)2 /I 1 — с

ф1 = 1 ф1 /

с г

1+/—

W

1 — с

,(1—с)2Ф1

(1—с)2 (ф2— фЗ +ФЗ с)—ф1е

,(1—с)2ч>1

л

(15)

Решение уравнения (14) представим в виде 1 * /

Ф1 =

-ln

(1 — с)2 (1 — с)((1 + /)(1 — с) —1)~ (1 — с)

-ln (1 — с)

1—

1+/.

1 ,1+/

ln 1+3(1 — с) —■

1

(

-ln

/

имеем:

(1 — с)2 / V 7 (1 — с)2 V1 + /

—с

0 <с<-

ф ~ —(/ +1) ln

— с

1

е

1

1

/

Бюллетень науки и практики / Bulletin of Science and Practice https://www.bulletennauki.com

Т. 5. №7. 2019 DOI: 10.33619/2414-2948/44

Заметим, что ф (о) = 0. Далее ф;=(ф-ф2 + ф2а)-

отсюда имеем ф = +

Ф1

(1 -a)((1 + Р)(1 -a)-1)'

Ф1

справедлива оценка:

(1 -a)((1 + Р)(1 -a)-1)

7 1

Ф2 ~ (Р + 1) -—1 + Р

+ ф2 (1 -a), ф2^0) =

1 - 2—

-ln

-a

1+ Р

-a

аналогично имеем

Фз ~ Ф2 =

(Р + 1)2

-ln

1+Р

-a V

1+ Р

-a

Подставляя найденные асимптотики в (13) имеем:

Ф(^) =

1

1 -a

—ф + s ф2

1

1-a

(

- sc ln

1 + —

-a

+s

2 Р-ln

1+Р

-a V

1 + Р

-a

1 + Р

ln

-a V

1 + Р

-a

Ф( 0 ) = 1 + s21ip + ...

Так как 1 — t = е~а/8, 0 <а< 1, то переменная I не может быть больше 1. Чтобы

построить асимптотическое решение при Т>1 введем еще одну новую переменную 5.

Пусть 5 = (7 — 1)Ёе[(1+р)8 , Т0:)=у(8). Тогда 8

[Э _Ё_

Ж = [ е(1+р)8^, *=1^5=0; 8^0; 8^0.

Уравнение (4) в новой переменной s примет вид:

%=(1+Р-У)

ds

у-(1+Р)

sV(1+—)

(16)

Асимптотическое решение уравнения (16) ищем в виде:

у = 1 + [ + 8(1 + [)2 у, + 82 (1 + [)2 у 2 +... Подставляя (17) в (16) получим:

у; +8у'2+... = —(у + 8у 2 +...)

(17)

(1+s(1+-)v1+(1+-)s2V2 +

Р

2

2

s

Бюллетень науки и практики / Bulletin of Science and Practice Т. 5. №7. 2019

https://www.bulletennauki.com DOI: 10.33619/2414-2948/44

Так как

V1+8V2 + ••

(1 + 8(1 + P)V1+(1 + P)82V2 + „

:eV1+s(v 2 — V (1 + /)) ev +

+ву в2(у3 - 2(1 + р)^у2 + (1 + Р)2 уЗ +

+ •••, 8 —^

Поэтому, после приравнивания коэффициентов при одинаковых степенях в, имеем

+1у 2—(1 + /)у 2^3+1(1 + /)2 v

v' = —vi <?V1

V'3 =—V3eV1 — V1 (Уз — V1 (1 + 3)) е

V

V'3 = —^3eV1 — V2 (V2 — V1 (1 + 3))е

—VeV1 f V3 — 3 (1 + /) vv 2 + (1 + /)2 v3 +1 V 2 — (1 + /) V 2^3 +1(1 + 3)2 v4 ^

(18)

(19)

(20)

Решение уравнения (18) представим в виде

s = sn

rv1 e , I —du, J1 u

s0 — const •

асимптотики имеют вид:

(

s = ■

1--+ O

1

VvfJJ, (s^—ад)

v V V1

s = — ln(—v)+s + lnK + O(v), V ^0 (s^ад),

где к =

K = exP IJ

1 ex — 1

dx

x

V'

J — ln (—s ) — ln (ln (—s )),

s ^ —ад.

отсюда имеем:

у1 = -св'0+ О ( в~2 ') ,5 — да; Аналогично, получим:

у 2 ~схв~\ ' — да; ~ с2в~' ,' — да;

Окончательно имеем

у (') = 1 + Р + в(1 + Р)2 с0в~" + в2 (1 + Р)2 с£~" + •••,' — да, в — 0. Теперь проведем сращивание.

Внутреннее решение ф(С) запишем через внешнюю переменную Р. = (1 + в 1п (1 -1 ))-1 -

\ 1 + Р , Л чЛЛ

1 + —^в 1п (1 -1)

— s(1 + 8 ln (1 — t))— ln (1 + 8 ln (1 — t))

3

JJ

при 8 >0 имеем: 1 — 8ln(1 — t). Внешнее решение T(t) = 1 + 8 ln 1 запишем через с: 1+с.

1 — t

e

Бюллетень науки и практики /Bulletin of Science and Practice Т. 5. №7. 2019

https://www.bulletennauki.com DOI: 10.33619/2414-2948/44

Внутреннее решение ф(С) запишем через переменную 5:

1

( 1

V1 + 0

f(

x In

VV

s In(-05 / s)

Y1 i

-s

V

V

1+ 0

s In (-05 / s)

Y

1 if 1 + 0 -+ s In (-05/ s) --s In (-05/ s)

1 + 0 V F 7 Jl 0

0

ii уУ

при в—>0 имеем:

1 + р-е(1 + 0)21п (-05 / е) - е(1 + 0)21п ((в / 0) 1п (-05 / в)). Упростив это выражение, получим:

1 + 0 - е(1 + 0)2 (1п (-5 ) + 1п (1п (-05 / е))).

Решение запишем через а:

1 + 0 + s(1 + 0)2 у

при в^0 имеем:

2( 1 1 + 0 + s(1 + 0) — In (s / 0)--

V s

И наконец, решение y(s) запишем через t:

с 1 (_Р_

0

0

1 + 0 1 + 0 + ЕЩ

-С

- In

W\

vsv

1 + 0

-С

'У

rt - 1Л

V 8 у

когда е—0 (1>1) имеем: 1+0. Чем быстрее затухают краевые эффекты, тем меньше ошибка асимптотического решения.

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

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

1. Зеленкин В. Г., Боровик С. И., Бабкин М. Ю. Теория горения и взрыва. Челябинск, 2011. 166 с.

2. Редкина Н. И., Ходаков Г. С. Механохимия и технологические свойства водных эмульсий высоковязких нефтепродуктов // Теоретические основы химической технологии. 2002. Т. 36. №4. С. 433-438.

3. Абдалиев У. К., Ташполотов Ы., Арзиев Ж. Физико-технические основы получения водоэмульсионного топливо // Вестник Ошского государственного университета. 2014. №3. С. 113-117.

4. Корягин В. А. Сжигание водотопливных эмульсий и снижение вредных выбросов. СПб.: Недра, 1995. 304 с.

5. Зельдович Я. Б., Баренбдатт Г. И., Либрович В. Б., Махвиладзе Г. М. Математическая теория горения и взрыва. М.: Наука, 1980. 478 с.

6. Померанцев В. В., Арефьев К. М., Ахмедов Д. Б. Основы практической теории горения. Л.: Энергоатомиздат, 1986. 309 с.

Бюллетень науки и практики / Bulletin of Science and Practice Т. 5. №7. 2019

https://www.bulletennauki.com DOI: 10.33619/2414-2948/44

References:

1. Zelenkin, V. G., Borovik, S. I., & Babkin, M. Yu. (2011). Teoriya goreniya i vzryva. Chelyabinsk, 166. (in Russian).

2. Redkina, N. I., & Khodakov, G. S. (2002). Mekhanokhimiya i tekhnologicheskie svoistva vodnykh emul'sii vysokovyazkikh nefteproduktov. Teoreticheskie osnovy khimicheskoi tekhnologii, 36(4), 433-438. (in Russian).

3. Abdaliev, U. K., Tashpolotov, Y., & Arziev, Zh. (2014). Fiziko-tekhnicheskie osnovy polucheniya vodoemul'sionnogo toplivo. Vestnik Oshskogo gosudarstvennogo universiteta, (3), 113-117. (in Russian).

4. Koryagin, V. A. (1995). Szhiganie vodotoplivnykh emul'sii i snizhenie vrednykh vybrosov. St. Petersburg, Nedra, 304. (in Russian).

5. Zeldovich, Ya. B., Barenbdatt, G. I., Librovich, V. B., & Makhviladze, G. M. (1980). Matematicheskaya teoriya goreniya i vzryva. Moscow, Nauka, 478. (in Russian).

6. Pomerantsev, V. V., Arefev, K. M., & Akhmedov, D. B. (1986). Osnovy prakticheskoi teorii goreniya. Leningrad, Energoatomizdat, 309. (in Russian).

Работа поступила Принята к публикации

в редакцию 19.05.2019 г. 22.05.2019 г.

Ссылка для цитирования:

Абдалиев У К., Ташполотов Ы., Садыков Э. Математическое моделирование процесса горения капли водоэмульсионного топлива // Бюллетень науки и практики. 2019. Т. 5. №7. С. 10-19. https://doi.org/10.33619/2414-2948/44/01

Cite as (APA):

Abdaliev, U., Tashpolotov, Y., & Sadykov, E. (2019). Mathematical Process Modeling of Burning a Drop of Aqueous Emulsion Fuel. Bulletin of Science and Practice, 5(7), 10-19. https://doi.org/10.33619/2414-2948/44/01 (in Russian).