Computer Modelling of Nonisothermal Transfer of Moisture during Drying of Wood with the Use of the Computing Environment Mathcad

Abstract

This study explores the features of process of drying internal heat and mass transfer in wood in the conditions of essentially nonisothermal transfer of moisture. The process of capillary conductivity of wood is analyzed on the model of colloidal capillary-porous solid. The density of moisture flows caused by the phenomenon of thermal and hydraulic conductivity is determined. Previously experimentally established fact of influence on the process of moisture removal of the current moisture content of wood is analytically confirmed. Dependences of the phase transition criterion taking into account the direction of the temperature gradient and wood moisture are analytically obtained. Recommendations are given to determine the thermogradient coefficient taking into account the asymmetry of moisture flows depending on the direction of the temperature gradient, and also for the correction in the entry form of the equation of moisture transfer with the influence of the effect of thermal and hydraulic conductivity.

Introduction

To describe the transfer of heat and mass of matter is usually used a generalized system of differential equations of transfer [ 1 ].

For the simplest case when the total pressure gradients are absent (heat and mass transfer processes in lowtemperature convective drying) A. V. Lykov [ 2 ] and M. S. Smirnov [ 3 ] give the following system of differential equations: For an unbounded plate, the initial and boundary conditions of the III kind have the form: ( ( ( ) ( ) )] ) (5) ( ) ( ) ( ) Symmetry condition: In criterion form of record [ 1 ] in (1) and (2) there are three dimensionless criterions: Lu (Lykov’), Ko (Kossovich’) and Pn (Posnov’), which, for example, N.I. Gamayunov [ 4 ] calls the criteria of non-isothermal transfer. On closer inspection (1) and (2), it is obvious that in non-isothermal moisture exchange, the latter is significantly influenced not only by the thermogradient coefficient δ, but also by the phase transformation criterion ε. It is especially necessary to take into account when calculating the drying parameters based on using the phenomenon of thermal and hydraulic conductivity [ 1, 5 ].

In non-isothermal conditions, in the presence of heat and mass transfer potentials, heat and mass transfer in gas mixtures is caused by thermodynamic driving forces. In gas blends moisture transfer occurs under the influence of temperature gradient (thermal diffusion of moisture) [ 1 ]. In this case, there is an additional flow of moisture [ 2 ]: δ - the coefficient of thermal diffusion, 0К-1 (sometimes called the coefficient of Soré);

t – the temperature gradient, 0К.

The presence of a temperature gradient causes a diffusion transfer of vapor under the influence of the gradient of its partial pressure in wet solids in addition to the Soré effect. In addition, the movement of the liquid in the capillary-porous solid in the direction of the heat flow can be caused by the presence of "pinched" air. When the temperature rises, the pressure of the "pinched" air increases and the air bubbles expand. As a result, the liquid in the capillary pore moves in the direction of the heat flow. Therefore, at P = const the total moisture flow is equal to: ( ) (9) The first member of equation (9) shows the isothermal mass transfer, the second is the process of thermal and hydraulic conductivity.

From the theory of drying it is known that such a phenomenon as thermal and hydraulic conductivity [ 5 ] can create an additional flow of moisture, if the temperature gradient is negative value, and thus accelerate the process of moisture removal (drying). But if the temperature gradient is positive value, thermal and hydraulic conductivity can significantly complicate moisture removal. Thus it is absolutely obvious that such inhibiting effect can extend only to process of removal of a liquid moisture (up to its full stop). At the same time, the molecular and molar transport of vaporous moisture inside the capillary-porous solid is not extended due to the different physical nature of these phenomena [ 6 – 8 ].

The inhibiting effect of thermal and hydraulic conductivity (non-isothermal flow) can reduce and even stop the liquid flow that has arisen due to the sufficient amount of the jamming pressure, but it can not make it negative. This is due to the following objective phenomena [ 9, 10 ]. The motion of the wetting fluid in a single through cylindrical capillary under the action of surface tension forces under laminar regime is determined by the equation [ 10 ].

( ) In the linear approximation, ignoring in equation (10) the first two terms due to their smallness, we get for the horizontal capillary (φ = 0) for a vertical capillary (φ = 900)

- the length of the liquid column in the capillary; t – the driving time;

– the liquid density; - the radius of the capillary; φ – the angle of the capillary;

– the viscosity; ζ – the surface tension. Equations of fluid motion (10 – 12) do not take into account the effects of clamped air [ 5 ]. The appearance of air bubbles when fluid flow in wood a significant role is played by the Jamen’ effect (figure 1). It consists in the fact that with a slow filling of capillaries, the formed closed air inclusions cause a sharp increase in the hydrodynamic resistance of the medium. Such a picture is observed at movement of the aerated liquid, for example oil, through the porous medium [ 11 ]. Classical works of academician E. F. Votchal [ 12 ] were devoted to the study of the effect of Jamen on the movement of fluid through the vessels of living wood, unfortunately not received in the future the necessary development. In essence, The effect of Jamen is an integral expression of the influence of internal boundary conditions on the motion of a liquid through a complex capillary-porous medium.

We can describe the filtration of a liquid as it moves through a porous medium by Darcy's law: ̅ or when it reaches the capillary Poiseuille's formula: с – the permeability coefficient; – the pressure gradient. (13) (14) At a sufficiently stable temperature, the viscosity coefficient can not change its value depending on the structure of the medium. The change in the filtration rate of the liquid ̅ can be attributed to the appearance of a foreign phase – gas bubbles, causing a change in the internal boundary conditions of the medium. This is most clearly seen from the law of Darcy: at constant and

[ 11 ] there is a decrease in u by reducing the effective value of the permeability coefficient.

Poiseuille's formula in its standard form does not apply to the case of non-continuous fluid flow and the discontinuity can be represented as [ 9 ] ̅ , (15) where α – the correction factor depending on the number of gas bubbles, the coefficient of surface tension in relation to the material of the capillary walls.

It should be noted that formally reducing ̅ can also be attributed by increasing the effective value of the viscosity coefficient of the liquid, due to its "carbonation", considering the permeability coefficient with the same value. However, according to N.I. Osnach [ 9, 15 ] reduction of filtration rate in the presence rosaries of Jamen in wood capillaries occurs due to the joint effect of both reducing conductivity and increasing the effective viscosity value, but with a significant predominance of the first factor. The gas bubble formed for one reason or another in the capillary will have a significant resistance to the movement of the liquid. The chain of bubbles - rosaries of Jamen - can completely close the capillary.

Research methods

During the study, a specially planned computational experiment was implemented. At the same time, the simulation of wood heat and mass transfer processes was carried out in the Mathcad-14 computing environment using special software [ 13 ]. Analysis of moisture transfer processes was carried out using a physical model of colloidalcapillary-porous solid for coniferous wood [ 14 ]. In addition, the appropriate mathematical apparatus was used to analyze the integral-differential equations.

Research results

Analysis of the effect of Jamen in wood capillaries showed the following.

For example, the capillary formed a chain of rosaries of Jamen containing n bubbles. Consider it i-th gas bubble. In dynamic mode, the left and right meniscus will have different curvature (figure 1). We determine the value of the additional resistance generated at the boundaries of the i-th liquid drop.

The pressure Pi, which in the literature is usually called Laplace’ at each point of the curved surface will be directed towards the center of the curvature of the corresponding elementary platform .

The total force of Laplace pressure on the surface of the meniscus А1А1 is where – accordingly, the radii of the smallest and largest curvature of the surface element.. For the surface of the meniscus А2А2 power Laplace’ pressure equal Thus, the backpressure force associated with the i-th drop in the promotion of the latter To calculate (18), it is necessary to know the subintegral expression as a function of independent variables, which will be very difficult to change during the deformation of the meniscus when it moves along with the liquid along the capillary. As a result, the expression (18) becomes virtually non-numbered.

However, it is possible to estimate the amount of backpressure. Однако оценить величину противодавления возможно.

Let's assume in the first approximation that each bubble introduces some additional resistance where - the value of the mean radius of curvature of the meniscus drops.

When the whole chain starts to move, the total capillary back pressure will be .

For wood capillaries (model colloidal capillary-porous solids [ 16 ], the equivalent radius of early tracheids pine) characteristic size the surface tension of the liquid at Т = 353 0К: Then the back-pressure produced by one bubble .

Researchers [ 9, 10 ] cautiously suggest that the presence of 50 ... 100 air bubbles in the capillary makes it impervious.

However, we can argue that the presence of already one bubble makes the capillary impervious (when the liquid moves inside), since taking into account the area of the capillary, the backpressure is Taking into account the impossibility of moisture movement through the capillaries into the wood under the influence of a temperature gradient, the following results were obtained in the computational experiment. Table 1 shows the calculated total density of moisture flow in different directions of the temperature gradient (figure 2).

∬ ( ∬ ( ∬ {( An analysis of the data in table 1 and figure 2 shows that: 1. The total flow of moisture also depends significantly on the moisture content of the wood. 2. The most significant thermal and hydraulic conductivity affects the removal of free moisture. If the moisture content of the wood is less than 20%, its influence becomes insignificant. It is necessary to consider it at construction of conditions of drying of wood.

A very important parameter in the calculation of heat and mass transfer processes is the criterion of phase transition ε. Analysis of the system of differential equations of heat and mass transfer (1-2) shows that the special importance of knowledge of the exact magnitude ε becomes in solving the problems of non-isothermal transfer, as in (1) ε largely determines the temperature of the wood. Moreover, the decrease ε, other things being equal, increases the temperature of the wood.

Table 2 shows the calculated values of the phase transition criterion. Also of note are the following: 1. The phenomenon of thermal and hydraulic conductivity significantly affects the value of the phase transition criterion ε. This is due to the fact that the change in the direction of the temperature gradient leads to a significant decrease (or increase) in the flow of liquid moisture, which changes the vapor-liquid ratio, and hence the value ε. 2. The free moisture in the positive direction of the temperature gradient is removed in the form of steam (ε = 1,0), that is, the inhibitory effect of the positive temperature gradient is so high that it overlaps the flow of moisture caused by the wedging pressure. The value of ε with a wood moisture of, %

10 0,616 0,487 0,544

20 0,728 0,555 0,630

30 0,805 0,614 0,677 > 30 1,0 0,31 0,487 The obtained information about the magnitude ε will be used in the future. When solving differential equations of heat and mass transfer of the form (1-2), it is necessary to present the value of the phase transformation criterion as where ( ) ( ) – coefficients depending on the direction of the temperature gradient;

( ) ( ) – the value of the phase transition criterion (table 2).

Figure 3 shows the values of the coefficients ( ) ( ).

( ) ( ) ( ) ( ) (20)

( ) ( ) ( ) ( ) The actual non-isothermal transfer is calculated by solving the second equation of the system of differential heat and mass transfer equations (2). The following circumstances should be taken into account (table 3).

When solving the system, there is some asymmetry of the moisture flows when the direction of the temperature gradient vector changes.

Equation (2) takes the form of (25) where – the coefficient taking into account the asymmetry of moisture flows.

( ) ( ) ( ) ( ) ( ) ( ) (26) (27) (28)

Conclusion

Strictly mathematically confirmed the hypothesis about the impossibility of the movement of moisture through the capillaries into the wood under the action of the temperature gradient.

The refined values of the phase transition criterion ε and the temperature coefficient δ, which can be used in modeling the processes of non-isothermal moisture transfer in the process of drying wood by solving a system of differential equations of heat and mass transfer, are obtained. At the same time, it is necessary to take into account some asymmetry of moisture flows when changing the direction of the temperature gradient vector.

References

1. A.V. Lykov . Theory of heat and mass transfer . Gosenergoizdat , Moscow, 1963 , 535 p.

2. A.V. Lykov . On systems of differential equations of heat and mass transfer in capillary-porous solids . Journal of physics and engineering , XX VI, 1 , 1974 . Pages 18 - 25 .

3. M.S. Smirnov . On the system of differential equations of the drying process . Journal of physics and engineering , IV, № 9 , 1961 . Pages 40 - 44 .

4. H.I. Gamayunov. About criteria for non-isothermal transfer . Journal of engineering physics , III, 12 , 1960 . Pages 58 - 62 .

5. N.R. Galyautdinov , S.R. Mukhametzyanov , I.F. Khakimzyanov . Investigation of the processes of vacuum drying of lumber in oscillating mode , Bulletin of Kazan technological University, 18 , 8 , 2015 . Pages 185 - 187 .

6. R.R. Safin , R.R. Khasanshin , I.F. Khakimzyanov , S.R. Mukhametzyanov , P.A. Kainov . Increasing the Energy Efficiency of the Process of Oscillating Vacuum-Conductive Drying of Wood by Means of a Heat Pump // Journal of Engineering Physics and Thermophysics , Vol. 90 , No. 2, March , 2017 . Pages 310 - 317 .

7. R.R. Safin , R.R. Khasanshin , I.F. Khakimzyanov , S.R. Mukhametzyanov . Non-volatile Facility for Vacuum Drying of Thermolabile Materials // Procedia Engineering . Volume 206 , 2017 , Pages 1063 - 1068 .

8. R.R. Safin , R.R. Khasanshin , R.V. Salimgaraeva , E.A. Beliakova . Heat treatment of crushed wood in rotary drum dryers // Key Engineering Materials. 2017 , Volume743, Pages 378 - 382 .

9. N.A. Osnach . The permeability and the conductivity of the wood . Forest industry , Moscow, 1964 . 182 p.

10. V.I. Patyakin , Y.G. Tishin , S.M. Bazarov . Technical hydrodynamics of wood. Forest industry , Moscow, 1990 , 304 p.

11. F.A. Trebin . Neftprominvest sand reservoirs . Gostoptekhizdat , Moscow, 1945 , 139 p.

12. E.F. Votchal . On the movement of pasoki in plants . Kushnerev & Co., Moscow , 1897 , 390 p.

13. E.E. Shishkina. Energy-Saving Technology of Convective Lumber Drying Based on Controlled Moisture Transfer in Wood: Dr. Eng. Sci. Diss . Yekaterinburg, 2016 . 336 р.

14. A.G. Gorokhovsky , E.E. Shishkina . Models of moisture transfer in colloidal capillary-porous structure of wood. Woodworking: technologies, equipment, management of the XXI century . Proceedings of the II international Eurasian Symposium , Ekaterinburg, 2008 . Pages 41- 46 .

15. N.A. Osnach . About the permeability of wood // Woodworking Industry, 1961 , № 3 , Pages 11 - 13 .

16. A.G. Gorokhovsky . Technology of Lumber Drying Based on Modeling and Optimization of Heat and Mass Transfer Processes in Wood: Dr. Eng. Sci. Diss . Yekaterinburg, 2008 , 263 p.