SYSTEM MODELLING OF NON-STATIONARY DRYING PROCESSES
Abstract and keywords
Abstract (English):
The kinetics of the drying process in continuous drum dryers differs from the drying of single objects in a batch mode. Drying process is affected by too many factors; hence, it is practically impossible to obtain an analyt- ical solution from the initial equations of heat and mass transfer, since the duration of drying depends on the opera- ting parameters. Therefore, it is of high theoretical and practical importance to create a highly efficient rotary drum dryer. Its design should be based on an integrated research of non-stationary processes of heat and mass transfer, hydrodynamics of fluidized beds, and drying kinetics in the convective heat supply. The experiment described in the present paper featured sunflower seeds. It was based on a systematic approach to modelling rotary convective drying processes. The approach allowed the authors to link together separate idealized models. Each model characterized a process of heat and mass transfer in a fluidized bed of wet solids that moved on a cylindrical surface. The experiment provided the following theoretical results: 1) a multimodel system for the continuous drying process of bulky mate- rials in a fluidized bed; 2) an effective coefficient of continuous drying, based on the mechanics of the fluidized bed and its continuous dehydration. The multimodel system makes it possible to optimize the drying process according to its material, heat-exchanger, and technological parameters, as well as to the technical and economic characteristics of the dryer.

Keywords:
System modelling, continuous drying, heat and mass transfer, drum dryer, fluidized bed of wet solids
Text
Publication text (PDF): Read Download

The general theory of heat and mass transfer in capil- lary porous and other dispersed media belonged to Prof. Lykov. It was based on the thermodynamics of irrever- sible processes and the theory of generalization of vari- ables. According to Prof. Lykov, sets of equations are to be solved as a single complex process [1].

The theory of heat and mass transfer is based on solving sets of linear equations with boundary condi- tions, which corresponded to constant and variable po- tentials in a medium that varied according to established laws. Thus, it was intended for stationary material and medium [2].

To intensify and improve technological processes, one needs reliable, physically valid modelling metho- ds [3]. This is especially important for energy-inten- sive drying processes of wet solids in a fluidized bed that is moving along a cylindrical surface, e.g. drum dryers [4].

 

In this paper, modelling means a physical  ana- lysis of heat and mass transfer, as well as the hydrody- namics of the processes that occur a rotary drum dryer, their mathematical description,  and  possible  solutions by analytical or numerical methods involving various software [5]. The analytical and numerical methods were based on preliminary data on the kinetics of drying and heating of individual particles. Such information was obtained either from the available model representations or from experimental data. In most cases, preference was given to direct experimental data, which took into account possible effects of anisotropy of heat and mass transfer properties and irregular geometric shape of par- ticles of real solids.

An adequate description of the continuous drying technology of wet solids in a fluidized bed requires an integrated approach to the problem. Such an approach requires a system analysis of hydrodynamic, diffu- sion, and thermal processes complicated by an overlap

 

 

Copyright © 2019, Antipov et al. This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material for any purpose, even commercially, provided the original work is properly cited and states its license.

 

 

 

of various phenomena. A complete theoretical picture of a continuous drying process should be based on a mathematical model that would link a set of typical structures, or idealized models, each of which reflects a particular type of transfer or transformation. The opti- mal way to develop a new drying technique is to com- bine the multimodel system of drying processes with experimental studies on the kinetics of moisture remo- val in a fluidized bed [6].

The research objective was to develop a multimo- del system for the continuous process of drying and heat transfer in a fluidized bed of wet solids.

 

STUDY OBJECTS AND METHODS

The study was based on a complex of general and specific scientific methods. The general scientific me- thods involved analysis and synthesis, testing a theory with practice, interpretation of the results obtained, etc. The specific scientific methods included the abstract logical method, the method of modelling, the empiri- cal method, the method of statistical  probability,  etc. The theoretical and methodological foundation of the research included studies conducted by Russian and fo-

 

authentic trial equipment were used to conduct the ex- periments and test the physical and mathematical models of the drying and steam treatment processes.

S

 

F

T

FS

ST

The systematic approach in modelling the convec- tive drying processes. The process of continuous dry- ing is an original research subject that can be marked as a certain system S, while the model of the drying pro- cess M represents a different system. According to the systematic approach principle, the drying process inte- racts with the external environment E. Depending on the research objective, the study may feature different ratios between the subject and the external environment. These relations represent a model of relations of the external environment with the subject. Since our ideas of the su- bject and the external environment are also models, the following models can be proposed: M is the model of the subject S, M is the model of external environment at the inlet, M is the model of external environment at the output, M is the model of the connections between the external environment and the subject at the inlet, and M is the model of the connections between the external environment and the subject at the output.

The combination of all these models

 

reign experts in the field of drying, such as Ginzburg,

 

M   = (M , M

 

, M , M

 

, M ),                    (1)

 

Frolov, and Lykov.

 

RS               F         FS          S

 

ST         T

 

The research featured sunflower oil seeds.

The thermophysical characteristics of the vege- table raw material were determined according to the non-stationary thermal mode method of two tempera- ture-time points developed by Volkenstein. The method of differential thermal analysis was used to identify the intervals of temperature zones of moisture evaporation with different forms and energy of moisture  binding with the material. It was accompanied by the method of differential scanning calorimetry, which was used for quantitative measurement of heat flows that occurred when the trial sample and the control sample were un- dergoing a simultaneous programmed heating. The methods of high-performance gas chromatography, atomic absorption spectroscopy, infrared spectroscopy, capillary electrophoresis, and acid hydrolysis were used to determine the content of vitamins, amino acids, and other quality indicators of the wet solids. The measure- ment errors did not exceed the values established in the current standards for quantitative analysis of the quality indicators of the wet solids. The main part of the theo- retical and experimental research was carried out on the premises of the Voronezh State University of Engine- ering Technologies (Voronezh, Russia) and the Bo- brovsky Vegetable Oil Plant (Bobrov, Russia).

The research objective was achieved by the synthesis and analysis of classical and novel analytical and empi- rical methods in the sphere of heat and mass transfer and food dehydration studies. The obtained relations, the ap- proximating equations, and the simulation results corre- sponded with the experimental data. The measurement results underwent statistical processing. The procedures and design solutions did not contradict the established methods  of  rational  design  and  engineering.  Modern

computer mathematical programmes, instruments, and

 

forms  a  multimodel  system  of  transfer  phenomena

(Fig. 1).

A drum dryer as a subject of system modelling of non-stationary processes. Fig. 2 shows the experimen- tal drum dryer used in the study of the continuous dry- ing process.

The pressure type fan (2) was fixed on the angular steel frame (1). The fan was designed to supply atmo- spheric air through the duct (4) into the heater (6). The heater consisted of several sections. The flow rate of the air supplied to the heater was regulated by means of the gate valve (5). The temperature and the relative humidi- ty of the air entering the drying chamber were measured with the dry (7) and the wet (8) thermometers.

At the core of the whole design there was a steel drum (11) with a diameter of 0.3 m and a length of

1.3 m. On the inner surface of the drum there was the channel nozzle (12). The nozzle contained longitudinal slots to feed the drying agent to the bed of wet granu- lar material. The drying drum (11) was supported by two

 

 

Fig. 1. Multimodel system of transfer phenomena.

 

 

view B

 

 

Fig. 2. Experimental drum dryer with a channel nozzle.

 

 

pairs of rollers (13). It was driven by the electric motor

(26) and the gearbox (24), which were kinematically in- terconnected by the chain transmission (25).

During the drying process, the drying agent was sup- plied directly to the zone of the channel nozzles under the dryable material.

The pipe (9) was assembled on the body of the loa- ding chamber (10). The pipe fed the wet bulk material to the drying drum. The spiral rotation of the drum moved the material to the discharge chamber (14).

The retaining ring (15) was assembled on the end surface of the drum. The ring slid in the groove of the fixed flange (16), through which the bulk material was continuously unloaded. During the drying process, the interior of the drum (11) was under a slight negative pressure due to the fact that the performance of the fan (21), which took away the waste drying agent, was seve- ral times higher than that of the delivery fan (2).

The temperature control of the drying agent was car- ried out directly in the bed of the wet solids. The ther- mocouples (17, 18, and 19) were installed on the bracket

(20) (Fig. 2, view B).

The temperature was recorded with the electronic au- tomatic self-recording potentiometer (23).

The relative humidity of the waste drying agent was measured with the hygrometer (22), a sorption-frequency single-channel digital device. The semiconductor ther- moanemometer (27) measured the speed of the drying agent at the inlet and the outlet of the drum, as well as above and below the surface of the bed.

The angle of the drum, the frequency of its rotation, and the speed and temperature of the drying agent were set at the initial stage of the experiments. After that, there began a continuous supply of proportioned  wet bulk material. Over the next 60–70 min, the material was sampled to measure its thermophysical characteris- tics, while the process of wet material feeding continued in the same mode.

95

 

After the sampling, the feeding of the wet bulk ma- terial to the drying drum and the movement of the drum stopped simultaneously. The speed and the temperature of the drying agent remained constant, according to the experimental conditions. After that, all the material was discharged from the drum into the receiving container, and its volume and weight were measured.

Taking into consideration the design feature of rota- ry drum dryers and the methods of continuous drying of wet solids in a fluidized bed, the model of continuous drying process can be represented as a multimodel sys- tem of transfer phenomena (Fig. 3).

The continuous drying process as a system of transfer phenomena. The model of the continuous dryi- ng process of wet solids in a fluidized bed displays the properties of the material and the coolant and the nature of non-steady processes that occur in the fluidized bed. It also makes it possible to improve the drying at various initial parameters of the material and the coolant.

The model of the product subjected to drying and the coolant connects the model of the system with the exter- nal environment and control actions. It also calculates the properties of the material and the coolant.

The following models of transfer phenomena cor- respond with the continuous process in a rotary drum dryer (Fig. 3):

  • the model of the product subjected to drying;
  • the coolant model;
  • the model of the movement of the wet solids on the cy- lindrical surface;
  • the model of the coolant fed to the fluidized bed;
  • the hydrodynamic model of the fluidized bed;

– the model of the complex heat and mass transfer;

– the model of the drying process in the fluidized bed;

and

the model of the technical and economic characteristics.

The model of the product subjected to drying. If we consider the wet product as a subject of drying with

 

 

 

Fig. 3. Multimodel system of the continuous drying process of wet solids in a fluidized bed.

 

 

its moisture content, temperature, density of individual particles, and average linear dimensions, we can con- struct a mathematical model of the technological prope- rties of the bulk material [8]. The model can describe the physical properties necessary for hydrodynamic and heat transfer processes. In this case, the linear dimen- sions of the bulk material determine its volume, surface area, equivalent diameter, sphericity (non-sphericity) co- efficient, and the specific surface area. The equilibrium moisture, angle of repose, and the bulk density of the wet solids are determined depending on their humidity and the moisture of the drying agent. The specific heat is de- termined depending on the temperature of the product.

The model of the product subjected to drying. Let us consider the wet solids as a subject of drying, chara- cterized by moisture content, temperature, density of individual particles of the product, and average linear dimensions. By doing so, we can constructmathe- matical model of the technological properties of the wet solids [8]. The model describes  the  physical  proper- ties necessary for conducting hydrodynamic and heat transfer processes. In this case, the linear dimensions of the bulk material are determined by its volume, sur- face area, equivalent diameter, sphericity coefficient (non-sphericity), and the value of specific surface. Equi- librium moisture, angle of repose, and the bulk density of the wet solids are determined depending on their hu- midity and the moisture of the drying agent. The specific heat depends on the temperature of the product.

The coolant model. The wet and heated atmospheric air is considered as a binary mixture of the components of dry air and steam [9]. The parameters of the drying are  calculated  by  the  method  of  superposition  of  the

component  parameters.  Tables  and  various  functional

 

and empirical relations are used to calculate the parame- ters of the components.

The model of the movement of wet solids on the cylindrical surface. In drum dryers, the material is dried in a fluidized bed while it is purged with a drying agent [5]. Fig.4 shows the flow chart of the bulk material in the transverse direction.

The angle of the bed in the cross section is determined by the angle of repose of the bulk material. It is assumed that the trajectory of the particles passes along a radial arc, whose radius equals the distance from the axis of the drum. When the particle gets to the end of the arc and sur- faces, it falls down along the chordal surface of the flow.

The particle moves in the axial (translational) direc-

tion only if the flow surface has a certain axial angle, i.e.

 

 

Direction of rotation

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fig. 4. Scheme of the movement of the wet solids in the rota- ting drum (cross section).

 

 

 

when there is a difference in the levels of the bed at the inlet and the outlet holes of the drum.

During translational motion, the trajectory of the particle follows the chordal plane. In this case, the tra- jectory is not perpendicular to the axis of the drum: it forms a < 90° angle with it and is perpendicular to the plane passing through the middle of the chords.

The model of the coolant fed to the fluidized bed. The drying agent enters the drying chamber of  the drum dryer through the channels formed by the channel nozzles [6] and the outer cylinder of the drum through the side cavity of the drum (Fig. 5).

The side cavity of the drum was divided by a par- tition in such a way that the drying agent entered the channels under the nozzles that are situated under the material while the drum is rotating. The channel nozzle has a gap of constant width. The drying agent enters the drying chamber through this gap and contacts the layer of dispersed material.

Thus, the task is reduced to the calculation of gas distribution with an outflow through the lateral perme- able surface, i.e. the layer of the dispersed material. It is possible to make an assumption that the parameters of

 

area of the product or the specific surface area of the bed and the speed with which the drying agent flows around the particles.

When solving the problem of the hydrodynamics of a bed, one has to determine the hydraulic resistance of the bed of the dispersed material. It is necessary to know the characteristics of the pore channel of the bed, its coeffi- cient of the hydraulic resistance, and the flow rate of the drying agent in the channel.

For generality, it is assumed that the flow velocity around the particle and the flow speed in the channel are equal and related to the flow velocity in the direction of filtration in the whole liquidized bed by the ratio:

w

 

0

= ξw.                                        (2)

 

The tortuosity coefficient ξ is calculated as follows:

 

ξ = 1 + (π/21)(1ε)2/3.                       (3)

 

The bed voidage ε is the ratio of voids between the particles in the layer and the volume of the bed. It can be expressed as the following ratio:

 

the drying agent and the height of the bed are constant.

Thus, one can calculate the hydrodynamics of the flow

 

ε = 1 – ρ

 

/ ρ .                                (4)

 

p

 

d

of the drying agent. In this case, the bed height may be its height in the middle section of the material flow. The parameters of the drying agent, such as density and ki- nematic viscosity, are assumed to be constant.

Two models were considered when choosing the mathematical model for the coolant supply from the channel nozzle slit through the bed of dispersed material. The first was the model for calculating the distribution of velocity and pressures along the z-shaped collectors; the


The specific surface of the material particles in the

bed can be calculated through the specific surface of the

s

particles S and layer porosity ε:

 

sg            s

S   = S (1 ε).                                 (5)

 

The loss of pressure during the movement of the dry- ing agent through the granular layer can be calculated similarly to the pressure loss in pipelines:

 

second was the model of the constant-section air distribu- tor with a longitudinal slit of constant width.

 

Δ = λ S

 

hρw2 / (2ε2).                       (6)

 

g

 

sg

The hydrodynamic model of the fluidized bed. Ac- cording to the hydrodynamics of the fluidized bed of the wet solids in rotating drum dryer, the drying agent flows around the particles in the fluidized bed of the granular product and between the particles of the material in the channels [6].

In this connection, we can point out external (flow-a- round), internal (filtration), and mixed hydrodynamics problem.

When solving the problem of heat and mass trans- fer,  it  is  necessary  to  determine  the  active  surface

 

 

 

Fig. 5. Scheme of the coolant flow in the rotating drum dryer.


The model of the complex heat and mass trans- fer. In the general case,  when  organizing  continu- ous drying of wet solids in a fluidized bed, the coolant is supplied to the drying chamber at the expense of an external heat source (mechanical energy). As  a  re- sult, there is a continuous forced flow of the heat transfer agent particles around the product.  The  pro- cess of continuous evaporation of moisture  from  the free surface of the moving material occurs,  in  this case, in the  boundary layer and  is caused  by veloci- ty and temperature gradients. As a  result,  there  oc- curs a continuous diffusion flow of the medium. This continuous flow is constantly directed from  the  sur- face of the product into the depth of the coolant flow. The difference in the concentration of the vapor-gas mixture near the evaporation surface and in the main heat-wave flow of the coolant leads to a density diffe- rence in the vapor-gas mixture. It results in a continuous free natural heat and mass transfer.

Fig. 6 demonstrates the external heat and mass trans- fer under conditions of continuous dehydration as a com- bination of simple (marginal) modes of motion: forced, free, and diffusive.

To calculate the simultaneous continuous heat and mass transfer processes, we used the method based on the superposition of the absolute values of Nusselt num-

 

 

                                           

 

 

er

 
Fig. 6. Motion modes during continuous drying: Re – forced; Arfree natural; PWdiffusive.

Fig. 7. The motion of the medium near the heat exchange surface in mixed convection.

 

 

 

ber [9]. In a complex continuous process, the absolute

 

The effective heat transfer coefficient α

 

is deter-

 

value of the transfer rate is determined according to the projection values of the transfer rate of simple processes onto two mutually perpendicular planes.

The velocity attitude of the forced motion of the me-

 

mined by a nonlinear differential equation, presented in dimensionless form:

dy/dx = A(y/x) + Rb,                          (13)

 

Re

 
dium U

makes an arbitrary angle φ with the planes, in

 

where y = t T ; x = u u ; A = αS / c

 

K; Rb = r(u u  /

 

Ar

 
which the velocity vectors of the free motion U

lie.

 

i                            e                       s        m                                 i         e

 

Fig. 7 shows that the calculated dependence of the forced and the free motion with an arbitrary mutual ori- entation is

 

[c (t T)] is Rehbinder number.

m

 
The solution for the equation (13) is determined by the following power function:

 

Nu = [(Nu4

 

 

 

ORe

 

cos4φ ± Nu4

 

)0,5) + Nu2

 

 

 

ORe

 

sin2φ]0.5.  (7)

 

Rb = c[(t T ) / (T T )]m

 

S

 

/ c K)n.         (14)

 

OAr

 

i             f           i

t

s

When φ = 0, the forced and the free motions happen

in the same plane:


The effective heat transfer coefficient for the flu- idized bed of wet solids under continuous dehydration conditions is determined by the following relation:

 

Nu = (Nu4

 

 

ORe

 

± Nu4

 

)0.25.                     (8)

 

α  = (c

 

K / S )[lnRblnc + mln(t T ) / (T T )] / n. (15)

 

OAr

 
When φ = π/2, the forced and the free motions hap-

er             m              s

 

i             f           i

 

pen in mutually perpendicular planes:

 
The model of the drying process in the fluidized

bed. Let us consider the general case of an approximate

 

Nu = (Nu2

 

 

ORe

 

± Nu2

 

)0.5.                       (9)

 

mathematical  description  of  a  continuous  process  of

 

OAr

 
If the critical geometrical dimensions of the forced

moisture transfer:

 

and the free motions are the same (L = L

 

= L  ), formu-

 

dx /dτ = Kx (τ),

 

Re            Ar                                                                                                    1                          1

la (8) coincides with the following relation:

dx /[G / U(τ)]Kx (τ),                         (16)

4             4                      4                                                                                                                           2                                   1

 

Nu =

 

NuORe  + NuOAr .                       (10)

 

where, in terms of optimal control methods [8], the pa-

 

In accordance with the nature of the flow of the me-

dium in the liquidized bed, the calculated relation of the

 

rameters x

1

 
control act

and x

2

 
.

are phase variables, and U(τ) is the

 

intensity of the continuous transfer of the complex pro- cess meets conditions (8) and (10).

The absolute value of Nusselt number in forced mo-

 

ion

opt

 
The system of equations (16) is a combination of the kinetics of the continuous drying and the balance rela- tion of moisture in the material and coolant.

 

tion conditions is calculated by the formula for the ball:

 

Let us define the optimal control U

 

(τ) which trans-

 

Nu

 
OAr

= 2 + 0.56(Ar × Pr)0.25[Pr / (0.846 + Pr)]0.25,        (11)

 

fers the continuous process from

1                 1

 

2                  2

– the given initial state:

 

where 1 < Ar Pr < 105. In conditions of free motion it is calculated by the formula for the sphere:

 

x (0) = x

 

Nu

 
ORe

= 2 + 0.03(Re0.54Pr0.33  + 0,35Re0.58Pr0.36,       (12)

 

0 and x (0) = x 0                              (17)

 

– and in the specified final state:

 

1    f             1

 
f and x ) = x f                              (18)

 

where 0,6 < Pr < 8×103 and Re < 3×105.

In formulae (11) and (12), the critical value is the

 

x ) = x

 

2    f             2

 

equivalent particle diameter of the wet solids.

 

so that the functional assumes the minimum value:

t к

 

er

 
The effective heat transfer coefficient α

is used as a

 

I = òU (t )dt .                            (19)

 

characteristic of heat transfer intensity. It takes into ac-

count the total amount of heat spent on continuous dry-

 

0

Let us establish a function with auxiliary variables

 

ing in a rotating drum dryer.

 

λ , λ

 

and λ :

 

0      1                 2

 

 

 

H = λ U – λ Kx  + λ (GK/U )x

 

(20)

 

In rotary drum dryers, the warm-up periods are short,

 

0               1        1          2                  t     1

 

and the drying rate is constant. Hence, the drying rate

 

1

 
and write the system of equations for the function λ :

dl1 / dt = H / x1 = l1K - l2 (GK / UТ );

dl2  / dt = H / x2  = 0.

 

 

 

 

(21)

 

can be represented as a function of moisture content [10]:

p

 
du/dτ = – K(u u ).                           (27) For rotary drum dryers, u = dz/dt. In this case, if we

 

The optimal control of the continuous process results from the condition that the function H has the maximum value:

2

 

take into account Eq. (25) and Eq. (26), the drying rate can be represented as follows:

p

 
du/dz = – (KSρ/q)(u u ).                    (28)

 

 

 

Thus,

 

H / UТ = l0 - l2 (GK/ UТ ) x1 = 0,

 

(22)

 

 

The  heat  is  continuously  supplied  to  the  dryable product. It heats the material and evaporates the mois-

 

Uopt (t ) =

 

(l2  / l0 )GKx1 .

 

(23)

 

ture:

 

The following function solves the equations of the system (7), which describe the kinetics of continuous drying with the boundary condition (8):

 

dT/dz = [(αS (tT) / v) + r(du/dz)] / (c  + c u),  (29)

 

s                                                                    m           w

 
S = S (1ε) / (V ρ ).                         (30)

 

0    - Kt

 

s           p                            p   m

 

x1 (t ) = x1 e     .

 

(24)

 

In Eq. (29), we point out, first, the heat supplied to

 

In the rotating drum dryer, the dryable product is, as a rule, continuously loaded into the drying chamber on one side and unloaded on the other. The height of

 

the material,

 

 

 

q

 

s

1

= αS (t T),                                (31)

 

the bed at the input exceeds the height of the bed at the output, due to the rotation frequency and the horizontal angle of the drum. In general, the coolant is fed to the drying chamber and passes through the fluidized bed of the bulk material [5].

To study the dynamics of the continuous dehydra- tion mechanism [9], the original scheme of interaction of the coolant flow and the fluidized bed of the wet solids moving along the cylindrical surface can be presented as (Fig. 8):

For practical calculations, the surface of the fluidized (sliding) bed is represented as a plane, and the cross sec- tion of the material flow is represented by a segment of a

 

second, the heat spent on heating the mass of absolutely dry material and moisture contained in it,

2               m           w

 
q  = v(c  + c u)(dT/dz),                       (32)

 

third, the heat spent on the evaporation of moisture from the material,

3

 
q  = vr(du/dz).                             (33)

 

Using the continuity equation (25), we express the equation for the temperature of the material through its moisture content:

 

circle. To describe the continuous drying mode, the con-

tinuity equation for the flow of the dryable material can

 

dT/du = – S (t T)] / c

 

K(u u )] + r / c

 

,    (34)

 

s

 

e

am

am

be represented as

 

where c   = c


+ c u.

 

w

 

am

m

q = ρvS,                                      (25)

 

Note that the density of the dryable material is the function of its moisture content u:

ρ(u) = a bu.                                (26)

 

 

 

Fig. 8. Diagramme of interaction of the coolant flow and the fluidised bed of the dryable product moving along the cylindri- cal surface.


To formulate the balance of the coolant, we divide the fluidized bed of the material of length l into a bed of infinitely small sections of length Δz. Let us assume that, within each of such moving sections, the material is perfectly mixed and has a constant temperature and moisture content. When moving from section to section, the temperature and the moisture contents of the mate- rial change to infinitely small values.

Let us single out one of these sections and consider the material flow in the moving element of the wet solids (Fig. 9).

The flow G(1 + u) of the wet solids enters the volume of the moving element, and the flow G(1 + u + du) comes out with its moisture content being u + du. The coolant supplied to the fluidized bed is determined by the expen- diture function L(z). When they enter the fluidized bed, the moisture content and the temperature of the coolant are constant throughout the whole layer. They are func- tions of the z coordinate at the output. The amount of moisture evaporated from the elementary volume of the

material equals the amount of moisture absorbed by the

 

 

 

s

 

am

dT/du = – S (t T) / c

 

where dT/du = ν(dT/dτ).


K(u u )] + r / c


,   (43)

 

e

 

am

By  integrating  Eq.  (43)  with  the  initial  condition

i

T(0) = T , we get:

 

T(u) = T φ(u) + {t + [rK(u u ) / αS ]}[1φ(u)],  (44)

i                                                       e               s

 

where φ(u) = (u u ) / (u u ).

e             i          e

 

 

 

 

 

 

 

 

Fig. 9. Material flows in the moving element of the wet solids.

 

coolant that enters the elementary volume:

 

s

 
dL(z) / V = – Gdu.                                (35)

 

The heat balance equation is similar to the moisture balance equation (26) for the elementary volume:

d[L(I I )] / V = G(di + rdu),                  (36)

 

Formula (44) gives us values that are close to those

we obtain after integrating Eq. (37), while the ratio error does not exceed 2.0%. This allows us to declare Eq. (43) and its analytical solution (44) applicable for analytical studies of continuous dehydration in a fluidized bed of bulk material.

Eqs. (37), (38), (40), and (41) make it possible to cal- culate the moisture content and the temperature of the dryable product in the fluidized bed along the length of the drying drum. We know the values of moisture con- tent and coolant temperature only at the surface of the material. To calculate the parameters of the coolant at the outlet of the drying drum, we use the formulae for coolant mixing [11].

The equations for the moisture content and the tem-

perature of the coolant at the outlet from the elementary

 

i             s

volume of the fluidized bed (Fig. 9) can be formulated

 

w

 

, c

m

is

where i = (c

 

+ c u) is the enthalpy of the material;


on the basis of  Eqs. (31) and (32), if we assume that the

 

I = c

 

t + (r + c t)x is the enthalpy of the coolant; c

 

 

hc      s

 

temperature tj  and the moisture content xj  of the coolant

 

s

 

hc

specific heat capacity of the coolant and the steam.

Thus, the continuous drying process on a cylindrical


at the outlet from the elemental volume are constant, and the coolant rate through this elementary volume is the

 

surface can be approximately described by a system of

 

difference in rates dL = L

 

L , where L is the coolant

 

ordinary differential equations:

 

j

rate when it is suppli

 

j + 1           j                          j

to         flu dized be  in the section

 

 

p

 
du/dz = – (KSρ / q) (u u ),                (37)

 

j

 
[0, z ]:

ed      the        i               d

 

 

s

 

m

dT / dz = [(αS (t T) / ν) + r(du/dz)] / (c with initial conditions

u(0) = u , T(0) = T


 

w

+ c u), (38)

 

 

(39)


x = x + (GV / dL )(u       u ),                  (45)

 

j          i                  s            j        j + 1          j

i                 s             j        j         i                 j + 1           j

Ij = I (GV / dL )[(i i ) + r(u        u )],      (46)

 

j            j             j             hc           s   j

t = (I rx ) / (c   + c x ),                        (47)

 

i                         i

 

and  a  system  of  balance   relations   obtained   after Eq. (35) and Eq. (36) were integrated with initial condi- tions x(0) = x and t(0) = t :

i                                i

 

x = x + (GV / L)(u u),                      (40)

i                 s                 i

 

where i = c T ; j = c T ; I = c  t + (r + c t )x .

i          m   i                m   j     s           hc  i                      s  j     j

 
The values of moisture content and  the  tempera- ture of the coolant at the outlet of the drying chamber are calculated from the parameter values of the coolant flows that are coming out of all the elementary volumes of the fluidized bed, according to using the mix formu-

 

t = {I rx (GV / L) [(i i ) + r(u u )]} / (c

 

+ c x). (41)

 

i                              s                           i

 

i                  hc          s

 

lae of the flows:

n

 

Upon integrating Eq. (37), we receive an analytical

expression for the moisture content of the material in the coordinates of the length of the rotating drum:

 

xê  = (1 / L)

 

å

j=1

 

n

 

xjdLj;

 

(48)

 

 

p

 
u(z) = (u

 

+ me–a) / (1 + ue–a),                   (42)

 

Iк  = (1/ L)å I jdLj;

j=1  

 

(48)

 

 

where   m   =   a(u

 

 

–   u )/ρ(u );   n   =   b(u

 

 

–   u )/ρ(u );

 

t = (I rx ) / (c   + c x ).                          (49)

 

f            f             f             hc           s   f

 
i                 e             i

i                 e             i

 

α  =  Kρ(u )V(z)/q ;  q

 

=  G(1  +  u );  ρ(u )  =  a  –  bu ;

 

e                  i          i

 

e                 e                                    p

 

If we use Eq. (42) and take into account that the ki-

 

ρ(u ) = a bu ; V(z) is the volume of the fluidized bed

e                              i                                                                                                                              netics of continuous drying is described by Eq. (27),

 

from the upload point to the coordinate z.

Unlike Eq. (37). Eq. (38) makes it possible to obtain an analytical solution only with the assumption that the heat capacity of the product remains average. If we agree

 

while the flow of the bulk material moving on the cylin- drical surface is presented as an ideal extrusion model (25), we can determine the equation of the effective con- tinuous drying coefficient:

 

m

 
that c = c

+ c u

 

, where u

 

is the average moisture

 

w

 

am

am

content of the material, we obtain the following equation:

 

= {q / [V ρ(u )]}ln[ρ(u )(u u ) / ρ(u ) (u u )]. (50)

 

e                         c         e

 

f        i          e

 

i          f           e

 

 

 

Relation (50) differs from others that determine the

coefficient of drying. It is based on the analysis of the

 

For a convective drying process, the energy efficien- cy can be expressed by the following relation:

 

mechanics of a fluidized bed of wet solids on a cylindri-

 

η = (t

 

t ) / (t

 

t ).                          (52)

 

cal surface with regard to the continuous dehydration.

The model of the technical and economic charac- teristics. To assess the energy performance of a drying unit with the convective method of heat supply, one can assess the use of the drying agent. The energy losses are determined by the difference between the amount of the supplied and the usable energy [11, 12–15].

Various coefficients of efficiency are used as energy

criteria. In the general case, they are defined as the ratio

 

t           1         2           1         0

 

The thermoeconomic analysis combines exergy ana- lysis and economic optimization. The criterion for the thermoelectric optimization is a composition of additive functions. These functions should quantify the exergy, the equipment costs, etc.

The most general formula for the so-called thermo- electric criterion is

ì      éå cei ei + å kn ùü

 

of the usable energy E

 

to the expended energy E :

 

{min C}= ïmin ê j               n

 

úï.               (53)

 

1

 

1

 

/ E

2.

η = E

 

2

 

(51)


í      ê

î

ï      êë

 

å enk

k


úý

þ

úûï

 

 

Table 1. Baseline input

 

Input parameter                                                                                                   Variants of the computational experiment

1                 2                  3                4                 5                 6

1.1. Rotary drum dryer

Dryer drum length, m

1.2

1.2

1.2

1.2

1.2

1.2

Outer radius of the drum, m

0.15

0.15

0.15

0.15

0.15

0.15

Channel nozzle radius, m

0.115

0.115

0.115

0.115

0.115

0.115

Width of the slits of the channel nozzle, m

0.02

0.02

0.02

0.02

0.02

0.02

Number of slits channel nozzles, pcs.

12

12

12

12

12

12

Number of slits in the channel nozzle through which the coolant is supplied, pcs.

4

4

4

4

4

4

Coefficient of local resistance

0.16

0.16

0.16

0.16

0.16

0.16

 

1.2. Product

 

 

 

 

 

Density of the material particles, kg/m3

770

770

770

770

770

770

Geometrical dimensions of the material particles, m:

 

 

 

 

 

 

length

10.7

10.7

10.7

10.7

10.7

10.7

width

5.0

5.0

5.0

5.0

5.0

5.0

thickness

3.3

3.3

3.3

3.3

3.3

3.3

Initial moisture content, kg/kg

0.105

0.105

0.105

0.105

0.105

0.105

Initial temperature, °C

19

20

14

17

13

14

Specified final moisture content, kg/kg

0.0547

0.0705

0.0546

0.0828

0.0574

0.0387

Set final temperature, °C

56

55

59

42

55

66

 

1.3. Coolant

 

 

 

 

 

Barometric pressure, kPa

100.5

100.5

100.4

100.3

100.3

100.3

Outside temperature, °C

16

15

14

16

14

16

Outside air humidity, %

73.0

73.0

81.8

82.5

81.9

72.0

Temperature of the drying agent at the inlet of the drying chamber, °C

180

210

240

180

210

240

Moisture content of the drying agent at the inlet of the drying chamber, kg/kg

0.008413

0.007886

0.008295

0.009550

0.008314

0.008314

Temperature of the drying agent at the outlet of the drying chamber, °C

149

167

187

134

165

185

Specified moisture content of the drying agent at the outlet of

the drying chamber, %

0.0243

0.0275

0.0332

0.0305

0.0294

0.0342

1.4. Process parameters

Dum angle, rad

 

0.03490

0.03490

0.05236

0.05236

0.01745

0.01745

Drum rotation frequency, 1/sec

 

0.0250

0.0583

0.0250

0.0583

0.0417

0.0417

Product feed rate, kg/sec

 

0.0106

0.0174

0.0167

0.0348

0.0159

0.0121

Drying agent consumption, m3/sec

 

0.04010

0.03868

0.04559

0.04396

0.04529

0.04199

Radius of the circle touching the bed at the input point, m

 

0.03

0.01

0.02

0.03

0.02

0.01

Consumption coefficient of the drying agent

 

0.6

0.6

0.6

0.6

0.6

0.6

Absolute roughness of the air duct wall, m

 

0.1×10–3

0.1×10–3

0.1×10–3

0.1×10–3

0.1×10–3

0.1×10–3

Filling rate of the drum, %

 

25

35

30

25

30

35

 

 

Table 2. The results of the computational experiment (output)

 

Input parameter                                                                                                                                                                   Variants of the computational experiment

1                           2                           3                           4                           5                             6

2.1. Calculation results for the product model subjected to drying

 

Volume of the particle, m3 Surface area of the particle, m2 Sphericity coefficient

0.7506×10–7

0.1608×10–3

0.53507

0.7506×10–7

0.1608×10–3

0.53507

0.7506×10–7

0.1608×10–3

0.53507

0.7506×10–7

0.1608×10–3

0.53507

0.7506×10–7

0.1608×10–3

0.53507

0.7506×10–7

0.1608×10–3

0.53507

Aspheric coefficient

1.87

1.87

1.87

1.87

1.87

1.87

Equivalent particle diameter, m

0.5234×10–2

0.5234×10–2

0.5234×10–2

0.5234×10–2

0.5234×10–2

0.5234×10–2

Equilibrium moisture content in the material, kg/kg

0.02183

0.02180

0.02183

0.02190

0.02193

0.02193

Angle of friction, rad

0.7345

0.7345

0.7345

0.7345

0.7345

0.7345

Loose weight density, kg/m3

398.5

398.5

398.5

398.5

398.5

398.5

Specific heat capacity of the material, kJ/(kg×K)

1.5160

1.5164

1.5140

1.5152

1.5136

1.5140

Specific surface of the particles, m2/m3

2142.62

2142.62

2142.62

2142.62

2142.62

2142.62

2.2. Calculation results for the coolant model

Moisture, %

1.33448

1.25207

1.31597

1.51116

1.31889

1.31894

Wet thermometer temperature, °C

45.2109

47.8192

50.2885

45.4979

47.9154

50.2922

Specific heat of dry air, kJ/(kg×K)

1.022

1.028

1.035

1.022

1.028

1.035

Specific heat capacity, kJ/(kg×K)

2.710

3.200

3.880

2.710

3.200

3.880

Specific evaporation heat, kJ/K

2015.20

1900.50

1766.00

2015.20

1900.50

1766.00

Coefficient of kinematic viscosity of the drying agent, m2/sec

0.3029×10–4

0.3136×10–4

0.3032×10–4

0.3000×10–4

0.3111×10–4

0.3029×10–4

Drying agent density, kg/m3

0.84373

0.84722

0.90496

0.85300

0.85459

0.90661

Specific volume of wet air, m3/(kg·sec)

1.31526

1.40113

1.49053

1.32022

1.40488

1.49211

Heat conductivity coefficient of the drying agent, W/(m·K)

0.3773×10–1

0.3998×10–1

0.4207×10–1

0.377×10–1

0.3998×10–1

0.4207×10–1

Prandtl number

0.7017

0.6946

0.6906

0.7041

0.6956

0.6907

Schmidt number

0.5638

0.5202

0.4512

0.5584

0.5160

0.4507

2.3. Calculation results for the model of the movement of wet solids along the cylindrical surface

Specific consumption of the drying agent, kg/kg

62.9435

50.9856

40.1518

47.7182

47.4229

38.6292

Minimum design airflow per drying, kg/sec

3.0371×10–2

2.7698×10–2

3.0584×10–2

3.3362×10–2

3.2481×10–2

2.8045×10–2

Minimum estimated volume flow rate of the drying agent m3/sec

0.03995

0.03881

0.04559

0.04405

0.04563

0.04185

Fictitious speed of the drying agent through the material bed, m/sec

0.2205

0.1941

0.2385

0.2417

0.2369

0.2108

Speed of the drying agent, reduced to the full cross section of the bed, m/sec

0.4570

0.4024

0.4943

0.5010

0.4910

0.4368

Porosity of the bed

0.4825

0.4825

0.4825

0.4825

0.4825

0.4825

Specific surface of the material in the bed, m2/kg

1.5913

1.5913

1.5913

1.5913

1.5913

1.5913

Equivalent pore channel diameter, m

0.1740×10–2

0.1740×10–2

0.1740×10–2

0.1740×10–2

0.1740×10–2

0.1740×10–2

Tortuosity coefficient of the channels

1.3679

1.3679

1.3679

1.3679

1.3679

1.3679

Length of the pore channels, m

0.09375

0.1193

0.1070

0.09375

0.1070

0.1198

Equivalent Reynolds number

78.9714

67.1395

85.3087

87.4081

82.6012

75.4855

Hydraulic resistance coefficient of the bed

1.0080

1.1418

0.9516

0.9347

0.9746

1.0431

Bed resistance, Pa

507.731

447.649

601.354

572.002

573.984

515.488

Heat transfer coefficient, kW/(m2 × K)

0.206419

0.186365

0.224280

0.223457

0.220131

0.201593

Residence time in the drying chamber, min

7.899

6.710

6.016

2.389

6.315

9.723

Drying coefficient, 1/sec

0.1959×10–2

0.1330×10–2

0.2580×10–2

0.2168×10–2

0.2242×10–2

0.2734×10–2

2.4. Results of the calculation for the model of the coolant supplied to the fluidized bed

Cross-sectional area of the air distributor, m2

0.9712×10–2

0.9712×10–2

0.9712×10–2

0.9712×10–2

0.9712×10–2

0.9712×10–2

Perimeter of the air distributor, m

0.8350

0.8350

0.8350

0.8350

0.8350

0.8350

Equivalent diameter of the channel nozzle, m

0.04653

0.04653

0.04653

0.04653

0.04653

0.04653

Speed of the drying agent at the beginning of the drum, m/sec

4.1287

3.9852

4.6940

4.5261

4.6631

4.3233

Reynolds number

6342.55

5907.68

7202.19

7020.14

6973.78

6641.71

Average flow rate of the drying agent from the channel nozzle, m/sec

0.4177

0.4029

0.4749

0.4579

0.4718

0.4374

Coefficient of friction of the air nozzle

0.03705

0.03761

0.03609

0.03628

0.03633

0.03670

Coefficient of friction of the air nozzle friction

1.14102

1.14166

1.13992

1.14014

1.14019

1.14061

Total resistance of the air nozzle, Pa

8.2

76

11.4

9.6

10.6

9.7

Slit parameter

5.931

5.931

5.931

5.931

5.931

5.931

Air duct parameter

5.082

5.097

5.058

5.062

5.064

5.073

2.5. Results of the calculation for the hydrodynamics model of the fluidized bed

Volumetric capacity of the dryer for wet material, m3/sec

0.266×10–4

0.4366×10–4

0.4191×10–4

0.8733×10–4

0.3990×10–4

0.3036×10–4

Dryer productivity according to absolutely dry material, kg/sec

0.9593×10–2

0.1575×10–1

0.1511×10–1

0.3149×10–1

0.1439×10–1

0.1095×10–1

Dryer productivity according to the evaporated moisture, kg/sec

0.483×10–3

0.543×10–3

0.762×10–3

0.6991×10–3

0.6849×10–3

0.7260×10–3

Angle between the surface and the axis of the drum, rad

0.02775

0.02911

0.02819

0.02775

0.02818

0.02913

Radius of the circle touching the bed at the output of the product, m

0.06331

0.04494

0.05383

0.06331

0.05382

0.04497

Volume of the drum occupied by the bed, m3

0.01246

0.01745

0.01496

0.01246

0.01496

0.01745

Area of the middle section of the bed, m2

0.01039

0.01454

0.01246

0.01039

0.01246

0.01454

Radius of the circle touching the bed in the middle section, m

0.04646

0.02742

0.03677

0.04646

0.03676

0.02744

Thickness of the bed in the middle section, m

0.06855

0.08759

0.07823

0.06853

0.07824

0.08756

Width of the bed in the middle section, m

0.1516

0.1660

0.1593

0.1516

0.1593

0.1660

 

 

The rest Table 2

 

Effective area of the bed, m2

0.1819

0.1932

0.1912

0.1819

0.1912

0.1992

Distance between the beginning of the drum and the middle section of the bed, m

0.5931

0.5982

0.5945

0.5930

0.5945

0.5982

Influence coefficient of the flow rate of the drying agent in the dense blown

bed on the performance of the dryer

2.24200

1.41040

2.52000

2.40090

2.74770

1.98300

2.6. Results of the calculation for the model of complex heat and mass transfer

Equivalent Reynolds number

71.1793

67.2305

81.9623

79.8905

79.3629

75.5839

Archimedes number

455.991

479.825

565.384

463.839

487.412

455.778

Limit value of Nusselt number for natural convection conditions

2.4107

2.4355

2.5338

2.4231

2.4459

2.5354

Nusselt number for simultaneous processes

5.5678

5.4071

5.8927

5.8149

5.7896

5.7035

Heat transfer coefficient, kW/(m2×K)

0.04014

0.04131

0.04737

0.04191

0.04423

0.04585

Specific heat flow, kW/m2

5.4107

6.6998

8.9866

5.6374

7.1691

8.6980

Specific mass flow, kg/(m2×sec)

0.2685×10–2

0.3525×10–2

0.5089×10–2

0.2797×10–2

0.3772×10–2

0.4925×10–2

2.7. Results of the calculation for the model of drying in a fluidized bed

The final design value of the moisture content in the material, kg/kg

0.0546

0.0704

0.0545

0.0828

0.0573

0.03856

The final design value of the material temperature, °C

56

55

59

42

55

66

The final calculated value of the moisture content in the drying agent,

kg/kg

0.0243

0.0276

0.0332

0.0306

0.0296

0.0342

The final design value of the temperature in the drying agent, °C

156.664

171.018

191.130

138.289

171.987

193.168

Effective heat transfer coefficient, kJ/(m2×K)

1.4338×10–3

1.0687×10–3

1.3248×10–3

2.4044×10–3

1.4310×10–3

1.0144×10–3

2.8. Results of the calculation for the model of technical and economic characteristics

Dryer productivity for moisture removal, kg/sec

0.4834×10–3

0.5441×10–3

0.7631×10–3

0.7001×10–3

0.6862×10–3

0.7276×10–3

Exergy of the drying agent at the inlet to the drying chamber, kJ

35.32

48.15

62.78

35.39

48.74

61.56

Specific exergy, kJ/kg

2227.66

2443.20

2515.93

1683.07

2289.72

2380.89

Energy efficiency

7.2

6.0

5.0

10.5

6.4

5.3

Capacity for evaporation according to evaporated moisture, kg/m3

0.9696×10–2

1.0913×10–2

1.5307×10–2

1.4041×10–2

1.3763×10–2

1.4593×10–2

 

 

RESULTS AND DISCUSSION

The research  was  based  on  the  informational and structural scheme of a convective drying unit model. For all its components, we developed the mathematical mod- els in accordance with the analytical multimodel system for the continuous drying process of wet solids in a flu- idized layer. As a result, we constructed an automated calculation system for the continuous process of convec- tive drying, which can be applied to a rotary drum dryer (Tables 1 and 2).

The practical result of the study consisted in assess- ing and comparing the quality indicators of sunflower seeds of natural moisture. First, one sample of sunflower seeds was dried on an experimental drum dryer with a channel nozzle. Second, another sample was dried on an industrial drum dryer with a lifting vane system at the vegetable oil plant ZAO ZRM Bobrovsky. Finally, the experimentally obtained results were compared with the

 

The effect of the drying mode on the change in the quality of oil in the sunflower seeds (Table 4) was mea- sured by changing the acid, peroxide, and iodine num- bers at a different initial seed moisture. The heating temperature did not exceed the maximum permissible temperature for the particular humidity. It ensured the inactivation of enzymes, i.e. lipase and lipoxygenase.

Table 4 shows that the acid values of the oils in the studied modes were somewhat reduced. This can be ex- plained by the fact that low molecular organic acids were distilled together with the water steam during the dry- ing process. The peroxide numbers somewhat increased with increasing temperature, which can be explained by the catalytic effect of temperature on fat oxidation due

 

Table 4. Effect of drying process of sunflower seeds in a drum

dryer with a channel nozzle on the quality of vegetable oil

 

computed results obtained from the mathematical model.

 

Drying agent

 

Acid number,

 

Peroxide

 

Iodine

 

The experimental data show (Table 3) that after the

sunflower seeds were dried in a drum dryer with a chan-

 

 temperature, °C     mg KOH                number, % I     number, g I      

2                                 2

10.56% moisture

 

 

nel nozzle, the difference between the maximum and

0

1.80

0.016

151.6

minimum  humidity  of  individual  seeds  decreased  by

130

1.71

0.021

148.1

2.34 times. This can be explained by the same residence

150

1.68

0.024

145.6

 

time in the drying zone and the uniform distribution of

the coolant flow in the fluidized bed.

 

 170                        1.65                   0.035                142.2           

14.45% moisture

 

 

0

1.82

0.015

149.8

130

1.71

0.017

149.2

150

1.65

0.019

147.9

170

1.63

0.030

144.5

 

 
Table 3. Change in humidity of single sunflower seeds during

the drying process in the drum dryer with a channel nozzle, %

 

 

Before drying                               After drying

 

18.42% moisture

 

min

max

min

max

 

0

1.85

0.014

146.3

14.02

14.69

5.19

5.51

 

130

1.79

0.017

143.8

12.28

12.68

4.19

4.34

 

150

1.76

0.018

141.2

10.04

11.08

3.65

4.10

 

170

1.72

0.025

140.5

 

 

 

Table 5. Physical and chemical indicators of sunflower seeds

and oil

 

rial in a rotary drum dryer, most researches determine

the average values of heat transfer coefficients. The pro-

 

                                                                                                                posed approach for calculating the effective heat transfer

 

Parameters                                            Drying method

 

coefficient in a fluidized bed provides the required repro-

 

In a drum dryer with

 

In an industri-

 

ducibility and differs from the experimental data by no

 

                                         a channel nozzle          al drum dryer   

 

more than 2.0% (Table 2).

 

Drying agent tempera-

ture, °C

Seed moisture, %

 

150                              280

 

The energy performance of rotary drum dryers with a convective method of heat supply can be assessed ac-

 

initial

final

 

14.35

7.15

 

14.35

7.20

 

cording to the degree of the coolant use. The energy

losses  are  determined  by  the  difference  between  the

 

Oil content on abso- lutely dry matter, %

 

55.82                           56.36

 

amount of supplied and usable energy. It is more diffi-

cult to determine the optimal variant if it is necessary

 

Damaged seeds, %          4.15                             6.15

 

to satisfy several efficiency conditions. In this case, one

 

Phosphatides in the oil, %

Peroxide number, % I

2

0.019

0.025

 

Iodine number, g I

2

147.8

160.6

Nomenclature

 

 
Acid number, mg KOH

0.050                           0.049

1.65                             1.87

 

should use compromise criteria, e.g. capital and energy costs, capacity, quality of the finished product, reliabi- lity of the management system, level of environmental safety, etc.

 

to the presence of oxygen in the air. The iodine numbers decreased with increasing temperature. This resulted from the chemical reactions of breaking double bonds in

 

 

 

τ  –  residence  time  of  the  bulk  material  in  the  dryer

drum, sec;

G dryer capacity, kg/sec;

3

 
ρdb dry bulk density, kg/m ;

 

the carbon chain of the fats and the addition of organic compounds and radicals that were present in the air.

Table 5 features some results of the comparative pro- duction tests. They confirm the fact that the temperature of the drying agent destroys protein structure. The num-

 

V dryer volume, m3;

d

 

1

G  – amount of material in the drum, kg;

2

G  – the amount of the material leaving the drum per unit of time, kg/sec;

m

V   volume of material unloaded from the drum, m3;

x (τ), x (τ), – the moisture content of the product and the

 

ber of damaged sunflower seeds when dried in a dense             1             2

 

ventilated  bed  of  moving  seeds  is  significantly  lower than in the fluidized bed.

The analysis of the physicochemical parameters of the oil suggests that the structure of the drying agent largely determined the quality of the dried sunflower seeds: acid, peroxide, and iodine numbers decreased by 12, 24, and 9% respectively.

If we compare the data obtained from the practical tests and from the model (Table 1), we can conclude that the results are reproducible. The following optimal va- lues were also obtained while solving the problem of convective drying optimization: the initial moisture con- tent of sunflower seeds was 16–17%; the temperature of

the drying agent in the bed was 66–67%; the consump-

 

coolant, respectively, kg/kg;

K drying ratio, 1/sec;

G – consumption of the dryable product, kg/sec;

  1. (τ), L(z) coolant flow rate, kg/sec;

q – moving mass flow, kg/sec;

ρ – density of the dryable product, kg/m3;

u material speed, m/sec;

S – section area of the fluidized bed, m2;

a, b constants determined experimentally;

u

 

e

equilibrium moisture content of the material, kg/kg;

T – product temperature, °C;

α – heat transfer coefficient, kW/(m2×K);

t coolant temperature, °C;

S

s

– specific surface area of the material, m2/kg;

r   specific heat of vaporization, J/kg;

 

tion of drying agent was (3.2–3.4)×10–2  m3/sec; the angle of the drum was 0.61–0.70 rad; the drum rotation fre-

 

m

 

w

c , c

 

 

– specific heat capacity of dry material and water,

 

kJ/(kg×K);

 

quency was 3.6–4.2 min. These results agreed with the

data of the model presented in Table 2.

 

CONCLUSION

The proposed multimodel system of non-stationary drying processes for bulk materials has a number of ad- vantages. First, it leads to a block-modular construction and expedient aggregation of rotary drum dryers. Seco-

 

surface of the particle, m2; V  volume of the particle, m3;  ε – porosity of the fluidized bed;

p

 

s

S

p

– specific volume of the coolant, m3/kg;

x(z) moisture content of the coolant, kg/kg;

f

V – the volume of the fluidized bed of the bulk material

in the rotating drum, m3;

ρ(u ), ρ(u ), ρ(u ) – product density corresponding to the

 

e               i               f

 

nd, it optimizes the allowances on the inputs and out-

puts of technological operations and links them toge-

 

equilibrium, initial, and final  moisture content of the

material, kg/m3;

 

ther. Third, it develops requirements for the quality of

 

 

Nu

 
ORe

,  Nu

 

 

OAr

 

–  Nusselt  numbers  for  forced  and  free

 

raw materials and environmental conditions, in terms of

the high efficiency of the organization of its processing. Thus, when studying the specifics of heat transfer be-

tween the coolant and the solid particles of bulk mate-

 

movement forms, respectively;

c, m, n constants of the equation; n drum speed per minute, min-1; φ – drum angle, rad;

 

 

 

ψ – the angle between the surface of the bed and the axis

of the drum, rad;

Θ – friction angle of the material, rad;

R – the radius of the channel nozzle of the drum, m;

 

the outlet from the drying chamber and the temperature of the outside air, °C.

C unit exergy value;

i

 
e energy of the consumed raw materials and exergy;

 

r

 

0

– the radius of the circle touching the bed at the input, m;

A – coefficient that takes into account the effect of the average flow rate in the bed of the wet solids on the throughput of the dryer, A = f(Re) is determined exper- imentally;

ξ – tortuosity coefficient of the channel;

w flow rate in the direction of filtration, m/sec;


c  unit cost of exergy of the raw materials and energy;

ei

nk

e   – exergy of the products;

n

k   – capital and other associated expenditures for the

n-subsystem.

 

CONFLICT OF INTEREST

The authors declare that there are no conflicts of in- terest related to this article.

 

d

 

p

ρ , ρ

 

bulk density and particle density, kg/m3;

 

λ

 

g

– the resistance coefficient of the granular bed;

w

w

flow rate in the direction of filtration, m/sec;

h thickness of the bed, m;

ρ – density of the drying agent, kg/m3;


FUNDING

The research was conducted by the authors as a part of their work at the Voronezh State University of Engi- neering  Technologies  and  ZRM  Bobrovsky  vegetable

 

t , t , t

 

temperature of the drying agent at the inlet and

 

oil plant. 

References

1. Antipov S.T., Panfilov V.A., Urakov O.A., and Shakhov S.V. Sistemnoe razvitie tekhniki pishchevykh tekhnologiy[Systemic development of food technology]. Moscow: KolosS Publ., 2010. 762 p. (In Russ.).

2. Zueva G.A., Kokurina G.N., Padokhin V.A., and Zuev N.A. Issledovanie teplomassobmena v protsesse konvektivnoy sushki voloknistykh materialov [A research on heat and mass transfer in the process of convective drying of fibrous materials]. Russian Journal of Chemistry and Chemical Technology, 2010, vol. 53, no. 7, pp. 93-96. (In Russ.).

3. Akulich P.V., Temruk A.V., and Akulich A.V. Modeling and experimental investigation of the heat and moisture transfer in the process of microwave-convective drying of vegetable materials. Journal of Engineering Physics and Thermophysics, 2012, vol. 85, no. 5, pp. 951-958. (In Russ.).

4. Frolov V.F. Macrokinetic analysis of the drying of particulate materials. Theoretical Foundations of Chemical Engi- neering, 2004, vol. 38, no. 2, pp. 133-139. (In Russ.).

5. Bon J. and Kudra T. Enthalpy - Driven Optimization of Intermittent Drying. Drying Technology, 2007, vol. 25, no. 4, pp. 523-532. DOI: https://doi.org/10.1080/07373930701226880.

6. Vaquiro H.A., Clemente G., Garcia-Perez J.V., Mulet A., and Bond J. Enthalpy-driven optimization of intermit- tent drying of Mangifera indica L. Chemical Engineering Research and Design, 2009, vol. 87, no. 7, pp. 885-898. DOI: https://doi.org/10.1016/j.cherd.2008.12.002.

7. Glouannec P., Salagnac P., Guézenoc H., and Allanic N. Experimental study of infrared-convective drying of hy- drous ferrous sulphate. Powder Technology, 2008, vol. 187, no. 3, pp. 280-288. DOI: https://doi.org/10.1016/j. powtec.2008.03.007.

8. Antipov S.T., Zhuravlev A.V., Kazartsev D.A., et al. Innovatsionnoe razvitie tekhniki pishchevykh tekhnologiy [Inno- vative development of food technology techniques]. St. Petersburg: Lan Publ., 2016. 660 p. (In Russ.).

9. Sazhin B.S., Otrubjannikov E.V., Kochetov L.M., and Sazhin V.B. Drying in active hydrodynamic regimes. Theoreti- cal Foundations of Chemical Engineering, 2008, vol. 42, no. 6, pp. 638-653. (In Russ.).

10. Padokhin V.A., Zueva G.A., Kokurina G.N., Kochkina N.E., and Fedosov S.V. Complex Mathematical Description of Heat and Mass Transfer in the Drying of an Infinite Cylindrical Body with Analytical Methods of Heat-Conduction Theory. Theoretical Foundations of Chemical Engineering, 2015, vol. 49, no. 1, pp. 54. (In Russ.)

11. Bobkov V.I. Issledovanie tekhnologicheskikh i teplo-massoobmennykh protsessov v plotnom sloe dispersnogo ma- teriala [A research on technological and heat and mass transfer processes in a dense layer of dispersed material]. Thermal Processes in Engineering, 2014, no. 3, pp. 139-144. (In Russ.).

12. Ol’shanskii A.I. Study of the heat transfer in the process of drying of moist materials from experimental data on moisture transfer. Journal of Engineering Physics and Thermophysics, 2014, vol. 87, no. 4, pp. 887-897. (In Russ.).

13. Ol’shanskii A.I. Heat transfer kinetics and experimental methods for calculating the material temperature in the drying process. Journal of Engineering Physics and Thermophysics, 2013, vol. 86, no. 3, pp. 584-594. (In Russ.).

14. Chin S.K. and Law C.L. Product quality and drying characteristics of intermittent heat pump drying of Ganoderma tsugae Murrill. Drying Technology, 2010, vol. 28, no. 12, pp. 1457-1465. DOI: https://doi.org/10.1080/07373937.20 10.482707.

15. Sokolowskyy Ya., Dendiuk M., and Bakaletz A. Mathematical modeling of the two-dimensional moistural and visco- elasticity states of wood in the process of drying. IAWS plenary meeting and conference “Forest as renewable source of vital values for changing world”. St. Petersburg - Moscow, 2009, p. 120.


Login or Create
* Forgot password?