|
1.IntroductionMeasurement of absorptivity is important for analysis and modeling of laser-material interactions. Absorption of laser light depends on the laser parameters: irradiance, wavelength, polarization, and angle of incidence; and material properties including composition, temperature, surface roughness, oxide layers, and contamination. Theoretical and experimental studies of absorptivity and reflectance of metals have concentrated on perfectly pure, clean, and flat surfaces, free of oxide layers, unlike real life material processing applications where metal surfaces are rough to some extent and are contaminated with impurities and oxide layers. Furthermore, little information is available on absorptivity of alloys, which are far more commonly used than pure metals in real processing applications. Thus the published values of absorptivity for pure polished surfaces are inadequate for laser-material interaction models. In addition, detailed analysis and modeling of laser-material interactions require material absorptivity over a wide range of temperatures. The standard collection of textbooks1 has data mainly for ideal materials at room temperature. In engineering applications, metals are oxidized and have impurities, surface irregularities, and defects; hence, their absorptivity can be considerably different from handbooks’ values.2 Lasers have been widely used for many years in various material processing applications. In situations that do not require narrow linewidth laser beams where the laser acts as an energy source only, the more energy-efficient diode-array lasers can be used. Such lasers are now commercially available in a wide variety of wavelengths and formats and were used as radiation sources in these experiments. Direct absorptivity measurements are difficult due to the need to change the sample temperature in controlled way. Absorptivity is usually determined indirectly as the residual of the measured reflectivity.3 However, reflectance is difficult to measure accurately due to diffuse scattering and thermal radiation, especially for measurements at high temperatures. To avoid these problems, direct calorimetry was used here in which a thin sample is heated uniformly by a diode array and the temperature evolution is measured using thermocouples and (or) an infrared camera. The small thickness provides a temperature uniformity across the sample. Heat losses due to convection and thermal radiation are derived directly from the experimental data and, by removing this contribution, the temperature-dependent absorptivity at a wavelength of 780 nm is obtained. In addition, the measurements of thermal loss as a function of temperature, containing information about both convective loss and thermal radiative loss, are also discussed. The simple and compact setup, together with new way of data processing, is a convenient tool for temperature-dependent absorptivity measurements for engineering applications. This study focuses on the absorptivity of aluminum, while some data for titanium and steel are also presented. Extensive previous theoretical and experimental works on aluminum allow comparison of absorptivity for the engineering grade samples and ideal metals. The various mechanisms that increase absorptivity in real materials are discussed. Measurement results for steel and Ti alloys are also presented. It is demonstrated that the oxidation process is nonuniform and lead to variations in absorptivity across the illuminated area. 2.Experimental Setup and Measurement ProcedureFor direct, calorimetric absorptivity measurements, a thin sample is heated by the diode laser array and temperature evolution is measured using thermocouples, as illustrated schematically in Fig. 1. The sample is thermally isolated and all heating is produced by the absorption of laser light. The irradiation is approximately uniform over the sample. The small sample holder was made from fused silica, which has minimal absorption of laser light. Thermal measurements indicate a high degree of uniformity over the sample. At some moment, the radiation is turned off and the sample cools via thermal radiation and convection cooling. For the thin samples used, the thermal diffusion time is much shorter than the heating time, so the temperature is expected to be uniform throughout the sample. The temperature-dependent absorptivity for a sample with uniform temperature distribution (uniform irradiation and thin sample) is determined by the relation where is the sample thickness, is the density, is the temperature-dependent specific heat, is the incident intensity, and is the thermal loss. Assuming thermal loss is a function only of sample temperature, i.e., the loss channels are identical with and without laser illumination, and the measured during cooling when the laser is off also applies during heating, then the absorptivity of the sample can be determined with the minimal ambiguity.2Square samples of industrial grade aluminum, 1-mm thick and wide, were used in these measurements. The diode-array laser produced up to 3-kW power at 780 nm. Radiation was transported to the target through a fused silica “duct” that spatially smoothed the irradiance distribution, as shown in Fig. 2, by multiple reflections from the walls to a uniformity of . Irradiance on the sample was limited to 10 to to heat the target slowly enough for convective cooling to be in a steady-state regime and to improve the temporal resolution of temperature recording. The sample temperature was measured by up to five thermocouples distributed over the surface to verify the uniformity of the temperature, as shown in Fig. 3 where typical data from five different thermocouples are presented. The thermocouples used were type K with an expected measurement accuracy within 3°C and a response time well under a second. The parameters and are calculated from the time derivative of the measured sample temperature, which is difficult to do accurately when the data are noisy. To improve the numerical accuracy, data smoothing techniques described in Appendix A were applied. Measurements were made over a sufficiently broad range of sample temperatures; hence, it was necessary to account for the temperature dependence of the density and heat capacity of the metal sample. Data from Ref. 4, shown graphically in Fig. 4, were used. A polynomial curve fit for the density gives with , , and a correlation coefficient . Similarly, the fitted curve for specific heat capacity gives with , , , , and . The resulting thermal loss curves, , calculated for several Al targets heated with laser irradiances of 12 and in various heating experiments, shown in Fig. 5, demonstrate high repeatability and show that the thermal loss is a function only of the temperature that is independent of the heating rate. Smoothing and processing of measured temperatures does not result in a significant deviation in the results. Small differences between each curve may be attributed to environmental fluctuations in the experimental area. One example is the position and orientation of the thermocouple wires. Because the thermocouple wires were not attached and positioned exactly in the same way in every experiment, their influence on thermal convection will vary from experiment to experiment. The increase in variation at higher temperatures is also consistent with this explanation. An approximate deviation on the order of for would result in an error of 5% or less for the absorptivity—comparable to the accuracy of measurement. 3.Absorptivity of AluminumA comparison of absorption of Al alloys before and after correction for thermal loss, shown in Fig. 6, illustrates the importance of accounting for thermal losses. The loss-corrected results are found to be only weakly dependent on heating rate, justifying our approach. The absorptivity for Al is about 35% and is nearly constant over this temperature range. Heating to 450°C did not produce a significant change in sample appearance or morphology, and the repetitive measurements produced similar results. An estimated energy measurement error of 5% is expected due to thermopile heat sensor sensitivity and laser beam nonuniformity. The scheme for the removal of thermal losses may result in another 5% numerical error in absorptivity as stated above. Hence, one should expect a relative accuracy better than 10% for the derived absorptivity. Absorption data for Al presented in textbooks1 are measured for ideal, vacuum-deposited films. Measurements here on as received, engineering grade, plate aluminum samples show much higher absorption, which is usually attributed to surface roughness, the presence of oxide films, and absorbing impurities.2 The comparative magnitude of these effects and their wavelength dependence is discussed below. Aluminum is one of the so-called near-free-electron metals, for which electrons in the metal can be treated as a free-electron gas, with the addition of interband transitions.2 The Drude model provides a good description of the dielectric constant. Figure 7 shows the pure Al absorption as a function of wavelength for irradiation at normal incidence based on the experimental measurements of the refractive index. Absorption typically increases with temperature, except for the resonance region, simply because the free electrons gain kinetic energy while the phonon population grows, which increases the electron-phonon collision frequency. Aluminum has an absorption peak at 820 nm due to an interband transition, which is clearly shown by the blue line in Fig. 7. Absorption at 795 nm () is about 13%, and 2.5 times higher than the absorption at 1 μm (), while for 1.3-μm light () absorption falls to 2.7%. Absorption data for liquid Al (red dots in Fig. 7) show that the elimination of the metal structure makes absorptivity insensitive to wavelength. Experimental data presented for 1-μm radiation absorption of bulk Al2 demonstrates very different results, with absorption more than four times higher than that measured for the ideal sample.1 Also, the absorption is nearly temperature independent. In real situations, aluminum is coated with an oxide layer that is transparent up to the UV range but has a high refractive index, for a wavelength of 0.78 μm,2 that can affect the reflectivity. The complex amplitude of a wave reflected from the oxide film of thickness on the aluminum substrate is given by7 Here, is the oxide refractive index and is the complex refractive index of Al; is the amplitude of the wave reflected from the bulk aluminum oxide, given by , is the complex reflected amplitude from the bulk aluminum material, and is the angle of incidence. The complex value is . Reflectivity from bulk Al at this wavelength is and absorptivity . The value for is given by . For the p-polarized wave, all data are related to the magnetic field amplitude and the electric field amplitude for the s-polarized wave.In the case of normal incidence, the calculated absorptivity as a function of the oxide layer thickness is shown in Fig. 8, which shows that the presence of an oxide layer can substantially increase absorption. For a typical thickness of about 50 nm, absorption will be above 20%. Variability of the thickness will result in some average absorption value, but it will be higher than that predicted for pure aluminum. The reduction in absorptivity for a thicker oxide film, shown in Fig. 8, is the result of wave interference in the film. Thus, the presence of an oxide film increases the absorptivity but not sufficiently to explain the observed results. The effect of the oxide layer is more noticeable for shorter wavelength light. The additional absorption is explained by the surface roughness resulting in enhanced absorption and conversion of incident radiation to surface electromagnetic waves (plasmons). General calculations of absorption on a rough metal surface can be found, e.g., in Ref. 2, but are not helpful because detailed information about the surface structure is not usually available. Other possibilities are the surface defects and inclusions.2 The cold-rolled samples used for these measurements have parallel groves in the roll direction. Figure 9 shows the absorption measured with the electric field of the light parallel or orthogonal to the grooves. The higher absorption observed for the orthogonal field orientation is due to the fact that it is p-polarized with respect to the corrugated surface, thus having higher absorption than the parallel field which is s-polarized. These data directly demonstrate the effect of surface corrugations. Heating the samples to 400°C does not cause irreversible material changes in Al, as repetitive measurements on the same sample gave similar results and no visible changes were observed on the surface. However, when samples were heated close to their melting temperature, or , slight surface modifications were noticeable and an increase in absorptivity was found after cooling the sample (cf. Fig. 10). A probable source for this increase is the modification of the oxide layer. The increase in absorptivity after each thermal cycle is evidence for the modification of the oxide layer in Al. Figure 8 shows that an increase of the oxide layer thickness from 50 to 70 nm is sufficient to explain the results. Modification of surface roughness and defect generation also may be factors. 4.Absorptivity of Titanium and SteelMeasurements on Ti-6Al-4V samples were performed in a similar manner. The temperature distribution over the samples was uniform, and readings of all thermocouples were very similar. Uniform changes in the metal color caused by irradiation were observed, but these did not affect the measured absorptivity of approximately 0.5 and was only weakly dependent on the temperature as shown in Fig. 11. Absorptivity measured for 4130 steel on as-received samples was close to 0.5, nearly independent of temperature up to , and was generally increasing, though not uniformly, above this temperature. After the initial test, discolored spots could be seen on the sample surface, which may be attributed to oxidation (cf. Fig. 12). After cooling to room temperature, absorption measured on the already oxidized surface was about 0.8, higher than the as-received sample up to about 400°C and then was approximately the same at higher temperatures. Temperature variations over the sample’s surface resulting from changes in absorptivity at various locations on the surface due to oxidation were observed, as shown in Fig. 13. As a result, absorptivity is calculated using the mean temperature measured by the thermocouples and is shown in Fig. 14. Let us mention that the absorptivity increase measured during the initial illumination may be partially explained by the heat released during the oxidation process. As a comparison, we also tested stainless steel 316, which possesses higher resistance to corrosion and oxidation, under similar conditions and observed little absorptivity change in multiple heating and cooling cycles up to 500°C, as shown in Fig. 15. In particular, the sudden onset of surface modification at observed in steel 4130 is not present in stainless steel 316. Visible surface discoloration was observed in stainless steel 316 when samples were heated to over 600°C, likely as a result of surface oxidation. More detailed studies at higher temperatures are desired. Note that, to the best of our knowledge, most of the previous measurements were done using relatively small laser spots where the effect of nonuniform oxidation may be overlooked. The importance of oxidation and its nonuniformity makes the absorptivity sensitive to the environment. For example, airflow, which affects oxidation,2,8 can increase absorptivity. The oxidation process also releases additional energy which, in our experiments, can be counted as enhanced absorptivity. 5.Thermal LossesThermal loss consists of two primary components: convective cooling and thermal radiation. At low temperatures, radiative loss is small and can be disregarded. The heated surface induces convective flow near it, which cools the sample. Convective flow near a vertical wall is characterized by the Grashof number, Gr, Here is the gravitational acceleration, is the air kinematic viscosity, is the sample size, and is the difference between the sample temperature and the ambient temperature. For the developed convection, . In our case, for , , and , . The cooling flux is given by the expression where is the air heat conduction and Nu is the Nusselt number. For our range of Gr, .8 Here, Pr is the Prandtl number, the ratio of viscosity to the thermal diffusivity, which for air is . Note that viscosity and thermal conduction vary with temperatureThe above value of the Nusselt number was obtained for experiments with the temperature independent viscosity and thermal conduction. Our data are fitted using the same Eq. (4) with viscosity and thermal conduction taken at the sample temperature and an adjustable numerical coefficient . Thus, the thermal flux, after averaging over the plate, can be represented as Here, the factor 2 accounts for convection losses from each side of the sample. At room temperature, and . ForComparison of Eq. (7) with and experimental data for Al, Ti, and steel samples are shown in Fig. 16. At low temperatures, Eq. (7) provides a good description of convective losses and the losses for Al and Ti are similar. The convection losses are independent of materials. At higher temperatures, losses are expected to exceed those given by Eq. (7) due to radiative losses and will be discussed later. The above estimates were made for steady-state convection. It can be different in a nonstationary regime. Let us estimate the time for a steady-state flow onset. Heat diffuses from the sample surface to the air over the distance , inducing a convective flow with the velocity . In a convective flow, buoyancy is compensated by viscosity In a steady state, the vertical temperature convection is compensated by the lateral thermal transport Combining the last two equations, one gets for the onset time For , , and in air, the onset time is small, . Practically, the convection is expected to be in a stationary regime. This conclusion is consistent with the absorptivity results that were measured for different pump intensities. Now, with an estimate for the convective loss, one can subtract it from and obtain the contribution from radiative losses. Convective losses are independent of the material but sensitive to the sample geometry. Radiative losses are independent of geometry but are sensitive to the material. Data for Al and Ti will be presented. The intensity of thermal radiation according to Kirchhoff’s law can be related to the blackbody radiation and reflectivity 8 Total emissivity is obtained after averaging over all angles and wavelengths The emissivity is related to the absorptivity and sensitive to the surface quality and the change from sample to sample. As a result, there is considerable variability in emissivity data.8 In our specific case, the emissivity can be estimated by The graph of emissivity for Al, Ti, and steel calculated from Eq. (13) is shown in Fig. 17. The data have been averaged over a few samples. It should be noted that the data variation from sample to sample is not small—larger than 20%. Nevertheless, the result is consistent with Fig. 16—steel radiates more than Ti and Ti more than Al. Let us estimate the radiative losses for Al. The blackbody radiation losses at 500°C are about . If we use an emissivity of about 0.25 from Fig. 17, the total radiative losses at 500°C would be of the total losses of from Fig. 5. 6.ConclusionTo the best of our knowledge, the first measurements of absorptivity of 0.8-μm light by engineering grade metals over a broad temperature range are presented. The measurements were performed utilizing a simple calorimetric setup. Use of laser diodes as a light source makes possible the uniform irradiation of large samples. The direct measurement approach simplifies the data processing and allows for accurate determination of absorptivity. Thermal losses due to convection and thermal radiation are shown to play an important role in energy balance and must be taken into account. A method to subtract thermal losses to derive absorptivity at elevated temperatures is demonstrated. It was found that up to a temperature of , the absorptivity of Al 6061 varies weakly with temperature, consistent with measurements at 1.3 μm (Ref. 3) and almost two times higher than absorptivity at . This result is expected due to the interband transition peak near 825 nm. The absorptivity of Ti-6Al-4V is also insensitive to the temperature and comparable with absorptivity at 1 μm. Measurements of steel 4130 show that oxidation can affect absorptivity. As oxidation occurs nonuniformly over the sample, the local absorptivity also varies and one can only talk about average absorptivity values. Evaluation of the thermal losses, which may be derived from our measurements, also provides useful information about convective and radiative losses. Finally, methods presented here can be applied in measurements of absorptivity for a variety of materials using laser diodes or other laser sources with different wavelengths. AppendicesAppendix ADifferentiation of noisy experimental data is typically unreliable, so we describe two smoothing techniques: least-squares polynomial fit and Savitzky–Golay derivative filter. Since the temperature signal is a well-behaved monotonic function during heating and cooling, finding least-squares polynomial fits for the recorded data over the heating and cooling periods is simple and produces fairly good approximations. We found tenth-order polynomial fits to the measured data followed by differentiation in some of our thermal calculations. In other cases, we utilized the method of least squares described by Savitzky and Golay,9 where a given number of data points are fitted to a low-order polynomial via a least-squares procedure to obtain the value of the smoothing interpolant at the center of the interval. The least-squares calculations may be carried out by convolution of the data points with properly chosen sets of integers. This method assumes that the data are evenly spaced and without breaks. Due to technical issues, our recorded time series data usually contain uneven spacing and short breaks where no useful data are recorded. This poses no problems for the polynomial fit algorithm but offers some complications when one tries to directly apply Savitzky–Golay filters. Because in most cases the simplistic polynomial fit agrees very well with the approximation found via the more sophisticated Savitzky–Golay method, as seen in Fig. 18, these methods were applied interchangeably and without distinction where appropriate. AcknowledgmentsThis work was performed under the auspices of the US Department of Energy by Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344. ReferencesE. D. Palik, Handbook of Optical Constants of Solids, Academic Press, San Diego, California
(1998). Google Scholar
A. M. Prokhorovet al., Laser Heating of Metals, Adam Higler, Bristol
(1990). Google Scholar
R. K. FreemanF. A. RigbyN. Morley,
“Temperature-dependent reflectance of plated metals and composite materials under laser irradiation,”
J. Thermophys. Heat Transfer, 14
(3), 305
–312
(2000). http://dx.doi.org/10.2514/2.6546 JTHTEO 0887-8722 Google Scholar
K. G. Mills, Recommended Values of Thermophysical Properties for Selected Commercial Alloys, Woodhead Publishing, Cambridge, England
(2002). Google Scholar
S. KrishnanP. C. Nordine,
“Optical properties of liquid aluminum in the energy range 1.2–3.5 eV,”
Phys. Rev. B, 47
(18), 11780
–11787
(1993). http://dx.doi.org/10.1103/PhysRevB.47.11780 PRBMDO 0163-1829 Google Scholar
S. KrishnanP. C. Nordine,
“Analysis of the optical properties of liquid aluminum,”
Phys. Rev. B, 48
(6), 4130
–4131
(1993). http://dx.doi.org/10.1103/PhysRevB.48.4130 PRBMDO 0163-1829 Google Scholar
L. LandauE. Lifshitz, Electrodynamics of Continuous Media, Pergamon Press, Oxford, UK
(1960). Google Scholar
E. EckertR. Drake, Analysis of Heat and Mass Transfer, McGraw Hill, New York
(1972). Google Scholar
A. SavitzkyM. Golay,
“Smoothing and differentiation of data by simplified least squares procedures,”
Anal. Chem., 36
(8), 1627
–1639
(1964). http://dx.doi.org/10.1021/ac60214a047 ANCHAM 0003-2700 Google Scholar
BiographyAlexander M. Rubenchik received his PhD degree in theoretical physics from the Institute of Nuclear Physics, Novosibirsk, in 1974, and the Dr. Sci. degree from the Space Research Institute, Moscow, in 1983. Since 1992, he is a physicist at Lawrence Livermore National Laboratory at Livermore. His main scientific interests are in physics of laser-matter interactions, nonlinear optics including the studies of optical damage, laser physics and plasma physics. He has over 350 publications, 4 100 R&D awards. He is a member of OSA and DEPS. Sheldon S. Q. Wu received his PhD degree in electrical engineering with specialization in photonics from the University of California at San Diego in 2009. Since then, he has conducted ongoing research in the NIF & Photon Science Directorate at Lawrence Livermore National Laboratory. He has made technical contributions in the fields of diode-pumped alkali lasers and precision Compton scattering light sources for nuclear photonics applications. His main scientific interests are laser physics, laser-matter interactions, classical and quantum electrodynamics, high field physics, and quantum optics. |