За допомогою числовог моделi, розробленог на базi запропонованих фiзичног та математичног моделей, виконано моделювання процесу сушт-ня флютингу тфрачервоним випромтюванням. Верифжащя числовог моделi показала збiжнiсть результатiв числового моделювання темпера-тури поверхт флютингу, тривалостi та швид-костi суштня з даними фiзичного експерименту. За маси квадратного метра сухого флютингу вiд 0,112 кг/м2 до 0,2 кг/м2 вгдхилення результатiв становило до 5 %
Ключовi слова: пашр для гофрування, флютинг, тфрачервоне випромтювання, вологовм^т, три-
вал^ть суштня, числове моделювання
□-□
С помощью численной модели, разработанной на основе предложенных физической и математической моделей, выполнено моделирование процесса сушки флютинга инфракрасным излучением. Верификация числовой модели показала сходимость результатов числового моделирования температуры поверхности флютинга, продолжительности и скорости сушки данным физического эксперимента. При массе квадратного метра сухого флютинга от 0,112 кг/м2 до 0,2 кг/м2 отклонения результатов составило до 5 %
Ключевые слова: бумага для гофрирования, флютинг, инфракрасное излучение, влагосодержа-ние, продолжительность сушки, численное моделирование
UDC 676.056.5:66.085.1(.43.3)
|DOI: 10.15587/1729-4061.2017.967411
NUMERICAL MODELING OF PHYSICAL FIELDS IN THE PROCESS OF DRYING OF PAPER FOR CORRUGATING BY THE INFRARED RADIATION
A. Karvatskii
Doctor of Technical Sciences, Professor* E-mail: anton@rst.kpi.ua V. Marchevsky PhD, Professor** E-mail: marchevsky@mail.ru O. Novokhat PhD, Senior Lecturer** E-mail: novokhatoleh@gmail.com *Department of chemical, polymer and silica engineering*** **Department of Machines and Apparatus of Chemical and Oil Refinery Productions*** ***National Technical University of Ukraine «Igor Sikorsky Kyiv Polytechnic Institute» Peremohy ave., 37, Kyiv, Ukraine, 03056
1. Introduction
In connection with the growing volume of production and demand for corrugated cardboard [1], there has been a significant increase in the production of paper for corrugating - fluting.
The most energy-intensive and metal-intensive stage in the production of fluting is a drying process on the drying cylinders of paper machine [2]. However, this method has several disadvantages. First of all, it is the high cost and metal consumption of drying equipment. There may also occur the sticking of glue on the surface of drying cylinders after gluing the fluting. Moreover, in the course of water vapor production for the drying cylinders, there are emissions of fuel combustion products. In addition, there are design limitations on the way to intensify the process of drying.
The specified shortcomings are a substantial reason to search for alternative techniques of drying and to develop new, or modernize the existing, designs of drying plants that would ensure energy efficiency, reduce the cost of equipment and decrease anthropogenic impact on the environment.
One of the promising ways of drying is the radiation method using infrared (IR) radiation that has a capacity to penetrate the volume of material, thereby improving efficiency of the drying process [3]. The fluting is rapidly
heated from the heat that is released in its volume when exposed to the IR radiation, reducing the duration of its heating and drying in general. This makes it possible to reduce the number of drying cylinders in a paper machine and, consequently, the cost of water vapor. In addition, when compared to drying cylinders, the equipment that employs IR emitters differs by lower metal consumption and cost. The lack of direct contact between a fluting and IR emitters (heaters) makes it possible to apply the drying process immediately after its gluing. In this case, under condition of using "clean" electricity (for example, nuclear power), the IR heaters do not emit combustion products of fossil fuel into the environment.
Development of industrial plants for drying the fluting requires the knowledge of kinetic regularities in this process, first of all, changes in moisture content and temperature of fluting during its drying. This requires examining the kinetics of drying process of fluting by the radiation method and developing a procedure for calculating its basic parameters required for designing and constructing the appropriate drying equipment.
Therefore, a relevant task is to develop mathematical and numerical models for the process of IR drying of fluting to determine kinetic parameters of the process, which are required to design new equipment with the minimization of the number of physical experiments.
©
2. Literature review and problem statement
Theoretical foundations of the process of radiation-convection drying of colloidal capillary-porous bodies, including paper, in the form of a system of integral-differential equations are described in [4]. However, for the practical application of fundamental provisions of the theory of radiation-convection drying, each particular case requires additional identification and substantiation of radiation model of the environment. It is also necessary to determine the character of borders reflection, conditions of uniqueness, in particular boundary conditions, and numerical method for solving the stated problem, etc. This significantly complicates the practical application of this theory.
The existing mathematical models for the drying process of fibrous material, including paper or cardboard, for the most part describe a contact or a contact-convective heat transfer method.
Thus, for example, [5] describes a mathematical model for drying the fluting by a contact method on the drying cylinders. The model is based on the balance of interrelation between mass and energy.
Similar mathematical model for drying the paper on a multi-cylinder paper machine by the contact method is described in [6].
Common way of drying the paper is a contact-convec-tive one, in particularly sanitary-hygienic paper on a tissue cylinder. Mathematical model that describes this method of drying is given in [7].
Less common are other types of drying. Thus, authors of [8] describe modeling the filtration method for drying of tissue paper with a decrease in water consumption at a paper processing plant.
The described mathematical models of the paper or cardboard drying process with the use of IR radiation do not make it possible to compute the drying process of fluting due to the absence of required kinetic laws.
Thus, authors of [9] describe modeling the drying process by a radiation method. However, the model lacks the distribution law of IR radiation by thickness of the material.
[10] described radiation, convection, and radiation-con-vective drying methods for coated paper. However, there were many empirical dependences employed for the examined paper during modeling. This complicates the application of the presented model to describe the drying for other types of paper. It was not specified how temperature of paper is calculated during drying.
There are a large number of mathematical models that describe the radiation method of drying the products in food industry.
There is a known model of the radiation drying of carrot slices [11]. The article defines the impact of emitters' power on the speed of drying. It is shown that the duration of drying is determined by using the resulting regression equation. However, this technique does not allow defining the kinetic patterns in the drying of other materials, in particular fluting.
[12] described modeling the drying process of potato slices. A comparison of the results of calculation and experiments is presented. A change in moisture content of the material during drying is proposed to define by logarithmic formulas. However, their diversity indicates the need to confirm the possibility of applying these formulas for other types of material.
A mathematical description of kinetic regularities in the drying process of pomegranate seeds is given in [13]. Nevertheless, the method for calculating the radiation distribution by thickness of the material is not specified. There is no a procedure for determining the coefficients of heat and mass transfer, either.
Therefore, it is a relevant task to develop a mathematical model for the radiation method of fluting drying to define important parameters of this process. First of all, these are the temperature and moisture content of fluting, velocity and duration of drying. The accuracy of the obtained data must be of precision sufficient for the engineering calculations of equipment with IR heaters.
3. The aim and tasks of the study
The aim of present work is the numerical modeling of physical fields in the process of radiation drying of fluting and the verification of adequacy of results of the numerical analysis to the data of physical experiments.
To achieve the set aim, the following tasks had to be solved:
- to formulate a physical model based on the analysis of physical processes that occur during radiation drying of fluting;
- to develop, based on the physical representations of the process, a mathematical and a numerical model of the process of drying the fluting by IR radiation;
- to verify the developed numerical model of the process of radiation drying of fluting by the data of physical experiments.
4. Physical model of the process of radiation drying of fluting
The medium between an IR emitter and fluting is accepted as absorbing and emitting and comprises air and water vapor.
Wet fluting is a selective or "grey" body that has reflective, absorbing and transmitting capacity. The specified physical properties change with the change in the ratio of a square meter of the fluting mass and the amount of moisture in it. At the thickness of standard fluting (up to 0.2 mm), almost all of the IR radiation is absorbed in it. Therefore, we can assume that the lower side of fluting is opaque to radiation.
When passing through a layer of fluting, the IR radiation causes the thermal energy release in it, the amount of which decreases in the direction of heat flux. In this connection, there is a gradient of temperature that is directed away from the more heated surface layers of fluting to those less heated inside it. Due to the presence of a temperature gradient, some amount of heat is transferred by thermal conductivity in the direction of the IR radiation passage. However, under conditions of low temperature difference between the surfaces and small thickness of fluting, the amount of this heat is negligible.
Heat from the surfaces of fluting is released into the environment due to convective heat transfer from the surfaces and is removed along with the formed vapor.
When the wet fluting is heated, the process of drying starts. Its intensity is the largest in the surface layers, where
the amount of the released heat is maximal, and it decreases into the depth of a layer of the material. At the surface of fluting, there forms a layer of saturated water vapor.
The driving force behind the process of drying in the first period is the difference between the pressure of the saturated water vapor in the boundary layer at the surface of fluting and the partial pressure of vapor in the environment. The driving force in the second period of drying is the difference between the magnitudes in the current and resulting equilibrium moisture content [14].
During drying, evaporation of moisture can occur not only from the surfaces of fluting, but also from its deeper layers.
Fluting as a capillary-porous body contains free moisture, moisture in micro capillaries and fibers, and adsorption-associated moisture [4]. In the first period of drying, the free moisture evaporates, and the rate of the process is limited by the intensity of heat flow. In the second period of drying, the rate is initially limited by the diffusion of moisture into the area of evaporation (removal of moisture in micro capillaries), and by the end of drying it is limited by the destruction of adsorption bond between moisture and the fluting and its evaporation.
Prior to the start of drying, moisture content and temperature are the same throughout the entire thickness of fluting. During drying, in the surface layers where there is the largest amount of adsorbed heat from the IR radiation, the temperature value is the highest while the moisture content is the lowest. Since, during the IR radiation passage through fluting, the amount of the released heat decreases, then the temperature, accordingly, falls while the moisture content increases. At the end of the drying process, the moisture content reaches an equilibrium value and levels along thickness of the layer of fluting.
In a paper machine, IR heaters are expedient to install on a free run of the paper web between drying cylinders, or replacing them. Under these conditions, inside the paper web and at its boundary with the surrounding environment, the heat is transferred due to the radiation, conductive and convective constituents. The radiation component occurs under the action of IR heaters, conductive - due to heat transfer inside the material by thermal conductivity and convective - due to the contact between a surface of the paper and a gaseous heat carrier. Thus, the heat exchange in a paper web (fluting) is complex - radiation-conductive (RCH). That is why the material of fluting can be assigned to the class of media, semi-transparent to the IR radiation.
5. Mathematical model of the process of radiation drying of fluting
In a general case, the equation of RCH for a semi-transparent medium, by which we accept fluting and a gaseous environment over it in the process of IR drying, is in the divergent form [15, 16]:
q^ = -IgradT = -1VT
is the density vector of heat flux, which is defined by the Fourier law, W/m2; l is the coefficient of thermal conductivity, W/(m-K); V is the Hamiltonian, m-1; qr is the density vector of radiation heat flux, W/m2; qv is the internal source of heat, for example, due to the evaporation and moisture transportation, W/m3.
Divergence of the radiation heat flux is determined by equation [16]:
divqr =J Kv J TvdQ-4nnJ;I0v(T)
dv,
(2)
where Kv and nv are the spectral absorption coefficient (m-1) and index of refraction, respectively; fl is the solid angle, sr; Iov is the Planck function, W-s/(m2- sr; Iv is the vector of spectral intensity of radiation (W-s/(m2-sr)) for direction S in the solid angle dfl; So corresponds to the boundary of the region; v is the radiation frequency, Hz.
The vector of spectral intensity of radiation is determined by the equation of energy transfer [16]:
Iv (s)= Iv (s0 ) exp -J Kvds
v s0
s œ s
+J nVl0vKvexp -J Kvds'' ds'.
(3)
For unambiguous equation of RCH (1), it is necessary to record the appropriate initial and boundary conditions.
Equation (1) refers to the class of integral-differential equations, which is why its solution in a three-dimensional statement by numerical methods is a non-trivial task [15, 16]. However, technological modes and equipment for the radiation drying of fluting allow us to consider this complex process as a one-dimensional one [3, 14]. In this case, the accuracy of such approximation is sufficient to perform engineering calculations when designing the equipment.
Let us consider main assumptions in the one-dimensional problem to formulate the basic interrelations of non-stationary heat exchange in a layer of fluting based on precise equations for determining the radiation component of heat exchange [17]: - temperature field and other fields that are derivatives from it, are one dimensional and vary only in the direction of the Oz axis by thickness of the layer of fluting (Fig. 1);
cp pft=-div (q i+q r )+qv-
(1)
where cp is the specific mass isobar heat capacity, J/(kg-K); p is the density, kg/m3; T is the absolute temperature, K; t is the time, s;
Fig. 1. Schematic of estimated region for the formulation of one-dimensional non-stationary mathematical model of RCH regarding the drying of fluting by IRR: 1 — a layer of air with water vapor; 2 — a layer of wet fluting
- the modeling region consists of two layers - the upper and the lower. The upper layer is a mixture of air with water vapor (gas), and the lower one is a wet fluting. Wet fluting is a capillary-porous body (from here on denoted as the system) that undergoes the process of radiation drying. That is, the system is a two-phase one;
- thermal-physical properties of wet fluting (mass isobar heat capacity, density and coefficient of thermal conductivity) are effective values and depend on temperature;
- the layer of gas is considered to be immobile because temperature of fluting is higher or equal to the surrounding gaseous heat carrier (hence the convective component as the source of heat supply of heat is not taken into account);
- a two-phase system from the point of view of radiation properties is selective, semi-transparent, non-scattering, as well as absorbing and radiating. The upper boundary of the system (IR emitter) is opaque, and a boundary between the layers is translucent. The lower side of fluting is assumed to be opaque for the simplification;
- all the boundaries of the system are diffuse, that is, a reflection coefficient does not depend on the direction;
- radiation properties (absorption rate, index of refraction, the degree of blackness and reflection coefficients of the system's boundaries) are selective (depend on the radiation frequency).
Taking into account the described assumptions, equation of RCH of type (1) for the system gas - porous body for the direction of axis Oz takes the form:
C„„ (T)P„, <T f = l..(T)0+,,„
qvi =
0, i = 1;
-P„ <T)r£, i =2,
T = T(z), 0<z<Sair + Sf = SS,
- on the upper opaque boundary IR heater - wet air (2=0) we accept boundary conditions of the first kind:
T|z=0 = T1 = T(W(4
(6)
- where W(t) is the dependence of change in the electric power of radiation heater on the time, W. At this boundary we also assign the degree of blackness e1v = fe1 (v) or reflection coefficient R1 v = 1 - e1 v to determine the radiation component of equation (2);
- at the translucent boundary between the layer of gas and fluting we consider boundary conditions of absolute contact or of the fourth kind:
f(T} = 0; lin ■ q i} = 0,
(7)
(4)
where i=1, 2 is the phase index (for i=1 - a mixture of air with water vapor, for i=2 - wet fluting (Fig. 1)); cpieff <T) cpi eff (T) is the effective mass isobar heat capacity of the ith phase, J/(kg-K); p i eff (T) is the effective density of the i-th phase, kg/m3; X i eff (T) is the effective coefficient of thermal
conductivity of the i-th phase, W/(kg-K); -jqZ" is the divergence of density of radiation heat flux density, W/m3; qr is the density vector of radiation heat flux, W/m2; T is the absolute temperature, K; t is the time, s; z is the coordinate, m;
where {T} = T + - T is the temperature on the left and on the right in the vicinity of translucent boundary;
{n ■ q +}=n --q +-n + ■ q -;
q+i = -i2eff(T)VT, qi=-i1eff(T)VT
is the density vector of heat flux in the vicinity of translucent boundary form the sides of wet fluting and air, respectively, W/m2; V = — is the Hamiltonian in the one-dimension-9z
al case, m-1; n-, n+ are the vectors of external normal to the translucent boundary from the side of wet fluting and air, respectively. At this boundary, it is also necessary to define and assign the reflection coefficient to the layer of gas Rairv = fairv(v) and the reflection coefficient to the layer of fluting Rfv = ffv (v);
- at the bottom opaque boundary, which is the lower side of fluting - environment (z=SE), we consider boundary conditions of the third kind:
n ■<-l2eff <T)VT)| =a2 <T - Ts ),
(8)
where a2 is the coefficient of heat transfer, W/(m2-K); Ts is the ambient temperature, K. At this boundary, we also assign the degree of blackness e2v = fe2 (v) or reflection coefficient
r2,v 1 e2,v.
is the volumetric density of source of heat due to the evaporation and moisture transport of fluting, W/m3; r is the mass
heat of the water vaporization, J/kg; — is the relative rate
dt
of evaporation and moisture transport, s-1 (kg/(kg-s)) [4].
Equation (4) in the effective values of thermal-physical magnitudes takes into account the content of water vapor in the air or moisture in a porous fluting by the rule of additiv-ity according to their relative humidity. Initial conditions for (4), t=0:
The density of radiation flow in the one-dimensional system, which consists of two layers, different in their properties, of selective-emitting and absorptive media, is described by a system of transport equations [17]:
i ^=-Ki,v pv+(z, Hi )-n2Av (T <z))];
Xi^M = Ki,v [I;(z, Hi )-n2Av <T <z))],
(9)
(5)
where Sair, Sf, SE are the thickness of layers of wet air, wet fluting and the total thickness of layers, respectively, m. Boundary conditions for (4), t>0:
where i=1, 2 refers to the wet air and wet fluting, respectively; Iv+ ,Iv- are the vectors of spectral intensity of radiation in the positive and negative directions, respectively, W-s/(m2-sr); Ibv(T(z)) is the Planck function, W-s/(m2-sr); |=|cos 9|; j is the angle between a vector of the intensity of radiation and axis z, rad; nv is the spectral index of refraction; Kv is the spectral absorption coefficient, m-1; v is the radiation frequency, Hz.
In (9) and hereinafter, the absorption coefficients of media are also effective values. For example, for a wet fluting this coefficient is determined by formula:
K2,v = Kfv(1 - w) + Kwvw,
(10)
where Kf,v, Kwv are the absorption coefficients of dry fluting and water, respectively, m-1; w is the volumetric share of moisture in fluting.
The density vector of radiation flow in any part of the system of wet air or fluting is determined by relation [18]:
Vtr [ 1 ]
qr (z) = 2nNJ[I; (z, n,) -1; (z,n,)] n, Idv, (11)
where vtr is the limit of thermal radiation transmittance in the IR region of the spectrum, Hz.
Upon differentiating the integrand expression (11 by parameter z, we obtain a formula for determining dqr, which
9z
is included in equation (4) through partial derivatives: dqr (z)
9z
' = 2nJU
91+
(z ^ i )
91-
iz h )
9z
9z
i !
lv.
(12)
Since the density vector of radiation flow qr (z) coincides with the Oz axis, then in (12) and hereafter in the text, instead of the vector we write its component by z.
To exclude partial derivatives in the integrand expression (12), it is necessary to substitute the appropriate expressions from the system of transport equations (9) that yields as a result:
dqr
9z
-2nKi,v NJ[T;(z,h)-^v-(z,h)]dh -2n2vIbv(T(z))idv. (13)
translucent boundary between the layers of wet air and fluting. Spectral values of magnitude Rair,v are defined from the Fresnel formulas [18]. In this case, it is assumed that n1v=1, similar for the gas medium. Spectral magnitude Rf,v is found from the condition of preserving the balance of radiation energy on the translucent boundary:
(1 - Rair,v ) n2,v = (1 - Rf,v)n2,v ^ 1 - Rair,v = (1 - Rf,v )
(16)
where nj,v=1.
The value of degree of blackness Sf,v=1-Rf,v on the translucent boundary is determined by using Fresnel formulas [18] under conditions n2v>ni,v, n2v>nf,v n2v = nf v, ni,v=1: - perpendicular component of the reflection coefficient:
Rs»=
n2,v-(1 V )
^ + n2,v-(1 )
(17)
- parallel component of the reflection coefficient:
n2,v-(1 )
Rpv(^)=r I-7-r
n2,v^^ n2,v-(1 )
- mean value of the reflection coefficient:
R v(^) = 2 (Rs»+ Rp,v(^)).
(18)
(19)
Then the selective hemispherical degree of blackness of the surface of wet air - fluting in the direction of wet air will be determined by formula:
1
eairv = 2j(1 - Rv(|i))id|a, Rair,v = 1 - eair,v.
The degree of blackness of the surface of wet air in the direction of fluting is defined based on (16):
(20) fluting
To determine magnitudes Iv+, Iv-, which are included in (12), (13), it is necessary first to record boundary conditions for the radiation on the opaque and translucent boundaries. Under conditions of the accepted model of diffuse reflection, boundary conditions take the form [17]:
I; (0 ) = n2,v (1 - R1,v ) Ibv (Tz=o ) + 2R1,v J Î; (o, ) ^ ;
0
I- (s-ir )= 2Rair,vjlv+(S-ir, ^ + 2 (1 - Rf,v)jl;(s +ir, ^ )M^2,
0 0
I+ (s+ir ) = 2Rfji;(S +ir, ^ + 2 (1 - Rair,v)jlv+(s-ir, ^
0 0 1
Iv (Ss) = n2,v (1 - R2,v )Ibv (Tz=Ss )+ 2R2,v J Iv+ (Ss, ^
where Tz_o, Tz_^ are the temperatures of respective opaque boundaries of the system gas - capillary-porous body, K.
In the systems of equations (14) and (15), I+(0) and I-(8S) are related to the non-transparent boundaries of the system gas - capillary-porous body, and Iv (S±ir) - to the
1 -Rair:
Hence
- = (1 -Rf,v) n
(21)
e f,^^. (22)
, n2,v
That is, provided n2,v>1, inequality 8f,v<8air,v holds.
In order to obtain final equations for deter-
(14)
mining the density of heat flux of radiation and its divergence, it is necessary to solve the appropriate system of equations. This system includes (14), (15), a formal solution (9) for the boundary surfaces of the system and the current values of coordinate z [17].
(15) Thus, the system of equations (4)-(22) is a complete mathematical formulation of the selective problem on radiation-conductive heat exchange that describes the drying process of fluting by the
infrared radiation. In this case, expressions for qr and dqL are
dz
accurate analytical formulas for determining the density of radiation flow and its divergence in the one-dimensional case.
The range of waves from the industrial emitters is close to the range of the largest absorbing capacity of water.
That is why it is possible to simplify a mathematical model (4)-(22) by using a transition from a selective to the "grey" model of environment.
For the transition from a model of selective environment to a simpler "gray" environment, it is necessary to perform the following transformations. All radiation properties of phases and boundaries of the system are averaged, for example, by Roseland [17, 19]. The Planck function in (9), (13)-(24) is
replaced with expression —. The value of a in the expression n
is the Stefan-Boltzmann constant, W/(m2-K4).
Based on the scientific literature review, we found that in order to model the radiation drying of materials with physical properties that are close to fluting, it is possible to apply a simplified model based on the Burgers' equation [12, 20]. In this case, the results obtained describe the drying process with accurately sufficient for the engineering calculations.
Considering data [12, 20], it is possible to perform further simplification of the mathematical model (4)-(22) based on the Burgers' equation [4]. In this case, heat transfer from an IR emitter to the surface of fluting can be described by the Stefan-Boltzmann law, disregarding the layer of wet air above it. It is, however, necessary to employ the approximation of optically thin layer with regard to absorption [19].
The medium, which is a wet fluting, is assumed to be translucent "grey" and absorbing; the boundaries of its layer and the surface of an IR heater - diffuse. It is possible in this case to confine ourselves to considering only the layer of fluting's material, but under more complex conditions for external heat exchange than for the two-phase selective medium. It is also required to complement the model with a kinetic equation that describes a change in moisture content over time. At this approximation, mathematical statement of the problem, considering the kinetics of drying [4], includes the following system of equations:
v- <T)r- <t)|t=i.. <t>S-^-p. (T)
qr (z)=qr (Tz=0 )e0 ; qr (Tz=0 ) = epr (Ka-v )o(t,4r-Tz=0 );
dt = b [ps (Tz=0 )-pp (Tz=0 )]+ k (u - um ),
(23)
the drying speed, s-1 (kg/(kg^s));
dt
e pr (Ka-v ) = -
1
1_ -1+ 3KaA
is the resulting degree of blackness of surfaces of the IR heater of the layer of fluting, which takes into account the absorption of layer of wet air [16]; sIR is the degree of blackness of radiating surface of an IR heater; sz=0 is the degree of blackness of the surface of fluting; FIR, Ff is the area of surfaces of IR-heater and fluting, respectively, m2; Ka-v is the "grey" value of absorption coefficient of wet air, m-1; Sair is the thickness of layer of wet air, m; a is the Stefan-Boltzmann constant, W/(m2-K4); TIR, Tz=0 is the absolute temperature of surfaces of IR heater and fluting, respectively, m; p is the moisture exchange coefficient, s/m; k is the drying constant, s-1; g is the mass of square meter of dry fluting, kg/m2; ps, pp are the pressure of vapor saturation at the surface of a layer of fluting and the partial pressure of water vapor in the adjacent layer of environment (wet air), respectively, Pa; u, um are the mean relative moisture content in the current moment and relative equilibrium in the layer of fluting, respectively.
A formula for determining the divergence of radiation flow can be received by differentiating, by parameter z, the integrand expression, which is an indicator of the power of exponent in the second equation of system (24).
9qr (r ) d
qr (Tz=
-|K(T(z))dz
9z 9z
= Ik (T (z)) = const = -qr (Tz=0 )Ke-Kz.
(24)
The equation of thermal conductivity in the system of equations (23) can also be written through a thermal diffu-sivity coefficient. Then, with regard to the latter and (24), the system of equations (23) takes the final form:
dT=a (T)dT + qr (Tz=0)
9x= af (i)9z2 + cpe. (T )peff (T) qr (Tz=0 ) = epr (Ka-v )o(t,rr-Tz4=0 );
^ [ps (Tz=0 )-pp (Tz=0 )]+ k (u
Ke-Kz -
Pf (T)r
du
Pef- (T)cpe- (T) dt'
(25)
where
l.
where cp e--(T) is the temperature dependence of mass isobar heat capacity of wet fluting, J/(kg-K); p e--(T) is the temperature dependence of effective density of wet fluting, kg/m3; p .(T) is the temperature dependence of density of absolutely dry fluting, kg/m3; T is the absolute temperature, K; t is the time, s; X eff(T) is the temperature dependence of effective coefficient of thermal conductivity of wet fluting, W/(m-K); z is the coordinate, m; qr(z) is the density vector of radiation heat flux, W/m2; r is the mass heat of water vaporization, J/kg; K(T(z)) is the temperature dependence of "gray" value
of effective coefficient of absorption of wet fluting, m-1; ^ is
cp eff peff
is the effective thermal diffusivity of fluting, m2/s.
For a wet fluting, effective "gray" value of absorption coefficient is expressed by the formula according to the known law of additivity:
K = Kf (1 - w) + Kww,
(26)
where Kf, Kw are the absorption coefficients of dry fluting and water, m-1; w is the volumetric share of moisture in fluting.
The rest of effective properties, which are included in the system of equations (25), are received similarly using the additivity property.
The system of equations (25) is a nonlinear one by the thermal-physical properties and source terms. Thus, for example, qr depends on the temperature at the surface of fluting Tz=0 in fourth power.
a=
Initial conditions for (25), t = 0 : T = T(z), 0 <z <Sf,
Tir = T (W),
n-(-1eff (T)VT)| z=0 = at (T - T„ ); n (T)VT)| =a2 (T - Ts2 ),
and 8uk+1, i = 1,n, respectively, and the desired temperature and relative humidity in the nodes of computational grid are determined by formulas:
where s. is the thickness of a layer of wet fluting, m. Boundary conditions for (25), t>0: - on the opaque boundary of an IR emitter, we consider boundary conditions of the first kind:
T,k+1 = Tk + STk+\ i = 1, n; uk+1 = uk + Suk+1, i = íñ
(31)
(28)
where W is the dependence of change in the electric power of a radiation heater on the time, W. At this boundary, we also assign a «gray» value of the degree of blackness of the working surface of IR emitter sIR;
- at the upper and lower boundaries of a layer of fluting, we assign boundary conditions of the third kind:
(29)
(30)
where k is the number of iteration at each step of the integration by time.
In this case, the system of discrete equations for relative humidity is solved first in the iteration loop of solving the problem (25)-(30). Thus, a system of equations of thermal conductivity and kinetics is resolved by the method of successive approximations. In order to improve convergence in the solution of a system of equations, it is possible to employ the method of top relaxation.
A criterion for obtaining the numerical solution of the initial problem at each integration step by time is the fulfillment of condition:
|STk| <eT, i = 1,n; ISukl <eu, i = ü,
(32)
where n is the vector of external normal to the surfaces
of sides of wet fluting; V = — is the Hamiltonian in the
dz
one-dimensional case, m-1; aj, a2 are the coefficients of heat transfer from the upper and lower sides of a layer of fluting, respectively, W/(m2-K); TS1, Ts2 are the ambient temperatures, above a layer of fluting from its upper and lower sides, respectively, K. At the upper boundary of a layer of fluting, we also assign the degree of blackness of the surface of fluting Sz=0.
Mathematical statement of the problem on complex heat transfer taking into account the kinetics of the process (25)-(30) is a complete mathematical formulation of the nonlinear one-dimensional problem on drying the wet fluting by the radiation method.
6. Technique for the numerical solution of the problem on drying by the radiation method
In order to numerically solve the stated problem, described by equations (25)-(30), we employed the finite difference method [23].
Uniform grid is used for the discretization of computational domain
z = h(i-1), i = Ü, h = -f,
where n is the number of grid nodes. To obtain finite-difference analogs of the original differential equations (25), we use an absolutely stable implicit three-point scheme with the first order approximation by time and the second -by coordinate. Since the systems of obtained discrete equations for determining the temperature and relative moisture content are non-linear, then, in order to solve them, we apply the linearization by Newton's method by temperature.
Systems of linearized equations at each iteration step are resolved successively by economical method of three-point sweep relative to discrepancy vectors 8Tk+1, i = 1,n
where sT, su are the errors, assigned in advance, in determining the temperature (K) and relative humidity, respectively.
Software realization of the numerical technique is implemented in the Mathcad programming language [24]. A description of the algorithm in the process of its development and application is given in [25]. By using the formulated mathematical model, whose numerical problem is described by equations (25)-(30), it is also possible to develop a software code in other programming languages.
7. Verification of numerical model for the process of radiation drying of fluting
A comparative analysis of results of the numerical modeling with data on physical experiments to the study of kinetics of drying of fluting by the IR radiation was conducted based on such parameters as the duration and speed of the drying process (at attaining the same values of moisture content) and the fluting surface temperature. The values of fluting drying speed were determined by the results of numerical and physical experiments (25).
To compare the results, we used the samples of fluting of mark B-1 with a square meter weight of 0.112, 0.125 and 0.140 kg/m2, as well as the samples from non-standard fluting with increased specific weight of 0.2 kg/m2.
As the source of radiation in the experiments, we used ceramic electric IR emitter of type ECH4 with nominal capacity of 1 kW and the temperature of working surface about 833 K [22]. The length of electromagnetic waves generated by the IR emitter is 3.3...3.7 ^m. Results of comparison are given in Table 1.
Designations in Table 1: g - square meter weight of fluting, kg/m2; q - density of heat flux at the surface of fluting, W/m2; if - drying time at the end of the first period, s; xE - drying time to resulting equilibrium moisture content, s; T1, T2 -absolute temperatures of the surface of fluting in the first period of drying and at the end of the second period
(resulting temperature), K;--speed of drying in the first
9t 1
period of drying, s-1 (kg/(kg-s)).
Table 1
Comparison of experimental data with the data of numerical analysis
g, kg/m2 q, W/m2 T1, s* T2, s* T1, K* T2, K* — , s-1* 3t 1 '
0.112 6835 56/58 189/187 363/367 496/495 0.0192/0.0190
0.125 6835 61/64 203/201 365/368 496/496 0.0178/0.0175
0.125 4300 116/117 244/246 342/345 460/462 0.0120/0.0119
0.125 3665 144/145 318/317 339/343 434/432 0.0090/0.0089
0.140 6835 65/67 212/214 367/370 497/499 0.0166/0.0163
0.200 6835 90/92 252/253 371/384 504/505 0.0112/0.0099
Note: * - data of numerical analysis are given through a fraction
As a result of comparing the data of numerical modeling with experimental data, we established that the discrepancy between them does not exceed 5 % on the interval of change in the weight of square meter of dry fluting from 0.112 kg/m2 to 0.2 kg/m2 and density of heat flux at the surface of fluting from 3665 W/m2 to 6835 W/m2.
8. Discussion of results of numerical modeling of physical fields in the process of drying of fluting
Numerical modeling of the drying process of fluting by the radiation method allows us to determine kinetic patterns and basic parameters required to intensify the process and design of drying equipment.
The numerical model developed makes it possible to determine duration of the process of drying the fluting to the degree of its dehydration assigned in advance. In this case, we also define the current degree of fluting dehydration, maximum value of fluting temperature and its thermal-physical parameters considering the moisture content. In addition, by using a numerical model, it is possible to calculate the capacity of IR heaters and their quantity. This makes it possible to design new drying equipment of the radiation type.
However, the mathematical and numeric models presented do not take into account the phenomenon of fluting shrinkage during drying that affects its thermal-physical properties. A change in the structure of fluting during the crystallization of glue that cements it is also not considered. But the results of verification of the numerical model indicate that the mentioned shortcomings do not essentially influence its accuracy.
The stated mathematical model and the developed numerical model of the radiation drying of capillary-porous materials of the fluting type can also be used for modeling the drying process of other types of paper, for example sanitary-hygienic [21].
Under high-temperature drying mode of fluting, deterioration in its physical-mechanical properties is possible. That is why we plan to complement the mathematical model developed by the kinetic patterns in the process of fluting destruction.
9. Conclusions
1. Based on the analysis of physical processes that occur during drying of fluting by the IR radiation, we formulated a physical model of the process. The physical model takes into account the absorbing and emissive capacity of the environment translucent for the IR radiation, diffuseness in the properties of medium boundaries, temperature dependence of physical properties and the kinetics of drying process.
2. We formulated a mathematical model and developed a numerical model of the process of radiation drying of fluting, which considers the kinetics of drying process, the translucency to the IR radiation of the material using the approximations of the Burgers' models, "gray" medium and diffuse reflection of boundaries. Software realization of the numerical method is implemented in the MathCAD programming language. The formulated models of drying process will also make it possible to develop a software code in other programming languages.
3. Verification of the numerical model of the radiation drying of fluting revealed a discrepancy between the results of numerical modeling and the data of physical experiment within the range of 5 %. A comparison of the results was conducted by the following parameters: surface temperature of fluting, duration and speed of fluting drying. The weight of square meter of dry fluting varied between 0.112 kg/m2 and 0.2 kg/m2. Heat flux density at the surface of fluting changed from 3665 W/m2 to 6835 W/m2.
References
1. Kolchyna, Y. A. Rynok kartona v Ukrayne (sostoyanye y problemy) [Text] / Y. A. Kolchyna // Upakovka. - 2013. - Issue 2. -P. 22-26.
2. Lu, T. Numerical and experimental investigation of paper drying: Heat and mass transfer with phase change in porous media [Text] / T. Lu, S. Q. Shen // Applied Thermal Engineering. - 2007. - Vol. 27, Issue 8-9. - P. 1248-1258. doi: 10.1016/j.applther-maleng.2006.11.005
3. Marchevskiy, V. M. Infrachervone nagrivannya i sushinnya paperu y kartonu [Text] / V. M. Marchevskiy, O. A. Novohat // Himich-na inzheneriya, ekologiya ta resursozberezhennya. - 2011. - Issue 2 (8). - P. 42-44.
4. Lykov, A. V. Teoriya sushki [Text] / A. V. Lykov. - Moscow: Energiya, 1968. - 472 p.
5. Ghodbanan, S. Steady-State Modeling of Multi-Cylinder Dryers in a Corrugating Paper Machine [Text] / S. Ghodbanan, R. Aliza-deh, S. Shafiei // Drying Technology. - 2015. - Vol. 33, Issue 12. - P. 1474-1490. doi: 10.1080/07373937.2015.1020161
6. Heo, C. H. Dynamic modeling of paper drying processes [Text] / C. H. Heo, H. Cho, Y.-K. Yeo // Korean Journal of Chemical Engineering. - 2011. - Vol. 28, Issue 8. - P. 1651-1657. doi: 10.1007/s11814-011-0046-0
7. Ottosson, A. A mathematical model of heat and mass transfer in Yankee drying of tissue [Text] / A. Ottosson, L. Nilsson, J. Berg-hel // Drying Technology. - 2016. - Vol. 35, Issue 3. - P. 323-334. doi: 10.1080/07373937.2016.1170697
8. Weineisen, H. Modeling Drying and Energy Performance of Industrial Through-Dryers [Text] / H. Weineisen, S. Stenstrom // Drying Technology. - 2008. - Vol. 26, Issue 6. - P. 776-785. doi: 10.1080/07373930802046443
9. Dhib, R. Infrared Drying: From Process Modeling to Advanced Process Control [Text] / R. Dhib // Drying Technology. - 2007. -Vol. 25, Issue 1. - P. 97-105. doi: 10.1080/07373930601160908
10. Khansary, M. A. Modeling drying of a coated paper [Text] / M. A. Khansary, F. K. Q. Joogh, A. Hosseini, J. Safari, E. Allahyari, N. S. Zadeh, A. H. Sani // International Journal of Modeling, Simulation, and Scientific Computing. - 2014. - Vol. 05, Issue 01. -P. 1350019. doi: 10.1142/s1793962313500190
11. Kocabiyik, H. Drying of carrot slices using infrared radiation [Text] / H. Kocabiyik, D. Tezer // International Journal of Food Science & Technology. - 2009. - Vol. 44, Issue 5. - P. 953-959. doi: 10.1111/j.1365-2621.2008.01767.x
12. Doymaz, I. Infrared drying of sweet potato (Ipomoea batatas L.) slices [Text] / I. Doymaz // Journal of Food Science and Technology. - 2011. - Vol. 49, Issue 6. - P. 760-766. doi: 10.1007/s13197-010-0217-8
13. Doymaz, I. Drying of pomegranate seeds using infrared radiation [Text] / I. Doymaz // Food Science and Biotechnology. - 2012. -Vol. 21, Issue 5. - P. 1269-1275. doi: 10.1007/s10068-012-0167-1
14. Marchevsky, V. Kinetics of Corrugated Board Flute Drying with the Use of Infrared Radiation [Text] / V. Marchevsky // The Advanced Science Journal. - 2015. - Vol. 2015, Issue 6. - P. 69-72. doi: 10.15550/asj.2015.06.069
15. Kolesnikov, A. V. The complex heat exchange model at growing of large alkali halide crystals [Text] / A. V. Kolesnikov, V. I. Deshko, Yu. V. Lokhmanets, A. Ya. Karvatskii, I. K. Kirichenko // Functional Materials. - 2010. - Vol. 17, Issue 4. - P. 483-487.
16. Zigel', R. Teploobmen izlucheniem [Text] / R. Zigel', Dzh. Hauell; B. A. Hrustalev (Ed). - Moscow: Mir, 1975. - 934 p.
17. Deshko, V. Y. Modelyrovanye nestatsyonarnoho protsessa vytyahyvanyya krystallov yz rasplava [Text] / V. Y. Deshko, A. Ya. Kar-vatskyy, Yu. V. Lokhmanets, A. O. Hulenko // Matematychne modelyuvannya. - 2011. - Issue 2 (25). - P. 75-79.
18. Sergeev, O. A. Teplofizicheskie svoystva poluprozrachnyih materialov [Text] / O. A. Sergeev, A. A. Men. - Moscow: Izd-vo stand-artov, 1977. - 288 p.
19. Ocisik, M. N. Slozhnyj teploobmen [Text] / M. N. Ocisik; N. A. Anfimov (Ed.). - Moscow: Mir, 1976. - 616 p.
20. Seyed-Yagoobi, J. Heating/drying of uncoated paper with gas-fired and electric infrared emitters - fundamental understanding [Text] / J. Seyed-Yagoobi, H. Noboa // Drying 2004 - Proceedings of the 14th International Drying Symposium (IDS 2004). -Sao Paulo, 2004.
21. Marchevskiy, V. M. Sushinnya sanitarno-gigienichnogo paperu z vikoristannyam infrachervonogo viprominyuvannya [Text] / V. M. Marchevskiy, O. A. Novohat, L. G. Voronin, O. O. Tatarchuk // Himichna inzheneriya, ekologiya ta resursozberezhennya. -2015. - Issue 1 (14). - P. 29-31.
22. Marchevsky, V. Paper drying process for corrugation (fluting) using radiant energy [Text] / V. Marchevsky, O. Novokhat, O. Tsep-kalo // Ukrainian Journal of Food Science. - 2015. - Issue 2. - P. 310-321.
23. Deshko, V. I. Control of radiation-conductive heat exchange at crystal growth from melt [Text] / V. I. Deshko, A. Ya. Karvatskii, A. V. Lenkin, Yu. V. Lokhmanets // Functional Materials. - 2008. - Vol. 15, Issue 2. - P. 229-234.
24. Mathcad. Engineering math software that allows perform, analyze, and share your most vital calculations [Electronic resource]. -Available at: http://www.ptc.com/engineering-math-software/mathcad/
25. Karvatskyi, A. Ya. Rozv'iazannia neliniinykh nestatsionarnykh zadach teploprovidnosti z vykorystanniam CAD-system [Text] / A. Ya. Karvatskyi, A. Yu. Pedchenko // Matematychne ta komp'iuterne modeliuvannia. Seria: Fizyko-matematychni. - 2016. -Issue 13. - P. 67-77.