CC BY
28
УДК 664.7 ББК 36.821 М-54
Подгорный Сергей Александрович, кандидат технических наук, доцент, доцент кафедры автоматизации производственных процессов ФГБОУ ВПО «Кубанский государственный технологический университет», 350072, г. Краснодар, ул. Московская, 2;
Косачев Вячеслав Степанович, доктор технических наук, профессор, профессор кафедры технологического оборудования и систем жизнеобеспечения ФГБОУ ВПО «Кубанский государственный технологический университет», 350072, г. Краснодар, ул. Московская, 2, тел.: 8(861)275-22-79;
Кошевой Евгений Пантелеевич, доктор технических наук, профессор, заведующий кафедрой технологического оборудования и систем жизнеобеспечения ФГБОУ ВПО «Кубанский государственный технологический университет», 350072, г. Краснодар, ул. Московская, 2, тел.: 8(861) 275-22-79;
Схаляхов Анзаур Адамович, доктор технических наук, доцент, профессор кафедры технологий, машин и оборудования пищевых производств, декан технологического факультета ФГБОУ ВПО «Майкопский государственный технологический университет», 385000, Республика Адыгея, г. Майкоп, ул. Первомайская, 191, тел.: 8(8772)57-04-12.
МЕТОДЫ РЕШЕНИЯ ЗАДАЧИ ПОТЕНЦИАЛОПРОВОДНОСТИ ДЛЯ ЗЕРНА ЭЛЛИПСОВИДНОЙ ФОРМЫ
(рецензирована)
В работе проведено моделирование процесса использование потенциала переноса процесса потенциалопроводности при сушке зерна с учетом особенностей их геометрии с целью определения оптимальных режимов.
Ключевые слова: потенциалопроводность, моделирование, эллипсовидная форма, изопотенциальные поверхности, метод сеток
Podgorny Sergey Alexandrovich, Candidate of Technical Sciences, assistant professor of the Department of Production Process Automation of FSBEIHPE «Kuban State Technological Technology», 350072, Krasnodar, 2 Moscow Str.;
Kosachev Vyacheslav Stepanovich, Doctor of Technical Sciences, professor, professor of the Department of Technological Equipment and Life Support Systems of FSBEI HPE «Kuban State Technological University», 350072, Krasnodar, 2Moscow Str., tel.: 8 (861) 275-22-79;
Koshevoy Eugenii Panteleevich, Doctor of Technical Sciences, professor, head of the Department of Technological Equipment and Life Support Systems of FSBEI HPE «Kuban State Technological University», 350072, Krasnodar, 2Moscow Str., tel.: 8 (861) 275-22-79;
Skhalyakhov Anzaur Adamovich, Doctor of Technical Sciences, associate professor, professor of the Department of Technologies, Machinery and Equipment for Food Production, Dean of the Technological Faculty of FSBEI HPE «Maikop State Technological University», 385000, Republic of Adyghea, Maikop, 191 Pervomayskaya Str., 191, tel.: 8 (8772) 57-04-12.
METHODS FOR SOLVING THE PROBLEM OF POTENTIAL CONDUCTIVITY FOR ELLIPSOID GRAIN
(reviewed)
The process of the use of the transfer potential of the potential conductivity process when drying grain taking into account the characteristics of their geometry in order to determine the optimum conditions has been simulated.
Keywords: potential conductivity, modeling, elliptical shape, isopotential surface, grid method.
Моделирование процесса потенциалопроводности при сушке зерна необходимо для определения механизма и рациональных режимов сушки. Особенностью зерна как объекта сушки является эллипсовидная форма и необходима разработка методов решения учитывающих эту особенность. В работе [1] рассмотрен метод сеток, который позволяет решить эту задачу. В данной работе проводится сравнительная оценка различных методов.
Полученные значения методом сеток сравнивались с аналитическим решением, предложенным Лыковым А.В. [2] для сферической частицы, которое может быть представлено программным методом:
Sph
if Fo . т < 0.0035
1 + 3 ' Fo 1 т - б' ч/Fo
VFo ^ Т'' ‘I —+ 2' Z
I'F
Hi
( ( i- ... ^
I ierfc I --------- II. if т ф 0
l l J Jj
1 otherwise
(1)
r
z r 6
I
-ip E 1 n
I 2 2
! L П -П
2 2 i . , 11
'71 ' Fo Iх J j otherwise
где @sph(^)~ среднее значение потенциала в шаре для момента времени т;
)
erf с 14! d?
ierfc (х) •
n
X
2
Fo(t) = а©-т/Яэф - критерий гомохронности Фурье; - интегральная функция ошибок.
Используя (1) и результаты расчетов методом сеток [1] были проведены сравнительные расчеты изменения поля во времени для шара. Сравнение результатов расчетов представлено на рисунке 1.
Рисунок 1. Сравнение результатов расчета по сетке (MaxNum = 6, MaxTime = 30, Rph = 1,8467 мм, а© = 2 мм2/с, Ат =0.05 с, Ah = 0,30778 мм)
с аналитическим решением Лыкова
Из представленных данных на рисунке 1 видно, что при даже небольшом числе узлов сетки возможно получить достаточно точные решения в широком диапазоне критерия гомохронности Фурье.
Важным параметром при расчете является обобщенный радиус (Яэф), используемый при расчете критерия Фурье, который для шарообразной частицы может быть определен по формуле:
V f • ж • 7? , R rh
шара _ 3 эф _ эф
v
S 4 ■ ж ■ R \
шара эф
R . = з
эф
S
3
(2)
В связи с тем, что коэффициент формы в этом случае равен единице. При деформации осей сфера преобразуется в трехосный эллипсоид, определяющий размер которой можно оценить по параметрам относительной деформации этих осей.
В таком случае удобно использовать соотношения объемов и поверхностей коаксиальных эллипсоидов. Формула расчета объема эллипсоида:
Формулу площади поверхности эллипсоида нельзя выразить при помощи простейших функций, хотя выражение известно:
где LegendreF(x,k), LegendreE(x,k) - эллиптические интегралы Лежандра:
Г 1
LegendreF (x’k) |
/ 2 2
dt LegendreE (x’k) \j 1 k 'sin (t) dt
2 2 k ' sin (t)
(5)
параметры m и 9 определяются следующими соотношениями:
2 2 _ 2
а Ь о
m ( а ’ Ь ’ с) = --'----------- и
2 2 _ 2
b а с
Г
( а * с ) = Arc Sin
2 I с I
(6)
0
0
2
a
Полученные результаты показывают, что в частицах равного объема чем ближе форма частицы к сферической, тем скорость переноса потенциала медленнее.
Для проведения сравнительного анализа представленного сеточного метода решения был использован метод Галеркина [3], который заключается в
поиске пробных функций, удовлетворяющих исходному дифференциальному уравнению и для одномерных задач, позволяет получать полуаналитические решения.
В качестве пробных функций использованы полиномы следующего вида:
Pol c : = 1
P2l 4
3X - 1
P4.S! :=
s4 c2
35' <= ~30-$ + 3
(7)
P6 4 :=
„6 e4 c 2
231 “ 315+ 105“ 5
16
Из этих полиномов образованы пробные функции, которые удовлетворяют граничным условиям задачи, составляем из них ряд, являющийся пробной функцией решения данной краевой задачи:
Do ( X ) := А ’ ( РО ( X ) “ Р2 ( X ) ) + А ’ ( РО ( X ) “ Р4 ( X ) ) + А ’ ( РО ( X ) “ Р6 ( X ) ) (8)
где параметры Aj - функции во времени, а Pi (x) - пробные координатные функции.
Решая матричное уравнение:
1
А := М 'В
где
(РО ( X) - Р2( X))' (РО ( X) - Р2 (X)) dX (РО ( X) " Р4(Х))(Р0(Х) " Р2(Х)) dX (РО ( X) " Р6(Х))(Р0(Х) " Р2(Х)) dX
(РО ( X) " Р2( X)) ‘ (РО ( X) - Р4 (X)) dX (РО ( X) " Р4(Х))(Р0(Х) " Р4(Х)) dX (РО ( X) " Р6(Х))(Р0(Х) " Р4(Х)) dX
(РО ( X) " Р2( X))' (РО ( X) “ Р6(Х)) dX (РО ( X) " Р4(Х))(Р0(Х) " Р6(Х)) dX (РО ( X) " Р6(Х))(Р0(Х) " Р6(Х)) dX
(9)
в
0 Г 1
J
L О
Со - с РО с х ) — Р2 с х ) ) ах
со • с РО ( х ) — Р4 с х ) ) ах
со (РО(Х) — Р6(Х)) ах
а Со - приравнивается к единице, исходя из значения безразмерного избыточного потенциала в начальный момент времени.
В результате были получены значения параметров Ai, которые соответствуют начальному моменту времени.
На рисунке 2 представлен вид пробной функции, описывающий начальное распределение избыточного потенциала от безразмерной координаты ^.
2
2
8
1
1
0
• 1
M
0
0
0
Рисунок 2. Пробная функция, описывающая начальное распределение потенциала
Подставляя ряд из пробных функций в исходное дифференциальное уравнение, и производя интегрирование полученных уравнений в интервале от 0 до 1 (область эллипсоида в безразмерной координате) получаем матрицы, представляющие собой коэффициенты системы линейных дифференциальных уравнений, используемые для расчетов параметра А; во времени.
rrrd d
L; <р<
-dXdX d
1 1
(X) " P2( X)) |'S(X)
-S(X)---(P 0( X) P 2( X))
dX dX
(P0( X) - P 2( X))
rrrd d
L; <p<
-dXdX d
1 1
(X) " P2(X)) Г S(X)
-S(X)---(P 0( X) P 2( X))
dX dX
(P0(X) - P4(X))
rrrd d
L; (p<
-dXdX d
1 1
(X) " P2( X)) |'S(X)
-S(X)---(P 0( X) P 2( X))
dX dX
(P0( X) - P6(X))
rrrd d 1 1
'-------(P 0( X) - P 4 (X)) |'S(X) ■
.d X d X J
-—S(X)'-—(P0(X) P 4( X))
dX dX
1 I rrrd d 1
I II I-----(P0( X) - P6( X)) I S (X ) •
I N LdXdX J
(P 0(X ) - P 2(X))
— S(X) —(P0(X) P6( X))
d X d X
Г (Pо( X) - P2( X)) I d X
rrrd d 1 1
-----(P 0( X) - P4(X)) Г S( X) ■
-d X d X J
-—S(X)'-—(P0(X) P 4( X)) dX dX
1 I rrrd d 1
I II I-----(P0( X) - P6( X)) |' S (X ) •
I || LdXdX J
(P 0(X ) - P4(X))
— S(X) —(P0(X) P6(X))
d X d X
Г (P о ( X) - P4 ( X)) I d X
rrrd d 1 1
------(P 0( X) - P 4 (X)) |'S(X) ■
.d X d X J
-—S(X)'-—(P0(X) P 4( X))
dX dX
1 I rrrd d 1
I M |------(P0( X) - P6( X)) | S (X ) •
I II LdXdX J
(P 0(X ) - P 6(X))
— S(X) —(P0(X) P6( X))
d X d X
I" (P 0 ( X) - P6(X))|dX
dX
dX
V(X)
V(X)
V(X)
dX
dX
dX
dX
dX
V(X)
V( X)
V( X)
dX
dX
dX
dX
dX
V(X)
V( X)
V( X)
dX
dX
dX
Решение системы дифференциальных уравнений проводили методом Рунге Кутта с изменяемым шагом по времени, используя в качестве условия Коши вектор начальных параметров А;, а в качестве матрицы правых частей
( MIC 'x + MIC 'x + MIC 'x 'l
I 0 > О 0 0 >1 1 0 ’ 2 2 |
DC tori := ! MIC -X + MIC -X + MIC -X \
x> | 1-0 0 1-1 1 1 - 2 2 |
MIC -x + MIC 'X + MIC 'X
V 2-00 2-11 2-2 2 J
~ 1
где Mic := E-ч» -M -с
Система решена с помощью стандартной функции Rkadapt в среде MathCAD, что дало возможность получить коэффициенты ряда
Z := Rkadapt (х> 0 > 0Л > 200 > D)
Далее находится среднее значение потенциала по формуле:
г1
I ITz (РО(Х) - Р2 ( X ) ) + Z (РО(Х) - Р4 ( X ) ) + Z -(РО(Х) _ Р6 ( X ) ) Т V( X ) 1 dX
I LL П’1 П’2 П’З JJ
V(X) dX
Используя расчеты моделирования диффузии в эллипсоиде [4] и полученное среднее значение потенциала, как методом сеток, так и методом Галеркина, построили график, сравнивающий изменение избыточного потенциала эллипсоидального тела во времени (рисунок 3).
■ метод сеток[6]
Обобщенное время, Fo ■метод Галеркина
потенциал, [1]
3. Зависимость среднего избыточного потенциала от времени, полученная разными методами
Рисунок
0
1 -
n
0
Как видно из рисунка метод сеток наиболее близко приближен к аналитическому решению потенциалопроводности в эллипсоидальном теле. Так же заметно, что решение методом Галеркина имеет наибольшее отклонение, но, не смотря полученные несовпадения, отклонения считаются допустимыми.
Проведен анализ потенциалопереноса при сушке зерен риса с учетом особенностей их геометрии. Форма зерен различных сортов эллипсовидная с максимальным габаритом между 5 и 9 мм по длине.
В процессе сушки зерен наиболее важными параметрами являются, объем семени (V) и площадь его поверхности (S) [4]. Для их определения использовали модель эллипсоида, оси которого соответствуют габаритам семени. Как установлено ошибка аппроксимации составляет от 3,6 до 4,3 % от экспериментальных величин.
Выводы
Использование изопотенциальных поверхностей для снижения размерности задачи основано на введении размерной координаты нормированной от 0 до 1, покрывающей весь объем эллипсоида.
Установлено, что метод сеток наиболее близко приближен к аналитическому решению потенциалопроводности в эллипсоидальном теле. При этом решение методом Галеркина имеет несколько большее отклонение, которые можно считать допустимыми.
Разработанная методология моделирования позволит определять развития градиентов потенциала в процессе сушки, что необходимо для прогнозирования развития трещиноватости зерна при сушке.
Литература:
1. Лыков А.В. Теория теплопроводности. М.: Высш. школа, 1967. 600 с.
2. Постановка задачи потенциалопроводности для зерна эллипсовидной формы и обоснование решения численными методами / С.А. Подгорный [и др.] // Новые технологии. 2015. Вып. 2. С. 20-28.
3. Флетчер К. Численные методы на основе метода Галёркина / пер. с англ. М.: Мир, 1988. 352 с.
4. Li Z., Ye J., Kobayashi N. Modeling of Diffusion in Ellipsoidal Solids: а Simplified Approach to Solving Some Drying Problems // 16th International Congress of Chemical and Process Engineering (22-26 August 2004). Prague: Czech Republic, 2004.
References:
1. Lykov A.V. The theory of heat conduction / A.V. Lykov. M.: Higher School, 1967. 600 p.
2. Statement of the problem of potential conductivity for ellipsoid grain and reasons for the solving by numerical methods / S.A.Podgorny [and oth.]//New Technologies. 2015. Iss. 2. P. 20-28.
3. Fletcher K. Numerical methods based on the Galerkin method/ Tr. From English. M.: Mir, 1988. 352 p.
4. Li Z., Ye J., Kobayashi N. Modeling of Diffusion in Ellipsoidal Solids: a Simplified Approach to Solving Some Drying Problems // 16th International Congress of Chemical and Process Engineering (22-26 August 2004). Prague: Czech Republic, 2004.
