# Thermal Blooming

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

A high powered laser going through the air will cause the air to heat up due to absorption. If the laser beam lasts longer than a small fraction of a second, the hot air will expand and become less dense. Less dense air has a lower index of refraction, so the density gradient across the beam acts like a lens which in this case makes the beam expand. This expansion is called thermal blooming. The phenomenon is non-linear - if you try to compensate for not enough intensity on target by increasing the beam power without also making the beam wider to compensate, you can just end up heating the air more and actually reducing the beam intensity.

Bending of the beam from an unperturbed circular shape (dashed circle) into a crescent shape due to thermal blooming in wind[1]

If there is wind, or if you are moving quickly through the air (like if you are flying an airplane), or your target is (like if you are shooting at missiles or artillery shells), you are continually getting fresh air into your beam. This makes the air that has been in your beam for longer hotter than the newer air, forming a lens which bends your beam in addition to expanding it.

Adaptive optics methods can help to solve this, but would-be laser engineers need to be careful because a straightforward implementation just ends up making the problem worse (to counteract the expansion, you focus the beam more which heats up the air more which makes your beam diverge even more). Alternatively, you can use high power laser pulses that are so short that the air does not have time to expand.

Does your beam need to worry about thermal blooming? For relatively long duration beams moving relative to the air you can compute the thermal distortion factor, ${\displaystyle N}$, to find out. If ${\displaystyle N}$ is less than one, you can ignore thermal blooming. If ${\displaystyle N}$ is greater than this, you will need to correct for blooming if you want to keep a tight spot size. The larger ${\displaystyle N}$, the more difficult the blooming will be to correct.

To compute ${\displaystyle N}$ [1] (Warning! Math ahead!)

${\displaystyle {\begin{array}{lll}N_{c}&=&{\dfrac {16{\sqrt {2}}}{\pi }}{\dfrac {\beta \,P\,R^{2}}{R_{a}\,v\,D^{3}}}\\N_{E}&=&R/R_{0}\\f&=&\left\{{\begin{array}{ll}1&{\mbox{if }}N_{E}=0\\{\dfrac {2}{N_{E}^{2}}}\,\left(N_{E}-1+\exp[-N_{E}]\right)&{\mbox{otherwise}}\end{array}}\right.\\N_{F}&=&{\dfrac {\pi }{4}}\,n\,{\dfrac {D}{S}}\\q&=&\left\{{\begin{array}{ll}1&{\mbox{if }}N_{F}=1\\{\dfrac {2\,N_{F}^{2}}{N_{F}-1}}\,\left(1-{\dfrac {\ln[N_{F}]}{N_{F}-1}}\right)&{\mbox{otherwise}}\end{array}}\right.\\N_{V}&=&v_{T}/v\\s&=&\left\{{\begin{array}{ll}1&{\mbox{if }}N_{V}=1\\{\dfrac {2}{(N_{V}-1)^{2}}}\,\left(N_{V}\ln[N_{V}]-(N_{V}-1)\right)&{\mbox{otherwise}}\end{array}}\right.\\N&=&N_{c}\,f\,q\,s\end{array}}}$

where ${\displaystyle S}$ is the beam spot size diameter at the target, ${\displaystyle D}$ is the beam diameter at the aperture, and ${\displaystyle R}$ is the range to the target. ${\displaystyle P}$ is the beam power. The relative wind speed to the beam is ${\displaystyle v}$ at the aperture and ${\displaystyle v_{T}}$ at the target, and ${\displaystyle n}$ is the refractive index of the air. ${\displaystyle R_{a}}$ is the characteristic absorption length of the wavelength of light you are using for the given atmospheric conditions, and ${\displaystyle R_{0}}$ is the net attenuation length including both absorption and scattering. Remember that ${\displaystyle R_{a}}$ is the absorption length, not the total attenuation length, so primarily scattering phenomena like clouds or fog or mist won’t change it much.

The parameter ${\displaystyle \beta }$ is the atmospheric thermal coefficient

${\displaystyle \beta ={\frac {(-dn/dT)}{\rho \,C_{P}}}}$

where ${\displaystyle dn/dT}$ is the rate of change of the refractive index with increasing temperature, ${\displaystyle \rho }$ is the air’s mass density, ${\displaystyle C_{P}}$ is the specific heat capacity of the air at constant pressure, and ${\displaystyle \beta }$ in sea level air on earth under typical conditions at visible and near visible wavelengths is around ${\displaystyle \beta =8.3\times 10^{-10}{\mbox{m}}^{3}/{\mbox{J}}}$.

If there is no wind speed, you can take the speed at which the heated air rises out of the beam path (a process known as convection) instead.

${\displaystyle v_{\mbox{convection}}=\left({\frac {2\,P\,g}{\rho \,C_{P}\,R_{a}\,T}}\right)^{1/3}}$

where ${\displaystyle g=9.8\,{\mbox{m}}/{\mbox{s}}^{2}}$ is the acceleration due to gravity. Normally, however, if you are outside this will be small compared to any natural wind in the environment. And if you are inside the range is probably too short to matter.

For pulsed beams, thermal blooming can be less of an issue. For one thing, the pulse is usually so fast that wind has a negligible effect. If your pulse length is shorter than your minimum beam diameter divided by the speed of sound, thermal blooming is much less of an issue because the air will not have had enough time to expand to make a lens. The equivalent of the thermal distortion factor ${\displaystyle N}$ for pulsed beams are the parameters ${\displaystyle N_{l}}$ and ${\displaystyle N_{s}}$ for pulses long and short compared to the time for sound to cross their beam diameters, respectively. [2]

${\displaystyle {\begin{array}{lll}N_{l}&=&{\dfrac {1}{6\pi }}\,{\dfrac {\beta \,E\,R^{2}}{R_{a}\,D^{2}\,S^{2}}}\\N_{s}&=&{\dfrac {1}{30\pi }}\,{\dfrac {\beta \,E\,R^{2}\,(t\,c_{s})^{2}}{R_{a}\,D^{2}\,S^{4}}}\end{array}}}$

where ${\displaystyle E}$ is the pulse energy, ${\displaystyle t}$ is the pulse duration, ${\displaystyle c_{s}}$ is the speed of sound (approximately 330 m/s in sea level air at shirt-sleeve temperatures). If ${\displaystyle D/c_{s}, use ${\displaystyle N_{l}}$. If ${\displaystyle S/c_{s}>t}$, use ${\displaystyle N_{s}}$. If you just happen to be somewhere between these, either take your best guess or get a grant, hire several grad students and post-docs, and perform the research to figure out the proper formula yourself. Again, if ${\displaystyle N_{l}}$ or ${\displaystyle N_{s}}$, as appropriate, is less than one thermal blooming is not much of an issue. If they are, either use sophisticated non-linear adaptive optics beam controls or adjust your beam parameters so it is less than one.

It is interesting to note that for diffraction-limited beams, ${\displaystyle {\frac {R}{S\,D}}}$ does not depend on distance so long as the beam is focused at the target. So a pulsed beam governed by ${\displaystyle N_{l}}$ has the same amount of thermal blooming at any distance.

## Credit

Author: Luke Campbell

## References

1. Frederick Gebhardt, “Twenty-five years of thermal blooming: an overview”, Proc. SPIE 1221, Propagation of High-Energy Laser Beams Through the Earth’s Atmosphere, (1 May 1990); doi: 10.1117/12.18326
2. H. Kleiman and R. W. O'Neil "Thermal blooming of pulsed laser radiation", Appl. Phys. Lett. 23, 43 (1973). ${\displaystyle N_{l}}$ differs by a factor of 2 from ${\displaystyle T_{lc}}$ and ${\displaystyle N_{s}}$ by a factor of 4 from ${\displaystyle T_{s}c}$ because those give the time-dependent change in the laser spot and we are interested in the total fluence dumped on the target during the time the spot is on. So the presented results are for time-averaging the ${\displaystyle T}$ parameters.