tjhei/aspect — pull request #3 by gassmoeller (2014-09-20 09:19:15+00:00)
Title: New zhong version a1.

Hi Timo,

Hier sind unsere Versuche für den Zhong Benchmark. Leider hat auch die Verringerung der Dichtevariation im Model die Ergebnisse des Papers nicht völlig reproduziert. Das vielversprechendste Modell (a1j.prm) hat eine rms Geschwindigkeit von 34,2 statt 32,6 m/s (im Paper sind andere Studien angegeben, die alle um 32,6 +- 0.5 liegen). Außerdem ist das Ergebnis auch optisch nicht ganz richtig (siehe Bild), wir bekommen zwar tetrahedrische Upwellings, aber die sind paarweise mit kleinen "ridges" verbunden, was eigentlich nicht der Fall sein sollte. Eine Idee, die ich noch nicht ausprobieren konnte, wäre, das Modell noch mal mit "realistischen" Werten zu rechnen (Standardwerte für alles), und die Rayleighzahl über die Viskosität einzustellen, damit sollten dann eigentlich keinerlei Probleme mit der Skalierung auftreten, und die Dichteänderungen im Modell sollten auch relativ klein sein (im Bereich von wenigen Prozent). Bei a1j.prm habe ich das so ähnlich gemacht, allerdings mit nichtdimensionalen Einheiten. Dort kann ich die Skalierung zwischen alpha und g nicht viel über den Faktor 100 legen (alpha=1e-2, g=100), bei mehr beginnt sich die Lösung zu verändern.

![zhong-a1-viscosity-scaling](https://cloud.githubusercontent.com/assets/7582930/4344796/603bfcae-40a6-11e4-97f0-e36a0c62c5de.png)

Add CMakeLists for zhong benchmark.

Update zhong.cc.

Add parameter file that leads to same result as a1g but withouth zhong material model.

Add a model setup with simple material model, reduced density changes and rayleigh number scaling via viscosity.