STRUCTURAL AND MATHEMATICAL MODEL OF THE THERMAL CONDUCTIVITY OF CONCRETE
Vladimir Gryzlov1, Alla Kaptyushina2
12Cherepovets State University Lunacharsky Prospect, 5, Cherepovets, Russia
Corresponding author: gryvs@mail.ru
Abstract
Introduction: The significance of energy saving in residential construction is associated with the reliable determination of design thermal conductivity characteristics of construction materials. The authors describe a concept of the structural and mathematical modeling of concrete thermal conductivity. The concept is based on methods of the structural approach and generalized conductivity theory. Purpose of the study: The study is aimed at developing an adequate structural and mathematical model to determine the thermal conductivity of concrete. Methods: The authors used the statistical method as a representation of the statistical homogeneity of a multi-component composite closely related to thermal homogeneity, which is understood as a medium with effective thermal resistance constant in space. Results: The authors developed a structural and mathematical model to determine the thermal conductivity of concrete. The model accounts for concrete structural factors existing at the time when the structure formation of concrete is, for the most part, complete. The model also accounts for the thermal and physical properties of concrete components as well as macro- and mesostructural features of concrete. The paper presents calculated data for potential macro- and mesostructural factors, which makes it possible to determine the thermal conductivity of concrete characterized by low thermal conductivity and use the results derived for the approximation of prediction trends related to the thermal properties of such concretes in the operation and adaptation period during the use of building envelopes.
Keywords
Thermal conductivity, structural and technological factor, mathematical model, prediction of thermal and physical properties.
Introduction
The issue of the reliable determination of concrete thermal conductivity is one of the most important in modern construction material engineering. The potential of improving heat-protection efficiency and resource-saving in the manufacturing and operation of building envelopes depends on its successful solution.
The following model of concrete thermal conductivity was suggested in the earlier paper of the author (Gryzlov, 2008):
h(t) = ho - Aho(t) e - (t - to)/T (1)
where l(t) is concrete thermal conductivity at fixed time t; \0 is concrete thermal conductivity when the structure formation process is, for the most part, complete (28 days of hardening); A \0(t) is an increment that occurs as a result of operations during time t; t0 is time of intense structure formation; t is time of thermal conductivity relaxation.
The minus sign before \ 0(t) indicates that thermal conductivity adaptation is related to balancing constructive and destructive processes occurring in concrete that, in general, lead to a decrease in internal stress and transition to a stable equilibrium state of exponential nature. Equation (1) can be considered a general equation for the formation and prediction of concrete thermal and physical properties
during the operation of building envelopes. The solution to the prediction problem reduces to the optimization of \0 and t. The \0 parameter is a structural and technological aspect of the property, which can be found at the stage of structure selection and formation; the t value depends not only on the internal parameters of \ but also on the concrete adaptation mode.
According to many researchers (Buzhevich, 1970; Ivanov, 1974; Karamyan, 1976; etc.), the thermal conductivity of concrete is closely related to its density since pores in concrete are filled with air to a large extent. Such an approach is determinative in selecting the design values of thermal conductivity coefficients, and it is captured in regulatory documents (Ministry of Regional Development of Russia, 2012). Table 1 provides empirical equations to determine the thermal conductivity coefficient for inorganic particulate materials, recommended for a wide range of materials, including concrete with a density of 400 < y < 1800 kg/m3. These equations support the concept of the predominant influence of the material density on its thermal conductivity. However, they also indicate a significant variation of data on the thermal conductivity coefficient at the same density (AA > 50%), which is even more significant in the case of extreme values (Figure 1).
Table 1. Empirical equations to determine the thermal conductivity of concrete
No. Author Equation
1 Buzhevich G. A. X = 0.38r10-3 - 0.12
2 Vlasov O. E. X = 0.2 Y + 0.05 Y2
3 Ivanov I. A. X = 0.0005r - 0.25
4 Kaufman B. N. X = 0.11Y 111.68 y + 0.022
5 Karamyan K. O. X = 0.046 + 0.16Y2
6 Nekrasov V. P. X = 1.16 il(0.0196 + 0.22y2) - 0.16
7 Spektor B. V. X = 0.029 + 2.19 • 10-4 y
X is the thermal conductivity coefficient; X is density
0 8 0.7 u 0.6
0
J 05
> 0.4
Ï3
u
1 0.3
0
u
1 0.2
QJ
E 0.1 0
,6
/ y / /
1 I S 3
/ /2
M y ^ 7
r - ' 4;5
v_ _
0
500 1000 1500
Concrete density, kg/m3
2000
Figure 1 Graphical representation of the equations given in Table 1.
It should be noted that this variation is a consequence of experimental studies conducted by the authors with regard to a certain group of concretes, which confirms that, in the generalized analysis of various types of concrete composites, the structural concept of thermal conductivity exists. The equations given are expressly linear, which indicates a certain bias in the authors' reasoning.
When developing the concept of the influence of concrete porosity on its thermal conductivity, a number of researchers suggested determining thermal conductivity as models of particulates with account for their porosity (Chudnovsky, 1962; Dulnev and Zarechnyak, 1974;
Krischer, 1934). An analysis of these models has shown that the models by Krischer and Bernshtein are the most appropriate for concrete. To a certain extent, these models interpret the dependence of thermal conductivity on porosity and can be transformed into linear functions of the following type:
A = A - a Pr
c m gen
(2)
In the case of fixed values of matrix and inclusion thermal conductivity (Table 2), the equations become identical and have the same value of the dimensionless factor a.
Table 2. Models of particulates' thermal conductivity with account for porosity
Author
Model
Equation
* Equation transformed into the following form: 1 = f(P)
Krischer O.
100
(100 - PrgJ + Prgen
100 - Pr Pr
X = Xl-—+ X
100
100
A = A - 0.23Pr
c m
A = A - 0.23Pr
cm
Bernshtein R. S.
X = X-
4Pr X
_sen + _Ln _ 2pr )
1 + ^ X2
A = A - 0.23Pr
cm
Note. A , A 2 are matrix and inclusion thermal conductivity values; Pr n is general porosity, fractions of volume;
*example: A, = 0.54 W/(m • °C; A 2 = oJw/(m • °C).
at 1 — const
The most important thing arising out of the analysis of these models is the confirmation of the statement that, for concretes with continuous matrix structure and minimum porosity, the matrix, i.e. the cement paste (sand-cement mortar), is the main heat conductor.
With regard to the recipe and technological factors as well as macrostructure, mathematical models
describing the dependence of the thermal conductivity coefficient on the inclusion volume are quite interesting (Table 3). While having different mathematical structures, these models, at fixed values, also transform into the linear form, which confirms their qualitative and not quantitative nature.
Table 3. Mathematical models of thermal conductivity with account for the inclusion volume
Author Equations 'Equation transformed into the following form: 1 = /(P) at 1 m / 1 _ = const
Odelevsky V. I. ^c = ^ / \ 1 - p 1 p 1 - p I P-i 3 J A = A - 0.22 P c m
Dovzhik V. G. ^ = ^ + - 2P(Xm - 2A.m + + P(K - ^J \ / A = A - 0.22 P c m
Missenard, A. ^c = ^ I 1 -3VP - (1 -ft) A = A - 0.22 P c m
p is the relative inclusion volume; p = A m/ A in; example: A m = 0.54 W/ (m • °C); A inc = 0.3 W/ (m • °C).
The equations given in Tables 2 and 3, with relevant assumptions, can be transformed into the following form: \ c =f(\ m, j), and they have a more information-bearing structure than the equation (Table 1).
kc =f (Prgen) = k„, ~a (1 ~ J / Jctd)
(3)
where jai is concrete true density, a is a dimensioniess factor.
Ac = f(P) = In - a(y - 1.23C)/yai
(4)
keff = 27/1,- vr + 27 4h kj Vi Vj / (ki + kj)
(5)
will be characterized by the maximum thermal resistance, while the second one will be characterized by the minimum one.
where jd is aggregate density in a lump; C is cement consumption per 1 m3 of concrete.
However, each of the equations, taken separately, does not provide sufficient insight into the influence of the differential distribution of pores, defectiveness of the structure, granulometric, phase and mineral composition of the components on the thermal conductivity of concrete. In this context, the purpose of the study was to:
- develop a structural and mathematical model to determine the thermal conductivity of concrete that would reflect the influence of its structural factors existing at the time when the structure formation of concrete is, for the most part, complete.
Methods
Actual multi-component composites consist of particles of no particular form. Therefore, in order to determine the effective thermal conductivity of such media, it is recommended to use the statistical method. This method can be applied to statistically homogeneous media. Visual interpretation of the statistical homogeneity of a medium means that grains of heterogeneous parts shall be sufficiently small as compared to the characteristic dimensions of the body and shall be distributed in a random manner but uniformly in space. Statistical homogeneity is closely related to thermal homogeneity, which is understood as a medium with effective thermal resistance constant in space. When applying the basic concepts of the generalized thermal conductivity theory (Odelevsky, 1951) to a multi-component matrix system, the equation for the calculation of the effective thermal conductivity, provided that the conditions of structural and statistical homogeneity are met, can be transformed into the following form:
a)
b)
Figure 2. Regarding the calculation of the maximum and minimum thermal resistance of the model: a — the layers are perpendicular to the heat flow; b — the layers are in parallel to the heat flow
The maximum (R ) and minimum (R ) thermal
s max' x miw
resistance of this model are described below. Let us assume that A , A are thermal conductivity coefficients
m mc
with regard to the matrix (binding agent) and inclusion (aggregate); P is the relative inclusion volume, P = l -P is the matrix volume; f} = Xm/X. is the ratio of the thermal conductivity coefficients.
Scheme a. Heat flow Q does not depend on the point on the model surface. In this case, the heat flow is expressed clearly through the difference of temperatures on the opposite surfaces.
AT
Q = —
R
where r = 1 ~ /'"1C + ^
m
^inc
For a unit layer:
/
R =
h
\
mc
1 -
V L J
h
mc
Àinc
(6)
Since (h.JL) = iAit, equation (6) can be transformed into the following form:
^ = [1 +P( p-l)]/V
(7)
where A eff is the thermal conductivity coefficient of the system, W/(m degree); A , A , v., v. are thermal conductivity coefficients and relative volumes of components, respectively.
which will correspond to the maximum thermal resistance of the model.
Scheme b. Heat flow Q depends on the point on the model surface. In points related to the inclusion:
0mc = A77(IAinc)
(8)
This equation is the first approximation to the actual situation. If we present a two-component model of concrete in the form of a layer structure with layers located perpendicular (Figure 2 a) and in parallel (Figure 2 b) to the heat flow (Krischer's model), it is obvious that the first case
In points related to the matrix:
Om = AT/(L/lm)
(9)
To calculate the average value of Q on the model surface, let us denote the inclusion surface area as S. ,
and the matrix surface area — as S . The total area is
m
equal to: S=S +S .
^ inc m
The average density of the heat flow on the surface S equals to:
ß = ß
inc
inc
5 m ß —m .
(10)
Let us apply the values of Qnc and Qm to equation (10). Then, we will obtain the following:
Q = AT
s x
inc inc
S m X m
\
V
SL
SL J
(11)
Hence, the thermal resistance of the model is as follows:
L
R =
S
inc
inc
For unit thickness, we will obtain:
1
R =-
S
inc
Xinc
m
(12)
Taking into account that S /S corresponds to P, we will
& inc ~ '
transform equation (12) into the following form:
Rmin = [1+ P (ß - 1) / ( P + (1 - P)ß] / Im
(13)
which will correspond to the minimum thermal resistance of the model.
By comparing equations (7) and (13), a general structural model of the thermal conductivity of concrete can be found:
Xeff = lm / [1 + P (ß - 1)x]
(14)
Figure 3. Nomographic chart of the calculated values of x'
In the model under consideration, p is a boundary condition, a peculiar kind of the criterion of sensitivity that
divides concretes into ordinary ones (p < 1) and those with low thermal conductivity (p >1), which is adopted as one of the control factors at the microstructural level.
Results and discussion
The model presented (14) can be considered universal since, through algebraic transformations, a number of known mathematical models of thermal conductivity based on the principles of the generalized conductivity theory (Table 4) are transformed into the form (14), which shows a good approximation of this model with regard to the structural interpretation of the thermal conductivity of concrete. The analysis of the models transformed points at their intermediate value between the values of R and
max
Rmin in equations (7) and (13), i.e. each of these models reflects a certain principle of inclusion location between the maximum and minimum values of the thermal resistance indicators and, therefore, the thermal conductivity of concrete. Being structural and mechanical, they do not reflect the influence of physical and chemical processes occurring in the contact area, i.e. they do not take into account the level of the mesostructure.
Table 4. Transformation of known mathematical models of thermal conductivity into the form (14)
In this equation, x is a dimensionless control factor, which is a product of x' (macrostructural level factor) and x'' (mesostructural level factor). To determine x', the expression 1 / [p + (1 - p)p] from equation (13) is used. It should guarantee the minimum error when determining the thermal conductivity of concrete. A nomographic chart of the calculated values of x' is given in Figure 3.
Author, source
Odelevsky V. I.
Dovzhik V. G.
Missenard, A.
Khlevchuk V. R. Kim L. N.
Dependence type, transformation into the form (14)
i i
1 + (p- 1)P-
P 1 - P p-1 3
3
2(p + l-2P(p-1)
(14.1)
^ -
2V + - 2P(^m - W
1 1
am + ^inc + 2-
1 + (P- 1) P-
m '*mc 3
2ß + l-2P(P- 1)
(14.1)
-■ 1 ~ß i -3 4P (i -p)
i+(p-1) p-
i+(p-1)
(14.2)
1 -3 4P(1 -p)
1 1
Xm Xm
1 + (P- 1) P-
1 + (P- 1)1 P3 - P
(14.2)
= À
c m
Àm m
X. = À
c m
The assessment of this level was conducted experimentally and made it possible to derive a condition for the assignment of x". First of all, it was established that layered models consisting of the cement paste and the aggregate, imitating a contact area, have an integrative indicator of thermal conductivity that is less than the additive one by 10-13% (Kirpichev and Konakov, 1949); at (fi > 1), the following condition is in place: Aca <A,nc<Am(A ca is the thermal conductivity coefficient of the c ontact area).
To assess the impact of the contact area, we will represent the inclusion as a set n of the fractions of the aggregate. Each fraction consists of cubes with the edge (d1 < d2 < ...< d). The cubes are oriented in such a way that the heat flow is perpendicular to a surface. The relative share of the ith fraction of the aggregate is = 1. Then, with a known relative aggregate volume p, the values of Pca (contact area volume) and Pm (matrix volume) can be det ermined by the following equations:
Pc.a. = 2d P (3S ai di + 6d S ai di + 4d2) / S ai di
3.
Pm 1 Pc.a. P
(15)
With the thickness values known, we can find the thermal resistance of the parallelepiped R by the following equation:
Rp=L[Plhnc + 2 dP / Àc.a. S cadi + + (1 -P-2dP/I. ai d,)/ V] / S
(16)
To transform equation (16) into the form corresponding to the general structure of x"(14), we will introduce following parameters: D = d/d,, i.e. the ratio between the thickness of the contact layer an d the size of the aggregate grain, and p = im/ii.e. the ratio between the thermal conductivity coefficient of the matrix and the thermal conductivity coefficient of the contact area.
After the transformation, we will obtain:
Xeff = xm /[1 + P (fi - 1) x' (1 + 2D (p - 1) / (fi - 1)] (17)
Let us adopt x" = 1 + 2Dp - 1) / (P - 1),
By substituting values of x' and x" in equation (13), we will obtain:
where d is the thickness of the contact layer that depends not on the size of the aggregate particles but on the aggregate technology and material.
Let us consider a parallelepiped with thickness L consisting of the layers of the matrix, contact layer and the aggregate, joined up in series, in a concrete wall. The cumulative thickness values will be as follows:
for the aggregate layers — p L; for the contact area layers — 2dPL / z a. d; for the matrix layers — (L - LP -2dP L)/Z a. d..
Xeff = lm [P + P (1 - P)] / [p + 2D P (p - 1)]
(18)
where [p + fi (1 - p)] / [fi +2D P (p - 1)] can be generally adopted as the structural and technological factor fy) at the stage of the formation of integrative quality, in this case — the thermal conductivity of concrete.
The results of calculating x for theoretical and practically possible values of p fi, D, p in concretes with low thermal conductivity are given in Table 5.
Table 5. Results of calculating the macro- and mesostructural factor x in concretes with low thermal conductivity at
fixed values of v , S, D, p
p p D P X P P D P X P P D P X
0.8 2.0 0.2 2.0 0.517 0.75 2.0 0.2 2.0 0.543 0.7 2.0 0.2 2.0 0.614
1.5 0.555 1.5 0.581 1.5 0.654
1.2 0.581 1.2 0.606 1.2 0.680
0.15 2.0 0.535 0.15 2.0 0.584 0.15 2.0 0.633
1.5 0.566 1.5 0.615 1.5 0.665
1.2 0.585 1.2 0.635 1.2 0.685
0.12 2.0 0.548 0.12 2.0 0.596 0.12 2.0 0.645
1.5 0.572 1.5 0.622 1.5 0.671
1.2 0.589 1.2 0.638 1.2 0.688
0.8 1.5 0.2 2.0 0.604 0.75 1.5 0.2 2.0 0.622 0.7 1.5 0.2 2.0 0.646
1.5 0.662 1.5 0.678 1.5 0.701
1.2 0.703 1.2 0.718 1.2 0.739
0.15 2.0 0.632 0.15 2.0 0.649 0.15 2.0 0.672
1.5 0.679 1.5 0.695 1.5 0.716
1.2 0.710 1.2 0.725 1.2 0.745
0.12 2.0 0.650 0.12 2.0 0.666 0.12 2.0 0.693
1.5 0.689 1.5 0.704 1.5 0.726
1.2 0.715 1.2 0.729 1.2 0.750
continuation of the Table 5.
p ß D P X P ß D P X P ß D P X
2.0 0.684 2.0 0.700 2.0 0.716
0.2 1.5 0.764 0.2 1.5 0.777 0.2 1.5 0.834
1.2 0.822 1.2 0.833 1.2 0.843
0.8 2.0 0.722 0.75 2.0 0.726 0.7 2.0 0.752
0.15 1.5 0.788 0.15 1.5 0.800 0.15 1.5 0.811
1.2 1.2 0.833 1.2 1.2 0.843 1.2 1.2 0.853
2.0 0.747 2.0 0.760 2.0 0.774
0.12 1.5 0.802 0.12 1.5 0.813 0.12 1.5 0.825
1.2 0.840 1.2 0.849 1.2 0.859
Thus, knowing or assuming the structural and technological characteristics of concretes with low thermal conductivity, using equations (1) and (18), it is possible to approximate the predictive trends of heat-protective properties of these concretes in the operation and adaptation period during the use of building envelopes.
Conclusions
1. The authors developed an adequate structural and mathematical model to determine the thermal conductivity coefficient of concrete. The model accounts for concrete structural factors existing at the time when the structure formation of concrete is, for the most part, complete.
2. This model, albeit a bit complicated, makes it possible to assign the value of the thermal conductivity
coefficient of concrete on a case-by-case basis, taking into account its structural and technological characteristics.
3. It has been confirmed that in concretes with continuous structure, the matrix is the main heat conductor and thermal conductivity factor; the effectiveness of the matrix can be improved or reduced through its modifications or application of the binder with a lower thermal conductivity coefficient.
4. For concretes with low thermal conductivity, the aggregate with the thermal conductivity that is significantly lower than the thermal conductivity of the matrix shall be selected and used, and standard operating procedures, ensuring an increase in the degree of airborne dispersion of the mesostructure and, therefore, an increase in its thickness, shall be applied.
References
Bernshtein, R. S. (1948). Thermal conductivity of a layer. In: G. G. Knorre (ed.) Studies of natural fuel combustion processes. Moscow-Leningrad: Gosenergoizdat, pp. 82-88.
Buzhevich, G. A. (1970). Lightweight concrete with porous aggregates. Moscow: Stroyizdat, 272 p.
Chudnovsky, A. F. (1962). Thermal and physical characteristics of disperse materials. Moscow: Fizmatgiz, 456 p.
Dovzhik, V. G. (1972). Thermal conductivity of keramzit. Stroitel'nye materialy, 3, pp. 21-23.
Dulnev, G. N. and Zarechnyak, Yu. P. (1974). Thermal conductivity of mixtures and composite materials. Leningrad: Energiya, 264 p.
Gryzlov, V. S. (2008). Universal model of the thermal conductivity of lightweight concrete. Stroitel'nye Materialy, 1, pp. 44-46.
Ivanov, I. A. (1974). Technology of lightweight concrete with artificial porous aggregates. Moscow: Stroyizdat, 287 p.
Itskovich, S. M. (1977). No-fine concrete. Moscow: Stroyizdat, 119 p.
Karamyan, K. O. (1976). Temperature- and moisture-induced deformations of the building envelope made of lightweight concrete with natural porous aggregates. Yerevan: Ayastan, 308 p.
Kaufman, B. N. (1955). Thermal conductivity of construction materials. Moscow: State Publishing House of Literature on Construction and Architecture, 160 p.
Khlevchuk, V. R., Kim, L. N. and Shteinman, B. I. (1983). Calculation of the thermal conductivity of lightweight concrete depending on structural and technological factors. Lightweight concrete house construction. Moscow: Central Research and Design Institute of Residential Buildings, 121 p.
Kirpichev, P. V. and Konakov, P. K. (1949). Mathematical basics of the similarity theory. Moscow-Leningrad: Publishing House of the USSR Academy of Scienses, 104 p.
Krischer, O. (1934). Der Einfluß von Feuchtigkeit, Körnung und Temperatur auf die Wärmeleitfähigkeit körniger Stoffe (Die Leitfähigkeit des Erdbodens). Beihefte zum Gesundheits-Ingenieur, Series I, No. 33.
Ministry of Regional Development of Russia (2012). Regulations SP 50.13330.2012. Thermal performance of the buildings. Revised edition of SNiP 23-02-2003 (Construction Rules and Regulations). Moscow: Federal Autonomous Organization "Federal Center for Standardization", 95 p.
Missenard, A. (1968). Thermal conductivity of solids, liquids, gases and their compositions. Moscow: Mir, 464 p.
Nekrasov, A. A. (1940). Revisiting the dependence of the thermal conductivity coefficient of disperse systems on the external pressure. Technical Physics, 10 (2), pp. 168-173.
Odelevsky, V. I. (1951). Calculation of the generalized conductivity of heterogeneous systems. Technical Physics, 21 (6), pp. 667-677.
Spektor, B. V., Svintsitskaya, L. G. and Gutgarts, G. L. (1962). Thermal and physical characteristics of keramzit perlite concrete. Kiev: Academy of Construction and Architecture of the Ukrainian Soviet Socialist Republic.
Starostin, F. D. (1935). Calculation of the thermal conductivity of porous solids. Otopleniye i Ventilyatsiya, 3, pp. 24-25.
Vlasov, O. E. (1938). Fundamentals of construction heat engineering. Moscow: Military and Engineering Academy of the Workers' and Peasants' Red Army, 94 p.
СТРУКТУРНО - МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТЕПЛОПРОВОДНОСТИ БЕТОНА
Владимир Сергеевич Грызлов1, Алла Германовна Каптюшина2
12Череповецкий государственный университет
Проспект Луначарского, 5, Череповец, Вологодская область, Россия
1Е-та11: gryvs@mail.ru
Аннотация
Актуальность энергосбережения в строительстве жилых зданий непосредственно связана с достоверным определением расчетных теплотехнических характеристик строительных материалов. Изложена концепция структурно - математического моделирования теплопроводности бетона, в основу которой, положены методологии структурного подхода и обобщенной теории проводимости. Цель исследования. Разработка адекватной структурно - математической модели по расчету теплопроводности бетона. Метод. Использован статистический метод, как интерпретация статистической однородности многокомпонентного композиционного материала, которая тесно связана с термической однородностью, под которой понимается среда с постоянным по пространству эффективным термическим сопротивлением. Результаты. Разработана структурно - математическая модель по расчету теплопроводности бетона, учитывающая влияние его структурных факторов, сформированных на момент завершения, в основном, процесса структурообразования. В разработанной модели учитываются теплофизические свойства компонентов бетона, его макро- и мезо-структурные особенности. Приведены расчетные данные практически возможных макро - мезо - структурных факторов, что позволяет определять величину теплопроводности малотеплопроводного бетона и использовать полученные результаты для аппроксимации прогнозных трендов теплозащитных свойств этих бетонов в эксплуатационно - адаптационном периоде работы ограждающих конструкций.
Ключевые слова
Теплопроводность, структурно - технологический фактор, математическая модель, прогноз теплофизических свойств.