Research Article  Open Access
Extracting Concrete Thermal Characteristics from Temperature Time History of RC Column Exposed to Standard Fire
Abstract
A numerical method to identify thermal conductivity from time history of onedimensional temperature variations in thermal unsteadystate is proposed. The numerical method considers the change of specific heat and thermal conductivity with respect to temperature. Fire test of reinforced concrete (RC) columns was conducted using a standard fire to obtain time history of temperature variations in the column section. A thermal equilibrium model in unsteadystate condition was developed. The thermal conductivity of concrete was then determined by optimizing the numerical solution of the model to meet the observed time history of temperature variations. The determined thermal conductivity with respect to temperature was then verified against standard thermal conductivity measurements of concrete bricks. It is concluded that the proposed method can be used to conservatively estimate thermal conductivity of concrete for design purpose. Finally, the thermal radiation properties of concrete for the RC column were estimated from the thermal equilibrium at the surface of the column. The radiant heat transfer ratio of concrete representing absorptivity to emissivity ratio of concrete during fire was evaluated and is suggested as a concrete criterion that can be used in fire safety assessment.
1. Introduction
Fire safety for structural members in buildings is necessary to be assessed as failure of structural members due to fire can cause extensive property damage or loss of life. Many research efforts have been conducted to assess fire safety for structural members [1–3]. For a reinforced concrete (RC) member exposed to fire, the increase in temperature in both steel and concrete causes a decrease in strength and stiffness of the member. Therefore, most design codes for RC structures limit temperatures in RC member subject to fire [4–8]. Alternatively, effective sections can be considered for fire safety assessment of RC members such that the strength of RC member under fire is evaluated with a reduced strength for reinforcing steel based on its temperature while the strength of concrete over 500°C is ignored [7, 8]. Temperature variations in the RC member exposed to fire depends on the thermal properties of concrete and steel and the shape of the crosssections. The thermal properties of concrete and steel are obtained from fire test or using prediction methods. While many methods were proposed to predict temperature variations in concrete members subject to fire [9–13], there exists difficulty to choose the thermal characteristics, especially thermal conductivity of concrete. This work suggests a simple method to identify thermal characteristics of concrete under fire.
The thermal conductivity of concrete is a major property that affects the temperature variations in concrete members subject to fire. In normal concrete, the thermal conductivity varies from 0.5 W/(mK) to 3.3 W/(mK) depending on the type of aggregate and the concrete mixture used [15–19]. As the thermal conductivity of concrete is relatively lower than that of other construction materials such as 53.3 W/(mK) at 20°C for structural steel [20], RC members exhibit good performance in fire if extensive spalling is prevented [14, 21]. The thermal conductivity of concrete has also been reported to be affected by the measuring techniques [19, 22–24]. To measure the thermal conductivity of concrete at a test temperature, it requires the concrete specimen to reach the status of thermal steadystate for the test temperature [23, 24]. As the thermal conductivity of concrete is relatively low, it takes significantly long time for a concrete brick to reach the status of thermal steadystate for a test temperature which is typically a high temperature. For example, the soak period of concrete bricks (230 mm × 114 mm × 65 mm) for a test temperature of 600°C is four days. This relatively long timeperiod is necessary to measure the thermal conductivity of concrete following ASTM C1113 [14]. If a test time period shorter than the ASTM standard is to be used, it would be necessary to develop a method to measure thermal conductivity from temperature variations in RC under unsteadystates condition.
The heat flux is transferred from one medium to the other by conduction, convection, and radiation [25]. In general, the heat flux is transferred from fires to structures by both convection and radiation [26, 27]. The convective heat flux is proportional to the difference in temperatures between the convection currents due to fire and the structural surface [26, 27]. The radiant heat flux of fire will be transferred to the surface of a structure by absorptivity ratio of the surface while the structure emits the radiant heat flux in proportion to its emissivity [28]. The transferred heat flux from a fire to the structure is then transferred into the structure by conduction. In fire tests of RC structures, the thermal steadystate condition when the temperatures of gas in fire test furnaces and the whole body of test structure become equal has been rarely observed [11]. However, if a decreasing temperature change with time is observed, at the surface of RC structure during a fire test using a standard fire, the thermal radiation properties of RC structure might be estimated by considering thermal equilibrium of the surface.
In this study, we suggest an inverse methodology to identify thermal conductivity of concrete from time history of temperature variations during fire test. The developed methodology was applied to find the thermal conductivity of RC column exposed to a standard fire. Fire test of an RC column was conducted following ISO834 standard fire curve [29]. The profiles of time history of temperatures at six different locations across the RC column section were obtained. The thermal conductivity of RC columns was then determined by minimizing errors between the numerical solution and the temperature observations with time. The determined thermal conductivity was compared to that observed from experiments following ASTM C1113 [24] for verification of the proposed method. Moreover, the thermal radiation properties of the RC column were estimated considering the thermal equilibrium at the surface of the RC column in fire. The use of the thermal radiation properties of RC columns as a fire safety index is discussed.
2. Experimental Methods
In this work, two experimental programs were conducted as shown in Figure 1. First, fire test of an RC column was prepared to obtain the time history of temperature variations in a section of the column. The temperature variations will be used to extract the thermal conductivity of the column using a proposed numerical method. To verify the results of the thermal conductivity extracted from RC column testing, the second experimental program to measure the thermal conductivity of concrete bricks was conducted following ASTM C1113 [24].
2.1. Fire Test of RC Column
An RC column was tested to evaluate the time history of temperature variations following ISO834 standard fire curve [29] for 180 minutes. The ISO834 standard fire curve is presented as where is the gas temperature in fire test furnace at time in minutes. is the ambient temperature. The RC column was constructed following the design specification proposed by the Korean Concrete Institute (KCI) [30]. The RC column had the square section of 350 mm by 350 mm and the height of 1500 mm. The column was reinforced with 822 mm steel rebars longitudinally and 10 mm hoop bars placed at 300 mm spacing along the column height. Grade 420 steel rebars, which have the yield strength of 420 MPa, were used to reinforce the column. Details of column reinforcement are shown in Figure 2. A concrete having the 28day characteristic compressive strength of 39.9 MPa was used for the column. The maximum aggregate size of 25 mm pea gravel was selected for the concrete. The slump of fresh concrete was controlled in the range from 180 mm to 220 mm. The mixture proportion of the concrete is presented in Table 1. The temperatures in concrete following the standard fire curve were measured at the center section of the column located at 750 mm from the bottom of the column. Across the section of the column, temperature was measured at the surface, 20 mm, 40 mm, 60 mm, 80 mm, and 175 mm from both sides, respectively. The temperatures in longitudinal steel rebars were also measured for all 8 rebars at the center section where the temperatures of concrete were measured. The locations of the thermocouples are shown in Figure 2. The column was placed in fire test furnace as shown in Figure 3 and the time histories of temperatures following the standard fire curve for 180 minutes were recorded.

2.2. Thermal Conductivity Measurement of Concrete Bricks
To verify the proposed inverse method to identify thermal conductivity, a standard test method to measure thermal conductivity of concrete was used as shown in Figure 4. In this standard method, the thermal conductivity of concrete bricks constructed with the mixture proportion presented in Table 1 was measured following ASTM C1113 [24]. For the main part of the harness wire, platinum wire of 0.33 mm diameter was used. The specimens consisted of three 230 mm × 114 mm × 65 mm straight concrete bricks. The thermal conductivity of concrete bricks was measured at nine test temperatures of 50°C, 100°C, 200°C, 300°C, 400°C, 500°C, 600°C, 670°C, and 780°C. The thermal conductivity of concrete bricks was measured three times at each test temperature. Heating rate for the furnace temperature was 55°C/hour. Soak periods to make concrete bricks reach thermal steadystate with the test temperature were different for each test temperature.
3. Numerical Methods
In a thermal system without the heat generation, onedimensional unsteadystate heat equation considering both specific heat and thermal conductivity as functions of temperature can be formulated as where is the temperature of concrete at the distance from the surface at time and is the density of concrete considered constant and equal to 2300 kg/m^{3} in all calculations. and are the specific heat and the thermal conductivity of concrete at temperature . Equation (2) can be expanded as For a finite difference approximation, all variables in (3) except the temperature are linearly transformed to be dimensionless by putting , , , and . Equation (3) is then rearranged as Putting , the constant term in the right side of (4) becomes one. A finite difference approximation of (4) will be expressed as where where is the number of time steps and is the number of distance steps. and are the step sizes for time and distance, respectively. and are the thermal conductivity and the corresponding derivative with respect to temperature evaluated at the temperature of while and are the specific heat and the corresponding derivative with respect to temperature evaluated at the temperature of . Using the quadratic forms of the thermal conductivity and the specific heat [6] as the coefficients and in (5) can be calculated as The normalizing factors and for the thermal conductivity and the specific heat are selected as those evaluated at 20°C using (7) and (8), respectively. While for the specific heat in (8) is fixed at 4.0 [6], for the thermal conductivity in (7) will be changed to find the calculated time history of temperature variations which gives the minimum error with the observed time history of temperature variations from the RC column exposed to the standard fire. The founded will be used to determine the thermal conductivity of concrete as a quadratic function of temperature using (7).
4. Results and Discussion
4.1. Thermal Conductivity of Concrete
After finishing the fire test of the RC column for 180 minutes, it was confirmed that there was no spalling of the RC column. The time histories of temperatures measured in concrete at the two locations for the same depth from the concrete surface are averaged and presented in Figure 5 for the depths of 0 mm (surface), 20 mm, 40 mm, 60 mm, 80 mm, and 175 mm, respectively. For the time histories of temperatures at 40 mm, 60 mm, 80 mm, and 175 mm in Figure 5, the temperature increases stop at the temperature of 100°C for a while and then continue. This can be attributed to the effect of capillary water on concrete thermal behavior. For normal concrete, it would have partially saturated capillary pores. This enables concrete to conduct heat rapidly up to the temperature of 100°C. For temperatures above 100°C, water in the capillary pores evaporates and the concrete with empty pores will conduct heat more slowly than that with partially saturated pores.
The time histories of temperatures measured in longitudinal steel bars located at the middle of a side and the corner of the square section are averaged, respectively. The average time histories of temperatures are presented in Figure 6. Because the longitudinal steel bars located at the corner of square section take heat flux from two sides, the temperature of the corner steel bars was always higher than the temperature of steel bars located at the middle of a side. This observation indicates that twodimensional heat transfer mechanism is necessary to predict steel temperature located at the corner of square section. However, the time history of temperature in the steel located at the middle of a side is almost identical to that in concrete at the depth of 60 mm from the surface as shown in Figure 7. Therefore, the temperatures in steel bars located at the middle of a side of the square section can be estimated to be equal to concrete temperatures at the same depth from surface. Using the time histories of temperatures at the six locations, a regression surface which gives the maximum value of 0.983 was achieved with value in (7) of 0.0115 as shown in Figure 8. The corresponding thermal conductivity at 20°C, is calculated as 2.26 W/(mK). In Figures 9(a), 9(b), and 9(c), the calculated time histories of temperature variations using the numerical analysis with of 0.008, 0.0115, and 0.015 are compared with the observed time histories of temperature variations, respectively. The time history of temperature variations in Figure 9(b) is obviously the optimal temperature variation. A monotonic sequence can be observed in Figure 9. The increase of thermal conductivity of concrete enables conducting more heat to the center of the square section of RC column. The thermal conductivity of concrete for RC column with respect to temperature is, then, determined using (7) as shown in Figure 10. The thermal conductivity of concrete bricks measured at the test temperatures 50°C, 100°C, 200°C, 300°C, 400°C, 500°C, 600°C, 670°C, and 780°C following ASTM C1113 [24] are also presented in Figure 10. At each test temperature, the thermal conductivity of concrete bricks was measured three times and the results were averaged. The thermal conductivity computation with respect to temperature of normal concrete proposed by the Euro code [6] () presented in (7) with of 0.008 is plotted in Figure 10. As shown in Figure 10, the standard thermal conductivity is less than until the temperature reaches 600°C. It can also be observed that the difference between the two measurements gets reduced as the temperature increases from 200°C to 600°C. This might be attributed to the fact that is the thermal conductivity for concrete with empty pores because water in the capillary pores evaporates during the soak period to enforce the thermal steadystate at test temperature. However, is the thermal conductivity for concrete with partially saturated capillary pores, which might be more preferable for practical use. The observations of over the temperature of 600°C agree well with . Moreover, the thermal conductivity proposed by the design code is less than . Therefore, it might be necessary to consider the effect of partially saturated capillary pores on the thermal conductivity of concrete for a conservative estimate of RC members’ thermal behavior in fire.
(a) , , W/(mK)
(b) , , W/(mK)
(c) , , W/(mK)
4.2. Radiant Heat Transfer Ratio of Concrete
Examining the temperature difference between the surface of RC column and ISO834 standard fire curve shown in Figure 11, it can be observed that the RC column is still in thermal unsteadystate. At the surface of RC column, the heat flux is transferred from fire to the column by both convection and radiation. The thermal equilibrium at the surface is formulated as where and are the heat flux components due to convection and radiation, respectively. and are the volume and the exposure area of the RC column, respectively. According to the amount of transferred heat flux and the surface temperature, the time rate of increase of temperature will be determined. The amount of heat flux by convective transfer is given by where is the coefficient of convective heat transfer. and are the gas temperature in fire test furnace and the surface temperature of the RC column in Kelvin. By assuming that the surrounding gas in fire test furnace is a black body and the RC column is a grey body, the amounts of heat flux by radiant transfer are formulated with the thermal properties of concrete as where is Stefan Boltzmann’s constant of 5.67 × 10^{−8} W/(m^{2 }K^{4}). is the thermal absorptivity of concrete for the gas in fire test furnace and is the thermal emissivity of concrete toward the gas. As indicated in (12), the radiant heat transfer keeps until is equal to . Therefore, when is equal to , the radiant heat transfer stops. We name the ratio, here, as the radiant heat transfer ratio. Considering this aspect, with time is formulated in Figure 12. As shown in Figure 12, the change of with time decreases as it tends to converge. An exponential equation can be used to approximate the convergence value of such as Using regression analysis of (13) with the formulation of with time, the radiant heat transfer ratio of 0.587 and of 59.2 minutes is determined as shown in Figure 12. Considering that the radiant heat transfer ratio represents the ratio of radiant heat flux inflow to outflow of the RC column, the radiant heat transfer ratio can be used as a reference value to assess fire resistance of RC elements. Concrete with a low the radiant heat transfer ratio will be capable of preventing thermal radiant heat transfer to the entire structure when exposed to fire. It is interesting to note that although the radiant heat transfer ratio cannot be directly compared with that for solar radiation because of the difference in wave lengths of thermal energy in fire test furnace from that of solar energy, the ratio determined in this study (0.587) is close to the value of 0.568 reported for concrete and based on solar absorptivity of concrete 0.5 over the emissivity of concrete as 0.88 [31]. There have been extensive studies in the literature about the ratio of the solar absorptivity over emissivity for selecting appropriate coating materials for spacecraft [28]. The regression analysis to find the radiant heat transfer ratio will be effective when the time rate of temperature increase in (10) significantly decreases such as in Figure 13. It is notable that the radiant heat transfer ratio is developed due to the black body assumption of the gas surrounding the RC column in fire test furnace. Moreover, the thermal radiation of compartment walls making up the fire test furnace is also neglected [32]. Therefore, if the information about the combined emissivity of gas and compartment walls in test fire furnace is available [32], the radiant heat transfer ratio can be divided by and can be used as a criterion for concrete for fire safety assessment.
5. Conclusions
An inverse methodology to find the thermal conductivity from time history of temperature variations in the thermal unsteadystate is proposed. A square RC column with cross section of 350 mm by 350 mm and the height of 1500 mm was tested in order to evaluate the time history of temperature variations following ISO834 standard fire curve. Using a numerical solution for thermal equilibrium model in the unsteadystate condition, the thermal conductivity of concrete for the RC column was determined in quadratic form as in W/(mK) for the temperature range from 20°C to 1200°C. Thermal conductivity of concrete was conservatively estimated using the proposed method when compared with the thermal conductivity of concrete bricks measured using ASTM standards. Moreover, another aspect of fire safety assessment was evaluated from the difference of the applied temperature (ISO834 standard fire curve) and the surface temperature in fire test furnace. The radiant heat transfer ratio representing concrete absorptivity over emissivity ratio was determined from regression analysis of formulation with time. It can be observed that concrete with relatively low radiant heat transfer ratio can protect structures during fire by preventing thermal radiant heat transfer to the structures from fire. It is suggested that the radiant heat transfer ratio shall be used as a concrete criterion for fire safety assessment.
Nomenclatures
:  Specific heat at the temperature of 
:  The coefficient of convective heat transfer 
:  The number of time steps 
:  The number of distance steps 
:  Thermal conductivity at the temperature of 
:  Heat flux components due to convection 
:  Heat flux components due to radiation 
:  The exposure area of RC column 
:  Specific heat at the temperature of 
:  Thermal conductivity at the temperature of 
:  Ambient temperature 
:  The surface temperature of RC column 
:  Gas temperature in fire test furnace 
:  Temperature of concrete at the distance from the surface at time 
:  The volume of RC column 
:  The variable for the thermal conductivity in(7) 
:  Thermal absorptivity of concrete for the gas in fire test furnace 
:  The coefficient for the specific heat in (8), which is fixed at 4.0 
:  Thermal emissivity of concrete toward the gas 
:  The step sizes for time 
:  The step sizes for distance 
:  Density of concrete 
:  Stefan Boltzmann constant, which is . 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The first author would like to acknowledge National Research Foundation (NRF) Grant funded by the Korean government (MOE) (no. 2013R1A1A2062784). The corresponding author would like to acknowledge a Grant (code no. ′09 R&D A01) from CuttingEdge Urban Development Program funded by the Ministry of land, Transport and Maritime Affairs of Korean Government.
References
 N. Pinoteau, P. Pimienta, T. Guillet, P. Rivillon, and S. Rémond, “Effect of heating rate on bond failure of rebars into concrete using polymer adhesives to simulate exposure to fire,” International Journal of Adhesion and Adhesives, vol. 31, no. 8, pp. 851–861, 2011. View at: Publisher Site  Google Scholar
 H. Sen, Z. Liaojun, Z. Hanyun, and L. Longzhong, “Analysis of thermal stress for bent column of concrete filled steel tube in exposure to fire,” Water Resources and Power, vol. 30, no. 1, pp. 177–179, 2012. View at: Google Scholar
 S. J. Choi, S. Kim, S. Lee, R. Won, and J. Won, “Mix proportion of ecofriendly fireproof highstrength concrete,” Construction and Building Materials, vol. 38, pp. 181–187, 2013. View at: Publisher Site  Google Scholar
 ACI Committee 216, Standard Method for Determining Fire Resistance of Concrete and Masonry Construction Assemblies, American Concrete Institute, Detroit, Mich, USA, 1997.
 Canadian Standards Association, Code for the Design of Concrete Structures for Buildings, CAN3A23.3M02, Canadian Standards Association, Toronto, Canada, 2002.
 CEN: DAFT ENV. Euro code 2: Design of Concrete Structures, 1993.
 Ministry of Land, Infrastructures and Transport, The legal proposal of fire safety of beams and columns made of high strength concrete, Notice of Ministry of Land, Infrastructures and Transport 2008334, Sejong, Republic of Korea, 2008.
 AS3600 Concrete Structures Standards Australia, 2009.
 T. T. Lie and R. J. Irwin, “Method to calculate the fire resistance of reinforced concrete columns with rectangular cross section,” ACI Structural Journal, vol. 90, no. 1, pp. 52–60, 1993. View at: Google Scholar
 Z. Huang, A. Platten, and J. Roberts, “Nonlinear finite element model to predict temperature histories within reinforced concrete in fires,” Building and Environment, vol. 31, no. 2, pp. 109–118, 1996. View at: Publisher Site  Google Scholar
 J.T. Yu, Z.D. Lu, and Q. Xie, “Nonlinear analysis of SRC columns subjected to fire,” Fire Safety Journal, vol. 42, no. 1, pp. 1–10, 2007. View at: Publisher Site  Google Scholar
 Y. Xu and B. Wu, “Fire resistance of reinforced concrete columns with L, T, and +shaped crosssections,” Fire Safety Journal, vol. 44, no. 6, pp. 869–880, 2009. View at: Publisher Site  Google Scholar
 V. K. R. Kodur, B. Yu, and M. M. S. Dwaikat, “A simplified approach for predicting temperature in reinforced concrete members exposed to standard fire,” Fire Safety Journal, vol. 56, pp. 39–51, 2013. View at: Publisher Site  Google Scholar
 M. B. Dwaikat and V. K. R. Kodur, “Hydrothermal model for predicting fireinduced spalling in concrete structural systems,” Fire Safety Journal, vol. 44, no. 3, pp. 425–434, 2009. View at: Publisher Site  Google Scholar
 K. Kim, S. Jeon, J. Kim, and S. Yang, “An experimental study on thermal conductivity of concrete,” Cement and Concrete Research, vol. 33, no. 3, pp. 363–371, 2003. View at: Publisher Site  Google Scholar
 S. B. Tatro, “Significance of tests and properties of concrete and concretemaking materials,” Tech. Rep. STP 169D, ASTM International, 2006. View at: Google Scholar
 T. D. Brown and M. Y. Javaid, “The thermal conductivity of fresh concrete,” Matériaux et Constructions, vol. 3, no. 6, pp. 411–416, 1970. View at: Publisher Site  Google Scholar
 K.Y. Shin, S.B. Kim, J.H. Kim, M. Chung, and P.S. Jung, “Thermophysical properties and transient heat transfer of concrete at elevated temperatures,” Nuclear Engineering and Design, vol. 212, no. 1–3, pp. 233–241, 2002. View at: Publisher Site  Google Scholar
 T. S. Yun, Y. J. Jeong, T.S. Han, and K.S. Youm, “Evaluation of thermal conductivity for thermally insulated concretes,” Energy and Buildings, vol. 61, pp. 125–132, 2013. View at: Publisher Site  Google Scholar
 CEN: DAFT ENV, Euro code 3: Design of Steel Structures, 1993.
 V. K. R. Kodur and L. Phan, “Critical factors governing the fire performance of high strength concrete systems,” Fire Safety Journal, vol. 42, no. 67, pp. 482–488, 2007. View at: Publisher Site  Google Scholar
 H. K. Kim, J. H. Jeon, and H. K. Lee, “Workability, and mechanical, acoustic and thermal properties of lightweight aggregate concrete with a high volume of entrained air,” Construction and Building Materials, vol. 29, pp. 193–200, 2012. View at: Publisher Site  Google Scholar
 ASTM C177 Standard Test Method for SteadyState Heat Flux Measurements and Thermal Transmission Properties by Means of the GuardedHotPlate Apparatus, American Society for Testing and Materials, 2010.
 “Standard test method for thermal conductivity of refractories by hot wire (platinum resistance thermometer technique),” ASTM C1113, American Society for Testing and Materials, 2009. View at: Google Scholar
 T. L. Bergman, F. P. Incropera, A. S. Lavine, and D. P. DeWitt, Fundamentals of Heat and Mass Transfer, John Wiley & Sons, New York, NY, USA, 7th edition, 2011.
 H. E. Anderson, “Heat transfer and fire spread,” USDA Forest Service, Research Paper INT69, International Forest and Range Experiment Station, Forest Service, U.S. Department of Agriculture, Ogden, Utah, USA, 1969. View at: Google Scholar
 J. I. Ghojel, “A new approach to modeling heat transfer in compartment fires,” Fire Safety Journal, vol. 31, no. 3, pp. 227–237, 1998. View at: Publisher Site  Google Scholar
 M. B. Wong and J. I. Ghojel, “Sensitivity analysis of heat transfer formulations for insulated structural steel components,” Fire Safety Journal, vol. 38, no. 2, pp. 187–201, 2003. View at: Publisher Site  Google Scholar
 “ISO834 Fire resistance testelements of building construction,” International Standard ISO384, Geneva, Switzerland, 1975. View at: Google Scholar
 Korean Concrete Institute, Standard Specification of Concrete, Ministry of Land, Infrastructures and Transport, 2009.
 M. M. Elbadry and A. Ghali, “Temperature variations in concrete bridges,” Journal of Structural Engineering, vol. 109, no. 10, pp. 2355–2374, 1983. View at: Publisher Site  Google Scholar
 B. Bøhm and S. Hadvig, “Calculation of heat transfer in fire test furnaces with specific interest in the exposure of steel structures. A technical note,” Fire Safety Journal, vol. 5, no. 34, pp. 281–286, 1983. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Jung J. Kim et al. 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.