Research Article  Open Access
Temperature Control in Spark Plasma Sintering: An FEM Approach
Abstract
Powder consolidation assisted by pulsed current and uniaxial pressure, namely, Spark Plasma Sintering (SPS), is increasingly popular. One limitation however lies in the difficulty of controlling the sample temperature during compaction. The aim of this work is to present a computational method for the assembly temperature based on the finite elements method (FEM). Computed temperatures have been compared with experimental data for three different dies filled with three materials with different electrical conductivities (TiAl, SiC, ). The results obtained are encouraging: the difference between computed and experimental values is less than 5%. This allows thinking about this FEM approach as a predictive tool for selecting the right control temperatures in the SPS machine.
1. Introduction
1.1. The SPS Process
SPS (Spark Plasma Sintering) consists in applying a pulsed current and a uniaxial pressure, and thus allows high compaction rates and level, particularly for materials that are difficult to densify. This also leads to the formation of original microstructures and materials with outstanding or new properties. In a recent review of SPS, Munir et al. [1] reported a dramatic increase in the number of papers on this subject or on materials sintered with it. SPS also raises numerous fundamental questions: which physical mass transport mechanisms control particle powder welding? Also investigated in the literature by Munir et al. [1] and Conrad et al. [2] are the issues of plastic deformation, electroplastic deformation, electromigration, spark plasma, or electric arcs at the particle. Surely, the right explanation will depend on the kind of powder involved, particularly electrical conductivity compared to that of the die: densification is expected to take a different form depending on whether the current flows through the powder or through the die. Note that it is the unconfirmed hypothesis of plasma formation which leads to the commercial name SPS. Pulsed Electric Current Sintering (PECS) which is sometimes found in the literature would seem more appropriate, but SPS has apparently been the accepted name.
1.2. Temperature Control
Figure 1 shows a schematic drawing of SPS apparatus used in the present study. Depending on die and sample geometries, and on electrical resistivities, current lines will spread in different ways. Therefore, highly different temperature fields can be generated in the setup. In particular, in the powder volume, temperature can be inhomogeneous, especially if the volume is high or complex. Can these internal temperatures be known? This paper addresses this issue. Pyrometric measurements yield the temperature of the external die wall. By placing thermocouples inside the holes at different depths, some local internal die temperatures can also be measured. However, the measurement inside the sample remains problematic: inserting a thermocouple is a local and destructive procedure. Indirect measurement could be considered by observing phase transformations. In this case the results would probably lack accuracy and it would only be possible to determine whether a temperature corresponding to a known phase transformation temperature has been exceeded. Three other limitations can be seen: (i) phase transformation temperatures are known for equilibrium configurations, whereas typical SPS sintering is a fast process occurring out of equilibrium; (ii) this approach requires materials with phase transformations; (iii) as with thermocouples, it is destructive. To address the key issue of sample temperature control, a solution might be to compute the temperatures of interest which could then be compared to those measured, thus leading to a predictive tool for selecting the right control temperature of the SPS machine, for any geometry or powder. Based on this approach, this paper aims at presenting our first results in the case of three different materials and geometries.
2. State of the Art
Determining the temperature field inside an SPS setup is not a new concern. Through analytical calculations based on the heat equation, Yucheng and Zhengyi [3] found a difference of 345°C between the edge and the centre of a TiB_{2} sample, compared with an experimental difference of 450°C. The order of magnitude is correct within the approximation made, and the difference between calculations and measurements can be accounted for the simplifications introduced and by the fact that unlike experimental conditions the calculated results correspond to stationary state. Matsugi et al. [4] developed a finite difference method (FDM) for coupling heat and electricity equations. From an experimental point of view, they placed thermocouples into punches, die, and samples with a conductive powder (Ti) and an insolating one (Al_{2}O_{3}). In these experiments, the stationary state was reached. The difference between calculations and experiments did not exceed 10°. Zavaliangos et al. [5] employed an FEM approach (ABAQUS code) for calculating the temperature of a graphite sample inside a graphite die. Taking into account thermal resistances between die and sample, a good fit was achieved between calculated and experimental differences of temperatures T between the external wall die and the centre of the sample. AnselmiTamburini et al. [6] also made use of FEM calculations (CFD code) disregarding in this case thermal resistances. With an Al_{2}O_{3} sample, the calculated T was 27% lower than the experimental temperature but absolute values were significantly different: the calculated ones were 200°C lower than the measured ones. Lastly, Vanmeensel et al. [7] used another FEM code (ANSYS) taking into account the contact resistances on the basis of many specific experiments. A good agreement was obtained between calculations and experiments for a graphite sample in the graphite die: the calculated T was 4% lower than the experimental one.
3. Methodology
3.1. The SPS Machine
A Dr. Sinter 2080 SPS apparatus was used in this work (Sumitomo Coal Mining Co., Japan), which is installed in the “Plateforme Nationale de Frittage Flash/CNRS” in Toulouse (PNF2/CNRSMHT, Paul Sabatier University, Toulouse, France). Pulse duration was 3.3 ms, and a pulse sequence consisting of twelve pulses followed by two periods of zero current was selected.
3.2. Temperature Measurements
To test our calculations in various configurations, specific SPS experiments were carried out with three cylindrical geometries and three materials (Figure 2): a conducting material TiAl (electrical resistance ρ_{e} = 610^{7} Ωm at 300 K), an insulating material Al_{2}O_{3} (ρ_{e} = 10^{8} Ωm at 300 K), and a material with an intermediate conductivity SiC (ρ_{e} = 2 Ωm at 300 K), referred to as “TiAl”, “Al_{2}O_{3}”, and “SiC” geometries, respectively. A point in space is defined by the radial distance r (e.g., the vertical edge of the TiAl sample is mm) and the longitudinal distance (height) , the origin of being taken as half the height of the assembly, that is, at the level of the axial symmetry plane. For the TiAl geometry, temperature measurements were obtained by microstructure characterization of the sample, using scanning electron microscope (SEM) observations: this intermetallic alloy exhibits a phase transformation which can be identified and correlated to transus temperature. Note that with this geometry, temperature was controlled during sintering by a pyrometer focused at at the back of a 3 mm depth hole, that is, mm. For the Al_{2}O_{3} geometry, local temperatures were measured by a Kthermocouple placed into similar holes. For , thermocouples were placed at , 11, 16, and 22 mm, respectively. For mm, measurements were also carried out above ( mm) and below ( mm). As far as the SiC geometry is concerned, measurements were performed at with thermocouples at and 22 mm, but also inside the sample with a specific platinum thermocouple at .
(a) TiAl sample
(b) Al2O3 sample
(c) SiC sample
3.3. Calculation
3.3.1. Calculation Principles
For each point of the solids and each time t, the local equations of electricity (a) and heat (b) had to be verified:
where is the electrical potential, T the temperature, and λ,ρ,C_{p} and σ are the thermal conductivity, mass density, thermal capacity, and electrical conductivity of the material, respectively. These four parameters, related to the three materials studied and to inconel and graphite, are given in Table 1 as a function of temperature.
(a) Materials data.  
 
(b)  

From a mathematical view point, these two secondorder differential equations have to be coupled each time with respect to the three coordinates. The heat source (right term of (b), does present an electrical origin as described by (a). Solving this equation system at every time and at every points of the assembly requires using a numerical approach based on finite differences or finite elements. Here, the finite elements method has been retained.
3.3.2. Simplifying Hypothesis
The general form of the local electricity equation is where ρ_{q} is the charge density. When time increases, temperature equally rises just like the charge density, especially in the SiC semiconductor. However, one has to consider that ρ_{q} is a constant value even in this case. Another simplification is that all data in Table 1 correspond to massive materials, not powders. This is not problematic for Al_{2}O_{3} and SiC since the corresponding temperature measurements have been achieved with massive samples placed in the die. For TiAl, the situation is different because powder has been compacted during experiment. Thus only those temperatures reached at the end of compaction have been taken into account for the purpose of comparison with calculations. As for TiAl, the material has been assumed to be stable without any phase transformation, and this is contrary to reality. Also transformation heats have been disregarded and the temperature dependence of λ,ρ,C_{p}, and σ has been determined as if temperature had always remained below the transus temperature. On the other hand, for both Al_{2}O_{3} and SiC, phase transitions do not exist in the temperature range of our experiments.
Electrical and thermal contact resistances, especially between die and punches, have been ignored. It has been shown by Zavaliangos et al. [5] and Vanmeensel et al. [7], that they can play an important role but they are highly dependent on surface states or applied pressures. For example, AnselmiTamburini et al. [6] have demonstrated that if the applied pressure is higher than 50 MPa (our situation), electrical resistance of the graphite assembly is similar to the one of a graphite block with the same geometry. Moreover, incorporating these resistances would require specific measurements for each experiment to be simulated. As our goal has been to get a simplistic predictive tool, we have disregarded them and preferred to investigate a posteriori the influence of this choice. As will be shown in Section 5 giving calculated and measured temperatures, this seems reasonable. The same goal prompted us to ignore the 0.2 mm width graphite paper located on the internal die wall to facilitate sample removal. In the particular case of Al_{2}O_{3} and SiC geometries, we ignored the influence of thermocouples and the holes drilled in the sample.
Voltage is made up of pulses but, as shown by AnselmiTamburini et al. [6], the period involved is long enough for modelling this signal by a continuous function. In our experiments, this function starts from 0 at time zero and increases up to a stage with permanent fluctuations imposed by the machine for following the temperature ramp imposed. In our calculations, a direct and constant voltage has been used as time starts with a value corresponding to the experimental value reached at the final stage of the experiment.
3.3.3. Limit Conditions
With respect to the heat flux, the following limit conditions have been imposed:
(i)for vertical walls: radiative heat flux with σ_{s} = the StefanBoltzmann’s constant, ε = graphite emissivity (ε= 0.8 from [7]), T_{e} = emission surface temperature (graphite), and T_{a} = absorption surface temperature (chamber walls) = 300 K(ii)for horizontal walls:(a)conductoconvective heat flux for the two Inconel surfaces in contact with cooling water, that is, with h_{c} = conductoconvective coefficient ( Wm^{2}K^{1} from [7]), T_{i} = Inconel horizontal wall temperature, and T_{w} = water temperature = 300 K.(b)Taking into account the assembly symmetry, it has been considered that for opposite horizontal surfaces, emitted and absorbed fluxes are counterbalanced.3.3.4. MEF Calculations
The COMSOL code featuring a tetrahedral mesh has been used. This kind of mesh is not as efficient as a cube one. However, it provided the required calculation convergence. The initial condition for temperature is obviously K (uniformly). Voltage is only known at the two Inconel extreme surfaces: 0 V on the bottom surface and the imposed voltage on the top surface. Thus, different voltage distributions within the assembly were tested:
(i)linear variation between bottom and top, (ii)linear variation of 0.5 V in the two inconel rams plus linear variation from 0.5 to 4.5 between them, (iii)linear variation of 0.5 V in the two inconel rams plus variations for each graphite pieces taking into account their size differences, in consistency with Ohm law.As a simple linear distribution was found to yield the same results as the other distribution, it was selected. The amplitude of voltage between the bottom and top imposed as an initial condition is derived from experiment as shown below.
To take into account time dependence (∂T/∂t term in heat equation), calculations have been made step by step, each step corresponding to a time increment. Therefore, a calculation corresponds to the progressive temperature rise at each point of the assembly. A sufficiently high step number permits to reach a stationary state.
4. Results
4.1. Measured temperatures
4.1.1. TiAl
A Ti_{49}Al_{47}Nb_{2}Cr_{2} powder prepared by gas atomisation has been used. The SPS conditions (temperature and pressure ramps) selected for compacting this alloy terminate in a 2minute holding time at 1190°C (this is the control temperature measured by the pyrometer) [8]. The corresponding voltage during this final stage was 4.2 V. Figure 3(b) shows the sample microstructure at mm, that is, at half the distance between the centre and the vertical edge of the TiAl billet ( mm). This microstructure is made of small dark contrasted grains with a structure and bigger grains with a clear streaked contrast, the socalled lamellar grains. This kind of lamellar grains corresponds to a coherent juxtaposition of very thin lamellae, either ordered phase or phase [9, 10]. Such a microstructure with dominant lamellar grains is obtained for temperatures slightly lower than the transus temperature, that is, 1330°C ± 10 [11]. Figure 3(a) shows the microstructure of the same sample at . It appears that grains are very large and lamellar. This microstructure could only be obtained if temperature is higher than the transus temperature. So, for a sample sintered at 1190°C, the temperature reached is lower than 1330°C at mm and higher than 1330°C at : 1330°C is attained somewhere between these two places, that is, around r = 4.5 mm. In this case, the difference between the sample temperature at mm and the temperature measured where the pyrometer is focused ( mm) is 140°C ± 10.
(a) 1190°C and 𝑟 = 0
(b) 1190°C and r = 9 mm
(c) 1175°C and 𝑟 = 0
In another experiment for which sintering temperature was 1175°C (Figure 3(c)), the microstructure at was the same as at 1190°C at mm. In other words, for a compaction with a temperature 15° lower, the temperature at the centre has not reached 1330°C. A 15° rise of the external temperature seems thus to produce a rise of the same order of the centre temperature. From these findings, it follows that the temperature difference between the centre of the sample and its vertical edges is approximately = 30°. This is consistent with what has been observed by Couret et al. [12] for homogeneous microstructures along the sample radius with a sintering temperature in the region.
4.1.2. Al_{2}O_{3}
Figure 4(a) shows the temperatures measured in the graphite die at and , 16, 11, and 6 mm, respectively, during an experiment where temperature has been increased in five stages. The temperature rises from the external wall of the die to the internal one. Figure 4(b) shows the same kind of measurements with thermocouples placed at the same radial distance ( mm) and at different heights: , +5 and –10 mm, respectively. The temperature decreases symmetrically with respect to the symmetrical plane of the assembly (). In both cases, the temperature ramp was the same, the controlling temperature was measured by the thermocouple situated at mm and the five voltages reached were successively 1.4 V, 2.1 V, 2.8 V, 3.5 V, and 4.4 V.
(a) At different radial distances ( ℎ = 0 )
(b) At different heights ( 𝑟 = 2 2 mm)
4.1.3. SiC
Measurements in SiC (Figure 5) were made in the graphite die at and at different depths ( and 22 mm) but also in the material itself (). The four short stages during temperature rise are associated with four voltages: 2.4, 3.0, 3.6, and 4.2 V, respectively. As in the Al_{2}O_{3} cases, the temperature rises from the external wall of the die to the internal one.
4.2. Calculated Temperatures
4.2.1. TiAl
The value of the voltage imposed as the initial condition in our calculations is 4.2 V, that is, the value reached at holding time during the sintering at 1190°C (Section 4.1.1). Figure 6 shows the radial variation of temperature in the horizontal symmetry plane for one over ten calculated steps, which correspond to the time and temperature increments (see Section 3.3.4). It shows that the number of steps was high enough to reach the stationary state. The step which yielded a temperature of 1190°C was selected. Note that 1190°C was reached before the last step, the stationary value being 1581°C. At this step, the temperature field is given in Figure 7(a). In Figure 8, plot (a) is derived from this field. It represents the radial variation of temperature at , where the pyrometer is placed. From this plot, the calculated temperature at and 9 mm can be compared to the experimental values given in Section 4.1.1.
(a) For a 4.2 V voltage, at the step giving 1190°C where the pyrometer focuses
(b) For 2.85 V which gives 1190°C as a stationary temperature, at the same place
On the other hand, by successive trials, the voltage needed in calculations to get 1190°C as a stationary temperature at the place of interest was determined to be 2.85 V. The corresponding temperature field is given in Figure 7(b), and plot (b) in Figure 8 gives the corresponding function at h = 0.
4.2.2. Al_{2}O_{3} and SiC
The same procedure described in the previous section was used for Al_{2}O_{3} and SiC geometries: for each voltage associated with each temperature stage (Sections 4.1.2 and 4.1.3), calculations were made up to the stationary state. For each one of these 9 calculations (5 for Al_{2}O_{3} and 4 for SiC), the step that yields a calculated temperature identical to the experimental temperature at and mm was selected. Tables 2(b) and 2(c) commented in next section list the temperatures obtained at this selected step.
(a)  
 
(b)  
 
(c)  

5. Comparison and Discussion
Tables 2(a), 2(b), and 2(c) give a comparison of calculated and measured temperatures. From these tables, it appears that calculated temperatures are relatively close to experimental ones: for TiAl, Al_{2}O_{3}, and SiC, the average differences between calculated and experimental values are %, % (standard deviation: %) and % (standard deviation: 2.8%), respectively. Therefore, the different simplifications assumed (Section 3.3.2) have a limited effect and can be justified a posteriori. It can be pointed out that as for AnselmiTamburini et al. [6] who did not take into account thermal resistances or Vanmeensel et al. [7] who factored them in, differences are systematically negative (except for one value among the 35 comparisons). Whatever the electrical resistivity of the material, lower (TiAl) or higher (Al_{2}O_{3}) than the graphite die, calculations maintain the same result quality. Besides, irrespective of the assembly shape (narrow part between punches and spacers or not, see Figures 2(a) and 2(b)), calculations give similar results in both cases.
Generally, SPS compactions are not achieved under thermodynamic equilibrium conditions. In the three situations investigated, temperature stages do not correspond to stationary states, that is, they are too short. For example, in the TiAl case, with the 4.2 V voltage maintained for 2minutes holding time, calculations give a stationary temperature of 1581°C, that is, 391°C above the experimental temperature (1190°C). On the contrary, a 2.85 V voltage is needed in the calculation to get a 1190°C stationary temperature. In this case, the temperature difference between the centre of the sample and the hole where the pyrometer is measering is unsurprisingly lower: 80° in comparison, 112° (Figure 8) and the sample temperature is less satisfactorily reflected in the calculation. The best agreement with experimental results is achieved with calculations performed under nonequilibrium conditions, a situation really reflecting what happens during the SPS process.
As a matter of fact, its nonequilibrium feature led to the choice of our method which consists in calculating step by step the temperature for a given voltage deduced from experiments and to select the step providing a calculated temperature equal to the measured temperature at a particular place. Subsequently, it is possible to compare the experimental and calculated temperatures where thermocouples were positioned, and to evaluate the temperature from the calculations at any place. Moreover, some other data, as the field of current intensity, can be extracted from modelling.
The present modelling is based on physics of continuous media, and obviously situated at a mesoscopic scale with a randomisation of the microscopic mechanisms. For example, this accounts for the “Al_{2}O_{3}” geometry, for which the discrepancies are lower in the symmetrical plane (), where the microscopic state is probably close to the average. This mesoscopic modelling allows a correlation at a macroscopic scale and thus its utilisation as a predictive tool for the industrial mastery of the process. It also provides comparative data for calculations at the microscopic and even nanoscopic scales (ab initio atomistic calculations…).
At least, it should be pointed out that this method does not allow for a direct prediction of the sample temperature without a first experiment allowing the determination of the tension. To achieve this goal, new calculations are being made using as input a time increasing voltage similar to the experimental one.
6. Conclusion
In this paper, finite elements calculations of the temperature field in samples and dies during a Spark Plasma Sintering experiment have been presented. These calculations were performed using a conducting material (TiAl), an insulating material (Al_{2}O_{3}), and a material with intermediate conductivity (SiC) and with various geometries of the encapsulating system. This method is based on the selection of the step reflecting the external temperature which is measured by an external pyrometer. The results were validated by comparison with experimental temperatures either measured by thermocouples or determined by the activation of a phase transition.
A good consistency, with a discrepancy lower than 5%, was found between the calculated and measured temperatures irrespective of place, assembly geometry, and conductivity of the material. The temperature difference between the external wall of the die and the sample is quite elevated, higher than 100°C. Since the SPS apparatus allows only the measurement of this external temperature, these calculations appear to be the only way to determine in advance the right external SPS temperature corresponding to the sample temperature desired. They should also support the evaluation of the temperature gradient in samples with a complex geometry.
Acknowledgments
The authors would like to thank Marc Thomas with DMMPONERA (Châtillon, France) for SEM images that allowed the temperature determination in TiAl samples and for fruitful discussions. They are also thankful for the PNF2 (Plateforme Nationale de Frittage Flash/CNRS in Toulouse, France) which provides SPS facilities.
References
 Z. A. Munir, U. AnselmiTamburini, and M. Ohyanagi, “The effect of electric field and pressure on the synthesis and consolidation of materials: a review of the spark plasma sintering method,” Journal of Materials Science, vol. 41, no. 3, pp. 763–777, 2006. View at: Publisher Site  Google Scholar
 H. Conrad, A. F. Sprecher, W. D. Cao, and X. P. Lu, “Electroplasticity—the effect of electricity on the mechanical properties of metals,” Journal of Materials Science, vol. 42, no. 9, pp. 28–33, 1990. View at: Google Scholar
 W. Yucheng and F. Zhengyi, “Study of temperature field in spark plasma sintering,” Materials Science and Engineering B, vol. 90, no. 12, pp. 34–37, 2002. View at: Publisher Site  Google Scholar
 K. Matsugi, H. Kuramoto, T. Hatayama, and O. Yanagisawa, “temperature distribution at steady state under constant current discharge in spark sintering process of Ti and ${\text{Al}}_{2}{\text{O}}_{3}$ powders,” Journal of Materials Processing Technology, vol. 146, no. 2, pp. 274–281, 2004. View at: Publisher Site  Google Scholar
 A. Zavaliangos, J. Zhang, M. Krammer, and J. R. Groza, “Temperature evolution during field activated sintering,” Materials Science and Engineering A, vol. 379, no. 12, pp. 218–228, 2004. View at: Publisher Site  Google Scholar
 U. AnselmiTamburini, S. Gennari, J. E. Garay, and Z. A. Munir, “Fundamental investigations on the spark plasma sintering/synthesis process: II. Modeling of current and temperature distributions,” Materials Science and Engineering A, vol. 394, no. 12, pp. 139–148, 2005. View at: Publisher Site  Google Scholar
 K. Vanmeensel, A. Laptev, J. Hennicke, J. Vleugels, and O. van der Biest, “Modelling of the temperature distribution during field assisted sintering,” Acta Materialia, vol. 53, no. 16, pp. 4379–4388, 2005. View at: Publisher Site  Google Scholar
 G. Molénat, M. Thomas, J. Galy, and A. Couret, “Application of spark plasma sintering to titanium aluminide alloys,” Advanced Engineering Materials, vol. 9, no. 8, pp. 667–669, 2007. View at: Publisher Site  Google Scholar
 H. Inui, M. H. Oh, A. Nakamura, and M. Yamaguchi, “Ordered domains in TiAl coexisting with ${\text{Ti}}_{3}\text{Al}$ in the lamellar structure of Tirich TiAl compounds,” Philosophical Magazine A, vol. 66, no. 4, pp. 539–555, 1992. View at: Google Scholar
 S. Zghal, S. Naka, and A. Couret, “Quantitative TEM analysis of the lamellar microstructure in TiAl based alloys,” Acta Materialia, vol. 45, no. 7, pp. 3005–3015, 1997. View at: Publisher Site  Google Scholar
 I. Ohnuma, Y. Fujita, H. Mitsui, K. Ishikawa, R. Kainuma, and K. Ishida, “Phase equilibria in the TiAl binary system,” Acta Materialia, vol. 48, no. 12, pp. 3113–3123, 2000. View at: Publisher Site  Google Scholar
 A. Couret, G. Molénat, J. Galy, and M. Thomas, “Microstructures and mechanical properties of TiAl alloys consolidated by spark plasma sintering,” Intermetallics, vol. 16, no. 9, pp. 1134–1141, 2008. View at: Publisher Site  Google Scholar
 http://www.hightempmetals.com/techdata/hitempInconel617data.php, 2008.
 http://www.hightempmetals.com/techdata/hitempInconel600data.php, 2008.
 M. PerezBravo, Nuevos intermetálicos gammaTiAl: estudio y optimización microestructural para su aplicación en aeronáutica, Ph.D. thesis, Universidad del Pais Vasco (UPV/EHU), 2005.
 http://www.matweb.com/search/SpecificMaterial.asp?bassnum=BCSC2A, 2008.
 http://www.matweb.com/search/SpecificMaterial.asp?bassnum=BA1A, 2008.
 D. Lundström, B. Karlsson, and M. Gustavsson, “Anisotropy in thermal transport properties of cast $\gamma $TiAl alloys,” Zeitschrift für Metallkunde, vol. 92, no. 11, pp. 1203–1212, 2001. View at: Google Scholar
 D. Veeraraghavan, U. Pilchowski, B. Natarajan, and V. K. Vasudevan, “Phase equilibria and transformations in Ti(25–52) at.% Al alloys studied by electrical resistivity measurements,” Acta Materialia, vol. 46, no. 2, pp. 405–421, 1998. View at: Google Scholar
Copyright
Copyright © 2010 G. Molénat 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.