Adaptively Refined Large-Eddy Simulations of Galaxy Clusters

Appendix D Fourier transform and structure functions

The continuous one-dimensional Fourier transform in k-space F(k) of some function in x-space f(x)is defined like

F(k) =12π-f(x)e-ikx𝑑x(Fourier transform). (D.1)

Using the Fourier transform on a function twice will produce the original function again, but mirrored at the origin. That’s why one conventionally defines an inverse Fourier transform5656Nature does not know about the inverse Fourier transform. If you have some optical device, which produces the Fourier transform of some image and you use it twice on your image you will get a mirrored image!, that will generate the not mirrored original function again, when used on the Fourier transform of a function

f(x) =12π-F(k)eikx𝑑x(inverse Fourier transform). (D.2)

In three dimensions one defines the Fourier transform like

F(𝒌)=1(2π)3/2-f(𝒙)e-i𝒌𝒙𝑑V, (D.3)
f(𝒙)=1(2π)3/2-F(𝒌)ei𝒌𝒙𝑑K. (D.4)

In cartesian coordinates the kernel of the Fourier transform e-i𝒌𝒙=e-i(kxx+kyy+kzz) separates and so the three-dimensional Fourier transform of a function which separates in cartesian coordinates f(𝒙)=a(x)b(y)c(z) is also separable

F(𝒌)=A(kx)B(ky)C(kz)=1(2π)3/2-a(x)e-ikxx𝑑x-b(y)e-ikyy𝑑y-c(z)e-ikzz𝑑z.

Thats why we like to use cartesian coordinates when we are using Fourier transforms.

D.1 Fourier transform of a delta function

An important result can be derived by computing the inverse Fourier transform of the Fourier transform of a delta function

δ(x-x0) =12π-12π-δ(x-x0)e-ikx𝑑xeikx𝑑k
=12π-e-ikx0eikx𝑑k=12π-eik(x-x0)𝑑k.

From this we get that the inverse Fourier transform of a constant is the delta function

12π-eik(x-x0)𝑑k=2πδ(x-x0).

Taking the complex conjugate of this equation and making use of the fact that δ*(x-x0)=δ(x-x0) we get as a definition for the delta function

δ(x-x0)=12π-e±ik(x-x0)𝑑k. (D.5)

Using this we can derive the astonishing result

-f(x)𝑑x=2πF(0) (D.6)

as can be seen from

-f(x)𝑑x =-12π-F(k)eikx𝑑k𝑑x=12π-F(k)-eikx2πδ(k)𝑑x𝑑k
=2π-F(k)δ(k)𝑑k=2πF(0).

D.2 Convolution theorem

The Fourier transform of the product of two function in k-space is

12π-F(k)G(k)eikx𝑑k=
=12π-12π-f(x)e-ikx𝑑x12π-g(x′′)e-ikx′′𝑑x′′eikx𝑑k
=1(2π)3/2-f(x)e-ikxg(x′′)e-ikx′′eikx𝑑x𝑑x′′𝑑k
=1(2π)3/2-f(x)g(x′′)-e-ik(x+x′′-x)𝑑k2πδ(x′′-(x-x))𝑑x𝑑x′′
=12π-f(x)-g(x′′)δ(x′′-(x-x))𝑑x′′𝑑x
=12π-f(x)g(x-x)𝑑x=h(x).

The integral h(x) is called convolution of the functions f(x) and g(x). So the convolution theorem says that

h(x)=12π-f(x)g(x-x)𝑑x=12π-F(k)G(k)eikx𝑑k. (D.7)

D.3 Autocorrelation and Wiener-Khinchin Theorem

The autocorrelation of a function is defined as5757Note that with our definition of the Fourier transform we cannot define the autocorrelation function as hAC(x)=12π-f(x)f*(x+x)𝑑x, because we could then not derive the Wiener-Khinchin theorem.

hAC(x)=12π-f*(x)f(x+x)𝑑x. (D.8)

The Wiener-Khinchin Theorem states that

12π-f*(x)f(x+x)𝑑x=12π-|F(k)|2eikx𝑑k, (D.9)

which can be proved in analogy to the convolution theorem

12π-f*(x)f(x+x)𝑑x=
=12π-f*(x)-f(x′′)δ(x′′-(x+x))𝑑x′′𝑑x
=12π-f*(x)-f(x′′)12π-e-ik(x+x′′-x)𝑑k𝑑x′′𝑑x
=12π-12π-f*(x)eikx𝑑x12π-f(x′′)e-ikx′′𝑑x′′eikx𝑑k
=12π-F*(k)F(k)eikx𝑑k=12π-|F(k)|2eikx𝑑k.

A special case of the Wiener-Khinchin theorem is Parseval’s theorem

-|f(x)|2𝑑x=-|F(k)|2𝑑k, (D.10)

which can be obtained from the Wiener-Khinchin theorem for x=0

hAC(0) =12π-f*(x)f(x)𝑑x=12π-|f(x)|2𝑑x
=12π-|F(k)|2eik0𝑑k=12π-|F(k)|2𝑑k.

D.4 Structure functions

A structure function of order p is defined as5858See Pope (2000).

Sp(f(x))=[f(x+x)-f(x)]p=12π-[f(x+x)-f(x)]p𝑑x. (D.11)

The second order structure function is related to the spectrum |F(k)| of the function f like

12π-[f(x+x)-f(x)]2𝑑x=22π-(1-eikx)|F(k)|2𝑑k, (D.12)

which can be proved5959A sketch of this prove can also be found in Pope (2000, Appendix G). by expanding the second order structure function

S2(f(x))=12π-[f(x+x)-f(x)]2𝑑x
=12π[-|f(x+x)|2𝑑x-2-f*(x)f(x+x)𝑑x+-|f(x)|2𝑑x].

Substituting x′′=x+x in the first term we get

S2(f(x)) =12π[-|f(x′′)|2𝑑x′′-2-f*(x)f(x+x)𝑑x+-|f(x)|2𝑑x]
=22π[-|f(x)|2𝑑x--f*(x)f(x+x)𝑑x].

Using Parseval’s and the Wiener-Khinchin theorem we obtain the final result

S2(f(x))=22π[-|F(k)|2𝑑k--|F(k)|2eikx𝑑k]
=22π-(1-eikx)|F(k)|2𝑑k.

The structure functions used in the theory of Kolmogorov are the so called longitudinal structure functions of the velocity, which are defined as

S2(v(l))=([𝒗(𝒙+𝒍)-𝒗(𝒙)]𝒍l)p=(v(𝒙+𝒍)-v(𝒙))p. (D.13)

They are related to the longitudinal velocity spectrum6060In the literature this is often called kinetic energy spectrum, but this is only true for incompressible flows. |V(k)|2 via equation (D.12). Sometimes also second order transverse structure functions are measured. These are defined as

S2(v(l))=(|[𝒗(𝒙+𝒍)-𝒗(𝒙)]×𝒍|l)p. (D.14)

The behavior of the second order transverse structure functions for homogeneous turbulence is uniquely determined by the longitudinal structure function (Pope, 2000, p. 192, Eqs. (6.28)). They also show the characteristic 2/3-slope as predicted for the longitudinal structure functions (Frisch, 1995, p.60).

In general structure functions of vectorial quantities like the velocity are tensors, e.g. the general second order structure function of the velocity can be defined as

Sij(𝒙,𝒍)=[vi(𝒙+𝒍)-vi(𝒙)][vj(𝒙+𝒍)-vj(𝒙)]. (D.15)

But it can be shown that for local isotropy only the longitudinal structure function S2(v(l))=S11 and the transversal structure S2(v(l))=S22=S33 are unequal zero (Pope, 2000). Since the transverse structure function is determined by the longitudinal structure function in case of local homogeneity, for homogeneous and isotropic turbulence Sij is determined by the single scalar function S11=S2(v(l)) (Pope, 2000).

The third order structure function used in Kolmogorov theory is defined as

S111(𝒙,𝒍)=[v1(𝒙+𝒍)-v1(𝒙)]3, (D.16)

which is often simply called S3(v(l)). So the famous four-fifths law of Kolmogorov is actually true only for one component of the third order structure function tensor, but again for homogeneous and isotropic turbulence the third order structure function tensor Sijk is uniquely determined by the single scalar function S111=S3(v(l)).