Turbulence is a non-linear process, and so the equations governing it are non-linear. In 1941, Kolmogorov attempted to model turbulence in the atmosphere using a statistical approach. The Kolomogorov theory of turbulence is based on the assumption that wind velocity fluctuations are approximately locally homogeneous and isotropic random fields for scales less than the largest wind eddies. (Beland, 1993; Beckers 1993).
How do phase variations due to atmospheric turbulence affect the intensity
of the incident electro-magnetic wave? Assuming the incoming electric field
portion of the electro-magnetic wave E is proportional to
with
the angular frequency and t time, the wave equation
reduces to
where k is the wavenumber for propagation in free space, n is the index
of refraction and the dispersion relation
holds. (Note that the wavenumber for propagation through the atmosphere
is
.)
This equation cannot be solved exactly, so we approximate the solution
using a perturbation expansion of n and E. This assumes that for the area
of the sky we are interested in, the index of refraction undergoes slow,
small perturbations due to the atmospheric turbulence. Let
where
is the incident electric field and
is the scattered field.
where
.
Dropping the second order
term,
Putting in the perturbation to E
Dropping the small
term and using
equation 22 becomes
This can be solved in terms of
using Green's function for the wave equation, giving
integrating over the scattering volume between source and telescope.
(Beland, 1993).
However, this solution is only valid for very weak perturbations of
E and the integral over the scattering volume between the source and telescope
is not very convenient since the effects of phase variations are not obvious.
A less limiting method is the Rytov method of smooth perturbations. Here
E is written as
where
and
,
A is the amplitude and S is the phase of E. Lutomirski and Yura, using
an extended Huygens-Fresnel theory find
here integrating over a surface.
The intensity of the incoming wavefront will tell us how the phase fluctuations
degrade the resolution.
The mutual coherence function (MCF) in the focal plane of the telescope
is given by
and is equal to
where
is the transverse coherence length. A spatially coherent wavefront will
have its spatial coherence degraded to
after passing through turbulence.
Most astronomical applications translate this spatial mutual coherence
function to the modulation transfer function (MTF) in the spatial frequency
domain. The modulation transfer function
where f is the focal length,
is the spatial frequency (cycles/length),
is the wavelength of light, and
is the atmospheric coherence length or Fried parameter.
is about 10 cm for a vertical path for visible wavelength and characterizes
the effect of turbulence on the incoming wavefront. Since optical turbulence
blurs out point source functions, there is a cutoff spatial frequency
and
gives the limit of resolution due to optical turbulence. (Beland, 1993;
Beckers, 1993.)
The diffraction limit of resolution for 4m telescope at 500 nm is given
by equation 17
Assuming ,
the angular resolution due to atmospheric turbulence R is
So turbulence in this case degrades the resolution by a factor of 100!