Adaptively Refined Large-Eddy Simulations of Galaxy Clusters

Appendix A Dimensional analysis

We want to write down the general equations of fluid dynamics in dimensionless form. Therefore we introduce the following dimensionless quantities

ri*=ril0 ri=ri*ri*ri=1l0ri*,
vi*=viv0,
t*=tl0 t=t*t*t=1t0t*,
ρ*=ρρ0.

Inserting these into the continuity equation (2.1) we get

ρ0t0t*ρ*+v0ρ0l0rj*(vj*ρ*) =0|l0ρ0v0,
l0v0t0Srt*ρ*+rj*(vj*ρ*)=0.

This derivation shows that solutions of the continuity equation are similar, if the Strouhal number Sr=l0v0t0 is the same. Flows with a Strouhal number Sr=0 are so called stationary flows. Nevertheless, the Strouhal number is most often set to one, by assuming vo=l0t0. Using the additional dimensionless quantities p*=pp0,σij*=σijσ0,g*=gg0 in the momentum equation (2.2) yields

ρ0v0t0t*(ρ*vi*)+ρ0v02l0rj*(vj*ρ*vi*) =-p0l0ri*p*+σ0l0rj*σij*+ρ0g0ρ*gi*|l0ρ0v02,
l0v0t0Srt*(ρ*vi*)+rj*(vj*ρ*vi*) =-p0ρ0v02Maiso-2ri*p*+σ0ρ0v02Re-1rj*σij*+ρ0g0l0ρ0v02Fr-1ρ*gi*.

The occurring dimensionless numbers are the isothermal Mach number Maiso, which is related to the Euler number Eu or the Ruark number Ru like Maiso2=Eu=Ru-1, the Froude number Fr, which is related to the Richardson number Ri like Fr=Ri-1 and the Reynolds number Re. All these numbers measure the importance of the term they are related to compared to the nonlinear advection term rj*(vj*ρ*vi*), e.g. for high Mach numbers, the pressure term ri*p* becomes less and less important compared to the advection term; for high Reynolds numbers the stress term rj*σij* becomes less and less important compared to the advection term, and the equation shows more and more nonlinear behavior. For a newtonian fluid with

σijηvirj=η0v0l0η*vi*rj*, (A.1)

where we introduced the dimensionless quantity η*=ηη0, we get σ0=η0v0l0 and therefore we can express the Reynolds number5252We neglected the second viscosity ζ. In principle there exists a second Reynolds number Re2=ρ0l0v0ζ0. like

Re=l0ρ0v02η0v0=ρ0l0v0η0. (A.2)

Playing the same game with the equation for the internal energy

tρeint+rjvjρeint=T(tρs+rjvjρs)-prjvj (A.3)

using eint*=eintu0 we get

l0v0t0Srt*ρ*eint*+rj*vj*ρ*eint* =-p0ρ0u0Ga1p*rj*vj*+σ0ρ0u0Ga2σij*rj*vi*.

The new dimensionless quantities that occur in the energy equation seem to have no name in the literature, but we will call them ”Gamma1” (Ga1) and ”Gamma2” (Ga2) for now, since they are related to the adiabatic coefficient. This can be seen by replacing p0p* according to equation

p0p*=(γ-1)ρ0u0ρ*eint*, (A.4)

which is valid for an an ideal, nonisothermal (γ1) gas. Doing this we get5353We cannot get rid of Ga2 in the same way, since therefore we would have to assume an equation relating σ0σij* to the internal energy. But this would only be possible, if we would assume that the internal energy is a tensorial quantity, which is not the way how internal energy is defined normally.

Srt*ρ*eint*+rj*vj*ρ*eint* =-(γ-1)ρ*eint*rj*vj*+Ga2σij*rj*vi*.

For a selfgravitating fluid we even have on more dimensionless quantity, which appears, when we write down the dimesionless form of the Poisson equation of gravity

g04πGρ0l0CGgi*ri*=ρ*. (A.5)

But this quantity CG also seems to have no name in the literature (e.g. Durst, 2007).