Научная статья на тему 'Исследование процессов формирования газо-паровой фазы в жидких смесях'

Исследование процессов формирования газо-паровой фазы в жидких смесях Текст научной статьи по специальности «Физика»

CC BY
38
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
VAPORIZATION / MOVEMENT OF PORE'S BOUNDARY / PORE'S HEAT TRANSFER / CONDITIONS OF CONTROLLED SWELLING

Аннотация научной статьи по физике, автор научной работы — Pavlenko А., Koshlak H., Cheilytko А., Nosov M.

Equilibrium conditions of the pore former agent in the material during the formation of the porous structure, pressure and temperature of the pore former agent gas inside the closed spherical pore, heat transfer between the growing pore and the surrounding mixture were researched. The calculation of the gas vapor area was made with obtained equations. This calculation needs for prediction of pore growth dynamics. Obtained dependences show that, in general, the peak values of the growth speed of vapor volume, movement of the pore boundary, mass flow, heat flow have smaller values under lower periods of oscillations. The character of changing of the calculation quantities under 130 oC and 180 oC are the same. This temperature range was chosen, because real swelling happens under the same range. Obtained equation of the overheat temperature of the pore former agent gas inside the pore clarifies the temperature (180 oC), under which the inertial period of formation of the gas microphase in the first heating stage is missing. The equilibrium conditions give a chance to assess energy parameters of the swelling process under stabilization of the predictable pore sizes.During the calculation of the pressure inside the closed spherical pore it was found that the bigger the difference between the chemical potentials of material-pore systems, the lower the gas pressure inside the pore. Since the convective heat transfer in a gas is directly proportional to the pressure, next statement can be made: to achieve minimum heat transfer of pore, it’s necessary to increase the difference between the chemical potentials of material-pore systems. Obtained methodology allows finding conditions of controlled swelling and conditions of controlled structure formation of the material with predicted thermophysical properties. It can be real only, because this methodology takes into account physical properties of the raw mixture, the chemical potential of the mixture components, levels of energy influence on the raw mixture and the impact of all above factors on the size of the gas-vapor area (pore). The differences of new methodology allow predicting the porosity of thermal insulatingmaterial and its thermophysical properties.These results are proposed to use in designing technological processes of production of porous materials for various purposes.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Study of the formation of gas-vapor in the liquid mixture

Equilibrium conditions of the pore former agent in the material during the formation of the porous structure, pressure and temperature of the pore former agent gas inside the closed spherical pore, heat transfer between the growing pore and the surrounding mixture were researched. The calculation of the gas vapor area was made with obtained equations. This calculation needs for prediction of pore growth dynamics. Obtained dependences show that, in general, the peak values of the growth speed of vapor volume, movement of the pore boundary, mass flow, heat flow have smaller values under lower periods of oscillations. The character of changing of the calculation quantities under 130 oC and 180 oC are the same. This temperature range was chosen, because real swelling happens under the same range. Obtained equation of the overheat temperature of the pore former agent gas inside the pore clarifies the temperature (180 oC), under which the inertial period of formation of the gas microphase in the first heating stage is missing. The equilibrium conditions give a chance to assess energy parameters of the swelling process under stabilization of the predictable pore sizes.During the calculation of the pressure inside the closed spherical pore it was found that the bigger the difference between the chemical potentials of material-pore systems, the lower the gas pressure inside the pore. Since the convective heat transfer in a gas is directly proportional to the pressure, next statement can be made: to achieve minimum heat transfer of pore, it’s necessary to increase the difference between the chemical potentials of material-pore systems. Obtained methodology allows finding conditions of controlled swelling and conditions of controlled structure formation of the material with predicted thermophysical properties. It can be real only, because this methodology takes into account physical properties of the raw mixture, the chemical potential of the mixture components, levels of energy influence on the raw mixture and the impact of all above factors on the size of the gas-vapor area (pore). The differences of new methodology allow predicting the porosity of thermal insulatingmaterial and its thermophysical properties.These results are proposed to use in designing technological processes of production of porous materials for various purposes.

Текст научной работы на тему «Исследование процессов формирования газо-паровой фазы в жидких смесях»

Дослиджено процеси формування газо-паровог фази в ргдких сумшах, як перебувають в стан термодинамiчног рiвноваги. Визначено термоди-намiчш параметри газо-паровог фази, на основi яких запропонована методика управлтня ттен-сивтстю фазового переходу в ргдких сумшах зi змтними реологiчними властивостями, зокре-ма в'язтстю i щтьтстю. Отримано матема-тичш залежностi, що використаш для анал^ зу та управлтня процесами спучування ргдких гелеоподШних сумшей в технологи виробництва пористих матерiалiв

Ключовi слова: пароутворення, реологiя редких сумшей, термодинамiчш параметри паро-

вог областi

□-□

Исследованы процессы формирования газо-па-ровой фазы в жидких смесях, находящихся в состоянии термодинамического равновесия. Определены термодинамические параметры газо-паровой фазы, на основе которых предложена методика управления интенсивностью фазового перехода в смесях взаимно не растворимых жидкостей с изменяющимися реологическими свойствами, в частности вязкостью и плотностью. Полученные математические зависимости реализованы для анализа и управления процессами вспучивания жидких гелеобразных смесей в технологии производства пористых материалов Ключевые слова: парообразование, реология жидких смесей, термодинамические параметры паровой области

UDC 536.775, 539.217.1

|DOI: 10.15587/1729-4061.2016J5428|

STUDY OF THE FORMATION OF GAS-VAPOR IN THE LIQUID

MIXTURE

А. Pavlenko

Doctor of Technical Sciences, Professor Department of Building Physics and Renewable Energy Kielce University of Technology Al. Tysi^clecia Panstwa Polskiego, 7, Kielce, Poland, 25-314 Е-mail: am.pavlenko@i.ua H. Koshlak PhD, Associate Professor Department of heat and gas supply, ventilation and heat power engineering Poltava national technical Yuri Kondratyuk university Pershotravnevyi ave., 24, Poltava, Ukraine, 36011

A. Cheilytko PhD, Associate Professor* Е-mail: cheylitko@yandex.ua M. Nosov* Е-mail: nosovmak@gmail.com *Department of heat and power engineering Zaporizhzhya state engineering academy Soborny ave., 226, Zaporizhzhya, Ukraine, 69006

1. Introduction

Significant part of the modern technological processes of structure formation of materials are based on the thermal impact on various raw mixtures, which implement in them physical and chemical processes that create conditions for structure formation. Such processes can be controlled, but it is necessary to find qualitative and quantitative dependencies of the changing thermal properties of materials in the process of their formation.

Investigation and control of the processes of structure formation of materials are difficult tasks, which are still unsolved. Clear understanding of the mechanism of structure formation makes it possible to develop a methodological basis of new technologies, including the technology of production of thermal insulating materials with predictable thermal properties.

However, the ways of formation of the porous structure are still little investigated, and a clear relationship between the porosity and the physical properties of the material were not found. A lot of experimental data show the relationship between the porosity of the material and its thermal properties [1-3]. It is obvious that the structure of the materials, in particular the porosity, determines their properties. But in the mentioned works this impact is shown in different ways. For example, in the [3] the thermal conductivity of

Fe (58.19 W/(m-K)) and of clay (3.26 W/(m-K)) differ by 18, but the thermal conductivity of insulation structures, which are made from granules of Fe and granules of clay with the same porosity are almost equal - 0.0403 W/(m-K) and 0.0402 W/(m-K).

Such results show that not only porosity, but also size and shape of the pores affect the properties of the material.

Since any material has an own characteristic distribution of the pore size, it is obvious that various researchers obtained conflicting information about the nature of the influence of pore size on the thermal properties of these materials. Current technologies of structure formation do not provide the prediction of the geometric structure, which means that there is no possibility to predict the properties of materials.

2. Literary analysis and the problem statement

In [1-4] the dependence of the thermal properties of porous materials on a structure was discussed, but recommendations about optimal structure were not given. In [5-7] the dependence of the mechanical properties on a structure of porous materials was analyzed and recommendations about structure formation with predicted properties were given, but there was no information about the thermal properties.

©

In [8] the influence of internal pore pressure on closed-cell elastomeric foams was explored. Changes of internal pore pressure with different hydrostatic loads were considered. Obtained results show that this pressure can significantly change macroscopic reaction and stability of closed-cell elastomeric foams. Also it shows that elastomeric foams with internal pore pressure have a higher stiffness, even with atmospheric pressure, than without it. But the method of calculating the internal pore pressure was not given. Also pressure was taken only as a function of density.

In [9] changes of the structure with a closed porosity under compression and extension, with different initial pressures in the pores, were researched. Experiment results show, that internal pressure has a positive effect under compression and negative under extension. The impact of deformation on the structures with closed porosity can lead to high initial pressure in the pores, which increases the total energy absorption and stiffness of the material under the process of deformation. But nothing was said about the methods of achieving certain pressure in the pores of the material and its calculation.

In [10] the micromechanical analysis of the porous material with internal pressure in the pores was made, the polymer BX-265 was taken as experimental material. The influence of initial pore pressure on the predicted elongation of the sample and the influence of the applied load (with pressure in the pores and without it) on the predicted break of the sample were shown. But nothing was said about the calculation methods of internal pore pressure.

In [11] the dependencies of thermal characteristics from the structure of materials were compiled and the task of the controlled pore formation by adjustable heat treatment of the raw gel-like mixture was formulated.

In [12-14] the main physicochemical formation processes of the gas-vapor area (pores) were analyzed, which were taken as a basis in this work.

All mentioned works have one main idea that materials structure affects their properties. But there are no practical recommendations about formation of specific material structure. Above information cannot be generalized, because it's contradictory. That's why existing technologies of the thermal swelling (structure formation) can't give the required structure.

3. The purpose and objectives of the research

The purpose of this work is the research of regularities of gas-vapor phase formation in liquid mixtures, which are in a state of thermodynamic equilibrium.

To achieve this purpose, next objectives must be solved:

- to simulate the condition of thermodynamic phase equilibrium in the liquid mixture to control the process of pore formation;

- to determine the energy parameters of the swelling process (the pressure of the pore former agent gas inside the closed spherical pore, the conditions of thermal balance);

- to determine the overheating temperature of the liquid mixture to assess the growth dynamics of the gas-vapor area;

- to research the growth dynamics of the gas-vapor phase (pores).

4. Materials and methods of research

In this work, production technologies of porous heat-insulating materials were researched by the method of thermal swelling of the gel-like raw material mixture. When it is heated, inside the raw mass, are formed gas-vapor bubbles are formed, which after crystallization of the mixture form a porous space. Thermal conductivity of the material will depend from the size and shape of the resulting pores.

To obtain material as the active silica material, tuff, rotten stone, silicate - lump were used. The chemical composition of the materials is given in Table 1.

Table 1

The chemical composition of active silica materials

Name of material Chemical composition, %

SiO2 AI2O3 Fe2O3 CaO MgO

Tuff 93 2.8 0.5 0.6 0.1

Rotten stone 83 5.8 4.5 1.4 0.8

Silicate - lump 70 1.4 1.0 0.6 0.5

To obtain the raw mixture, quartz sand with SiO2 -97.7 % was used. The calculation of the amount of sodium hydroxide in the raw mixture is done in terms of Na2O. Clay DN-1 has the following chemical composition, pts. wt.: Al2O3 - 33; Fe2O3 - 1.17. Refractoriness is 1690oC.

The ratio of the raw mixture components for the porous material production is shown in Table 2.

Table 2

The ratio of the raw mixture components

Component pts. wt.

Siliceous rock 100

Mineral filler 1-75

Clay 1-75

Sodium hydroxide 1-15

Sodium bicarbonate 1-8

Water 50-125

Active silica material, mineral filler and clay were ground to the next specific surface areas: silica material - 1000 cm/g, filler - 4000 cm/g, clay - 1000 cm/g. Then ground components, sodium hydroxide and sodium bicarbonate (with appropriate concentrations) were mixed according to the ratio. The resulting mixture was steamed and then granulated to a grains size =20 mm. The granules were dried with air for 10-15 minutes at a temperature of 100-120 °C and then were swelled with heat treatment in a drying drum at a temperature of 500-550 °C for 5-20 minutes.

The obtained material has a low thermal conductivity even at temperatures of 1000-1200 °C (Fig. 1). The thermal conductivity of the material was determined by thermal conductivity meter IT - X - 400. Samples with a cylindrical form (height 5 mm, diameter 15 mm) were put inside the meter and were processed by temperature influence from 200 to 700 °C. In this temperature range, the thermal conductivity of the material was determined according to the standard method, which is described in the device manual.

Fig. 1 shows the experimental dependence of thermal conductivity on the temperature of the material, for technologies of which the main technological stages were simulated.

Ê

S >

3

c

.e -

0,14 0.12 0,1 0,08 0,06 0,04 0,02 0

/

SL-

w(t) =

Fig.

100 200 300 400 500 600 700 800 900 1000 1100 1200

Temperature, °C

1. Dependence of thermal conductivity on the temperature

1 Pp(T) - Pg 1

V 1,5 p 1

pp(t) - Pgt

V 1,5 p

pgR

1,5 pt- C pgR

1,5 p

Pp(T) - P

1,5 p

; P„(T) - p=0.

c -

> 0, 1,5 pt

PgR

Pp(T) - Pg 1,5 p

< 0,(3)

With

(w-a)(w0 + «)_ "mT

5. Results of the research of regularities of phase transitions in liquid mixtures during heat treatment

5. 1 Equilibrium conditions of the pore agent in the material, during the formation of the porous structure

As the dynamic characteristic, which determines the direction of size changes in vapor pore, tension difference was taken, which was caused by the pressure in the vapor area and by the resistance of the boundary surface of the pore.

The equation, which characterizes the dynamics of growth or reduction of the vapor bubble [11]

dw

dT

1,5pw2 + Pg - Pp(T)_ 1,5pw2 Pp(T) - Pg

PgR

pgR

pgR

,(1)

(w + a)(w0-a) final equations will be the next

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

1,5 pt

(w„ -a)e pgR + w„ + a

w„ + a-(w„ -a) e

1,5pt

IgR

1,5 p Pp(T) - pg

w0 - tg

; a =

1,5 pt

_ ,pp(T)-pg.pp(T)-pg>0,

1,5 p

1,5 p

1,5 p 1,5 p

pp(t) - pg pp(t) - pg

w0tg

w0pgR

w = , _ 0 g—-; Pp(T) - Pg = 0.

1,5 pt

Pp(T) - Pg —g <0, 1,5 p

(4)

1,5 pw0t-pgR

where w - speed of growth of the vapor bubble; t - time of bubble growth; p - density; Pp - pressure inside the vapor area; Pg - pressure in the surrounding liquid; T - temperature; R - radius of the vapor bubble.

Increasing, decreasing and stabilization of the bubble sizes can be represented by three cases:

Pp(T) - pg

1,5 p

Pp(T) - pg

1,5p Pp(T) - pg

> 0,

< 0,

(Pp(T) - Pg

^ ln w- V 1,5 p

^Pp(T) - P w + |Pp(T) - pg

V 1,5 p

In the last case, when Pp(T) - Pg = 0, the gas-vapor area (pore) doesn't change in volume. In the technological aspect, relations between energetic parameters, which characterize the predicted pore size (average), were achieved. Therefore, thermophysical parameters also were achieved, including thermal conductivity. If w0 ^ 0, the latter equation loses the physical sense, because for this case w = 0.

The equation for finding the speed of size changing of the pore, can be written as

1,5 p

pp(t) - Pg

1,5 p

= 0, — = -

1,5 pt

arctg

C.

1,5 p

pp(t) - Pg

1,5 pt

1,5 pt

C;

dR dt '

w0pgR

1,5 pw0t-pgR

(5)

C; (2)

When Pp(T) - P > 0 - vapor bubble is increasing, Pp(T)- P < 0 - vapor bubble is decreasing, in case when Pp(T) - P = 0 - size of the vapor bubble is stabilized. As can be seen from the equations, Pp(T) depends on the external heat flow (temperature). Therefore, the process of pore formation can be controlled.

Speed of size changing of the vapor bubble can be found by the next equations

When solving the equation (5), duration of the swelling process can be found. Since the average value of the pore size is one of the main factors, which determine thermophysical properties of the sample, it this method gives a chance to predict the discussed properties. For example, the thermal conductivity of the material from section 4.

This method allows assessing the kinetic parameters of thermal swelling with sufficient accuracy. For chemical pore formation, the method should be supplemented by the energy estimate of the probability of such reaction, which can be assessed by the Gibbs energy equation.

The Gibbs energy of the pores embryo with radius r (under thermal or chemical swelling)

G(1) = 9(„1)M(1) + oQ,

(6)

a

w =

w =

where M(1) - mass of the pore former, which equals 4nr3/u'; Q - outside surface of the pore former, which equals 4nr2;

- chemical potential of the pore former gas (with infinitely large volume of this gas).

The Gibbs energy of the system, which consists of the pores embryo in the material

G = G(1) + G(2),

where G(2) - the Gibbs energy of the material.

Let's consider the equilibrium state of the system. If the external pressure p' equal to the pore former pressure and temperature of the system T' are constant, the equilibrium condition dG=0 will be the next

dU(1) + dU(2) - T'(dS(1) + dS(2)) +

+p'(dV(1) +dV(2) ) + ^!1)dM(1) + ^(2)dM(2) =0. (7)

Considering that

— | = T; ( —] = -p,

3SL ' UvJ p'

for gas bubble according to the Laplace's equation [15]

temperature and corresponding pressure, ^(1) (p(1),T) and 9(2)(p(2),T) will be obtained. Using analogy with the gas bubble in an infinite liquid volume (expansion 9(1) (p(1),t) in a series in p(1) - p(2), and of ^(2) (p(2),T) in a series in p(2) -p(1) [15]) 9(1) (p(1),T) = ^(2)(p(2),T) can be found.

From the resulting equation it follows that if boundary surface is not flat, chemical potentials of the pore former agent and material related to its own pressure and total temperature are equal. Therefore, equilibrium conditions of the pore former agent in the material under formation of the porous structure with surface tension can be written as

T(1) = t(2); p(1

9(1) (p(1),T) = 9(2) (p(2),T).

2o

-p(2) = —; r

(9)

5. 2. Pressure of the pore former agent gas inside a closed spherical pore

After differentiating the equation (8) by pressure p, with T=const and (9^ / 9p)T = u, the next equation will be obtained

r.(1) - -

2o

^ -9(2)-

dr = 0,

or, as dr^0,

9

(2) .

■9l1) =

2qu(1)

1

- r'=-i- p'- —h

find

dU(1) + dU(2) = T(1)dS(1) + +T(2)dS(2) p'- dV(1) - p'dV(2),

where p - external pressure (in this case, pressure of the material on the pore former agent).

Considering that the total mass of the system is constant

dM(2) = -dM(1),

from the equilibrium condition

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

^ dV(1) + (^1) -9(2) )dM(1) =0.

For spherical pore former agent V(1) = 4nr2dr;

dM(1) = ^dr.

That's why, the latter ratio can be written as 2qu(1) ^

P

As

(2)

2u(1)q(9rx

dp

2q r

du(1) dp

(10)

P

~(2)

is much smaller than u( ), and u( ) at relatively

low gas pressures can be equal to RT/p, then __^u(1) N

v9p j t

RT 2u(1)q(9r

(1)

2q r

dp

After integrating this equation by p>0

RT

J=J

2u(1)q (

dp

dp- J^

J r

dp

dp.

RT [ln (p(1) + c)] = -2q-

(1)

where c - integration constant. After transformations

2a uM

(8)

p(1) = e-c ■ e-RT r .

The constant e-c is found from equilibrium conditions, provided that when r , p = p0, i. e. the pressure of the pore former agent equals to pressure of gas under the same physical conditions, outside of the material

2a uW

p(1) = pc ■ e-RT r . (11)

From this equation it follows that in equilibrium the pressure of the gas pore former gas inside the pore is smaller than the pressure of this gas outside of the material by

2qu<1>

eRT r times.

If instead ^ and ^(2), take real values of the chemical potentials of the pore former agent and material, with the

Obtained equation also can be written as

E

__a_

p(1) = pc ■ e-RT,

(12)

r

where Ea = 2a r or like in (1), Ea = ^(2) -^(„1) - activation energy, which equals to the difference between the chemical potentials of the material and the pore former agent.

Equation (6) shows that the larger the difference between the chemical potentials of material-pore systems, the lower the pressure of gas inside the pores. Since the convec-tive heat transfer in a gas is directly proportional to the pressure, next statement can be made: to achieve minimum heat transfer of the pore, it's necessary to increase the difference between the chemical potentials of material-pore systems.

5. 3. Temperature of the pore former agent gas inside a closed spherical pore

From [14] it's known that the temperature of the phase transition of the spherical drop in vapor is higher than the temperature of the phase transition from the surface. This means that it is necessary to overheat the liquid inside the pores, so that it turned into vapor, at the temperature

AT =

2cu" T Lr ,

(13)

where u'' - volume of saturated vapor; L - heat of vaporization; q - coefficient of surface tension at the temperature TS (i. e. before the curvature of the surface was changed).

Equation (13) applies to vapor bubble, which is in equilibrium with the environment. It can be used to determine the degree of liquid overheat, which does not contain dissolved gases and which is boiling under constant external pressure on the liquid.

The temperature of the pore former gas, which is in equilibrium, can be found by the temperature of the material. After differentiating the equation (9) by T, with (9^/9T)p = -s and (9^ / dp)T = u, the equation will be

9f

9p

dp_ s(2) -s1

dp

dT .

(2) _ s(1)

9^'

dT

9^'

(2)

dp

dp .dT

9^'

(2)

9T

dT v(2) _ v(1)'

dT =

__. dp(2) _ v(1) ■ dp(1)

T=

1

p(2) i i

P^ _ RTln p(1) _(2) 1 1

const.

(14)

After differentiating the equation (7) by temperature T, with (9^/9T) = -s and P=const, the equation will be

(1)

_ s(2) = _ fill +

9v(1) dT

At equilibrium conditions, the embryo radius of the pore former and its volume do not depend on the temperature, so if pressure is constant

2ov

7F.

With the equation (14) (when r the overheat

temperature of the pore former agent gas inside the pore will be

AT = -

1

Rln p(1)

2ov(1)

_(2)

— T

L0'

(15)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

where T0 - temperature of the pore former gas under the same physical conditions outside the pore.

5. 4. Heat transfer between the growing pore and the surrounding mixture

The equation of heat balance for a single vapor bubble

(T0-Ts M (f+1 )=-j(Ts )L (Ts)-q(Ts), (16)

where 8 - current thickness of the thermal boundary layer; j - density of the material flux; L - heat of condensation; 1(- coefficient of thermal conductivity of the liquid; q - heat flux density.

The equation for the pore inside the flux of mixture, will be

21

S r

dr

AT-l =_po eRT ~r

dt

-qp=,

(17)

where qP=const - heat flux density in the isobaric stage; t' -time of the first stage of pore growth.

In the main stage of pore growth (isobaric) heat, which is supplied to the growing pore from the surrounding mixture, will be

Nu ~ Ren Prm.

The average in time, pore radius will be - 1

17 J r(t)dt,

(18)

(19)

and with the dependence of time from pore radius in the isobaric stage of pore growth

Since the density of the material is almost independent of the pressure, at low gas pressures v = RT/p, so

r= t2 . Ip _ po

1,83 t_t' y _ . The Nusselt number

Nu = -° -r

l-AT

The Reynolds number

r dr

Re =---—.

v dt

(20)

(21)

(22)

If substitute these equations in the heat transfer equation, we will get

\n-1

(r)n

qp=const = constI ■ Prm

l-AT vn ^dt^

With average pore radius, will be

(23)

u

Ea

or

qp=const =const -at x

/ 2 |-X»-1

__p - pg

t-t' V p

dr ^n (^-cp dt,

(24)

Final heat balance equation of growing pore in the initial mixture with isothermal stage

AT-l,if+1 ) = -p„ e- ^ X

,dr

(dt

' (ucp )

-const-at p/

■N.n-1

'p - pg

t-t' M p

This equation will be used in the description of the growth process and stabilization of the gas-vapor area in the liquid raw mixture, under heat treatment and swelling. With the equation (25) the quantity of energy, required for the realization of the above process, can be found.

6. The analysis of the received results of calculation of the gas-vapor area of the pore

Calculation of the gas-vapor area (pore) will be made with using the equation (25). The calculation purpose is the prediction of growth dynamics of the pores. As dependences of pore size from time (5) with different levels of energy influence on the mixture (11), (15), (16) are known, the process of heat treatment and final pores size can be controlled. It means that thermophysical properties of the material can be predicted.

Calculations were made by the finite difference method (Fig. 2, 3).

The raw mixture (Table 2) was swollen under temperatures 110-190 oC. Temperatures were found by differential thermal analysis. In this temperature range the first exothermic effect, associated with evaporation of the moisture,

was observed.

12 10 8 6

4 -2 0

1200 1000

E

800 a- 600 -

¡3 400 -aL

200 0

180 c

1 °r

0,0 0,2 0,4 0,6 0,8 1,0 Time i-103, s

2 4 6 8 10 12 Time x-103, s

b

Fig. 2. Changing in the time of: a — vapor pressure; b — pore

radius under swelling temperatures 130 oC and 180 oC

Decreasing of the heat flow under increasing radius can be seen in Fig. 3, a. This happened because of reduction of the temperature difference between a raw mix and vapor. The radius changes in time (Fig. 2, b), but characteristic oscillations of the radius are not observed, because the mass flow is able to compensate the expansion effect of the vapor volume (reducing of the vapor density). The curves

of movement speed changes of the water-vapor boundary (Fig. 3, b) are similar in the character of changing in time. The initial jump of the mass flow and movement speed of the water-vapor boundary, are the result of the sharp decline of the temperature and vapor density, which also led to decreasing of vapor pressure (Fig. 2, a).

The oscillations of vapor pressure (Fig. 1, a) are the result of the combined effect of temperature and vapor density. And if pp = f(Tp) is used in calculations, results will be incorrect. Vapor pressure has an oscillatory character of changing in time with constantly decaying amplitude, which leads to some balance of vapor pressure.

The results of calculations for initial temperature to=130 °C can be seen in Fig. 2, 3. Reducing the

(25)

temperature to this value leads to a reduction of the oscillations of measured quantities, or to complete disappearance of them.

1,6 1,4 1,2 1,0 0,8 0,6 0,4 0,2 0.0

ISO T

—-,

130 <J

£ 4

0

ISO c

L /-

m

12 3 4 Time x-103, s

-0,2 0,0 0,2 0,4 0,6 Time T'103, s b

Fig. 3. Changing in the time of: a — heat flow to vapor; b — movement speed of the pore boundary under swelling temperatures 130 oC and 180 oC

The character of movement speed changes of the water-vapor boundary (Fig. 3, b) is consistent with the change of vapor temperature Tp, as with temperature 180 °C. In general, the peak values of the growth speed of vapor volume (Fig. 3, a), movement of the pore boundary, mass flow, heat flow have smaller values under lower periods of oscillations. The general character of changing of the calculated values under 130 oC and 180 oC is the same.

The character of changing of the measured quantities for the vapor pore radius 10 and 500 mkm is the same, as with 100 mkm, but amplitude quantities increase with decreasing radius. At a small value R, the curve changing shows a sharp decrease of the mass flow, which is caused by a decrease of heat gain from the water, as a result of its intensive cooling (rapid decline of temperature T). Increasing of the mass flow with R3=500 mkm is caused by the longer heat gain from the water.

Equations of the growth model of vapor phase were used for calculation of the distribution of pressure near the growing vapor value area, for certain radiuses. Results of calculation can be seen in Fig. 4.

Experimental and calculation results of pressure distribution in time have good coherence in the frequency of pressure change. The mismatch of the pressure change amplitude is due to the fact that in the calculation the drop with spherical shape was taken, but the real drop hasn't ideal spherical shape.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

In the calculation, the period of vapor phase formation until the size 1 mkm wasn't considered. This assumption is taken into account by the equation (15). But Fig. 4 shows that appearance of each individual bubble can be identified with the formation of a vapor layer with a certain size. Comparison

t=t

a

a

of the experimental and calculated data shows the relative error in the determination of the pulse pressure =15.8 %.

40

i 1 ->

1 1 1

1 _____1 I IHv\ /v —

T(1) = T(2).

-,(!) -

■p(2) =

2g

9(1) (p(1),T) = ^(2) (p(2),T).

These conditions allow estimating the energy parameters of the swelling process.

2. Pressure of the pore former agent gas inside the closed spherical pore in equilibrium state

-20

0 0,5 1 1,5 2 2,5 3 3,5 4 Time, ms

Fig. 4. Comparison of calculated (curve 1) and experimental (curve 2) data of growth of the vapor phase under boiling water

Verification of equations, based on indirect estimates (pressure, temperature) and comparison of the calculated pore size with real allows making a conclusion about the adequacy of the proposed methodology of assessing the main parameters of the swelling process.

The main advantages of the proposed method are the ability to assess the energy impact on the liquid raw mixture, determine the conditions of controlled swelling and controlled structure formation of the material with the predicted thermal properties. These results are proposed to use in designing technological processes of production of porous materials for various purposes.

The continuation of this work will be the study of physical and chemical phenomena in the material during its heat treatment, to optimize the composition of the raw liquid mixture.

7. Conclusions

1. Equilibrium conditions of the pore former agent in the

material during the formation of the porous structure, with account of surface tension

p(1) = P„ ■ e

_a

"rt

where Ea = cp(2) - cp^ - activation energy, equal to the difference between the chemical potentials of the material and pore former agent.

The bigger the difference between the chemical potentials of material-pore systems, the lower the gas pressure will be inside the pore. Since the convective heat transfer in a gas is directly proportional to the pressure, next statement can be made: to achieve minimum heat transfer of the pore, it's necessary to increase the difference between the chemical potentials of material-pore systems.

3. Overheat temperature of the pore former agent gas inside the closed spherical pore in equilibrium state

AT = -

1

Rln p1

,(i)

2av(1)

r(2)

- T

where T0 - temperature of the pore former gas under the same physical conditions outside the pore.

This equation was used for assessing the intensity of the material structure formation.

4. The gas-vapor area was calculated. Obtained equations give a chance to calculate size changing of the pore in time with different levels of energy influence on the mixture and give a chance to control structure formation of the material with predicted properties.

References

1. Bodnar'ova, L. Povedenija legkogo keramzita betona, podverzhennyh vozdejstviju vysokih temperatur [Text] / L. Bodnar'ova, R. Hela, M. Hubertova, I. Novakovu // Mezhdunarodnyj zhurnal grazhdanskoj, jekologicheskoj, strukturnoj, stroitel'stva i arhitek-turnogo proektirovanija. - 2014. - Issue 8 (12). - P. 1205-1208.

2. Nimmo, J. R. Porosity and Pore Size Distribution [Text] / J. R. Nimmo // Encyclopedia of Soils in the Environment. - London: Elsevie, 2004. - P. 295-303.

3. Shpac, A. Clasters nanostrukturnye materials. Vol. 3 [Text] / A. Shpac, P. Cheremskoj, Ju. Kunickij, O. Sobol // Poristost 'as a special state samoorganizovannoi structure in the solid state and materials. - Kyiv: Akademperiodika, 2005. - 516 p.

4. Freire-Gormaly, M. The Pore Structure of Indiana Limestone and Pink Dolomite for the Modeling of Carbon Dioxide in Geologic Carbonate Rock Formations [Text] / M. Freire-Gormaly // Department of Mechanical and Industrial Engineering University of Toronto. - 2013. - 85 p. - Available at: https://tspace.library.utoronto.ca/bitstream/1807/42840/1/Freire-Gormaly_Mari-na_201311_MASc_thesis.pdf

5. Eom, J.-H. Processing and properties of macroporous silicon carbide ceramics: A review [Text] / J.-H. Eom, Y.-W. Kim, S. Raju // Journal of Asian Ceramic Societies. - 2013. - Vol. 1, Issue 3. - P. 220-242. doi: 10.1016/j.jascer.2013.07.00

6. Komissarchuk, O. Pore structure and mechanical properties of directionally solidified porous aluminum alloys [Text] / O. Komis-sarchuk, Z. Xu, H. Hao // China Foundry. - 2014. - Vol. 11, Issue 1. - P. 1-7. - Available at: https://doaj.org/article/002c72e2e0 1345db8bf4fef190113057

7. Bajare, D. Lightweight Concrete with Aggregates Made by Using Industrial Waste [Text] / D. Bajare, J. Kazjonovs, A. Korjakins // Journal of Sustainable Architecture and Civil Engineering. - 2013. - Vol. 4, Issue 5. doi: 10.5755/j01.sace.4.5.4188

8. Lopez-Pamies, O. Effects of internal pore pressure on closed-cell elastomeric foams [Text] / O. Lopez-Pamies, P. Ponte Castañeda, M. I. Idiart // International Journal of Solids and Structures. - 2012. - Vol. 49, Issue 19-20. - P. 2793-2798. doi: 10.1016/j.ijsol-str.2012.02.024

9. Vesenjak, M. Influence of pore gas in closed-cell cellular structures under dynamic loading [Text] / M. Vesenjak, A. Öchsner, Z. Ren // German LS-DYNA Forum. - Bamberg, 2005. - Available at: https://www.dynamore.de/de/download/papers/forum04/ new-methods/influence-of-pore-gas-in-closed-cell-cellular

10. Aboudi, J. Micromechanics of Composite Materials: A Generalized Multiscale Analysis Approach [Text] / J. Aboudi, S. M. Arnold, B. A. Bednarcyk. - Elsevier, 2013. - 973 p.

11. Bratuta, E. Porous insulating materials [Text] / E. Bratuta, A. Pavlenko, H. Koshlak. - Kharkiv: "Eden", 2010. - 105 p.

12. Pavlenko, A. M. Peculiarities control the forming of the porous structure [Text] / A. M. Pavlenko, H. V. Koshlak, B. O. Usenko // Metallurgical and Mining Industry. - 2014. - Issue 6. - P. 50-55.

13. Pavlenko, A. M. Kinetics of evaporation in the homogenization [Text] / A. M. Pavlenko, R. A. Klimov, B. I. Basok // Prom. heat engineering. - 2006. - Vol. 28, Issue 6. - P. 14-20.

14. Vukalovich, M. P. Termodinamika [Text] / M. P. Vukalovich, I. I. Novikov. - Moscow: "Mashinostroenie", 1972. - 672 p.

i Надоели баннеры? Вы всегда можете отключить рекламу.