BemnuxßtyWT/Proceedings of VSUET, Т. 79, № 2, 201L.
Оригинальная статья/Original article_
УДК 51-74
DOI: http://doi.org/10.20914/2310-1202-2017-2-30-36
Моделирование процесса сушки древесного сырья для последующего использования его при копчении рыбы
Сергей Б. Луковкин 1 lukovkinsb@mstu.edu.ru Людмила И. Хохлова 1 xoxlovaluda@mstu.edu.ru Сергей А. Шиманский 1 shimanskiysa@mstu.edu.ru
1 Мурманский государственный технический университет, Морской институт, ул. Спортивная, 13, г. Мурманск, 183010, Россия Реферат. Модели многих технологических процессов (сушка древесины, копчение и сушка рыбы) описываются уравнениями в частных производных. К таким моделям относится и рассматриваемая здесь модель процесса сушки сырья на примере сушки древесины. Рассматриваемая в статье модель сушки сырья, основана на уравнении теплопроводности с переменными коэффициентами. Физические свойства древесины в процессе сушки изменяются, поэтому и коэффициенты уравнения изменяются по некоторому заданному закону. Модель становится существенно нелинейной. Как правило, получить решение для таких моделей в явном виде невозможно. В этом случае становится актуальной задача разработки численных методов решения модельных уравнений и написание соответствующих программ для проведения компьютерного моделирования технологических процессов. Наличие таких программ позволяет исследователям решать задачи по подбору оптимальных параметров и оптимальных условий протекания изучаемых технологических процессов, определять необходимое время достижения желаемых значений целевых функций. На основе полученных расчётных формул разработана программа моделирования теплового поля в пакете MatLab, которая может быть использована для моделирования процесса сушки древесины при различных условиях. В статье приводятся примеры модельных вычислений для различных значений параметров сырья. Технологические процессы, связанные с нагревом или сушкой исходного сырья, наблюдаются при копчении или при сушке рыбы. Предложенный здесь алгоритм может быть, после некоторой модификации, использован и для их моделирования.
Ключевые^лова:модели22^ании1™£ии£иими2ииинии1иехииии^
Modeling of drying process of wood raw material for its subsequent _use at fish smoking_
Sergei B. Lukovkin 1 lukovkinsb@mstu.edu.ru Ludmila I. Hohlova 1 xoxlovaluda@mstu.edu.ru Sergey A. Shimanskiy 1 shimanskiysa@mstu.edu.ru
1 Murmansk State Technical University, Maritime Institute, Sportivnaya Str., 13, Murmansk, 183010, Russia Summary. Models of many technological processes (drying, Smoking and drying fish) are described by equations. Such models and applies the model of drying process of raw materials for example wood drying. The article considered a model of drying of raw materials based on the heat equation with variable coefficients. Physical properties of wood during the drying process change, and therefore the equation coefficients are changed according to some given law. The model becomes essentially nonlinear. As a rule, to obtain the solution for such models explicitly impossible. In this case, it becomes the actual task of developing numerical methods for solving model equations and writing of appropriate programs for computer simulation of technological processes. The existence of such programs allows researchers to solve the problem of selecting the optimal parameters and the optimum conditions of the studied technological processes to determine the time required to achieve the desired values of the objective functions. Based on the calculation formulas developed simulation program of the thermal field in MatLab that can be used to simulate the drying process under different conditions. The article provides examples of model calculations for different values of parameters of raw materials. Processes associated with heating or drying of the feedstock, are observed when Smoking or drying fish. The suggested algorithm can be, after
some modifications, used for their modeling._
Keywords: modeling, programming, technology, process, drying, wood, equation, thermal conductivity_
Введение
Во многих странах мира уже достаточно давно ведутся работы по изучению использования альтернативных возобновляемых источников энергии, например таких, как биомасса (древесные отходы, щепа, гранулы из древесины). Процесс сушки древесного сырья включен в технологические процессы производства некоторых продуктов питания, например, при копчении и сушке рыбы [1].
Для цитирования Луковкин С.Б., Хохлова Л.И., Шиманский С.А. Моделирование процесса сушки древесного сырья для последующего использования его при копчении рыбы // Вестник ВГУИТ. 2017. Т. 79. № 2. С. 30-36. (М:10.20914/2310-1202-2017-2-30-3б
В России подобные исследования также набирают ход. На качественное и стабильное горение биомассы влияет её влажность. В настоящей работе предлагается модель для изучения процесса сушки древесины с целью определения его продолжительности, необходимой для достижения нужного уровня влажности материала.
For citation
Lukovkin S.B., Hohlova L.I., Shimanskiy S.A. Modeling of drying process of wood raw material for its subsequent use at fish smoking. Vestnik VGUIT [Proceedings of VSUET]. 2017. vol. 79. no. 2. pp. 30-36. (in Russian). doi:10.20914/2310-1202-2017-2-30-36
Метод решения задачи
Процесс сушки древесной коры происходит при следующих условиях: кора имеет форму бесконечного плоского слоя заданной толщины; осуществляется принудительный нагрев конвекцией высокотемпературного газа одной из сторон этого слоя. В начале процесса кора имеет повышенную влажность. Процесс заканчивается, когда кора становится «достаточно» сухой. Параметры древесины в процессе сушки изменяются, так как зависят от влажности материала [2-4].
Для изучения процесса сушки древесной коры можно предложить следующую математическую модель, основанную на уравнении теплопроводности:
ди д ди,
с ■ р--= — (Я--).
д( дх дх
(1)
где и = и(х. О - температура материала в точке с координатой х в момент времени х изменяется от 0 до S; время ( от 0 до некоторого момента (к, здесь S - толщина древесного слоя. Следовательно, область интегрирования уравнения (1) задаётся неравенствами:
0 < х < £; 0 < (< Ь
(2)
Зададим начальное условие (начальное распределение температуры вдоль оси х): при ( = 0
и (х,0) = и0, (3)
и граничные условия: при х = 0
при х = £
зди г тп
дх
яди=0,
дх
(4)
(5)
В формулах (1)-(5) используются следующие параметры: с - эффективная теплоёмкость; р - плотность материала, и = и(х, () -температура участка коры с координатой х в момент времени Я - коэффициент теплопроводности; а - коэффициент теплоотдачи; Аи = и, — и* > 0 - интервал температур, где ил - температура «сухой» зоны материала, и* - температура «влажной» зоны. Предполагается, что в начальный момент времени весь материал сырой и его температура не превосходит и*.
Затем начинается обдув одной из поверхностей коры горячим газом, что приводит к её нагреву и постепенному «высыханию». При этом параметры модели с, р. Я меняются по некоторым заданным законам. Процесс сушки должен завершиться при достижении некоторого заданного уровня содержания влаги в материале.
Обозначим через у долю влажного материала. Очевидно, что у может непрерывно изменяться от 1 (весь материал влажный) до 0 (весь материал сухой). Закон изменения у можно выбрать в виде:
У =
1, и < и* ил — и
ил — и1*
и„ < и < и,
(6)
0.
и > и,.
Те участки коры, где температура и < и* считаются полностью влажными, а там где и > и, - материал полностью сухой.
Эффективная теплоёмкость можно определить следующим образом:
с =
су, и < и*;
с* ■¥ + сл ■(1 — У).
с,. и > и„.
и* < и < и,; (7)
где с«, са - теплоёмкости влажного и сухого материала; g - доля влаги в элементарном объёме сухого материала; Ь - удельная теплота испарения влаги.
Коэффициент теплопроводности материала определяется по формуле:
Я =
Я*, и < и*;
я*-¥ + Ял ■(1 — У).
я,, и > и,.
и* < и < и,; (8)
где Яш, Я а - коэффициенты теплопроводности влажного и сухого материала.
Плотность материала определяется по формуле:
Р =
Р*
Ру, Ра
и < и*;
■¥ + Рл ■(1 — У).
и > и„.
и* < и < и,; (9)
где р ш, р а - плотности влажного и сухого материала.
Решение уравнения (1) при рассмотренных ограничениях (3)-(5) с учётом зависимостей (7)-(9) может быть получено только численно. Мы предлагаем использовать явную схему метода конечных разностей. Полученное решение при этом будет условно устойчивым.
(10)
Учитывая зависимость А = Д(и(х, г)) в виде (8) можно выразить
дЯ= и;
ёх йи дх Тогда уравнение (1) можно записать в виде:
ди дЯ ди „ д2и
с ■ р--=----+ Я--—,
дг дх дх дх
д 2и
дх2
(11)
или, с учётом (10), в виде:
ди дЯ .ди. 2 „
с ■ р--=--(-)2 + Я
дг ди дх
Для вычисления -дЯ будем использовать
ди
зависимость А = Д( у (и)):
дЯ _ дЯ ду _ (Я* — Я,) _ АЯ ди ~ ду ди ~ (и* — иа) = Аи' (12) АЯ = Яа — Я*
Проведём разбиение области интегрирования (2) следующим образом: равномерно разобьём отрезок [0, £] с помощью J узлов по переменной х; равномерно разобьём отрезок [0, (к] с помощью N узлов по переменной г. Вычислим шаг Ь по переменной х и шаг Аг по времени г:
к = £ /у — 1), А( = гк / (N — 1). (13)
Узлы сетки по г и по х:
гп = (п — 1)^ Аг, п = 1..N;
х. = (; — 1), ; = и.
Функцию и вычисленную в узле сетки х_]) обозначим и;
Заменим в уравнении (11) производные функции и(х, г) по формулам:
ди ~ и;+1 — иди „и*+1 — и^
дг ~
(14)
Аг ' дх
д2и ^ и;+1 — 2и;+и;_1
дх2 ~ к2
(15)
Подставим эти выражения в уравнение (11): с-р^
и;+1—и;_=я . (и;+1—и; )2+
Аг
ди
и; — 2и; + и;1 +Я^ --^)
(16)
Используя значения и;.+1, и;, и;_ 1 в слое г = ^ из соотношения (16) можно найти значение и".+1 в слое ^+1:
^ дЯ ^
и;+1 = и; +-
А
срк2
ди
+Я
и", — и"
;+1 ;
( и;, — 2и; + и; ^
V
V
к1
(17)
В формуле (17) индекс; изменяется от 2 до J-1. Для вычисления значений и в точках с индексами (; + 1, 1) и (; + 1, J) надо использовать граничные условия (4) и (5).
Рассмотрим граничное условие при х = 0. Разложим и(х, г) в ряд Тейлора в точке х = Ь:
и (к, г) = и (0, г) + к2 д2и
+к
ди
дх
2 дх2
о(к );
(18)
из формулы (18) выразим
ди
дх
и; — и; к д2и
2 дх2
Из уравнения (11) найдём
д 2и = 1 ( с!Ц —дЯ ди )2Ч
~дхГ "Я с 'р'~дТ ~~дй' ~дх
(19)
(20)
и подставим выражение (20), вычисленное при х = 0, в (19):
ди и; — и;
дх х=0 к
Г ; + 1= ; / ; Л (21)
и;+1 — и; дЯ | и; — иу ср-
к 1
2 Я
ди
дх
А ди I к )
\ 4 у /
Подставим выражение производной в краевое условие (4):
(
—Я
и; — и;
\
к 1 2 Я
и;+1 — и;
ср-
дЯ и; — и: ди
к
= а(итх —и;), (22)
Выразим из (22) и;+1:
и1;+1 = и;
Аг
с ■р
( \ 2а
—±1(и; — и ) +
7 V 1 тах'
к
2 Я
+2-{и; — и;) +
дЯ
ди
и; — и;
(23)
/ /
х=0
Формула (23) позволяет вычислить значение функции U в левом конце интервала при х = 0 в момент t = t+Дt, если известны значения U в момент времени t при х = 0 и х = Ь.
Рассмотрим теперь граничное условие при х = S. Поступим аналогичным образом, разложим Щ(х, 0 в ряд Тейлора в точке х = S при заданном t:
Щ(Я - к, t) = Щ(5, t) -
-к •
дЩ дх
+ -
к2 д Щ
Выразим из (24)
2 дх
дЩ дх
+о(к2)
(24)
дЩ дх
к д 2Щ
Щ(Я, t) - Щ(Я - к, t)
и; - щ;_1
к д Щ
2 дх2
Значение
2 дх2
(25)
д 2Щ
дх2
найдём из уравнения
(11) при х = Я, учитывая краевое условие (5): с • р дЩ
д 2Щ
дх2
X дt
_ ттп+1 ттп
(26)
дЩ дх
X At
Подставим (26) в (25), учитывая, что = 0:
и1 и1 -1
Выразим
_ ттп+1 ТТП с•р -
X
At
= 0;
2Х At
иг : и;+1 = и; + ——(Щ^ - и;). (27) к ср
Формулы (17), (23) и (27) задают вычислительный процесс для определения температуры Щ(х, t+Дt) при известных значениях Щ(х, 0.
Для компьютерного моделирования процесса сушки была разработана программа в среде МаАаЬ (см. приложение)
Программа работает в соответствии с алгоритмом, описанным выше. Основная программа формирует вектор параметров
модели, массивы времени t и координаты х. Функция nаgrеv2 проводит вычисления температуры в узлах сетки в соответствии с формулами (17), (23) и (27). Функция param вычисляет параметры модели в соответствии с формулами (6)-(9). Нужно отметить, что начальное распределение температуры вдоль оси х может быть нераномерным, т. е. оператор под номером 24 в функции nаgrеv2 можно задать так, что значения и(0, х), будут меняться по заданному закону. Например,
. х(1: J) - 0.25
и (1,1: J) = и 0-1--Г (28)
хи0 + 0.1ои 0;
(
и (1,1: J) = и 0abs
х(1: J) - 0.25
х( J)
(29)
хи0 + 0.1и0;
Функция (28) - «перевёрнутая» парабола, сдвинутая вверх на 1.1и0 и вправо по оси х на 0.25. Функция (29) - «модуль» х, сдвинутый вверх на 1.1и0 и вправо по оси х на 0.25.
Расчёты проводились при следующих значениях исходных данных: время t = 0.60 С, число дискретов 201; для переменной х толщина слоя древесины Я = 0.5, число дискретов 21; при х = 0 температура обдува 126 градусов; коэффициент теплоотдачи а = 20; начальное распределение температуры при t = 0: и0 = 30; температура влажной и сухой зоны: Щщ = 40 и Щ = 80; теплоёмкость влажной и сухой древесины Сщ = 3.6 и Са = 2.2; плотность влажной и сухой древесины р щ = 600 и р а = 400; теплопроводность влажной и сухой древесины Хк = 0.35 и Хй = 0.20.
На рисунке 1 представлены графики изменения температуры вдоль оси ОХ через равнее промежутки времени (0.3 с) от 0 до 60 при: а) постоянном и0 = 30 начальном значении температуры; Ь) при начальном распределении температуры в соответствии с формулой (28); с) при начальном распределении температуры в соответствии с формулой (29): а) при начальном значении, равном конечному распределению температуры, полученном в 1 случае, но «перевёрнутым» относительно оси ОХ. Это процесс сушки древесины, при котором древесина сушится сначала с одной, а затем с другой стороны.
Рисунок 1. Графики изменения температуры вдоль оси ОХ Figure 1. Graphs of temperature variation along the axis OX
Результаты и обсуждение
На основе полученных расчётных формул (17), (23) и (27) написана программа в среде МatLаb и проведено компьютерное моделирование процесса сушки древесины. Результаты работы программы для различных начальных условий распределения представлены в виде графиков распределения температуры u(t, x) (ось OY) вдоль оси ОХ при различных значениях времени t.
Важным моментом подобных расчётов является выбор шага дискретизации по времени t и пространственной переменной x [5-6]. В этом варианте использовалась явная схема сеточного метода, которая является условно устойчивой. В условиях классической постановки задачи для уравнения теплопроводности с постоянными коэффициентами шаг по x и по t рекомендуется выбирать так, чтобы значение выражения (1-2s) было положительным, где s = k ■ Dt / (Dx)2, а коэффициент k - удельная теплопроводность материала. Это означает, что s < 0.5 и шаги по t и по x должны быть согласованы. Для этой задачи выбор шага по t и по x делается путём подбора, так как коэффициенты модели меняются.
В качестве дальнейшего развития данной модели можно разработать алгоритм моделирования процесса сушки с применением метода прогонки - это будет неявная схема решения задачи, но более устойчивая. Кроме этого, можно аппроксимировать кусочно-линейные зависимости (7)-(9) с помощью гладких зависимостей типа «арктангенс», что также должно улучшить качество процесса моделирования процесса сушки.
Заключение
В результате проведенного исследования разработана модель процесса сушки древесины, которая, в частности, может быть использована для изготовления коптильных сред, используемых далее по назначению. Модель может быть использована в системах автоматического управления при альтернативных технологиях получения энергии с использованием биомассы.
Благодарности
Выражаем благодарность к.т.н., профессору А.М. Прохоренкову за постановку задачи и консультации в ходе выполнения работы и д.т.н., профессору Ю.В. Шокиной за консультации в области физических свойств древесины и обсуждение полученных результатов.
Приложение
1.
2.
3.
4.
5.
6. 7.
9.
10.
11. 12.
13.
14.
15.
16.
17.
18. 19.
20 21 22
23
24
25
26
27
28
29
30
31
32
33
34
% основная программа моделирования процесса сушки startnagrev2.m S=0.5; % толщина слоя по х tmax=60; % интервал времени по t kmax=201; % количество дискретов по t;
xv=linspace(0,S,21); tv=linspace(0,tmax,kmax); vpar(1)=126; % температура обдува vpar(2)=20; % alfa = коэфф. теплоотдачи
vpar(3)= 30; % начальное распределение температуры вдоль х
40;
80; 3.6
температура влажной
температура сухой зоны теплоёмкость влажной
теплоёмкость сухой
плотность влажной дре-
плотность сухой древе-
vpar(4) зоны vpar(5) vpar(6) древесины vpar(7)= 2.2 древесины vpar(8)= 600 весины vpar(9)= 400 сины
vpar(10)= 0.35 % теплопроводность влажной древесины
vpar(11)= 0.2 00 % теплопроводность сухой древесины uv=nagrev2(tv,xv,vpar); prn3d=input('введите 1 для формирования графика 3D = '); if prn3d==1 surf(xv,tv,uv) xlabel x; ylabel t; zlabel U(x,t); else k=1;
while k<kmax, plot(xv,uv(k,:),'k-') hold on k=k+10; end; end;
%plot(xv,uv(200,:),'r-') grid on
function u = nagrev2(t,x,vp) % здесь проводятся вычисления в соответствии с алгоритмом % см. формулы (17), (23) и (27) umax=vp(1); %температура обдува alf=vp(2); % коэфф. теплоотдачи u0=vp(3); % начальное распределение температуры вдоль х
u1=vp(4); % температура влажной зоны u2=vp(5); % температура сухой зоны
ЛИТЕРАТУРА
1 Kartha S., Larson E.D. Bureau for Development Policy, Energy and Atmosphere Programme // Bi-oenergy Primer: Modernised Biomass Energy for Sustainable Development, UNDP/BDP, 2000, Energy and Environment Group, New York, USA
2 Laurila J., Havimo M„ Lauhanen R. Compression drying of energy wood // Fuel Processing Technology. 2014. Т. 124. С. 286-289.
10. c1= vp(6); % теплоёмкость влажной древесины
11. c2=vp(7); % теплоёмкость сухой древесины
12. % плотность
13. r1=vp (8);
14. r2=vp (9);
15. % теплопроводность
16. L1= vp(10);
17. L2=vp (11);
18. J=length(x);
19. N=length(t);
20. h=mean(diff(x));
21. dt=mean(diff(t));
22. ss=2*dt/(hA2);
23.u=zeros(N,J);
2 4. u(1,:)=u0; % начальное распределение температуры вдоль х
25. for n=1:N-1
26. [c,L,r,pL]=param(u(n,1),u1,u2,c1,c2,L 1,L2,r1,r2);
27. % граничные условия слева: x=0
2 8. upp=alf*(-umax+u(n,1))/L-( u(n,2)-u(n,1) )/h;
29. u(n+1,1)=u(n,1) + (-2*L/h*upp+...
30. +pL*(u(n,2) -u(n,1))A2/(h*h))*dt/(c*r);
31. % граничные условия справа: x=S
32. [c,L,r,pL]=param(u(n,J),u1,u2,c1,c2,L 1,L2,r1,r2);
33. u(n+1,J)=u(n,J) + (dt*2*L)/(c*r*h*h)*(u (n,J-1)-u(n,J));
34. for j=2:J-1
35. [c,L,r,pL]=param(u(n,j),u1,u2,c1,c2,L 1,L2,r1,r2);
36. u(n+1,j)=u(n,j)+dt/(c*r*h*h)*(pL*(u(n ,j+1)-u(n,j))*2+...
37. +L*(u(n,j + 1)-2*u(n,j)+u(n,j-1)));
38. end
39.end
1. function [c,L,r,pL]=param(t,t1,t2,c1,c2,L1,L2, r1,r2)
2. % вычисления по формулам (6) - (8)
3. if t<t1 psi=1;
4. elseif t>t2 psi=0; else psi=(t2-t)/(t2-t1);
5. end;
6. if t<t1 c=c1; elseif t>t2 c=c2;
7. else c=c1*psi+c2*(1-psi);end;
8. if t<t1 L=L1; elseif t>t2 L=L2;
9. else L=L1*psi+L2*(1-psi);end;
10. if t<t1 r=r1; elseif t>t2 r=r2;
11. else r=r1*psi+r2*(1-psi);end;
12. if t<t1 | t>t2 pL=0; else pL=(L1-L2)/(t1-t2); end;
13. end
Zhao J. и др. Modeling conventional drying of wood: Inclusion of a moving evaporation interface // Drying Technology. 2016. T. 34. №. 5. C. 530-538.
4 Tkachov V. S., Yurkov S., Yurkova V. Definition of rational mode of wood diying with the use of adaptive, regression model // Bulletin of Piydniprovs' ka State Academy of Civil Engineering and Architecture. 2014. № 6. С. 12-15.
5 Галкин В. П. и др. Математическая модель системы контроля процесса сушки древесины в поле СВЧ //Вестник Московского государственного университета леса. 2015. Т. 19. №. 1.
'' Гороховский А. Г., Шишкина Е. Е., Чернышев О. Н. Оптимальное управление процессами тепломассообмена при конвективной сушке древесины // Современные проблемы науки и образования. 2014. N°. 6.
REFERENCES
1 Kartha S., Larson E.D. Bureau for Development Policy, Energy and Atmosphere Programme. Bio-energy Primer: Modernised Biomass Energy for Sustainable Development, UNDP/BDP, 2000. Energy and Environment Group, New York, USA
2 Laurila J., Havimo M., Lauhanen R. Compression drying of energy wood. Fuel Processing Technology. 2014. vol. 124. pp. 286-289.
3 Zhao J. et al. Modeling conventional drying of wood: Inclusion of a moving evaporation interface. Drying Technology. 2016. vol. 34. no. 5. pp. 530-538.
СВЕДЕНИЯ ОБ АВТОРАХ
Сергей Б. Луковкин к. филос. наук, доцент, кафедра автоматики и вычислительной техники, Мурманский государственный технический университет, Морской институт, ул. Спортивная, 13, г. Мурманск,
183010, Россия, lukovkinsb@mstu.edu.ru Людмила И. Хохлова к. филос. наук, доцент, кафедра математики, информационных систем и программного обеспечения, Мурманский государственный технический университет, Морской институт, ул. Спортивная, 13, г. Мурманск, 183010, Россия, xoxlovaluda@mstu.edu.ru Сергей А. Шиманский аспирант, кафедра математики, информационных систем и программного обеспечения, Мурманский государственный технический университет, Морской институт, ул. Спортивная, 13, г. Мурманск, 183010, Россия,
shimanskiysa@mstu.edu.ru
КРИТЕРИЙ АВТОРСТВА Сергей Б. Луковкин предложил методику проведения эксперимента
Людмила И. Хохлова выполнила расчёты
Сергей А. Шиманский написал рукопись, корректировал её
до подачи в редакцию и несёт ответственность за плагиат
КОНФЛИКТ ИНТЕРЕСОВ
Авторы заявляют об отсутствии конфликта интересов.
ПОСТУПИЛА 14.01.2017 ПРИНЯТА В ПЕЧАТЬ 10.02.2017
4 Tkachov V. S., Yurkov S., Yurkova V. Definition of rational mode of wood drying with the use of adaptive, regression model. Bulletin of Prydniprovs' ka State Academy of Civil Engineering and Architecture. 2014. no. 6. pp. 12-15.
5 Galkin V.P. et al. Mathematical model of the system of control of process of drying wood in the microwave. Vestnik Moskovskogo gosudarstvennogo uni-versiteta lesa [Proceedings of Moscow state forest University]. 2015. vol. 19. no. 1 (in Russian).
6 Gorokhovskii A.G., Shishkina E.E., Cher-nyshev O.N. Optimal control of heat and mass transfer processes in convective drying of wood. Sovremennye problemy nauki i obrazovaniya [Modern problems of science and education]. 2014. no. 1. (in Russian).
INFORMATION ABOUT AUTHORS
Sergei B. Lukovkin candidate of philosophical sciences, associate professor, automation and computer engineering department, Murmansk state technical university, Maritime institute, Sportivnaya Str., 13, Murmansk, 183010, Russia, lukovkinsb@mstu.edu.ru
Ludmila I. Hohlova candidate of philosophical sciences, associate professor, mathematics, information systems and software department, Murmansk state technical university, Maritime institute, Sportivnaya Str., 13, Murmansk, 183010, Russia, xoxlovaluda@mstu.edu.ru Sergey A. Shimanskiy graduate student, mathematics, information systems and software department, Murmansk state technical university, Maritime institute, Sportivnaya Str., 13, Murmansk, 183010, Russia, shimanskiysa@mstu.edu.ru
CONTRIBUTION Sergei B. Lukovkin proposed a scheme of the experiment
Ludmila I. Hohlova performed computations
Sergey A. Shimanskiy wrote the manuscript, correct it before
filing in editing and is responsible for plagiarism
CONFLICT OF INTEREST
The authors declare no conflict of interest. RECEIVED 1.14.2017 ACCEPTED 2.10.2017