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

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

CC BY
55
14
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВСПУЧИВАЮЩИЕСЯ ПОКРЫТИЯ / INTUMESCENT COATINGS / ТЕМПЕРАТУРА ПЕЧИ / OVEN TEMPERATURE / ПОТЕРЯ МАССЫ / WEIGHT LOSS / ОБРАБОТКИ ПОВЕРХНОСТИ / SURFACE TREATMENT / ЭФФЕКТИВНОСТЬ ПОКРЫТИЯ / PHASE TRANSFORMATION FRONT

Аннотация научной статьи по физике, автор научной работы — Tsapko Ju., Tsapko А.

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

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

Похожие темы научных работ по физике , автор научной работы — Tsapko Ju., Tsapko А.

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

Simulation of the phase transformation front advancement during the swelling of fire retardant coatings

Description of the intumescent coating behavior at the time of formation of the porous structure is a separate and challenging task, covering both stages of the thermal insulation process: swelling of the coating and subsequent heat transfer which is formed by the swelling. Therefore, there is a need to study the formation conditions of a barrier for heat conduction and reveal a mechanism of phase transition from the coating film to the coke layer. In this regard, the mathematical model of the phase transformation front advancement during the swelling of fire-proof coatings is developed. According to the dependencies, it is found that the front line of the phase transformations of the coating under high temperature passes instantly. It is foundexperimentally that under the action of the heat flux on the samples for a short time at 190÷200 °C, there is an intense swelling of the coating, the height of the expanded foam coke layer increased to 22÷38 mm. As a result of testing, it is revealed that the phase transformation front moves in the direction of high temperature to form foam coke. The foaming front boundary line in the form of a thin layer, which is slightly shifted towards the temperature, divides the coating into two parts. On the one side, there is a swollen coke layer, the outer part of which moves at a certain speed, on the other side the layer of the source material, where the temperature is not sufficient to start the foaming process and the speed of transformations is zero.

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

-□ □-

Виявлення мехатзму роботи покриття в буд^ вельних конструкциях дозволить оцнити умови експлуатаци виробiв, визначити ефективтсть застосування методiв i засобiв вогнезахисту. Розроблена математична модель та розрахова-но процес просування фронту фазових перетво-рень при спучуванн вогнезахисного покриття. Наведено результати експериментальних досл^ джень фронту фазових перетворень при спучу-ванн вогнезахисного покриття та встановлено вплив наповнювачiв на нього

Ключовi слова: спучуючi покриття, температура печi, втрата маси, оброблення поверхт,

фронт фазових перетворень

□-□

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

Ключевые слова: вспучивающиеся покрытия, температура печи, потеря массы, обработки

поверхности, эффективность покрытия -□ □-

UDC 614.842.5:349.211

|DOI: 10.15587/1729-4061.2017.73542|

SIMULATION OF THE PHASE TRANSFORMATION FRONT ADVANCEMENT DURING THE SWELLING OF FIRE RETARDANT COATINGS

J u. Tsapko

Doctor of Technical Sciences* V. D. Glukhovsky Scientific Research Institute for

Binders and Materials Kyiv National University of Construction and Architecture Povitroflotskyi ave., 31, Kyiv, Ukraine, 03680 E-mail: [email protected] А. Tsapko Postgraduate student* Е-mail: [email protected] *National University of Life and Environmental Sciences of Ukraine Heroiv Oborony str., 15, Kyiv, Ukraine, 03041

1. Introduction

Currently, the most common building material is wood, which belongs to the group of average-flammability combustible materials. To reduce these shortcomings, fireproofing provides the wood with the ability to withstand flame exposure and propagation.

One way of wood flameproofing is application of a fireproof coating, which for a while prevents the access of heat to wood.

Swelling of the coating is accompanied by softening of components while endothermic decomposition of fire retardants and blowing agents, leading to the formation of a dense foam coke layer, which causes the fireproof properties of the coating. In this case, the outer coating surface moves towards the influence of flame temperature, which provides the necessary insulation of the material, and the other coating layer, which has not undergone changes and is separated by a line of phase transition, is in reserve.

The main reactions that lead to the formation of a protective foam coke layer occur in the high-temperature region and are complicated to establish the role of the components in the swelling process and advancement of the phase transformation front of foam coke. So, there is the need to study the formation conditions of a barrier for heat conduction and reveal a mechanism of phase transition from the coating film to the coke layer.

2. Literature review and problem statement

In recent years, fire protection research works, aimed at synthesis of coatings using organic paints, refractory oxides and silicates have been known. These materials form temperature and heat-resistant ceramic phases during heating [1, 2]. The most common are crystal glass and enamel coatings [3], but they can't provide reliable protection of structures at long-term temperatures as they are unable to form a necessary coke layer, heated and destroyed.

Modern fire protection methods include the use of intumescent coatings, which are complex systems of organic and inorganic components [4] and have a high intumescent ability. The effectiveness of flame retardant coatings based on organic substances has been shown in [5]. Flame retardants based on polyphosphoric acids and blowing agents may significantly influence the porous layer of foam coke. The above works present the component compositions, efficiency, and thermal characteristics of the coating, but not the mechanism of coke formation, phase and thermal transitions of the coating in foam coke.

Description of the behavior of intumescent coatings, one of the objectives of which is integration of experimental data with available theoretical models, has been considered in [6]. It allows at least a general estimation of simplifications made. Therefore, a thermophysical model, whose solution is

©

given by polynomials that are not associated with a physical content has been examined.

The mathematical model of heating of fire-proof coatings, based on the laws of conservation of matter and energy has been considered. The models immediately assume a particular type of functional relationships with a set of uncertain factors and the problem is reduced to determining the numerical values of these factors, due to the high inaccuracy [7].

On the basis of the kinetics considered, the coating model is represented, which is a mixture of starting materials and products of their transformations, including gas [8]. An increase in the volume of gas when heated causes swelling, but the impact of temperature on the coating churning zone movement has not been determined. The coefficient of thermal conductivity of fire-resistant coatings in fire has been estimated on the basis of numerical simulation, where the phase transitions in the coating have not been taken into account [9]. The effect of inorganic fillers on fire-retardant water-based coatings has shown to be effective. However, the coating swelling mechanism has not been specified and the coating phase transitions in coke have not been found [10]. The authors have presented the analytical model to calculate the thermal conductivity of the porous structure of foam coke of the fireproof coating, which allows for the shapes of pores. However, the model does not account for what phase transitions of the coating occur. Therefore, simulation of the phase transformation front during the swelling of the fireproof coating, under the impact of the components is the unresolved element of ensuring fire resistance of building structures. This led to the need for research in this area.

3. The purpose and objectives of the research

The aim is to study the process of phase formation of foam coke during the swelling of fire-proof coatings.

To achieve this aim, the following objectives were accomplished:

- to conduct simulation of the process of advancement of the phase transformation front of the intumescent coating under thermal exposure;

- to determine timing features of full swelling of fireproof coatings under thermal exposure.

4. Materials and methods of the research of the phase transformation front during the swelling of fireproof coatings

4. 1. Experimental materials

The study was conducted using a system of ammonium polyphosphate (PFA), melamine pentaerythritol (PER) and a binder-based on PVA dispersion. Titanium dioxide, talc, aluminum and magnesium hydroxide and mixtures thereof were added to the above mixture.

Experimental samples of coatings were based on the system containing 18-20 % PFA, 12-14 % melamine, 10-12 % PER, 16 % PVA dispersion and water. The resulting mass was stirred with injecting the fillers in an amount of up to 10 % and applied on a 0.5-0.01 mm steel plate (Fig. 1).

The fillers do not significantly after the appearance and structure of the coating, slightly increasing its viscosity.

Fig. 1. Samples of the coating on a steel plate

4. 2. Method of determining the properties of samples

Simulation of the advancement process of the phase transformation front of the intumescent coating under thermal exposure was performed with the use of the main provisions of mathematical physics.

Experimental determination of the swelling process was performed according to the method given in [11]. The method consists in the fact that a sample with the flame retardant coating is placed in the sample holder and inserted into a test chamber. The sample was fixed so that the end of the control thermocouple touched the back surface of the sample. The heating of the test chamber was turned on, the temperature on the reverse surface and in front of the sample was measured with a thermocouple and the coating was expanded by a ruler. According to the measured values, the heat-insulating properties and changes of the coating were determined.

The criterion for determining the phase transformation front of the intumescent coating under thermal exposure is the foam coke formation of the coating thickness. Herewith, the phase transformation front is observed in the form of a thin layer between the foam coke layer and source material.

5. Simulation of the phase transformation front advancement during the swelling of fire-proof coatings

The model for obtaining the dependencies for calculating the rate of advancement of the phase transformation front of the intumescent coating under thermal exposure is proposed. The focus was on three areas (Fig. 2):

1 - environment, x<0;

2 - intumescent foam coke layer area, 0<x<Z(t) (Z - the coordinate of the front of phase transformation of the coating film in the swollen foam coke layer, m);

3 - material sample with solid material (coating film), Z<x<h (h - half the thickness of the sample, m).

The differential equation that describes the process:

where Ti is the temperature in the "i" region, oC; ^ = ^a"; ai is the thermal diffusivity in the "i" region, m2/s; t is the sample residence time in the high-temperature environment, s.

ifn-EUf

:i"pf!

Front of phase transformations

df

df

Building / materials design

Fig. 2. The scheme of the fire-proof coating swelling process: 1 — environment; 2 — foam coke layer; 3 — starting material

At the initial time (t=0), the temperature in the environment is maximum Tmax, and in the sample material - minimum Tm;n, thus, the initial and boundary conditions can be written as follows:

Tl(x,0) = f Ts(x,0) = T

Z(0) = 0.

(2)

At the border of the sample and the environment, the boundary condition of the fourth kind, which means maintaining the heat flow at the transition from the environment into the inner region of the sample (the law of conservation of heat flow) is given [12]:

1 dTi

= 1

df2 "dX

(3)

where li, l are the coefficients of thermal conductivity of the environment (1) and the foamed coke layer (2), W/(m-°C).

Given the continuity of the heat flow fields, the temperatures on the sample surface and the front of the phase transformation of the intumescent coating into the foam coke layer coincide:

Ti(0, t) = T2(0, t); f2[Z(t), t] = T3 [Z(t), t] = Tf.

(4)

(5)

I2 —1—=q+13 —1—

dx

! dT2

12 ^X"

x=Z(t)

dx

1 dT3

x=Z(t)

=jf.

(7)

Thus, the problem is reduced to the solution of the differential equation (1) with the boundary conditions (2)-(5), (7).

To solve the differential equation (1), we introduce di-mensionless variables:

t-f

t =

h2 '

Ti - If If - f

(8)

and then the problem is formulated as follows: - to solve the differential equation:

d2Ti = 1 dTL dt2 = ? dt,

with the boundary conditions:

Ti(t,0) - 0;

Ts( t,0) -1;

(9)

1i ^

1

= T dT2

x=0

x=0

Ti(0, t) = f2(0, t); f2(Z, t) = T3(Z, t) = Tf;

l2"aXT

X=z

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

T dT3

X=z

=JZ, dt

(10) (11)

(12)

(13)

(14)

where

+i =Tl,( ^ = 1); 1 i =^,( 13 = 1);

Due to the phase transformation (solid into in foam coke), the temperature is absorbed at a rate q, W/(m3-s), which is proportional to the speed of advancement of the phase transformation front of the coating [13], namely:

» dZ q = Jdt,

(6)

Z* _Tf - f. J = ^LJ.

Z = f ; Tf =

T - T

ß3f - T

To determine the limits of phase transformations and the speed of advancement of the phase transformation front, the solution of equation (9) takes the following form taking into account [14]:

where J is the proportionality factor showing the heat flux absorption rate per unit depth of advancement of the phase transformation front of the coating, W/m4.

The absorption rate of the heat flow q, coming from region 1 to region 2 is absorbed through the formation of the expanded layer of foam coke [6]. Therefore, on the border of phase transformation, the condition should be satisfied:

T = c + b O

where the error integral is

2

2 ? 2

0(?) = Je-? d?.

(15)

(16)

or

3

3

The coefficients c;, b; are determined from the boundary conditions (10)—(14). After substituting (15) to (10), we obtain the equation:

T1(X,0) = C1 + b1®(~) = C1 + b1 = 0;

T3(X,0) = c3 + b30(~) = c3 + b3 = 1.

Since

- x2

dX ^¡V nt

where the sign of the right side coincides with the sign X, then from (11) it follows that:

b1 = -l^b2.

c2 + b20

/ — \ Z

v ^V^

= c3 + b3®

2jt

= Tf.

Z = Z( t ) = 2PV?,

c2 + b20

2

= Tf;

c3 + b3o (P) = Tf.

Ï _

12e ®b2 -e-p2b3 = VJ ^2

VT2 7

Tf -Q(P); ®*(P) ;

b1 = --

nTf

(17)

(18)

(19)

b2 =-

Tf

P

VY2 7

P

VY2/

b =i=TL ;

3 o*(P)'

12 _f

P

v^2 y

1-TL t o*(P)

= Vn JP,

(27)

(28)

(20)

where

Given (15) and F(0)=0, the condition (12) is converted to the form:

(21)

After substituting (15) to (13), we get the equation for the swelling front boundary:

(22)

Since the right side of equation (22) does not change, the arguments of the function F are also constant. This is possible only if the coordinate of the coating swelling front is proportional to the square root of time, provided that:

(23)

where for convenience the proportionality factor is taken as being equal to the doubled, yet unknown constant p. As a result, (22) can be written in the form of two equations:

= 1

A

It follows that all the constants ci, bi depend on the parameter p, which in turn can be determined on the basis of the solution (28) only by approximate calculations [15].

Thus, the calculated dependences, allowing to get changes in the swelling dynamics of fire-retardant coatings are obtained. However, they provide an opportunity to directly calculate the coating swelling time and the phase transformation front movement depending on the effects of temperature.

6. The results of determining the phase transformation front during thermal coating swelling

Fig. 3 shows the calculation of P depending on the di-mensionless complex J, characterizing the heat-insulating properties of intumescent coatings on the border of the phase transformation front in a situation where half of the heat is absorbed by the coating:

Tf =n = 0,5.

(29)

(24)

After substituting (24), (23) in the heat flow continuity condition given Z = X, the equation (14) is converted to the following form:

(25)

From the system of equations (12), (17), (18), (20), (24), (25), with respect to the unknowns c;, bi, P, we find the solutions:

nTf

(26)

Fig. 3. Dependence of the coating swelling front coefficient on the heat flux absorption proportionality coefficient

The results correspond to the physical content of the problem solved, namely the rate of advancement of the phase transformation front decreases given the constant swelling rate of fire-proof coatings.

We solve graphically the equation (28) and obtain the value P=0.062 (Fig. 4).

0,05 0,1 0.1

Front swelling ratio, ¡5

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

Fig. 4. Graphical solution of equation (28): the dependence of the fire-proof coating swelling rate and thermal conductivity on the swelling coefficient.

So, the dependencies, which allow getting a picture of the swelling dynamics and provide an opportunity to directly determine the phase transformation front, are given.

7. Experimental studies of the phase transformation front during the swelling of fire-proof coatings and their results

To find the time of full swelling of fire-proof coatings (Z=h), we get the value of:

t = -

h2

4p2

(30)

Given h=1 mm and a=2.9-10-5 m2/s, the time of full swelling of fire-proof coatings is 2.2 s.

Experimental studies of the phase transformation front during the expansion of a flame retardant coating are shown in Fig. 5, 6.

Fig. 5. The process of coating swelling

Exposure of the coating samples to a heat flow at temperatures of about 190-200 °C in front of the sample caused intensive swelling and a slight temperature increase on the back sample surface, which lasted upon reaching a temperature of about 360-390 °C.

The height of the expanded layer of foam coke increased to 22-38 mm and the phase transformation front line appeared on the surface under the foam coke layer.

As a result of the tests, it is found that under thermal exposure the coating foams and moves towards the high temperature to form coke. The foaming front boundary line in the form of a thin layer, which is slightly shifted towards the temperature, divides the coating into two parts. On the one side, there is a swollen coke layer (dark), the outer part of which moves at a certain speed, on the other side - the layer of the source material, where the temperature is not sufficient to start the foaming process and the speed of transformations is zero (white), Fig. 6, a, b.

Fig. 6. Test results of the intumescent coating:

a — organic and inorganic + (10 %) Ti02; b — organic and inorganic + (10 %) Al(OH)3; c — organic and inorganic + (10 %) Mg(OH)2; d - organic and inorganic + (5 %) Ti02 + (10 %) Al(OH)3

However, in the organic-inorganic coating with the addition of TiO2 and Al(OH)3, a relatively low swelling rate and the absence of the starting material, as well as the phase transition boundary were observed (Fig. 6, d). A large amount of Al(OH)3 led to the suppression of the coating swelling process.

8. Discussion of the results of research of the advancement process of the phase transformation front

When swelling of fire-proof coatings under the thermal effect of high-temperature flame, as indicated by the research results (Fig. 3, 5, 6), there is a natural process of the coating swelling under the influence of temperature and the formation of the phase transition line, due to the mechanism of the foam coke layer formation, which slow down the heat transfer processes. It should be noted that the presence of more than 10 % of additives, such as aluminum hydroxide, reduces the swelling process and leads to the lack of the phase transition line. Obviously, such a mechanism of the effect of additives is a factor regulating the degree of coke destruction and thermal insulation efficiency. But at the same time, this mechanism adversely affects the intumescent ability of the coating. This agrees with the data of the known works [6, 9], the authors of which have also attributed the changes in the swelling process to the addition of mineral

fillers. In contrast to the results of the studies [6, 9], the data on the influence of fillers on the process of phase transition of the coating into the foam coke layer and changes in the insulation properties suggest the following:

- the main process regulator is not only the formation of the foam coke layer, but also heat resistance of the fire retardant coating;

- correct addition of mineral fillers has a significant influence on the phase transition of the fireproof coating into the vermiculite foam coke layer.

The results of determining the phase transition of the fireproof coating into the vermiculite coke layer (Fig. 6) indicate a mixed impact of fillers on the intumescent ability of the coating. Such uncertainty can't be resolved within the present study because this would need additional experiments to obtain more reliable data. In particular, this requires sufficient data for carrying out the effective swelling process and identifying the time of the beginning of heat resistance decline. Such identification will allow exploring the conversion of the coating, which moves in the direction of high temperature to form coke, and identifying the variables that significantly affect the beginning of the process transformation.

This work is a continuation of the research presented in [11], which fully shows the formation mechanism, movement and high-temperature heat insulation of foam coke.

9. Conclusions

As a result of the studies, the process of phase formation of foam coke during the swelling of fire-proof coatings was revealed, in particular:

- the simulation of the advancement process of the phase transformation front during the swelling of the fireproof coating was performed, the swelling front ratio was determined and settlement dependences allowing to get changes in the swelling dynamics of fire-proof coatings under the influence of temperature were obtained. According to the dependencies, the foam coke formation time under thermal decomposition of the coating, which makes up 2.2 s, subject to the double heat absorption, was calculated;

- the phase transformation front moves in the direction of high temperature to form coke, the foaming front boundary line in the form of a thin layer divides the coating into two parts. On the one side, there is a swollen coke layer, the outer part of which moves at a certain speed to the source of heat. On the other side - a layer of starting material, where the temperature is insufficient to start the foaming process and the speed of transformations is zero. Further research may be aimed at studying the process of establishing the relationship between the components and properties of coatings and their effect on the heat resistance of building structures.

References

1. Gyvlyoud, M. M. Temperaturostiyki sylikatni zahysni pokryttya dlya metaliv ta splaviv na osnovi napovnenogo polimetylfenilsy-loksanu [Text] / M. M. Gyvlyoud, O. I. Bashynskiy, S. Ya. Vovk // Zbirnyk naukovyh prazch Lvivskogo derzhavnogo universytetu BZHD. - 2011. - Issue 18. - P. 40-45.

2. Artemenko, V. V. Eksperementalni doslidzhennya vognezahysnyh pokryttiv metalevyh konstruktsiy na osnovi napovnenyh polialyu-mosyloksaniv [Text] / V. V. Artemenko // Zbirnyk naukovyh prazch LDU BZHD. Pozhezhna bezpeka. - 2014. - Issue 25. - P. 6-11.

3. Timofeeva, S. V. Materialy ponizhenoy pozharnoy opasnosti s pokrytiem na osnove zhidkih siloksanovyh kauchukov, otverzhdennyh metodom poliprisoedineniya [Text] / S. V. Timofeeva, A. S. Malyasova, O. G. Helevina // Pozharovzryvobezopasnost. - 2011. -Vol. 20, Issue 9. - P. 22-25.

4. Antsupov, E. V. Antipireny dlya poristyh materialov [Text] / E. V. Antsupov, S. M. Rodivilov // Pozharovzryvobezopasnost. -2011. - Vol. 20, Issue 5. - P. 25-32.

5. Gravit, M. V. Issledovanie vliyaniya razlichnyh faktotov nakoefitsient vspuchivaniya organorastvorimyh ognezaschitnyh pokrytiy [Text] / M. V. Gravit // Lakokrasochnye materialy i ih primenenie. - 2013. - Issue 6. - P. 12-16.

6. Nenahov, S. A. Fiziko-himiya vspenivayushchihsya ognezashchitnyh pokrytij na osnove polifosfata ammoniya (obzor literatury) [Text] / S. A. Nenahov, V. P. Pimenova // Pozharovzryvobezopasnost. - 2010. - Vol. 19, Issue 8. - P. 11-58.

7. Khalturinskiy, N. A. O mehanizme obrazovanija ognezaschitnyh vspuchivajuschihsja pokritijy [Text] / N. A. Khalturinskiy, T. A. Rudakova // Izvestija JuFU. Tehnicheskie nauki. - 2013. - Issue 8. - P. 220-232.

8. Sharshanov, A. Ja. Matematicheskaja model' vspuchivajuschihsja ognezschitnih pokritijy [Text] / A. Ja. Sharshanov // Problemi pozharnoy bezopasnosti. - 2011. - Issue 30. - P. 273-280.

9. Cirpici, B. K. Assessment of the thermal conductivity of intumescent coatings in fire [Text] / B. K. Cirpici, Y. C. Wang, B. Rogers // Fire Safety Journal. - 2016. - Vol. 81. - P. 74-84. doi: 10.1016/j.firesaf.2016.01.011

10. Fan, F. Effects of inorganic fillers on the shear viscosity and fire retardant performance of waterborne intumescent coatings [Text] / F. Fan, Z. Xia, Q. Li, Z. Li // Progress in Organic Coatings. - 2013. - Vol. 76, Issue 5. - P. 844-851. doi: 10.1016/j.porg-coat.2013.02.002

11. Kryvenko, P. Determination of the effect of fillers on the intumescent ability of the organic-inorganic coatings of building constructions [Text] / P. Kryvenko, Y. Tsapko, S. Guzii, A. Kravchenko // Eastern-European Journal of Enterprise Technologies. -2016. - Vol. 5, Issue 10 (83). - P. 26-31. doi: 10.15587/1729-4061.2016.79869

12. Taganov, I. N. Modelirovanie processov maso- I energoperenosa. Nelineynie sistemi [Text] / I. N. Taganov. - Leningrad: Himija, 1979. - 208 p.

13. Nenahov, S. A. Dinamika vspenivaniya ognezaschitnyh pokrytiy na osnovi organo-neorganicheskih sostavov [Text] / S. A. Nenahov, V. P. Pimenova // Pozharovzryvobezopasnost. - 2011. - Vol. 20, Issue 8. - P. 17-24.

14. Samarskiy, A. A. Vichislitelnaja teploperedacha [Text] / A. A. Samarskiy, V. P. Vabischevich. - Moscow: Editoril URSS, 2003. - 784 p.

15. Bahvalov, N. S. Chislennie metodi [Text]: ucheb. pos. / N. S. Bahvalov, N. P. Zhidkov, G. M. Kobel'kov. - Moscow: Nauka, 1987. - 600 p.

.................................................................................................................................................................................................................................KS

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