8 March 2013 Temperature rise in objects due to optical focused beam through atmospheric turbulence near ground and ocean surface
Author Affiliations +
Optical Engineering, 52(3), 036001 (2013). doi:10.1117/1.OE.52.3.036001
Abstract
We consider an optical beam propagated through the atmosphere and incident on an object causing a temperature rise. In clear air, the physical characteristics of the optical beam transmitted to the object surface are influenced primarily by the effect of atmospheric turbulence, which can be significant near the ground or ocean surface. We use a statistical model to quantify the expected power transfer through turbulent atmosphere and provide guidance toward the threshold of thermal blooming for the considered scenarios. The bulk thermal characteristics of the materials considered are used in a thermal diffusion model to determine the net temperature rise at the object surface due to the incident optical beam. These results of the study are presented in graphical form and are of particular interest to operators of high power laser systems operating over large distances through the atmosphere. Numerical examples include a CO 2 laser (λ=10.6  μm ) with: aperture size of 5 cm, varied pulse duration, and propagation distance of 0.5 km incident on 0.1-mm copper, 10-mm polyimide, 1-mm water, and 10-mm glass/resin composite targets. To assess the effect of near ground/ocean laser propagation, we compare turbulent (of varying degrees) and nonturbulent atmosphere.
Stoneback, Ishimaru, Reinhardt, and Kuga: Temperature rise in objects due to optical focused beam through atmospheric turbulence near ground and ocean surface

1.

Introduction

When a high-power optical beam is propagated through atmospheric turbulence, and is incident on an object, power is transferred, causing the temperature of the object to rise. Recent demonstrations by the Navy show a multikilowatt beam igniting the engine housing of a small boat from over 1-km distance in a maritime environment.1 Interest in this technology motivates the investigation of turbulence effects on beam-wave propagation and its impact on target heating. Laser beam propagation through the atmosphere in the presence of turbulence has been studied extensively.2,3 The Navy has developed high energy laser code for atmospheric propagation, a 3-D code to model high energy laser propagation, and has published many papers based on the results.45.6 The basic problem of laser-induced heating of a variety of solid objects has been studied.78.9 Many studies of laser heating until now have been performed in the context of laser machining, without regard to a large volume of turbulent media between the transmitter and solid object. While these studies give useful results for the laser power and dwell times common in laser machining, these parameters do not generally translate well to the problem of a laser beam propagated over a large distance through the atmosphere. In these cases, spot sizes are generally larger and dwell times generally longer. This paper provides a prediction of temperature rise due to a high power optical emitter operating over a large distance in a variety of plausible atmospheric conditions. To predict heating of solid objects due to the incident laser, beam-wave propagation, optical absorption, and heat diffusion are simultaneously considered. For a beam wave in turbulence, a statistical model is used to obtain the second moment, which is the long-term, average intensity. For short pulse durations, on the order of 10 ms or less, the fourth-order moment is also obtained. A beam of sufficiently high-power density causes heating of the atmosphere and nonlinear propagation, or thermal blooming, of the beam. The threshold transmit powers, beyond which thermal blooming has a discernible effect on the characteristics of the optical beam, are readily estimated for each considered scenario. A discussion of possible implications of exceeding this threshold follows.

For transmit powers below the threshold, thermal blooming effects are disregarded. The beam wave is incident on the object, and power is absorbed, causing a net heating of the material. To quantify the absorbed power the optical absorptivity of the object must be determined. We assume optical absorption in the form of intrinsic bulk absorption. For metallic objects, a Drude model is used to find the complex dielectric constant of the object at the desired optical frequency. For nonmetals, published optical absorptance values are used. The net temperature rise at the object surface is calculated using a thermal diffusion model, appropriate boundary conditions, and the bulk thermal characteristics of the target material. A flow chart of the process is provided in Fig. 1.

Fig. 1

Flow chart of laser heating of solid object through turbulence calculation.

OE_52_3_036001_f001.png

2.

Focused Beam in Atmospheric Turbulence

The expression for an optical focused beam in free space, as shown in Fig. 2, is given by Ref. 10.

(1)

U0(ρ,z)=A0(1+iξz)exp[ikz(kξ2)ρ2(1+iξz)],
where ξ is given by

(2)

ξ=ξ1+iξ2=1z0+i1R0.
Variable z0 is called the “Rayleigh range,” and gives an approximate dividing distance such that the field is considered in the “near field” for z<z0. The Rayleigh range is a function of wavenumber k and aperture radius W0:

(3)

z0=kW022.
R0 is the focal distance. The field at z=0 has Gaussian amplitude and a focused phase front. Figure 1 shows the focused beam in atmospheric turbulence and a target. We define the total input power at z=0, and amplitude constant A0 as follows:

(4)

Pt=(π2W02)A02

(5)

A02=|E|22η0.
The intensity distribution of a focused beam in the turbulence has been obtained, and in the range zizzc, where zi and zc depend on the inner scale l0 and outer scale L0 of the turbulence as follows:

(6)

zi=(0.39Cn2k2l05/3)1

(7)

zc=(0.39Cn2k2L05/3)1.
The intensity distribution is given by Refs. 10 and 11:

(8)

I(ρc,z)=A022W02W2exp(αEz)0tJ0(tρ')exp[t2(12)Dst5/3]dt,
where ρ'=(22)ρc/W; αE is the atmospheric attenuation constant. The structure function Ds and beam size W are given by

(9)

Ds=6.2Cn2k1/3z8/5W5/3.

(10)

W2=W02[(ξ1z)2+(1ξ2z)2].

Fig. 2

Focused beam and object.

OE_52_3_036001_f002.png

The atmospheric attenuation constant αE accounts for the absorption and scattering properties of the medium. Values for αE are widely known and reduction of laser transmittance due to atmospheric attenuation for distances <1km in relatively clear air is minimal.12 In the numerical calculations that follow, αE=0 is assumed unless otherwise specified.

The index of refraction structure constant Cn2 has been found in the range 1017 (very weak) to 1012 (very strong) over land, and 1016 (weak) to 1013 (strong) over the open ocean.11,13 The coherent intensity of the beam wave is given by Ref. 10:

(11)

|U|2=W02W2exp(2ρ2W22a0z),
where the constant of attenuation due to turbulence α0 is

(12)

α0=0.39Cn2k2L05/3.
These expressions can be used to model the intensity of an optical beam through turbulence at the location of a distant target.

3.

Short-Term Intensity

For a short term, on the order of 10 ms or less, the peak intensity is approximately close to the speckle intensity. In this context, speckle refers to the amplitude variations in the laser spot caused by constructive/destructive interference from incident rays of varying phase. When the laser dwell time is sufficiently short, the interference pattern is essentially static, and the use of the average intensity to determine peak intensity is not appropriate. The fourth-order intensity I2 can be approximately given by the following, using circular complex Gaussian assumption:

(13)

I22I2|U|4.
Under this approximation, we get short-term intensity [Eq. (14)], and when the time is not short, we use the average intensity [Eq. (15)].

(14)

Ishort[2I2|U|4]1/2.

(15)

I[I2|U|4]1/2.
Using Eqs. (14) and (15), we can calculate the intensity of the laser beam at the object surface. In the numerical examples, Eq. (14) is used for pulse durations less than 10 ms, and Eq. (15) is used for pulse durations greater than or equal to 10 ms. This arbitrary choice of a short-pulse cutoff results in a small discontinuity in the incident flux at 10 ms.

4.

Thermal Blooming in the Atmosphere

The air in the atmosphere has some absorption at optical wavelengths. The absorption of power produces local heating, thereby increasing the pressure and decreasing the density, which causes a decrease of the refractive index. The beam intensity profile is initially maximum on the axis, but as the heating occurs, the refractive index decreases on the beam axis, which causes a divergence of the beam from the beam axis. This is referred to as thermal blooming of the beam.14 Thermal blooming has been studied numerically and analytically using geometric optics, perturbation theory, and extended Huygens-Fresnel methods. If multiple pulses are used, the wind moves the minimum refractive index away from the beam axis, and the new pulse will be focused into the upstream direction of the wind and the center of the beam moves toward the axis giving enhancement of the peak power on the axis. The peak power can then become higher than the free space peak power. Numerical examples in Ref. 14 show that the peak power can be enhanced by 17%.

If we make a few assumptions about the conditions of operation, we can use an order-of-magnitude expression, following the lead of Ref. 14, to estimate the onset of thermal blooming effects. We first assume the beam wave is propagated through some transverse flow of velocity v, and that the air has some absorption coefficient α, which contributes to local heating of the propagation medium. From Ref. 14 we can estimate the minimum transmit power level at which thermal blooming effects may be significant:

(16)

Pminρ0λγπvW0αEL(γ1)(n01).
If we assume the following constants (typical values for sea level air): γ is the heat capacity ratio; γ=1.4 for air at 20°C; n0 is the initial index of refraction; (n01)=3×104 for air; ρ0 is the ambient pressure; ρ0=105Nm2 at sea level; αE is the atmospheric extinction coefficient; αE=4×104m1 for tropical air in rural environment;10 λ is the wavelength of radiation; λ=105m; v is the velocity of the transverse flow; v=13.6ms1 about 30 mph; W0 is the radius of the aperture; W0=0.05m; and L is the propagation distance to the target; L=0.5×103m, then we find that the threshold power is approximately
Pmin(2.9×107)πvW0L120kW.

The order-of-magnitude calculation shows that for the power densities considered in this analysis, it requires transmit power in excess of the upper limit of most CO2 laser systems to encounter thermal blooming effects. Additionally, this calculation did not consider the effect of atmospheric turbulence, which will, in general, further reduce the radiation power density along the path of the beam wave.

5.

Absorbed Power in the Object

An optical beam is incident on the object and power is dissipated into the medium as heat, which will raise the temperature. The power transmitted into the object La can be calculated as

(17)

La=ωϵ0ϵ2|Et|2=La0exp(α2z2).
A detailed proof comparing Eq. (17) to a widely used expression for total absorbed energy based on reflected power and conservation of energy is given in the Appendix.

The dielectric constant of the object is ε(ω)=ε+iε and the attenuation in the object is given by

(18)

α2=2Im(k2),
where k2=ωμε.

The transmitted electric field into the object at a distance z2 from the incident surface is given by

(19)

Et(z2)=TEiexp(ik2z2),
where T is the transmission coefficient. If we consider a metallic object, the dielectric constant at optical wavelengths is given by a Drude model:

(20)

ε(ω)=1ωp2ω(ω+i1τ),
where ωp is the plasma frequency and τ is the relaxation time. For copper, the plasma frequency and relaxation time were determined in Refs. 1516.17. Using the most recent number, the value of ε is computed as

(21)

εcopper=8.04×103+i1.85×103.
For many nonmetals, the optical constants in the mid-infrared have been directly observed. The relative permittivity of water at 10.6 µm has been determined:18

(22)

εwater=1.21+i0.02.
The complex index of refraction of polyimide at 10.6 µm has been published; from this the relative permittivity is computed as follows:19

(23)

εPI=3.23+i0.18.

Glass-fiber reinforced plastic (GFRP) is a mixture of glass fibers in a thermally cured plastic resin (e.g., polyesters) background material. Properties of GFRP, including volume fraction of glass fiber, melting temperature of resin, and typical thermal properties are described in Ref. 20. The complex index of refraction of silica glass at 10.6 µm has been published; from this, the relative permittivity is computed as follows:21

(24)

εglass=3.75+i2.0.
The complex index of refraction of polyethylene at 10.6μm has been published; from this the relative permittivity is computed as follows:22

(25)

εPE=2.25+i0.0025.
The dielectric properties of composite materials such as GRFP have been studied extensively at microwave frequencies using different mixing formulas.23 At 10.6 µm, the size of glass inclusions in polymer resin background material are usually much larger than the wavelength, and these mixing formulas may not be applicable. Therefore, in this paper, we use a simple average of Eqs. (24) and (25) to obtain a rough estimate of the GFRP permittivity.

6.

Heat Diffusion in Solid Objects

The absorbed power is given by Eq. (18), which will be the heat source for the temperature rise. The diffusion equation for the temperature is given by

(26)

Tt=D2T+LaρCs,
where T is the temperature, ρ(kg/m3) is the density, Cs(J/kg°C) is the specific heat of the solid, and D(m2/s) is the thermal diffusivity. The thermal properties for the four materials under consideration are summarized in Table 1.

Table 1

Thermal properties of various laser targets.

MaterialDensity (kg/m2)Specific heat capacity (J/kg°C)Thermal diffusivity (m2/s)Skin depth thermal time constant (s)Slab thickness thermal time constant (s)
Copper8.92×1030.39×1031.12×10−47.6×10−132.9×10−5
Water1.00×1034.18×1031.44×10−76.0×10−22.3
Polyimide241.42×1031.09×1037.75×10−89.0×10−41.1×102
Glass/resin composite201.90×1031.00×1033.20×10−72.8×10−51.1×102

Equation (26) needs to be solved with boundary conditions. A commonly used approximate boundary condition is the Neumann condition:

(27)

Tn=0.

This is equivalent to zero heat flow across the boundary. This boundary condition is appropriate when the heat transfer from the solid object to the surrounding media is low (i.e., a relatively good thermal conductor in still air). The thermal time constant τδ for a skin depth δs determines the approximate time required for heat generated by the conversion of laser energy in the skin depth region to diffuse a distance greater than the thickness of the skin-depth region.

(28)

δs=1α2

(29)

τδ=δs2D
If the laser is incident on the target for a short time relative to the skin depth of the material the heat diffusion will be approximately given as follows:25

When t is in the range  (tτδ):

(30)

ΔT(0,t)LatρCs.

Assuming heat has diffused out of the skin depth region, for relatively short times the heat will have not yet diffused far enough to reach the back surface of the metallic object located at z2=l. This time is determined by the thermal time constant as follows:

(31)

τl=l2D.
In this time interval heat diffuses as it would in a semi-infinite metallic object, and a Neumann boundary condition is assumed at z2=0 only. With this assumption, the temperature rise, ΔT, at z2=0 is given by the following expression.25

When t is in the range [τδt(τl/3)]:

(32)

ΔT(0,t)(2δsLaρCsl)(τltπ)1/2.
For relatively long times, the heat will have diffused through the thickness of the object, and a Neumann boundary condition assumed at both z2=0 and z2=l. With these assumptions, the temperature rise, ΔT, at z2=0 is given by the following expression.25

When t is in the range [(τl/3)t],

(33)

ΔT(0,t)(δsLaρCsl)(t+τl3).

7.

Numerical Examples

7.1.

Attenuation Due to Turbulence

We consider a beam focused on a solid copper object of thickness l=0.1mm. Figure 3 shows the predicted minimum laser dwell time required to heat the surface of the copper target to the melting point of copper using 10.6-µm radiation through a focused 5-cm aperture at 0.5 km. Although in practice it may be difficult to focus a laser at a distance of 0.5 km, we use the focused case to illustrate the worst-case degradation due to turbulence. The effect of turbulence on a collimated beam is less pronounced, and thus, the result will be closer to the free-space solution. The transmit power is varied from 10 to 100 kW, and various cases of atmospheric turbulence are considered: free space (no turbulence) and turbulence Cn2 equal to 1014, 1013, and 1012. For the thickness of copper and time range considered, Eq. (33) is used exclusively as heat has diffused a distance much greater than the thickness of the copper target.

Fig. 3

ΔT versus transmit power (Pt) at λ=10.6μm; L=0.5km; A0=5cm (focused); Pt=10 to 100 kW; copper thickness l=0.1mm.

OE_52_3_036001_f003.png

The plot in Fig. 3 can be used to predict the laser dwell time required to realize the onset of melting for a laser system of known transmit power. This result considers only when the absorbed power is sufficient to cause a sustained rise in temperature in excess of the melt temperature it applies generally to either pulsed or continuous wave (CW) operation, provided that the average power is held constant (e.g., 10 to 100 kW). In pulsed operation, the beam will have a higher power density, and thus, may cause a momentary spike in local temperature above the melt temperature, but during the off portion of the duty cycle the heat will diffuse out of the area and the temperature will decrease. Over many pulses the average temperature will rise at the same rate as a CW laser of equal average power.

Figure 4 shows the predicted minimum laser dwell time required to heat the surface of a l=10mm sample of Kapton (polyimide) target the glass transition temperature (400°C) using 10.6-µm radiation through a focused 5-cm aperture at 0.5 km. The transmit power is varied from 10 to 100 kW, and various cases of atmospheric turbulence are considered: free space (no turbulence) and turbulence Cn2 equal to 1014, 1013, and 1012. The discontinuity at dt=102s is due to the inclusion of the fourth-order moment for short pulse durations as described in Sec. 3. For the thickness of polyimide and time range considered, Eq. (30) is used for t<9×104s, and Eq. (32) is used for t<9×104s.

Fig. 4

ΔT versus transmit power (Pt) at λ=10.6μm; L=0.5km; A0=5cm (focused); Pt=10 to 100 kW; Kapton thickness l=10mm.

OE_52_3_036001_f004.png

Figure 5 shows the predicted minimum laser dwell time required to heat the surface of a l=1mm sample of water to its boiling point (100°C) using 10.6-µm radiation through a focused 5-cm aperture at 0.5 km. The transmit power is varied from 10 to 100 kW, and various cases of atmospheric turbulence are considered: free space (no turbulence) and turbulence Cn2 equal to 1014, 1013, and 1012. For the thickness of water and power ranges range considered, Eq. (30) is used exclusively as τδ=0.06s, which is greater than the maximum required laser dwell time to heat water to its boiling point.

Fig. 5

ΔT versus transmit power (Pt) at λ=10.6μm; L=0.5km; A0=5cm (focused); Pt=10 to 100 kW; water thickness l=1mm.

OE_52_3_036001_f005.png

Figure 6 shows the predicted minimum laser dwell time required to heat the surface of a l=10mm sample of GFRP to the melting point of polyester resin using 10.6-µm radiation through a focused 5-cm aperture at 0.5 km. The transmit power is varied from 10 to 100 kW, and various cases of atmospheric turbulence are considered: free space (no turbulence) and turbulence Cn2 equal to 1014, 1013, and 1012. For the thickness of GFRP and power ranges range considered, Eq. (32) is used exclusively as the required laser dwell times fall between τδ=2.8×105s and (τl/3)=1.1×102s.

Fig. 6

ΔT versus transmit power (Pt) at λ=10.6μm; L=0.5km; A0=5cm (focused); Pt=10 to 100 kW; glass-fiber reinforced plastic (GFRP) thickness l=10mm.

OE_52_3_036001_f006.png

7.2.

Attenuation Due to Turbulence and Atmospheric Extinction

The combined effect of atmospheric scattering and attenuation on a laser beam propagated through the atmosphere to a solid object target is a reduction in transmittance, which can, under the right conditions, become similar to or even greater in magnitude than the interference caused by atmospheric turbulence. Like the attenuation caused by atmospheric turbulence, the total extinction seen in a volume of atmosphere increases with the greater separation between the laser source and the target. This is easily deduced from Eq. (8), which contains the exponential attenuation term, exp(αEz). The product of extinction coefficient and linear distance is commonly called optical thickness: τE=αEz. An order-of-magnitude approximation is that the effect of atmospheric extinction should be considered when optical thickness τE1. For a distance of 0.5 km this corresponds to αE2km1. Figure 7 shows the effect of atmospheric extinction in this range on the laser dwell time required to heat a copper target to its melting point. Considering the free space (no turbulence) case, we can observe an increase in minimum laser dwell time when αE=2km1, which is in the range of the increase required in moderate turbulence. If optical thickness τE1, then the increase in laser dwell time will be minimal, and if τE1, then the effect of atmospheric extinction will dominate, and the relative effect of turbulence will be small.

Fig. 7

ΔT versus transmit power (Pt) at λ=10.6μm; L=0.5km; A0=5cm (focused); Pt=10 to 100 kW; copper thickness l=0.1mm.

OE_52_3_036001_f007.png

Figure 7 shows the effect of αE=2km1 on laser propagation over 500 m. In real-world conditions, this value of atmospheric extinction is not typical when visibility is relatively good. According to Ref. 13, for visibility ranging from hazy to clear (5 to 40 km) atmospheric extinction remains below αE=0.33km1 in a rural nonmaritime environment. According to Ref. 26, the atmospheric extinction coefficient for a maritime environment at 99% relative humidity is αE=0.45km1. In the presence of fog, however, when visibility can be reduced to 50 to 1000 m atmospheric extinction in the range αE=1100km1 can be expected. This is shown in a number of fog model calculations in Refs. 26 and 27.

8.

Conclusions

This paper presents a study on an optical beam propagated through turbulence and incident on various targets, generating heat and causing a temperature rise. The propagation of a beam wave in the atmosphere is discussed, including the fourth-order moment and the effect of thermal blooming. Transmitted power of 10 to 100 kW is absorbed in the target and converted into heat causing the temperature to rise. Numerical examples are given for 10.6 µm, free-space, and turbulence, aperture size W0=5cm, copper target thickness of 0.1 mm, Kapton target thickness of 0.1 mm, water thickness of 1 mm, and GFRP thickness of 10 mm. In each case, the beam is focused at 0.5 km, and the minimum laser dwell time required to heat the sample to its melting/boiling point is predicted. The effect of atmospheric extinction through absorption and scattering is discussed, and it is determined that the attenuation due to atmospheric extinction should be included in a dwell time calculation when τE1. A plot is presented which shows the predicted difference in required dwell time for a focused laser through the atmosphere in zero and moderate turbulence with τE=0 and τE=1.

Appendices

Appendix:

Proof for Absorbed Power Equation

The absorbed power-per-unit volume in the solid object target is given in Eq. (16). Likewise, the total absorbed power from z=0 to z=d is given by

(34)

0dLadz.
If the integration path is sufficiently large, the total power is equal to the incident power minus the reflected power. This approach has been used in several papers including Refs. 16 and 17. In the case of thin objects and/or short integration paths, Eq. (34) gives the correct absorbed power. In this paper, we consider solid copper objects, and because the skin depth is small relative to object thickness, we can use Eq. (34) with . To prove equivalence, we can show that the incident power Pinc minus the reflected power Pref is equal to Eq. (37), with d.

We start with the exact formulation:

(35)

Ei(1+R)=Et

(36)

EiZ0(1R)=EtZt,
where Ei and Et are the incident electric field and transmitted electric field, respectively, and where Z0= μ0/ε0 and Zt= μ0/εt are the characteristic impedances of free-space and object regions, respectively. From Eqs. (35) and (36), we can get the transmitted power:

(37)

|Ei|2Z0(1+R)(1R*)=|Et|2Zt*.
Taking the real part of Eq. (37), noting that (1+R)(1R*)=1+RR*RR*, and that RR* is purely imaginary, we get

(38)

Re{|Ei|22Z0(1|R|2)}=Re{|Et|22Zt*}.
We note that Zt=Z0/ε; ε=(n+in)2; ε=n2n2; ε=2nn; and we get

(39)

|Ei|22Z0(1|R|2)=Re{|Et|22Z0n*}=|Et|22Z0n.
Using the expression for La in Eq. (18), and evaluating Eq. (34) at, d we get

(40)

0ωϵ0ϵ2|Et|2exp(α2z)dz=ωϵ0ϵ2|Et|21α2=|Et|22Z0n,
where α2=k02n and ε=2nn.

Whereas the expression for transmitted power in Eq. (39) is derived from conservation of power and object interface boundary conditions, the expression for transmitted power in Eq. (40) is the total sum of the position dependent transmitted power given by Eq. (18). We can see that for objects where no energy transits the thickness of the object (i.e., object is thick relative to optical skin depth), these are equivalent.

Acknowledgments

This work was supported by the Office of Naval Research (N000141110358) and the National Science Foundation (ECCS0925034).

References

1. 

S. Ackerman, “Video: Navy Laser Sets Ship on Fire,” Wired, 8 April 2011, http://www.wired.com/dangerroom/2011/04/video-navy-laser-sets-ship-on-fire/ (5 March 2011).Google Scholar

2. 

L. AndrewsR. Phillips, Laser Beam Propagation Through Random Media, pp. 321–395, SPIE, Bellingham (2005).Google Scholar

3. 

A. Ishimaru, Wave Propagation and Scattering in Random Media, pp. 399–452, Academic Press, New York (1978).Google Scholar

4. 

P. SprangleJ. R. PenanoB. Hafizi, “Optimum wavelength and power for efficient laser propagation in various atmospheric environments,” Report NRL/MR/6790-05-8907, Naval Research Laboratory, Washington, DC (2005).Google Scholar

5. 

P. SprangleJ. R. PenanoB. Hafizi, “Propagation of intense short laser pulses in the atmosphere,” Phys. Rev. E 66(4), 046418 (2002).PLEEE81063-651Xhttp://dx.doi.org/10.1103/PhysRevE.66.046418Google Scholar

6. 

J. R. Penanoet al., “Propagation of ultra-short, intense laser pulses in air,” Phys. Plasmas 11(5), 2865 (2004).PHPAEN1070-664Xhttp://dx.doi.org/10.1063/1.1648020Google Scholar

7. 

Z. H. Shenet al., “Mathematical modeling of laser induced heating and melting in solids,” Opt. Las. Technol. 33(8), 533–537 (2001).OLTCAS0030-3992http://dx.doi.org/10.1016/S0030-3992(01)00005-6Google Scholar

8. 

J. H. Betchel, “Heating of solid targets with laser pulses,” J. Appl. Phys. 46(4), 1585–1593 (1975).JAPIAU0021-8979http://dx.doi.org/10.1063/1.321760Google Scholar

9. 

A. M. Prokhorovet al., Laser Heating of Metals, pp. 39–72, Adam Hilger/IOP, New York (1990).Google Scholar

10. 

A. Ishimaru, “The beam wave case and remote sensing,” in Laser Beam Propagation in the Atmosphere, J. W. Strohbehn, Ed., pp. 129–170, Springer-Verlag, New York (1978).Google Scholar

11. 

S.F. Clifford, “The classical theory of wave propagation in a turbulent medium,” in Laser Beam Propagation in the Atmosphere, J. W. Strohbehn, Ed., pp. 9–43, Springer-Verlag, New York (1978).Google Scholar

12. 

B. J. Lund, “Laser atmospheric attenuation tables for LTAS,” Report AL/OE-TR-1997-0063, USAF, Armstrong Laboratory (1997).Google Scholar

13. 

P. A. Fredericksonet al., “Estimating the refractive index structure parameter over the ocean using bulk methods,” J. Appl. Meteor. 39(10), 1770–1783 (2000).JAMOAX0894-8763http://dx.doi.org/10.1175/1520-0450-39.10.1770Google Scholar

14. 

J. L. WalshP. B. Ulrich, “Thermal blooming in the atmosphere,” in Laser Beam Propagation in the Atmosphere, J. W. Strohbehn, Ed., pp. 223–320, Springer-Verlag, New York (1978).Google Scholar

15. 

M. A. Ordalet al., “Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W,” Appl. Opt. 24(24), 4493–4499 (1985).APOPAI0003-6935http://dx.doi.org/10.1364/AO.24.004493Google Scholar

16. 

X. E Lin, “Laser pulse heating,” in Proc. Particle Accelerator Conf., New York, Vol. 2, pp. 1429–1431 (1999).Google Scholar

17. 

S. B. BoydenY. Zhang, “Temperature and wavelength-dependent spectral absorptivities of metallic materials in the infrared,” J. Thermophys. Heat Trans. 20(1), 9–15 (2006).JTHTEO0887-8722http://dx.doi.org/10.2514/1.15518Google Scholar

18. 

S. G. WarrenR. E. Brandt, “Optical constants of ice from the ultraviolet to the microwave: a revised compilation,” J. Geophys. Res. 113, D14220 (2008).JGREA20148-0227http://dx.doi.org/10.1029/2007JD009744Google Scholar

19. 

Z. M. Zhanget al., “Infrared refractive index and extinction coefficient of polyimide films,” Int. J. Thermophys. 19(3), 905–916 (1998).IJTHDY0195-928Xhttp://dx.doi.org/10.1023/A:1022655309574Google Scholar

20. 

R. Komanduri, “Machining of fiber-reinforced composites,” Mach. Sci. Technol. 1(1), 113–152 (1997).MSTEFW1091-0344http://dx.doi.org/10.1080/10940349708945641Google Scholar

21. 

R. KitamuraL. PilonM. Jonasz, “Optical constants of silica glass from extreme ultraviolet to far infrared at near room temperature,” Appl. Opt. 46(33), 8118–8133 (2007).APOPAI0003-6935http://dx.doi.org/10.1364/AO.46.008118Google Scholar

22. 

E. D. PalikG. Ghosh, Handbook of Optical Constants of Solids, pp. 957–987, Academic, New York (1998)Google Scholar

23. 

L. TsangJ. A. KongK. H. Ding, Scattering of Electromagnetic Waves, pp. 351–353 Wiley, New York (2001).Google Scholar

24. 

“DuPont Kapton Polyimide Film Datasheet,” Published 2011, DuPont, http://www.dupont.com/kapton/general/H-38479-4.pdf (5 March 2013).Google Scholar

25. 

M. Sparks, “Theory of laser heating of solids: metals,” J. Appl. Phys. 47(3), 837–849 (1976).JAPIAU0021-8979http://dx.doi.org/10.1063/1.322717Google Scholar

26. 

E. P. ShettleR. W. Fenn, “Models for the aerosols of the lower atmosphere and the effects of humidity variations on their optical properties,” Report AFGL-TR-79-0214, USAF, Air Force Geophysics Laboratory, Hanscom AFB, Massachusetts (1979).Google Scholar

27. 

C. Reinhardtet al., “Atmospheric channel transfer function estimation from experimental free-space optical communications data,” Opt. Eng. 51(3), 031205 (2012).OPEGAR0091-3286http://dx.doi.org/10.1117/1.OE.51.3.031205Google Scholar

Biography

Matthew Stoneback received his BS and MS in electrical engineering from the University of Washington in 2008 and 2009, respectively. From 2007 to 2011, he was an intern and engineer with Boeing Research and Technology in Seattle, WA. He is currently pursuing a PhD in electrical engineering at the University of Washington. His interests are in the areas of antennas, phased arrays, electromagnetic propagation, and optics.

Akira Ishimaru received his BS in 1951 from the University of Tokyo, and PhD in electrical engineering in 1958 from the University of Washington. From 1951 to 1952, he was with the Electrotechnical Laboratory, Tokyo, and in 1956 he was with Bell Laboratories, New Jersey. In 1958, he joined the faculty of the Department of Electrical Engineering of the University of Washington, where he was professor of electrical engineering and adjunct professor of applied mathematics, and is currently professor emeritus. He was also a visiting associate professor at the University of California, Berkeley. His current research includes waves in random media, remote sensing, object detection and imaging in clutter environment, inverse problems, millimeter wave and optical propagation and scattering in the atmosphere and the terrain, acoustic scattering in the ocean, ultrasound imaging, and optical diffusion in tissues. He is the author of Wave Propagation and Scattering in Random Media (Academic Press, 1978; IEEE Press-Oxford University Press Classic Reissue, 1997) and Electromagnetic Wave Propagation, Radiation, and Scattering (Prentice Hall, 1991).

Colin Reinhardt received his BSEE and BS in physics from the University of Washington in 2005, and PhD from University of Washington in 2010, as an American Society for Engineering Education and Department of Defence (ASEE-DoD) Science, Mathematics, and Research for Transformation (SMART) Fellow. He now works at SPAWAR Systems Center Pacific, in San Diego, California, with the Atmospheric Propagation Branch. His interests include electromagnetic wave propagation in random media, free-space optical communications, multiple-scattering theory, radiative transfer theory, optical turbulence, and communications theory.

Yasuo Kuga was received his BS, MS, and PhD degrees from the University of Washington, Seattle, in 1977, 1979, and 1983, respectively. He was an engineer at Matsushita Electric Industrial Corporation, Osaka, Japan, from 1971 to 1975. From 1979 to 1983, he was a research associate at the University of Washington. From 1983 to 1988, he was a research assistant professor of electrical engineering at the University of Washington. From 1988 to 1991, he was an assistant professor of electrical engineering and computer science at the University of Michigan. He is now professor of electrical engineering at the University of Washington. His research interests are in the areas of microwave and millimeter-wave remote sensing, high-frequency devices, and optics.

Matthew Stoneback, Akira Ishimaru, Colin Reinhardt, Yasuo Kuga, "Temperature rise in objects due to optical focused beam through atmospheric turbulence near ground and ocean surface," Optical Engineering 52(3), 036001 (8 March 2013). http://dx.doi.org/10.1117/1.OE.52.3.036001
Submission: Received ; Accepted
JOURNAL ARTICLE
9 PAGES


SHARE
KEYWORDS
Atmospheric propagation

Atmospheric optics

Turbulence

Atmospheric turbulence

Copper

Ocean optics

Thermal blooming

RELATED CONTENT

Variable-Atmosphere-Effects-On-High-Energy-Laser-Propagation
Proceedings of SPIE (December 19 1979)
Turbulent-thermal-blooming
Proceedings of SPIE (April 24 2008)
Vortex-laser-beam-in-a-randomly-inhomogeneous-medium
Proceedings of SPIE (November 25 2014)
Intermittency-of-laser-scintillation-over-sea
Proceedings of SPIE (October 13 2005)

Back to Top