Research Article  Open Access
Modeling of Viscosity and Thermal Expansion of Bioactive Glasses
Abstract
The behaviors of viscosity and thermal expansion for different compositions of bioactive glasses have been studied. The effect of phosphorous pentoxide as a second glass former in addition to silica was investigated. Consequently, the nonlinear behaviors of viscosity and thermal expansion with respect to the oxide composition have been modeled. The modeling uses published data on bioactive glass compositions with viscosity and thermal expansion. regression optimization technique has been utilized for analysis. Linear and nonlinear relations are shown to establish the viscosity and thermal expansion coefficients associated with oxide components of the glasses under study. The modeling allows the calculation of viscosity for a given temperature and, accordingly, the fusion temperature of these glasses along with the coefficient of thermal expansion. The established model relations also suggest first and secondorder phosphorusalkali and alkaline earth oxides interaction which is reflected on the model coefficient that calculates viscosity and thermal expansion.
1. Introduction
In a review of the importance of comprehensive information on glass and melt properties both for the glass science and technology by Mazurin [1], it is shown that the number of publications per year concerning glass property measurements have decreased quite considerably after 1975. Moreover, the frequency of some quite important but laborconsuming property measurements decreases steadily. Besides that, these measurements are quite time and money consuming; thus, the specialist on glass science and technology needs to predict the glass properties via calculation of it. This task must be based on modeling methods which use experimental property data and calculate glass properties including newly designed glasses.
The changes in the composition area of glass often lead to the change of the reliable modeling, that is, different models may be needed. Finding the most reliable model is possible only by comparing the results of property calculations with experimental data for glasses belonging to a system of interest. Next, the model with a minimal rootmeansquare error is selected. After that, the model is used for future calculations of glass properties. At present, this is the best way to obtain a property value that is as near to the true one as possible. This explains why modeling of glass properties is an enduring active field in glass science and technology [2].
Viscosity at given temperature and the coefficient of thermal expansion CTE for bioactive glasses are the glass properties which are focused on in this study. Bioactive glasses have the property of forming an apatite layer; thus, when the implants is coated with bioactive glass, it can chemically bond to bone [3]. However, if there is a large difference in coefficient of thermal expansions between the coating and the implant, the bioactive glass will crack resulting in incomplete coating. A second problem is that the fusion temperature is needed to be well estimated for efficient cast of the glass coating to the implant. The fusion temperature should be low enough to avoid phase transformation and degradation of the mechanical properties of the implant. Therefore, scientific efforts must be carried out to achieve glass compositions that are suitable in view of fusion temperature and coefficient of thermal expansion [4]. The glass designer intends to assemble bioactive glass with fusion temperature below 900°C to avoid possible phase transformation of for example, Tialloys; also, the value CTE should be around which should match that of Tialloys [5]. In addition, the trend in bioactive glass compositions is to set P_{2}O_{5} content to 6% with high CaO content to ensure the formation of apatite layer after implantation that is responsible for the biocompatibility of the glass.
In general, viscosity and thermal expansion for bioactive glass shows nonlinear behavior with glass composition. Therefore, any proposed model should express effective compositionrelated terms that describe these properties for bioactive glass. Consequently, the model calculation can be used in design technology in such glass category [6].
The topic of predicting the viscosity and thermal expansion of bioactive glass is still vigorous for different sets of bioactive glass compositions [6, 7], and the modeling study based on formerly published data is still on focus in the applied ceramics community [8–10]. This work presents modeling methodology, results, and discussion of the modeled bioactive class along with the comparison on the effect of CaO and MgO content on the properties of the bioactive glass.
2. CompositionViscosity and Thermal Expansion Model
The variation of viscosity with the absolute temperature is expressed as the following [11]: The parameter is constants, is the gas constant, and is the activation energy of the viscous process.
The natural logarithm of (1) readsas the following: The above equation shows that and − are linked linearly. Accordingly, values of ln can be calculated for a given set of viscosity at temperature and for a range of glass compositions.
is as function of glass composition as well as the coefficient of thermal expansion CTE [12]:
is the weight percentages wt% of the th component composing the glass and or is the coefficient associated with the th glass component. The terms () or () can also include nonlinear combination of more than one glass component. Finally, n represents the total number of terms.
The final step is to explicitly express the compositionviscosity or thermal expansion model that is specific to the bioactive glass composition set under study. In this study, the effort to present such a model is enlightened by the discussion of Kuppinger and Shelby [13], that is, inspecting possible interaction between ions composing the glass which may lead to the nonlinear relations of the bioactive glass composition with and CTE. In this study, phosphorousalkali and alkaline earth effect is found. So, the suggested compositionglass property model is stated as below:
, , , and are first and secondorder functions of phosphorousalkali and phosphorousalkaline earth content:
Equations (4a)–(6b) are the model equations presented in this study that calculate and CTE as function of the glass compositions.
3. Modeling Technique
Table 1 shows published bioactive glass composition data with their thermal properties [5]. In this table, Ts represent the temperature at constant viscosity level of μ = 10^{6.6} Pa·s, which corresponds to the Littleton softening temperature. Tg represents another constant viscosity level of μ = 10^{12} Pa·s, which is in the range of the expected glass transition temperature (μ = 10^{11}10^{12.3} Pa·s).

Utilizing (2), the two unknowns ( and ln) can be found by substitution of Ts and Tg in Table 1. Figure 1 is a plot of against –ln which shows linear relation connecting them. The figure also shows the equation of the leastsquares line passing through the data points. The reference values of and ln found utilizing Table 1 is used together with the glass compositions in regression [14] calculations. Accordingly, calculated glass properties ( and CTEc) are found as function of their composition.
The assumption of linear additive behavior is firstly examined. This assumption yields to large differences between reference experimental and calculated glass properties. These differences remain at large values even when several combinations of the wt% of oxide constituents are taken into account. The differences between the reference experimental and calculated glass properties are reduced only when phosphorousalkali and alkaline earth effect as shown in (4a)–(6b) is taken into account. The calculations astonishingly show that the nonadditive relations are equally needed for modeling of both and CTE. The final rootmeansquare error for and CTE calculations are 26.47 and 0.23, respectively. These are only around 2% of the reference experimental value for and CTE.
Table 2 shows the reference experimental , reference experimental ln , calculated via regression, calculated ln (utilizing calculated and the leastsquares relation shown in Figure 1), reference experimental CTE, and calculated CTE via regression. The final coefficients and that obtained by regression to calculate and CTE for bioactive glasses are shown in Table 3.


To find the temperature at a given viscosity value (an inverse problem), it needs to calculate and ln by the above methodology, then (2) can be solved for . This can be accomplished using any numerical technique such as the “finite step method.” Finally, setting to 20 Pa·sec; the fusion temperature for a bioactive glass can be calculated [15].
To conclude, Figure 2 is a way to show how and CTE respond to CaO and Na_{2}O content. and CTE are calculated by setting K_{2}O, MgO, and Na_{2}O to 3, 6, and 9 wt%, respectively, and let SiO_{2} varies with CaO to make total contents to 100 wt%. The same is done when varying Na_{2}O by setting K_{2}O, MgO, and CaO to 3, 6, and 11 wt%, respectively. The figure is presented to show an example of the nonlinear behavior of and CTE with CaO and Na_{2}O. Yet, and CTE are linear with other oxide contents as shown in (4a)–(6b)
(a)
(b)
(c)
(d)
4. Discussion
Table 3 shows that SiO_{2} has a trend to increase and to decrease CTE. At the same time, P_{2}O_{5} tends to decrease both and CTE. MgO shows similar behavior with SiO_{2}; also, K_{2}O behaves similar to P_{2}O_{5} with different values of their coefficients. Figure 2 describes the behavior of CaO and Na_{2}O more obviously. Na_{2}O shows a trend to increase both and CTE in nonlinear manner. The effect of Na_{2}O is distinguished from that of K_{2}O although they are both alkali oxides. Similarly, Figure 2 shows that CaO is decreasing and increasing CTE in nonlinear fashion. This effect of CaO is also distinguished from that of MgO although they are both alkaline earth oxides.
The aforementioned different type of contribution for both Na_{2}O and CaO to the viscosity and thermal expansion compared with that of K_{2}O and MgO, respectively, may be understood in view of that Na2O and CaO content is generally higher than that of K_{2}O and MgO in the bioactive glass as seen in Table 1. Consequently, large ion concentration of Na and Ca give a chance for PNa and PCa interactions. These interactions are seen in a way that P_{2}O_{5} is modifying the role of Na_{2}O and CaO to that observed in Figure 2 and reflected formally in (9)–(12). Nevertheless, those interactions on a relatively lowconcentration level of K_{2}O and MgO do not influence the properties apparently. This explanation is illuminated by the discussion of Kuppinger and Shelby [13] that suggest short range ionion interaction in the glass medium that yields nonlinear relations of glass viscosity and thermal expansion with the oxide content. They study sodium potassium borate glasses with high sodium and potassium oxides contents which led to ionion interaction between the alkali ions. Nevertheless, in this study the high sodium and calcium oxide contents do not led to mutual interaction between them, but phosphorousalkali and alkaline earth interaction takes place. The acidic nature of P_{2}O_{5} may explain the interaction with alkali and alkaline earth oxides of the basic nature.
It is worth to mention that multilayer coating is one of the experimental experiences to avoid crack generation which originates from CTE mismatch [7]. However, multilayer coating may lead to other problems such as thick coating. As a direct application of the modeling presented in this study, a bioactive glass can be designed with suitable fusion temperature and CTE passing up lengthy experiments and cost.
5. Conclusions
(1)Model relations of viscosity and thermal expansion with bioactive glass compositions have been obtained. The model allows the predictions of viscosity and thermal expansion for bioactive glass necessary in the field of glass science and technology.(2)The model relations are discussed for the nonlinear behavior of viscosity and thermal expansion (deviation from additivity) with Na_{2}O and CaO content. It is suggested that the deviation from additivity originates from shortrange interaction between alkali and alkaline earth ions with phosphorus as formally presented in the model.
References
 O. V. Mazurin, “History, perspectives, and problems of measurement and calculation of glass properties,” in Proceedings of the 3rd Balkan Conference on Glass Science and Technology, pp. 1–8, Varna, Bulgaria, 2005. View at: Google Scholar
 O. V. Mazurin, “Glass properties: compilation, evaluation, and prediction,” Journal of NonCrystalline Solids, vol. 351, no. 1213, pp. 1103–1112, 2005. View at: Publisher Site  Google Scholar
 L. L. Hench, “The story of bioglass,” Journal of Materials Science, vol. 17, no. 11, pp. 967–978, 2006. View at: Publisher Site  Google Scholar
 E. Saiz, M. Goldman, J. M. GomezVega, A. P. Tomsia, G. W. Marshall, and S. J. Marshall, “In vitro behavior of silicate glass coatings on Ti6Al4V,” Biomaterials, vol. 23, no. 17, pp. 3749–3756, 2002. View at: Publisher Site  Google Scholar
 J. M. GomezVega, E. Saiz, A. P. Tomsia et al., “Novel bioactive functionally graded coatings on Ti6Al4V,” Advanced Materials, vol. 12, no. 12, pp. 894–898, 2000. View at: Google Scholar
 E. Vedel, H. Arstila, H. Ylänen, L. Hupa, and M. Hupa, “Predicting physical and chemical properties of bioactive glasses from chemical composition—part 1: viscosity characteristics,” Glass Technology Part A, vol. 49, no. 6, pp. 251–259, 2008. View at: Google Scholar
 E. Vedel, D. Zhang, H. Arstila, L. Hupa, and M. Hupa, “Predicting physical and chemical properties of bioactive glasses from chemical composition—part 4: tailoring compositions with desired properties,” Glass Technology Part A, vol. 50, no. 1, pp. 9–16, 2009. View at: Google Scholar
 K. Yang, W. H. Zheng, and J. S. Cheng, “Effect of Li_{2}O on viscosity and thermal expansion of silicate glass,” Advanced Materials Research, vol. 403–408, p. 7074, 2011. View at: Google Scholar
 M. Wang, J. Cheng, M. Li, F. He, and W. Deng, “Viscosity and thermal expansion of sodalimesilica glass doped with Gd_{2}O_{3} and Y_{2}O_{3},” Solid State Sciences, vol. 14, no. 8, pp. 1233–1237, 2012. View at: Google Scholar
 D. Bellucci, V. Cannillo, and A. Sola, “Coefficient of thermal expansion of bioactive glasses: available literature data and analytical equation estimates,” Ceramics International, vol. 37, no. 8, pp. 2963–2972, 2011. View at: Publisher Site  Google Scholar
 G. Urbain, F. Cambier, M. Deletter, and M. R. Anseau, “Viscosity of silicate melts,” Transactions & Journal of the British Ceramic Society, vol. 80, pp. 139–141, 1981. View at: Google Scholar
 I. N'dala, F. Cambier, M. R. Anseau, and G. Urbain, “Viscosity of liquid feldspars—part I: viscosity measurements,” Journal of the British Ceramic Society, vol. 83, no. 4, pp. 105–107, 1984. View at: Google Scholar
 C. M. Kuppinger and J. E. Shelby, “Viscosity and thermal expansion of mixed alkali sodium potassium borate glasses,” Journal of the American Ceramic Society, vol. 68, no. 9, pp. 463–467, 1985. View at: Google Scholar
 R. J. Vanderbei, Linear Programming: Foundations and Extensions, Princeton University, Princeton, NJ, USA, 2001.
 J. R. Taylor and A. C. Bull, Ceramics Glaze Technology, Pergamon Press, New York, NY, USA, 1986.
Copyright
Copyright © 2012 Saad B. H. Farid. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.