UNSTEADY HEAT TRANSFER DURING ENCAPSULATION OF DISPERSED MATERIALS IN QUASI-LIQUEFIED STATE

Similar documents
1. Nusselt number and Biot number are computed in a similar manner (=hd/k). What are the differences between them? When and why are each of them used?

Examination Heat Transfer

ECE309 INTRODUCTION TO THERMODYNAMICS & HEAT TRANSFER. 3 August 2004

EXPERIMENTAL AND NUMERICAL STUDIES OF A SPIRAL PLATE HEAT EXCHANGER

Principles of Food and Bioprocess Engineering (FS 231) Exam 2 Part A -- Closed Book (50 points)

INVESTIGATION OF PHYSICAL DISSOLUTION OF BENZOIC ACID POLYDISPERSE MIXTURE

Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition

Chapter 11: Heat Exchangers. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

Time-Dependent Conduction :

Chapter 7: External Forced Convection. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

PHYSICAL MECHANISM OF CONVECTION

PREDICTION OF PHYSICAL PROPERTIES OF FOODS FOR UNIT OPERATIONS

THE EXPERIMENTAL STUDY OF THE EFFECT OF ADDING HIGH-MOLECULAR POLYMERS ON HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS

Forced Convection: Inside Pipe HANNA ILYANI ZULHAIMI

If there is convective heat transfer from outer surface to fluid maintained at T W.

Analysis, Design and Fabrication of Forced Convection Apparatus

Circle one: School of Mechanical Engineering Purdue University ME315 Heat and Mass Transfer. Exam #2. April 3, 2014

Transient Heat Transfer Experiment. ME 331 Introduction to Heat Transfer. June 1 st, 2017

University of Rome Tor Vergata

Theoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder

The Effect of Mass Flow Rate on the Effectiveness of Plate Heat Exchanger

External Forced Convection :

TankExampleNov2016. Table of contents. Layout

HEAT EXCHANGER. Objectives

Heat processes. Heat exchange

A population balance approach for continuous fluidized bed dryers

Keywords: Spiral plate heat exchanger, Heat transfer, Nusselt number

DESIGN OF A SHELL AND TUBE HEAT EXCHANGER

Theoretical and experimental studies on a latent heat thermal energy storage system (LHTES) containing flat slabs of phase change materials

Heat and Mass Transfer Unit-1 Conduction

Parallel Plate Heat Exchanger

FIND: (a) Sketch temperature distribution, T(x,t), (b) Sketch the heat flux at the outer surface, q L,t as a function of time.

Examination Heat Transfer

Energy. The influence of plate corrugations geometry on plate heat exchanger performance in specified process conditions

Heat Transfer Convection

Exergy Analysis of Solar Air Collector Having W Shaped Artificial Roughness

QUESTION ANSWER. . e. Fourier number:

8.1 Technically Feasible Design of a Heat Exchanger

Gas Technology Institute UTSR Fellow: Miad Karimi. > Summer 2017

Phone: , For Educational Use. SOFTbank E-Book Center, Tehran. Fundamentals of Heat Transfer. René Reyes Mazzoco

Heat Transfer Coefficient Solver for a Triple Concentric-tube Heat Exchanger in Transition Regime

Finite-Element Evaluation of Thermal Response Tests Performed on U-Tube Borehole Heat Exchangers

Transactions on Engineering Sciences vol 18, 1998 WIT Press, ISSN

ME 331 Homework Assignment #6

Simultaneous heat and mass transfer studies in drying ammonium chloride in a batch-fluidized bed dryer

Investigations Of Heat Transfer And Components Efficiencies In Two-Phase Isobutane Injector

Heat Exchangers for Condensation and Evaporation Applications Operating in a Low Pressure Atmosphere

The Modified Analogy of Heat and Momentum Transfers for Turbulent Flows in Channels of Plate Heat Exchangers

PROBLEM and from Eq. 3.28, The convection coefficients can be estimated from appropriate correlations. Continued...

Droplet behaviour in a Ranque-Hilsch vortex tube

Name: ME 315: Heat and Mass Transfer Spring 2008 EXAM 2 Tuesday, 18 March :00 to 8:00 PM

Entropy Generation Analysis for Various Cross-sectional Ducts in Fully Developed Laminar Convection with Constant Wall Heat Flux

CHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION

Introduction to Heat and Mass Transfer. Week 12

FORCED CONVECTION AND EFFECTIVE TRANSPORT PROPERTIES IN PACKED BEDS

Theoretical Study of Forced Convective Heat Transfer in Hexagonal Configuration with 7Rod Bundles Using Zirconia-water Nanofluid

Temperature distribution and heat flow across the combustion chamber wall.

ENTHALPY BALANCES WITH CHEMICAL REACTION

Experiment 1. Measurement of Thermal Conductivity of a Metal (Brass) Bar

INTERNATIONAL JOURNAL OF MECHANICAL ENGINEERING AND TECHNOLOGY (IJMET)

Calculation of Power and Flow Capacity of Rotor / Stator Devices in VisiMix RSD Program.

WTS Table of contents. Layout

Fluid Flow and Heat Transfer Characteristics in Helical Tubes Cooperating with Spiral Corrugation

Heat Transfer Analysis of Centric Borehole Heat Exchanger with Different Backfill Materials

123MEAN thermal properties KATEDRA MATERIÁLOVÉHO INŽENÝRSTVÍ A CHEMIE

Axial profiles of heat transfer coefficients in a liquid film evaporator

arxiv: v1 [physics.app-ph] 25 Mar 2018

Investigations of hot water temperature changes at the pipe outflow

WITPRESS WIT Press publishes leading books in Science and Technology. Visit our website for the current list of titles.

Heat Transfer Predictions for Carbon Dioxide in Boiling Through Fundamental Modelling Implementing a Combination of Nusselt Number Correlations

Scale-up problems are often perceived as difficult. Here the reaction calorimetry has proven to be

Forced Convection Around Obstacles

The Two-Phase Mathematical Model of. Dehydration and Granulation in a Fluidized Bed

CHAPTER 4 THERMAL CONDUCTIVITY AND VISCOSITY MEASUREMENTS

Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer

Transport processes. 7. Semester Chemical Engineering Civil Engineering

Experimental and Theoretical Investigation of Hydrodynamics Characteristics and Heat Transfer for Newtonian and Non-newtonian Fluids

S.E. (Chemical) (Second Semester) EXAMINATION, 2011 HEAT TRANSFER (2008 PATTERN) Time : Three Hours Maximum Marks : 100

The Analytical Solution of the Water-Rock Heat Transfer Coefficient and Sensitivity Analyses of Parameters

Countercurrent heat exchanger

Determination of heat losses of a concrete silo for sugar and a fan project

ELEC9712 High Voltage Systems. 1.2 Heat transfer from electrical equipment

Experimental Investigation of Single-Phase Friction Factor and Heat Transfer inside the Horizontal Internally Micro-Fin Tubes.

PHYSICAL AND CHEMICAL ANALYSIS OF FORMALDEHYDE ABSORPTION PROCESS

Level 7 Post Graduate Diploma in Engineering Heat and mass transfer

Convective Mass Transfer

An alternative turbulent heat flux modelling for gas turbine cooling application

Introduction to Heat and Mass Transfer. Week 14

MECHANISM BEHIND FREE /NATURAL CONVECTION

Takaki OHKOUCHI Hiroyuki OSAKABE Toshihide NINAGAWA Kiyoshi KAWAGUCHI

طراحی مبدل های حرارتی مهدي کریمی ترم بهار HEAT TRANSFER CALCULATIONS

Specific heat capacity. Convective heat transfer coefficient. Thermal diffusivity. Lc ft, m Characteristic length (r for cylinder or sphere; for slab)

Internal Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.

Chapter 7: External Forced Convection

Mixing Process of Binary Polymer Particles in Different Type of Mixers

PROBLEM 7.2 1/3. (b) The local convection coefficient, Eq. 7.23, and heat flux at x = L are 1/2 1/3

Basic Fluid Mechanics

Investigation of Heat Transfer on Smooth and Enhanced Tube in Heat Exchanger

PHY Final Exam Problem 1 (20 points)

CFD Study of the Turbulent Forced Convective Heat Transfer of Non-Newtonian Nanofluid

Transcription:

CHEMISTRY & CHEMICAL TECHNOLOGY Vol. 9, No. 4, 05 Chemical Technology Oleg Nagursky, Jaroslaw Gumnitsky and Victoria Vaschuk UNSTEADY HEAT TRANSFER DURING ENCAPSULATION OF DISPERSED MATERIALS IN QUASI-LIQUEFIED STATE Lviv National Polytechnic University, Bandera str., 7903 Lviv, Ukraine; nahurskyy@mail.ru Received: July 04, 04 / Revised: August 8, 04 / Accepted: March 8, 05 Nagursky O., Gumnitsky Ya., Vaschuk V., 05 Abstract. Experimental and analytical investigations of the heat transfer process during encapsulation of dispersed materials in quasi-liquefied state are presented. The heattransfer coefficients for different types of materials have been determined during their heating depending on air rate. Keywords: unsteady heat transfer, kinetics, quasiliquefaction, dispersed material.. Introduction During encapsulation of dispersed materials in quasi-liquefied state using film-forming solutions it is necessary to heat the particles till operating temperature. In the batch apparatus such heating is a separate technological stage []. In the continuous apparatus the heating may proceed in the separate area or directly in the area of coating growing in parallel with film-forming agent plating depending on design []. Irrespective of apparatus type in the certain periods of time the dispersed material is in the state of unsteady heat transfer. The intensity of this stage depends on heat carrier rate and temperature, as well as on material physical properties [3, 4]. Therefore, to produce and use the encapsulated materials the investigations of heat transfer kinetics during encapsulation of dispersed materials in pseudo-liquefied state are in demand. It will grant the possibility to employ a scientific approach to the choice of encapsulation technological parameters, provide the high productivity of the equipment and obtain the materials with predicted properties.. Experimental Heat transfer between the air and material layer during heating occurs at the moment of dispersed material loading into an apparatus and ends as particles are heated till the operating temperature. The experiments were carried out in the batch cylindrical apparatus of quasiliquefied state. During the experiments the working part of the apparatus was insulated to avoid heat losses throughout the walls. Before the material loading the apparatus was heated to the operating temperature, hence there was no heat exchange between the air and apparatus walls. At the beginning the dispersed material temperature was 93 K. The temperature of heat carrier was 348 K and it was measured using chromel-copel thermocouple and 7-channels smart convertor PVI-098 with computer recording. The change of heat carrier temperature was fixed from the moment of loading and till approaching of the air temperature at the apparatus outlet to the initial values. The used solid materials were: polydispersed layer granulated mineral fertilizers (ammonium nitrate, calcium nitrate, carbamide and nitroamophose) and layer of particles with irregular geometry seeds which may undergo the presowing encapsulation by chemical protectants for plants and chemical elements of additional fertilizing (fodder beet, spinach). 3. Results and Discussion The dependencies of heat carrier temperature on time at different values of liquefying air rates within the range between the first and the ssecond critical values are represented in Fig.. The analysis of the experimental results shows that the increase in air rate decreases the time necessary for the layer to be in temperature equilibrium with the heat carrier. The reason is that the higher air rate decreases the layer thickness at the boundary air particle surface and increases the amount of heat transferred from the particle at the same space of time. These observations are in an agreement with the results of heat-and-mass transfer in the quasi-liquefied state obtained by other authors [3, 7, 8].

498 Oleg Nagursky et al. ammonium nitrate nitroamophose calcium nitrate carbamide spinach seeds fodder beet seeds Fig.. Liquefying air temperatures at its different rates (w, m/s) in the layer during the heating of dispersed material in quasi-liquefied state

Unsteady Heat Transfer during Encapsulation of Dispersed Materials in Quasi-liquefied State 499 The heat-transfer coefficient α which determines the amount of heat donated or accepted by the surface unit for time unit is an important thermal characteristic of the heat transfer processes. The coefficient is determined experimentally and the theory of generalized variables is used for its generalization. According to the theory the criterial dependencies between dimensionless numbers are obtained. The coefficient α is determined on the basis of heat transfer equation [3] which includes the surface temperature of the solid matter. It is problematic to measure the surface temperature of the particle in quasiliquefied state and therefore to determine the coefficient α we used the method based on the theory of unsteady heat conductivity. On the basis of heat balance equations (Eq. ()) we determined the average temperature of the solid T = f τ. particles as a time function ( ) ( 0 ) ( ) Q= M C T T τ = Vρct t τ () p p in out where Q the heat transferred from the liquefying air to the particles, J/kg K; M p particles mass, kg; C p material heat capacity, J/kg K; V, ρ and c consumption (m 3 /s), density (kg/m 3 ) and heat capacity (J/kg K) of the liquefying air, respectively; t in and t out liquefying air temperature at the inlet and outlet of the dispersed material layer, respectively, K; τ time, s. The solution of the equation of unsteady heat conductivity with boundary conditions of the third order for the average temperature of the spherical solid particle in dimensionless units gives Eq. () [5]: T T0 Θ= = Bnexp( µ nfo) () ta T0 where Θ the function of Fo and Bi numbers; µ n roots of the characteristic equation; Fo = a τ Fourier R number; Bi = α R λ Bio number; t a temperature of the liquefying air in the material layer, K; T and T 0 average and initial temperatures of the material particles, respectively, K. The amount of heat Q transferred from the liquefying air to the material particles is calculated according to the heat balance equation for the heat carrier, t a is accepted as averaging temperature of the heat agent in the material layer for investigated space of time according to the measurements represented in Fig.. At the considerable space of time we may limit ourselves in Eq. () by the first term of series. After some transformations we obtain Eq. (3): ta T Θ= = B exp( µ n Fo) (3) t T a 0 where B = 6Bi ( + Bi Bi) µ µ To determine the root of the characteristic equation µ we find the logarithm of Eq. (3): aτ lnθ= lnb µ Fo= ln B µ (5) R Then we build the plot ln(θ) = f(τ), where a sloping angle of the corresponding kinetic curves is equal to: a λ tgβ = µ = µ (6) R с ρ R where λ the heat conductivity factor of the liquefying air, W/m K; c heat capacity of the liquefying air, J/kg K; ρ density of the liquefying air, kg/m 3. Using the value of Bi number, calculated according to Eq. (4) we calculate the values of the coefficient α. The obtained results are averaged relative to the layer height and represented in Fig.. One can see from Fig. that the coefficient α increases with the increase in air rate according to the linear law. It is explained by greater heat application with the increase of liquefying air amount and decrease of boundary heat layer around the particle due to the gas flow turbulization. The obtained results (Fig. ) are in the agreement with the data of other authors [6]. The experimental results are generalized in accordance with Eq. (7) [5]: n m Nue = AReePr (7) Taking into account that physical parameters of the air are changed within the narrow range, we assume Nu e ~ Pr 0.33 [3, 7]. To determine the coefficients A and n in Eq. (7) the experimental data were described as the dependence Nu e /Pr 0.33 = f(re e ) in the logarithmic coordinates (Fig. 3). One can see from Fig. 3 that the experimental points for all materials are approximated by practically parallel lines. Thus, the exponents near Reynolds number are equal for different materials. The difference in vertical placement of the lines is explained by the influence of particles size on the heat-transfer coefficient. The same as for the heat-transfer coefficient, the dependence of Nusselt number on the rate of heat agent flowed around the particles is correlated with the data from Ref. [8, 9] for different by size balls made from glass and plastics. The determined coefficients A and n are given in the Table. The value of the coefficient A is a function of particle size and exponent n is the same for all materials and equals approximately to 0.9. To take the effect of geometrical sizes of particles d and apparatus D a into consideration it is necessary to introduce the geometric simplex d/d a into Eq. (7). (4)

500 Oleg Nagursky et al. Fig.. Dependence of the heat-transfer coefficients on dummy rate of the liquefying air for the dispersed materials heating Fig. 3. Generalization of the experimental results concerning the heat-transfer coefficients during the dispersed material heating by liquefying air The values of A and n in Eq. (7) for a determination Material А n Calcium nitrate 0.06 0.9 Nitroamophose 0.07 0.9 Ammonium nitrate 0.00 0.9 Carbamide 0.09 0.9 Spinach seeds 0.053 0.9 Fodder beet seeds 0.059 0.9 Table The value A may be approximated by exponential Eq. (8) depending on the geometric simplex. d A= A (8) D a The experimental data are approximated in the best way by two curves (Fig. 4). The first curve is drawn k through the points that correspond to the coefficient A obtained for the seeds and model particles, the second for granulated fertilizers. The points on the graph are situated from the left to the right in ascending order of particles size. The exponent in Eq. (8) is the same for both curves and equals to 0.67. A for the first curve is 0.48; for the second 0.53.

Unsteady Heat Transfer during Encapsulation of Dispersed Materials in Quasi-liquefied State 50 essential influence on the streamline conditions of heat agent. Therefore, for the seeds the difference between theoretically calculated Nusselt number and experimental data is approximately 6 %. References Fig. 4. Coefficient A vs. geometric simplex d/d a for dry materials: seeds () and granulated fertilizers () Thus, Eq. (7) takes a view of Eqs. (9) and (0) for calcium nitrate, nitroamophose, carbamide, ammonium nitrate and model particles: 0.9 0.33 d Nue = 0.005Ree Pr Da for spinach and fodder beet seeds: Nu 4. Conclusions 0.9 0.33 e = 0.07Ree Pr Da d 0.67 0.67 (9) (0) The obtained generalized dependencies allow to determine the heat-transfer coefficients of the dispersed materials heating till the operating temperature. The difference between different dispersed materials may be explained by the different form of particles. Thus, the seeds are characterized by complex geometry having the [] Demchuk I.: PhD thesis, Lviv Polytechnic Nats. Univ., Lviv 99. [] Ovchinnikov L. and Lipin A.: Kapsulirovanie Mineralnykh Udobrenij vo Vzveshenom Sloe. Ivan. GKhTU, Ivanovo 0. [3] Tovazhnyansky L., Hotlinska A. and Leschenko V.: Procesy i Aparaty Khimichnoi Technologii. Nats. Univ. KhPI, Kharkiv 007. [4] Lipin A.: Intensyfikatsia Teplowykh i Massoobmennykh Procesov v Heterohennykh Sredakh. Ivan.GKhTU, Ivanovo 009. [5] Lykov A.: Teoria Teploprovodnosti. Vysh.shkola, Moskwa 967. [6] Gelperin N.: Osnovy Techniki Psevdoozyzeniya. Khimiya, Moskwa 964. [7] Atamanyuk V. and Gumnitskyy Ya.: Naukovi Osnovy Filtraciynogo Sushinnya Dispersnykh Materialiv. Vyd-vo Lviv. Politechn., Lviv 03. [8] Kunii D.: Promyshlennoe Psevdoozyzenie. Khimiya, Moskwa 976. [9] Flisjuk E.: D.Sc. thesis, St. Petersburg Gos. Khimpharm. Acad., St.Petersburg 006. НЕСТАЦІОНАРНИЙ ТЕПЛООБМІН ПІД ЧАС КАПСУЛЮВАННЯ ДИСПЕРСНИХ МАТЕРІАЛІВ У СТАНІ ПСЕВДОЗРІДЖЕННЯ Анотація. Наведено експериментально-аналітичні дослідження процесу теплообміну під час капсулювання дисперсних матеріалів у стані псевдозрідження. Визначено коефіцієнти тепловіддачі під час нагрівання дисперсного матеріалу у стані псевдозрідження у залежності від швидкості повітря для різних типів матеріалів. Ключові словa: нестаціонарний теплообмін, кінетика, псевдозрідження, дисперсний матеріал.