Die beiden Einheitsvektoren der Ebene werden mit
\[ \begin{align} \mathbf{i} &\coloneqq \left(\begin{array}{c} 1\\ 0 \end{array}\right),\\ \mathbf{j} &\coloneqq \left(\begin{array}{c} 0\\ 1 \end{array}\right)\\ \end{align} \]
bezeichnet. Nun definiert man weiter ein dreielementiges Erzeugendensystem $\left(\mathbf{i}_1, \mathbf{i}_2 \mathbf{i}_3\right)$ durch
\[ \begin{align} \mathbf{i}_1 &\coloneqq \left(\begin{array}{c} 1\\ 0 \end{array}\right) = \mathbf{i},\\ \mathbf{i}_2 &\coloneqq \left(\begin{array}{c} -\sin\left(30^\circ\right)\\ \cos\left(30^\circ\right) \end{array}\right) = \left(\begin{array}{c} -\frac{1}{2}\\ \frac{\sqrt{3}}{2} \end{array}\right) = -\frac{1}{2}\mathbf{i} + \frac{\sqrt{3}}{2}\mathbf{j},\\ \mathbf{i}_3 &\coloneqq \left(\begin{array}{c} -\frac{1}{2}\\ -\frac{\sqrt{3}}{2} \end{array}\right) = -\frac{1}{2}\mathbf{i} - \frac{\sqrt{3}}{2}\mathbf{j}, \end{align} \]
dessen Elemente jeweils um 120$^\circ$ gegeneinander rotiert sind. Weiter definiert man ein hiergegen um 90$^\circ$ rotiertes Erzeugendensystem $\left(\mathbf{j}_1, \mathbf{j}_2 \mathbf{j}_3\right)$ durch
\[ \begin{align} \mathbf{j}_1 &\coloneqq \mathbf{k}\times\mathbf{i}_1 = \left(\begin{array}{c} 0\\ 0\\ 1 \end{array}\right)\times\left(\begin{array}{c} 1\\ 0\\ 0 \end{array}\right) = \left(\begin{array}{c} 0\\ 1 \end{array}\right),\\ \mathbf{j}_2 &\coloneqq \mathbf{k}\times\mathbf{i}_2 = \left(\begin{array}{c} 0\\ 0\\ 1 \end{array}\right)\times\left(\begin{array}{c} -\frac{1}{2}\\ \frac{\sqrt{3}}{2}\\ 0 \end{array}\right) = \left(\begin{array}{c} -\frac{\sqrt{3}}{2}\\ -\frac{1}{2} \end{array}\right) = -\frac{\sqrt{3}}{2}\mathbf{i} - \frac{1}{2}\mathbf{j},\\ \mathbf{j}_3 &\coloneqq \mathbf{k}\times\mathbf{i}_3 = \left(\begin{array}{c} 0\\ 0\\ 1 \end{array}\right)\times\left(\begin{array}{c} -\frac{1}{2}\\ -\frac{\sqrt{3}}{2}\\ 0 \end{array}\right) = \left(\begin{array}{c} \frac{\sqrt{3}}{2}\\ -\frac{1}{2} \end{array}\right) = \frac{\sqrt{3}}{2}\mathbf{i} - \frac{1}{2}\mathbf{j}. \end{align} \]
Man beobachtet
\[ \begin{align} \mathbf{i}_1 + \mathbf{i}_2 + \mathbf{i}_3 = \mathbf{0}, & {} & \mathbf{j}_1 + \mathbf{j}_2 + \mathbf{j}_3 = \mathbf{0}. \end{align} \]
Für einen zweidimensionalen Vektor $\mathbf{v}$ kann man
\[ \begin{align} \mathbf{v} = u\mathbf{i} + v\mathbf{j} \end{align} \]
schreiben, hierbei gelten
\[ \begin{align} u = \mathbf{i}\cdot\mathbf{v}, & {} & v = \mathbf{j}\cdot\mathbf{v}. \end{align} \]
Da die $\mathbf{i}_k$, $\mathbf{j}_l$ jeweils paarweise linear unabhängig sind, kann man notieren
\[ \begin{align} \mathbf{v} = u_k'\mathbf{i}_k + u_l'\mathbf{i}_l = v_k'\mathbf{j}_k + v_l'\mathbf{j}_l. \end{align} \]
Da die Auswahl der $k$, $l$ nicht eindeutig ist und man auch alle drei Einheitsvektoren verwenden könnte, kann man eine weitere lineare Bedingung einführen. Man nutzt diese Freiheit, um die Schreibweise
\[ \begin{align} \mathbf{v} &= \frac{2}{3}\left(u_1\mathbf{i}_1 + u_2\mathbf{i}_2 + u_3\mathbf{i}_3\right) = \frac{2}{3}\left(v_1\mathbf{j}_1 + v_2\mathbf{j}_2 + v_3\mathbf{j}_3\right) \end{align} \]
zu fordern. Dies ist mit
\[ \begin{align} u_k = \mathbf{i}_k\cdot\mathbf{v}, & {} & v_k = \mathbf{j}_k\cdot\mathbf{v} \end{align} \]
erfüllt, denn hieraus folgt
\[ \begin{align} \mathbf{v} &= u\mathbf{i} + v\mathbf{j} = \frac{2}{3}\left(\frac{6}{4}u\mathbf{i} + \frac{6}{4}v\mathbf{j}\right)\nonumber\\ &= \frac{2}{3}\left(u\mathbf{i} + \frac{1}{4}u\mathbf{i} - \frac{\sqrt{3}}{4}u\mathbf{j} - \frac{\sqrt{3}}{4}v\mathbf{i} + \frac{3}{4}v\mathbf{j} + \frac{1}{4}u\mathbf{i} + \frac{\sqrt{3}}{4}u\mathbf{j} + \frac{\sqrt{3}}{4}v\mathbf{i} + \frac{3}{4}v\mathbf{j}\right)\nonumber\\ &= \frac{2}{3}\left(u\mathbf{i} + \left(-\frac{1}{2}u + \frac{\sqrt{3}}{2}v\right)\left(-\frac{1}{2}\mathbf{i} + \frac{\sqrt{3}}{2}\mathbf{j}\right) + \left(-\frac{1}{2}u - \frac{\sqrt{3}}{2}v\right)\left(-\frac{1}{2}\mathbf{i} - \frac{\sqrt{3}}{2}\mathbf{j}\right)\right)\nonumber\\ &= \frac{2}{3}\left(u_1\mathbf{i}_1 + u_2\mathbf{i}_2 + u_3\mathbf{i}_3\right), \end{align} \]
wobei
\[ \begin{align} u_1 &= u, & {} & u_2 = -\frac{1}{2}u + \frac{\sqrt{3}}{2}v, & {} & u_3 = -\frac{1}{2}u - \frac{\sqrt{3}}{2}v. \end{align} \]
verwendet wurde. Dies gilt analog für die $v_k$. Hieraus folgt
\[ \begin{align} 0 = \mathbf{v}\cdot\mathbf{0} = \mathbf{v}\cdot\left(\mathbf{i}_1 + \mathbf{i}_2 + \mathbf{i}_3\right) = u_1 + u_2 + u_3\tag{27.18}\label{eq:linear_condition_trivariate} \end{align} \]
und analog für die $v_k$. Für den Gradienten $\nabla\alpha$ eines Skalarfeldes $\alpha$ erhält man nun
\[ \begin{align} \nabla\alpha = \frac{2}{3}\left[\left(\mathbf{i}_1\cdot\nabla\alpha\right)\mathbf{i}_1 + \left(\mathbf{i}_2\cdot\nabla\alpha\right)\mathbf{i}_2 + \left(\mathbf{i}_3\cdot\nabla\alpha\right)\mathbf{i}_3\right]. \end{align} \]
Wegen
\[ \begin{align} \mathbf{i}_k\cdot\nabla\alpha = \frac{\partial\alpha}{\partial x_k} \end{align} \]
folgt
\[ \begin{align} \nabla\alpha = \frac{2}{3}\left(\frac{\partial\alpha}{\partial x_1}\mathbf{i}_1 + \frac{\partial\alpha}{\partial x_2}\mathbf{i}_2 + \frac{\partial\alpha}{\partial x_3}\mathbf{i}_3\right).\tag{27.21}\label{eq:grad_three_elements} \end{align} \]
Man stellt hiermit außerdem
\[ \begin{align} \frac{\partial\alpha}{\partial x_1} + \frac{\partial\alpha}{\partial x_2} + \frac{\partial\alpha}{\partial x_3} = 0\tag{27.22}\label{eq:gradient_linear_condition} \end{align} \]
fest. Es gilt
\[ \begin{align} \mathbf{i}_1 &= \frac{1}{\sqrt{3}}\left(\mathbf{j}_3 - \mathbf{j}_2\right) \end{align} \]
und zyklisch:
\[ \begin{align} \mathbf{i}_2 = \frac{1}{\sqrt{3}}\left(\mathbf{j}_1 - \mathbf{j}_3\right), & {} & \mathbf{i}_3 = \frac{1}{\sqrt{3}}\left(\mathbf{j}_2 - \mathbf{j}_1\right) \end{align} \]
Hieraus folgt
\[ \begin{align} \nabla\alpha = \frac{2}{3\sqrt{3}}\left(\frac{\partial\alpha}{\partial x_1}\left(\mathbf{j}_3 - \mathbf{j}_2\right) + \frac{\partial\alpha}{\partial x_2}\left(\mathbf{j}_1 - \mathbf{j}_3\right) + \frac{\partial\alpha}{\partial x_3}\left(\mathbf{j}_2 - \mathbf{j}_1\right)\right). \end{align} \]
Umsortierung ergibt
\[ \begin{align} \nabla\alpha = \frac{2}{3\sqrt{3}}\left(\left(\frac{\partial\alpha}{\partial x_2} - \frac{\partial\alpha}{\partial x_3}\right)\mathbf{j}_1 + \left(\frac{\partial\alpha}{\partial x_3} - \frac{\partial\alpha}{\partial x_1}\right)\mathbf{j}_2 + \left(\frac{\partial\alpha}{\partial x_1} - \frac{\partial\alpha}{\partial x_2}\right)\mathbf{j}_3\right). \end{align} \]
Es gilt
\[ \begin{align} \mathbf{v} &= \frac{2}{3}\left(u_1\mathbf{i}_1 + u_2\mathbf{i}_2 + u_3\mathbf{i}_3\right) = \frac{2}{3}\left[u_1\left(\begin{array}{c}1\\0\end{array}\right) + u_2\left(\begin{array}{c}-\frac{1}{2}\\\frac{\sqrt{3}}{2}\end{array}\right) + u_3\left(\begin{array}{c}-\frac{1}{2}\\-\frac{\sqrt{3}}{2}\end{array}\right)\right] = \frac{2}{3}\left(\begin{array}{c}u_1 - \frac{u_2}{2} - \frac{u_3}{2}\\\frac{\sqrt{3}}{2}u_2 - \frac{\sqrt{3}}{2}u_3\end{array}\right)\nonumber\\ &= \left(\begin{array}{c}\frac{2}{3}u_1 - \frac{u_2}{3} - \frac{u_3}{3}\\\frac{1}{\sqrt{3}}u_2 - \frac{1}{\sqrt{3}}u_3\end{array}\right). \end{align} \]
Hieraus folgt für die Divergenz
\[ \begin{align} D \coloneqq \nabla\cdot\mathbf{v} = \frac{2}{3}\frac{u_1}{\partial x} - \frac{1}{3}\frac{u_2}{\partial x} - \frac{1}{3}\frac{u_3}{\partial x} + \frac{1}{\sqrt{3}}\frac{u_2}{\partial y} - \frac{1}{\sqrt{3}}\frac{u_3}{\partial y}. \end{align} \]
Es gilt
\[ \begin{align} \frac{\partial}{\partial x} = \frac{\partial}{\partial x_1}, \end{align} \]
wegen
\[ \begin{align} \mathbf{j} = \frac{1}{\sqrt{3}}\left(\mathbf{i}_2 - \mathbf{i}_3\right) \end{align} \]
ist außerdem
\[ \begin{align} \frac{\partial}{\partial y} = \frac{1}{\sqrt{3}}\left(\frac{\partial}{\partial x_2} - \frac{\partial}{\partial x_3}\right).\tag{27.31}\label{eq:ddy_hex} \end{align} \]
Hieraus folgt
\[ \begin{align} D &= \frac{2}{3}\frac{\partial u_1}{\partial x} - \frac{1}{3}\frac{\partial u_2}{\partial x} - \frac{1}{3}\frac{\partial u_3}{\partial x} + \frac{1}{3}\left(\frac{\partial u_2}{\partial x_2} - \frac{\partial u_2}{\partial x_3}\right) - \frac{1}{3}\left(\frac{\partial u_3}{\partial x_2} - \frac{\partial u_3}{\partial x_3}\right)\nonumber\\ &= \frac{2}{3}\frac{\partial u_1}{\partial x_1} + \frac{1}{3}\frac{u_2}{\partial x_2} + \frac{1}{3}\frac{u_3}{\partial x_3} - \frac{1}{3}\frac{u_2}{\partial x_1} - \frac{1}{3}\frac{\partial u_3}{\partial x_1} - \frac{1}{3}\frac{\partial u_2}{\partial x_3} - \frac{1}{3}\frac{\partial u_3}{\partial x_2}\nonumber\\ &= \frac{2}{3}\left(\frac{\partial u_1}{\partial x_1} + \frac{\partial u_2}{\partial x_2} + \frac{\partial u_3}{\partial x_3}\right) - \frac{1}{3}\left(\frac{\partial u_2}{\partial x_1} + \frac{\partial u_2}{\partial x_2} + \frac{\partial u_2}{\partial x_3} + \frac{\partial u_3}{\partial x_1} + \frac{\partial u_3}{\partial x_2} + \frac{\partial u_3}{\partial x_3}\right). \end{align} \]
Mit Glg. (27.22) folgt
\[ \begin{align} D = \frac{2}{3}\left(\frac{\partial u_1}{\partial x_1} + \frac{\partial u_2}{\partial x_2} + \frac{\partial u_3}{\partial x_3}\right).\tag{27.33}\label{eq:div_3elements2d} \end{align} \]
Durch Kombination mit Glg. (27.21) erhält man weiterhin
\[ \begin{align} \Delta\alpha = \frac{2}{3}\left(\frac{\partial^2\alpha}{\partial x_1^2} + \frac{\partial^2\alpha}{\partial x_2^2} + \frac{\partial^2\alpha}{\partial x_3^2}\right). \end{align} \]
Es gilt
\[ \begin{align} u_1 &= \mathbf{v}\cdot\mathbf{i}_1 = \frac{2}{3}\left(v_1\mathbf{j}_1\cdot\mathbf{i}_1 + v_2\mathbf{j}_2\cdot\mathbf{i}_1 + v_3\mathbf{j}_3\cdot\mathbf{i}_1\right) = \frac{2}{3}\left(0\cdot v_1 - \frac{\sqrt{3}}{2}v_2 + \frac{\sqrt{3}}{2}v_3\right)\nonumber\\ &= \frac{1}{\sqrt{3}}\left(v_3 - v_2\right). \end{align} \]
Durch zyklisches Verschieben der Indizes erhält man
\[ \begin{align} u_2 = \frac{1}{\sqrt{3}}\left(v_1 - v_3\right), & {} & u_3 = \frac{1}{\sqrt{3}}\left(v_2 - v_1\right). \end{align} \]
Hiermit erhält man eine weitere Darstellung der Divergenz:
\[ \begin{align} D = \frac{2}{3\sqrt{3}}\left[\left(\frac{\partial v_2}{\partial x_3} - \frac{\partial v_3}{\partial x_2}\right) + \left(\frac{\partial v_3}{\partial x_1} - \frac{\partial v_1}{\partial x_3}\right) + \left(\frac{\partial v_1}{\partial x_2} - \frac{\partial v_2}{\partial x_1}\right)\right]\tag{27.37}\label{eq:div_3elements2d_mod} \end{align} \]
Es gilt
\[ \begin{align} \frac{\partial v_2}{\partial x_3} - \frac{\partial v_3}{\partial x_2} &= \frac{\partial v_2}{\partial x_3} + \frac{\partial v_2}{\partial x_2} - \frac{\partial v_2}{\partial x_2} - \frac{\partial v_3}{\partial x_2} = \left(\frac{\partial}{\partial x_3} + \frac{\partial}{\partial x_2}\right)v_2 - \frac{\partial}{\partial x_2}\left(v_2 + v_3\right)\nonumber\\ &= -\frac{\partial v_2}{\partial x_1} + \frac{\partial v_1}{\partial x_2} = \frac{\partial v_1}{\partial x_2} - \frac{\partial v_2}{\partial x_1}. \end{align} \]
Durch zyklische Vertauschung dieser Gleichung erkennt man, dass alle drei eingeklammerten Summanden in Glg. (27.37) gleich sind. Daraus folgen
\[ \begin{align} D &= \frac{2}{\sqrt{3}}\left(\frac{\partial v_2}{\partial x_3} - \frac{\partial v_3}{\partial x_2}\right),\\ D &= \frac{2}{\sqrt{3}}\left(\frac{\partial v_3}{\partial x_1} - \frac{\partial v_1}{\partial x_3}\right),\\ D &= \frac{2}{\sqrt{3}}\left(\frac{\partial v_1}{\partial x_2} - \frac{\partial v_2}{\partial x_1}\right). \end{align} \]
Es gilt
\[ \begin{align} \mathbf{v} &= \frac{2}{3}\left(v_1\mathbf{j}_1 + v_2\mathbf{j}_2 + v_3\mathbf{j}_3\right) = \frac{2}{3}\left[v_1\left(\begin{array}{c}0\\1\end{array}\right) + v_2\left(\begin{array}{c}-\frac{\sqrt{3}}{2}\\-\frac{1}{2}\end{array}\right) + v_3\left(\begin{array}{c}\frac{\sqrt{3}}{2}\\-\frac{1}{2}\end{array}\right)\right]\nonumber\\ &= \left(\begin{array}{c}\frac{1}{\sqrt{3}}v_3 - \frac{1}{\sqrt{3}}v_2\\\frac{2}{3}v_1 - \frac{v_2}{3} - \frac{v_3}{3}\end{array}\right). \end{align} \]
Für die Vorticity erhält man somit
\[ \begin{align} \zeta &= \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} = \frac{\partial}{\partial x}\left(\frac{2}{3}v_1 - \frac{v_2}{3} - \frac{v_3}{3}\right) - \frac{\partial}{\partial y}\left(\frac{1}{\sqrt{3}}v_3 - \frac{1}{\sqrt{3}}v_2\right)\nonumber\\ &= \frac{\partial}{\partial x_1}\left(\frac{2}{3}v_1 - \frac{v_2}{3} - \frac{v_3}{3}\right) - \frac{1}{\sqrt{3}}\left(\frac{\partial}{\partial x_2} - \frac{\partial}{\partial x_3}\right)\left(\frac{1}{\sqrt{3}}v_3 - \frac{1}{\sqrt{3}}v_2\right)\nonumber\\ &= \frac{2}{3}\frac{\partial v_1}{\partial x_1} + \frac{1}{3}\frac{\partial v_2}{\partial x_2} + \frac{1}{3}\frac{\partial v_3}{\partial x_3} - \frac{1}{3}\frac{\partial v_2}{\partial x_1} - \frac{1}{3}\frac{\partial v_3}{\partial x_1} - \frac{1}{3}\frac{\partial v_3}{\partial x_2} - \frac{1}{3}\frac{\partial v_2}{\partial x_3}\nonumber\\ &= \frac{2}{3}\left(\frac{\partial v_1}{\partial x_1} + \frac{\partial v_2}{\partial x_2} + \frac{\partial v_3}{\partial x_3}\right) - \frac{1}{3}\left(\frac{\partial v_2}{\partial x_1} + \frac{\partial v_2}{\partial x_2} + \frac{\partial v_2}{\partial x_3} + \frac{\partial v_3}{\partial x_1} + \frac{\partial v_3}{\partial x_2} + \frac{\partial v_3}{\partial x_3}\right). \end{align} \]
Mit Glg. (27.22) folgt
\[ \begin{align} \zeta = \frac{2}{3}\left(\frac{\partial v_1}{\partial x_1} + \frac{\partial v_2}{\partial x_2} + \frac{\partial v_3}{\partial x_3}\right). \end{align} \]
Es gilt
\[ \begin{align} v_1 &= \mathbf{v}\cdot\mathbf{j}_1 = \frac{2}{3}\left(u_1\mathbf{i}_1 + u_2\mathbf{i}_2 + u_3\mathbf{i}_3\right)\cdot\mathbf{j}_1 = \frac{2}{3}\left(u_1\mathbf{i}_1 + u_2\mathbf{i}_2 + u_3\mathbf{i}_3\right)\cdot\mathbf{j}\nonumber\\ &= \frac{2}{3}\left(\frac{\sqrt{3}}{2}u_2 - \frac{\sqrt{3}}{2}u_3\right) = \frac{1}{\sqrt{3}}\left(u_2 - u_3\right).\tag{27.45}\label{eq:turn_3elemets2d} \end{align} \]
Durch zyklisches Verschieben der Indizes erhält man
\[ \begin{align} v_2 = \frac{1}{\sqrt{3}}\left(u_3 - u_1\right), & {} & v_3 = \frac{1}{\sqrt{3}}\left(u_1 - u_2\right). \end{align} \]
Hieraus folgt eine weitere Darstellung der Vorticity:
\[ \begin{align} \zeta = \frac{2}{3\sqrt{3}}\left(\frac{\partial u_2}{\partial x_1} - \frac{\partial u_3}{\partial x_1} + \frac{\partial u_3}{\partial x_2} - \frac{\partial u_1}{\partial x_2} + \frac{\partial u_1}{\partial x_3} - \frac{\partial u_2}{\partial x_3}\right) \end{align} \]
Durch Umsortierung erhält man
\[ \begin{align} \zeta = \frac{2}{3\sqrt{3}}\left[\left(\frac{\partial u_3}{\partial x_2} - \frac{\partial u_2}{\partial x_3}\right) + \left(\frac{\partial u_1}{\partial x_3} - \frac{\partial u_3}{\partial x_1}\right) + \left(\frac{\partial u_2}{\partial x_1} - \frac{\partial u_1}{\partial x_2}\right)\right]. \end{align} \]
Analog zum Fall von Glg. (27.37) gelten auch hier
\[ \begin{align} \zeta &= \frac{2}{\sqrt{3}}\left(\frac{\partial u_3}{\partial x_2} - \frac{\partial u_2}{\partial x_3}\right),\tag{27.49}\label{eq:hex_curl_symm_0}\\ \zeta &= \frac{2}{\sqrt{3}}\left(\frac{\partial u_1}{\partial x_3} - \frac{\partial u_3}{\partial x_1}\right),\tag{27.50}\label{eq:hex_curl_symm_1}\\ \zeta &= \frac{2}{\sqrt{3}}\left(\frac{\partial u_2}{\partial x_1} - \frac{\partial u_1}{\partial x_2}\right).\tag{27.51}\label{eq:hex_curl_symm_2} \end{align} \]
Mit dem Hauptsatz der Vektoranalysis kann man das Vektorfeld $\mathbf{v}$ in der Form
\[ \begin{align} \mathbf{v} = \mathbf{k}\times\nabla\psi + \nabla\chi \end{align} \]
notieren. Dabei gelten
\[ \begin{align} \mathbf{k}\times\nabla\psi &= \frac{2}{3}\left(\frac{\partial\psi}{\partial x_1}\mathbf{j}_1 + \frac{\partial\psi}{\partial x_2}\mathbf{j}_2 + \frac{\partial\psi}{\partial x_3}\mathbf{j}_3\right),\\ \nabla\chi &= \frac{2}{3}\left(\frac{\partial\chi}{\partial x_1}\mathbf{i}_1 + \frac{\partial\chi}{\partial x_2}\mathbf{i}_2 + \frac{\partial\chi}{\partial x_3}\mathbf{i}_3\right). \end{align} \]
Hieraus folgen
\[ \begin{align} u_1 &= \frac{1}{\sqrt{3}}\left(\frac{\partial\psi}{\partial x_3} - \frac{\partial\psi}{\partial x_2}\right) + \frac{\partial\chi}{\partial x_1} = -\frac{1}{\sqrt{3}}\left(\frac{\partial\psi}{\partial x_2} - \frac{\partial\psi}{\partial x_3}\right) + \frac{\partial\chi}{\partial x_1}\tag{27.55}\label{eq:trivariate_helmholtz_eq_0},\\ u_2 &= -\frac{1}{\sqrt{3}}\left(\frac{\partial\psi}{\partial x_3} - \frac{\partial\psi}{\partial x_1}\right) + \frac{\partial\chi}{\partial x_2}\tag{27.56}\label{eq:trivariate_helmholtz_eq_1},\\ u_3 &= -\frac{1}{\sqrt{3}}\left(\frac{\partial\psi}{\partial x_1} - \frac{\partial\psi}{\partial x_2}\right) + \frac{\partial\chi}{\partial x_3}\tag{27.57}\label{eq:trivariate_helmholtz_eq_2}. \end{align} \]
und
\[ \begin{align} v_1 &= \frac{\partial\psi}{\partial x_1} + \frac{1}{\sqrt{3}}\left(\frac{\partial\chi}{\partial x_2} - \frac{\partial\chi}{\partial x_3}\right)\tag{27.58}\label{eq:trivariate_helmholtz_eq_3},\\ v_2 &= \frac{\partial\psi}{\partial x_2} + \frac{1}{\sqrt{3}}\left(\frac{\partial\chi}{\partial x_3} - \frac{\partial\chi}{\partial x_1}\right)\tag{27.59}\label{eq:trivariate_helmholtz_eq_4},\\ v_3 &= \frac{\partial\psi}{\partial x_3} + \frac{1}{\sqrt{3}}\left(\frac{\partial\chi}{\partial x_1} - \frac{\partial\chi}{\partial x_2}\right)\tag{27.60}\label{eq:trivariate_helmholtz_eq_5}. \end{align} \]
Wegen
\[ \begin{align} D = \Delta\chi, & {} & \zeta = \Delta\psi \end{align} \]
gelten weiter
\[ \begin{align} \Delta u_1 &= -\frac{1}{\sqrt{3}}\left(\frac{\partial\zeta}{\partial x_2} - \frac{\partial\zeta}{\partial x_3}\right) + \frac{\partial D}{\partial x_1},\tag{27.62}\label{eq:hex_laplace_u_1}\\ \Delta u_2 &= -\frac{1}{\sqrt{3}}\left(\frac{\partial\zeta}{\partial x_3} - \frac{\partial\zeta}{\partial x_1}\right) + \frac{\partial D}{\partial x_2},\tag{27.63}\label{eq:hex_laplace_u_2}\\ \Delta u_3 &= -\frac{1}{\sqrt{3}}\left(\frac{\partial\zeta}{\partial x_1} - \frac{\partial\zeta}{\partial x_2}\right) + \frac{\partial D}{\partial x_3}\tag{27.64}\label{eq:hex_laplace_u_3}, \end{align} \]
\[ \begin{align} \Delta v_1 &= \frac{\partial\zeta}{\partial x_1} + \frac{1}{\sqrt{3}}\left(\frac{\partial D}{\partial x_2} - \frac{\partial D}{\partial x_3}\right),\\ \Delta v_2 &= \frac{\partial\zeta}{\partial x_2} + \frac{1}{\sqrt{3}}\left(\frac{\partial D}{\partial x_3} - \frac{\partial D}{\partial x_1}\right),\\ \Delta v_3 &= \frac{\partial\zeta}{\partial x_3} + \frac{1}{\sqrt{3}}\left(\frac{\partial D}{\partial x_1} - \frac{\partial D}{\partial x_2}\right). \end{align} \]
Sei $\alpha$ ein Feld der Form
\[ \begin{align} \alpha\left(\mathbf{r}\right) = \exp\left(i\mathbf{k}\cdot\mathbf{r}\right). \end{align} \]
Man bezeichnet den Abstand zwischen den Zentren zweier Hexagone mit $d$ (diese Größe wird von nun an als Gitterkonstante bezeichnet). Man definiert nun für $1 \leq j \leq 3$ den zentralen Differenzenquotienten in j-Richtung durch
\[ \begin{align} \delta_j\alpha &= \frac{\alpha\left(\mathbf{r} + \mathbf{i}_j\frac{d}{2}\right) - \alpha\left(\mathbf{r} - \mathbf{i}_j\frac{d}{2}\right)}{d}. \end{align} \]
Hierfür gilt
\[ \begin{align} \delta_j\alpha &= \frac{\exp\left[i\mathbf{k}\cdot\left(\mathbf{r} + \mathbf{i}_j\frac{d}{2}\right)\right] - \exp\left[i\mathbf{k}\cdot\left(\mathbf{r} - \mathbf{i}_j\frac{d}{2}\right)\right]}{d}\nonumber\\ &= \exp\left(i\mathbf{k}\cdot\mathbf{r}\right)\frac{\exp\left(i\mathbf{k}\cdot\mathbf{i}_j\frac{d}{2}\right) - \exp\left(-i\mathbf{k}\cdot\mathbf{i}_j\frac{d}{2}\right)}{d} = \alpha\left(\mathbf{r}\right)\frac{\exp\left(i\mathbf{k}\cdot\mathbf{i}_j\frac{d}{2}\right) - \exp\left(-i\mathbf{k}\cdot\mathbf{i}_j\frac{d}{2}\right)}{d}\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{\exp\left(ik_j\frac{d}{2}\right) - \exp\left(-ik_j\frac{d}{2}\right)}{d} = \alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_jd}{2}\right). \end{align} \]
Definiert man
\[ \begin{align} s_j \coloneqq \sin\left(\frac{k_jd}{2}\right), \end{align} \]
kann man dies in der Form
\[ \begin{align} \delta_j\alpha = \frac{2i}{d}s_j\alpha \end{align} \]
notieren. Man definiert nun weiterhin einen einfachen Mittelungsoperator in j-Richtung durch
\[ \begin{align} \newoverline{\alpha}^{(j)} & \coloneqq \frac{\alpha\left(\mathbf{r} + \mathbf{i}_j\frac{d}{2}\right) + \alpha\left(\mathbf{r} - \mathbf{i}_j\frac{d}{2}\right)}{2} = \frac{\exp\left[i\mathbf{k}\cdot\left(\mathbf{r} + \mathbf{i}_j\frac{d}{2}\right)\right] + \exp\left[i\mathbf{k}\cdot\left(\mathbf{r} - \mathbf{i}_j\frac{d}{2}\right)\right]}{2}\nonumber\\ &= \exp\left(i\mathbf{k}\cdot\mathbf{r}\right)\frac{\exp\left(i\mathbf{k}\cdot\mathbf{i}_j\frac{d}{2}\right) + \exp\left(-i\mathbf{k}\cdot\mathbf{i}_j\frac{d}{2}\right)}{2}\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{\exp\left(ik_j\frac{d}{2}\right) + \exp\left(-ik_j\frac{d}{2}\right)}{2} = \alpha\left(\mathbf{r}\right)\cos\left(\frac{k_jd}{2}\right).\tag{27.73}\label{eq:hex_simple_average} \end{align} \]
Definiert man
\[ \begin{align} c_j \coloneqq \cos\left(\frac{k_jd}{2}\right), \end{align} \]
kann man dies in der Form
\[ \begin{align} \newoverline{\alpha}^{(j)} = c_j\alpha \end{align} \]
notieren.
Die linearisierten shallow water equations Glg.en (13.173) - (13.174) auf der f-Ebene kann man mit dem Geopotential $\Phi \coloneqq gh$ in der Form
\[ \begin{align} \frac{\partial\Phi}{\partial t} + \Phi_0\nabla\cdot\mathbf{v} = 0, & {} & \frac{\partial\mathbf{v}}{\partial t} + f_0\mathbf{k}\times\mathbf{v} - \nabla\Phi = 0 \end{align} \]
notieren. Hierbei ist $\Phi_0 \coloneqq gD$. Mit den Glg.en (27.33) und (27.45) führt dies auf die Evolutionsgleichungen für die Komponenten bezüglich des dreielementigen Erzeugendensystems:
\[ \begin{align} \frac{\partial\Phi}{\partial t} + \frac{2}{3}\Phi_0\left(\frac{\partial u_1}{\partial x_1} + \frac{\partial u_2}{\partial x_2} + \frac{\partial u_3}{\partial x_3}\right) &= 0,\\ \frac{\partial u_1}{\partial t} - \frac{f_0}{\sqrt{3}}\left(u_2 - u_3\right) + \frac{\partial\Phi}{\partial x_1} &= 0,\\ \frac{\partial u_2}{\partial t} - \frac{f_0}{\sqrt{3}}\left(u_3 - u_1\right) + \frac{\partial\Phi}{\partial x_2} &= 0,\\ \frac{\partial u_3}{\partial t} - \frac{f_0}{\sqrt{3}}\left(u_1 - u_2\right) + \frac{\partial\Phi}{\partial x_3} &= 0. \end{align} \]
Räumliche Diskretisierung ergibt
\[ \begin{align} \frac{\partial\Phi}{\partial t} + \frac{2}{3}\Phi_0\left(\delta_1u_1 + \delta_2u_2 + \delta_3u_3\right) &= 0,\tag{27.81}\label{eq:swe_lin_hex_0}\\ \frac{\partial u_1}{\partial t} - \frac{f_0}{\sqrt{3}}\left(\newoverline{u_2}^{3} - \newoverline{u_3}^{(2)}\right) + \delta_1\Phi &= 0,\tag{27.82}\label{eq:swe_lin_hex_1}\\ \frac{\partial u_2}{\partial t} - \frac{f_0}{\sqrt{3}}\left(\newoverline{u_3}^{(1)} - \newoverline{u_1}^{3}\right) + \delta_2\Phi &= 0,\tag{27.83}\label{eq:swe_lin_hex_2}\\ \frac{\partial u_3}{\partial t} - \frac{f_0}{\sqrt{3}}\left(\newoverline{u_1}^{(2)} - \newoverline{u_2}^{(1)}\right) + \delta_3\Phi &= 0.\tag{27.84}\label{eq:swe_lin_hex_3} \end{align} \]
Dabei wurde für die Rekonstruktion der Windkomponenten für die Berechnung der Coriolisbeschleunigung das arithmetische Mittel der zwei nächsten Geschwindigkeitskomponenten in der gesuchten Richtung verwendet. Setzt man für alle Felder eine monochromatische ebene Welle an, also
\[ \begin{align} \psi\left(\mathbf{r}, t\right) = \newhat{\psi}\exp\left(i\mathbf{k}\cdot\mathbf{r} - i\omega t\right) \end{align} \]
für jedes Feld $\psi$ mit einer komplexen Amplitude $\newhat{\psi}$, erhält man
\[ \begin{align} -i\omega\newhat{\Phi} + \frac{4i}{3d}\Phi_0\left(s_1\newhat{u}_1 + s_2\newhat{u}_2 + s_3\newhat{u}_3\right) &= 0,\\ -i\omega\newhat{u}_1 - \frac{f_0}{\sqrt{3}}\left(c_3\newhat{u}_2 - c_2\newhat{u}_3\right) + \frac{2is_1}{d}\newhat{\Phi} &= 0,\\ -i\omega\newhat{u}_2 - \frac{f_0}{\sqrt{3}}\left(c_1\newhat{u}_3 - c_3\newhat{u}_1\right) + \frac{2is_1}{d}\newhat{\Phi} &= 0,\\ -i\omega\newhat{u}_3 - \frac{f_0}{\sqrt{3}}\left(c_2\newhat{u}_1 - c_1\newhat{u}_2\right) + \frac{2is_1}{d}\newhat{\Phi} &= 0. \end{align} \]
Dabei ist zu beachten, dass zwischen den Gleichungen ein Phasenversatz besteht, welcher daraus entsteht, dass die einzelnen Gleichungen an unterschiedlichen Orten gelten. Dieser wurde herausgekürzt. In Matrixform erhält man
\[ \begin{align} \left(\begin{array}{cccc} -i\omega & \frac{4i}{3d}\Phi_0s_1 & \frac{4i}{3d}\Phi_0s_2 & \frac{4i}{3d}\Phi_0s_3 \\ \frac{2is_1}{d} & -i\omega & -\frac{f_0}{\sqrt{3}}c_3 & \frac{f_0}{\sqrt{3}}c_2\\ \frac{2is_2}{d} & \frac{f_0}{\sqrt{3}}c_3 & -i\omega & -\frac{f_0}{\sqrt{3}}c_1\\ \frac{2is_3}{d} & -\frac{f_0}{\sqrt{3}}c_2 & \frac{f_0}{\sqrt{3}}c_1 & -i\omega \end{array}\right)\left(\begin{array}{c} \newhat{\Phi}\\ \newhat{u}_1\\ \newhat{u}_2\\ \newhat{u}_3 \end{array}\right) = \mathbf{0}. \end{align} \]
Multipliziert man dies mit der imaginären Einheit $i$, erhält man
\[ \begin{align} \left(\begin{array}{cccc} \omega & -\frac{4}{3d}\Phi_0s_1 & -\frac{4}{3d}\Phi_0s_2 & -\frac{4}{3d}\Phi_0s_3 \\ -\frac{2s_1}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_3 & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & \omega & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \frac{if_0}{\sqrt{3}}c_1 & \omega \end{array}\right)\left(\begin{array}{c} \newhat{\Phi}\\ \newhat{u}_1\\ \newhat{u}_2\\ \newhat{u}_3 \end{array}\right) = \mathbf{0}. \end{align} \]
Nichttriviale Lösungen exisitieren genau dann, wenn die Determinante der Koeffizientenmatrix dieses linearen Gleichungssystem verschwindet:
\[ \begin{align} \left|\begin{array}{cccc} \omega & -\frac{4}{3d}\Phi_0s_1 & -\frac{4}{3d}\Phi_0s_2 & -\frac{4}{3d}\Phi_0s_3 \\ -\frac{2s_1}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_3 & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & \omega & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \frac{if_0}{\sqrt{3}}c_1 & \omega \end{array}\right| &\stackrel{!}{=} 0\nonumber\\ \Leftrightarrow\omega\left|\begin{array}{ccc} \omega & -\frac{if_0}{\sqrt{3}}c_3 & \frac{if_0}{\sqrt{3}}c_2\\ \frac{if_0}{\sqrt{3}}c_3 & \omega & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{if_0}{\sqrt{3}}c_2 & \frac{if_0}{\sqrt{3}}c_1 & \omega \end{array} \right| + \frac{4}{3d}\Phi_0s_1\left|\begin{array}{ccc} -\frac{2s_1}{d} & -\frac{if_0}{\sqrt{3}}c_3 & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & \frac{if_0}{\sqrt{3}}c_1 & \omega \end{array}\right| & \nonumber\\ - \frac{4}{3d}\Phi_0s_2\left|\begin{array}{ccc} -\frac{2s_1}{d} & \omega & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \omega \end{array} \right| + \frac{4}{3d}\Phi_0s_3\left|\begin{array}{ccc} -\frac{2s_1}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_3\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & \omega\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \frac{if_0}{\sqrt{3}}c_1 \end{array} \right| &= 0 \end{align} \] \[ \begin{align} \Leftrightarrow \omega^4 - \frac{\omega^2f_0^2}{3}\left(c_1^2 + c_2^2 + c_3^2\right) + \frac{4}{3d}\Phi_0s_1\left|\begin{array}{ccc} -\frac{2s_1}{d} & -\frac{if_0}{\sqrt{3}}c_3 & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & \frac{if_0}{\sqrt{3}}c_1 & \omega \end{array}\right| & \nonumber\\ -\frac{4}{3d}\Phi_0s_2\left|\begin{array}{ccc} -\frac{2s_1}{d} & \omega & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \omega \end{array}\right| + \frac{4}{3d}\Phi_0s_3\left|\begin{array}{ccc} -\frac{2s_1}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_3\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & \omega\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \frac{if_0}{\sqrt{3}}c_1 \end{array} \right| &= 0 \end{align} \]
Es gilt
\[ \begin{align} \frac{4}{3d}\Phi_0s_1\left|\begin{array}{ccc} -\frac{2s_1}{d} & -\frac{if_0}{\sqrt{3}}c_3 & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & \frac{if_0}{\sqrt{3}}c_1 & \omega \end{array}\right| &= -\frac{8\omega^2}{3d^2}\Phi_0s_1^2 + \frac{8\Phi_0f_0^2}{9d^2}\left(s_1c_1s_3c_3 + s_1c_1s_2c_2 + s_1c_1s_1c_1\right) + \frac{8if_0}{3\sqrt{3}d}s_1\left(s_3c_2 - c_3s_2\right)\nonumber\\ &= -\frac{8\omega^2}{3d^2}\Phi_0s_1^2 + \frac{8\Phi_0f_0^2}{9d^2}s_1c_1\left(s_1c_1 + s_2c_2 + s_3c_3\right) - \frac{8if_0}{3\sqrt{3}d}s_1\left(s_2c_3 - s_3c_2\right). \end{align} \]
Die verbleibenden zwei Teildeterminanten ergeben sich analog durch zyklische Vertauschung der Indizes:
\[ \begin{align} -\frac{4}{3d}\Phi_0s_2\left|\begin{array}{ccc} -\frac{2s_1}{d} & \omega & \frac{if_0}{\sqrt{3}}c_2\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & -\frac{if_0}{\sqrt{3}}c_1\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \omega \end{array}\right| &= -\frac{8\omega^2}{3d^2}\Phi_0s_2^2 + \frac{8\Phi_0f_0^2}{9d^2}s_2c_2\left(s_1c_1 + s_2c_2 + s_3c_3\right) - \frac{8if_0}{3\sqrt{3}d}s_2\left(s_3c_1 - s_1c_3\right),\nonumber\\ \frac{4}{3d}\Phi_0s_3\left|\begin{array}{ccc} -\frac{2s_1}{d} & \omega & -\frac{if_0}{\sqrt{3}}c_3\\ -\frac{2s_2}{d} & \frac{if_0}{\sqrt{3}}c_3 & \omega\\ -\frac{2s_3}{d} & -\frac{if_0}{\sqrt{3}}c_2 & \frac{if_0}{\sqrt{3}}c_1 \end{array} \right| &= -\frac{8\omega^2}{3d^2}\Phi_0s_3^2 + \frac{8\Phi_0f_0^2}{9d^2}s_3c_3\left(s_1c_1 + s_2c_2 + s_3c_3\right) - \frac{8if_0}{3\sqrt{3}d}s_3\left(s_1c_2 - s_2c_1\right). \end{align} \]
Somit lautet das charakteristische Polynom schlussendlich
\[ \begin{align} \omega^4 - \omega^2\left[\frac{f_0^2}{3}\left(c_1^2 + c_2^2 + c_3^2\right) + \frac{8\Phi_0}{3d^2}\left(s_1^2 + s_2^2 + s_3^2\right)\right] + \frac{8\Phi_0f_0^2}{9d^2}\left(s_1c_1 + s_2c_2 + s_3c_3\right)^2 = 0.\tag{27.96}\label{eq:disp_rel_hex_c} \end{align} \]
Für den Fall gut aufgelöster Wellen $\left|\mathbf{k}\right| \ll \frac{1}{d}$ gelten
\[ \begin{align} c_j \approx 1, & {} & s_j \approx \frac{k_jd}{2}. \end{align} \]
Setzt man dies in Glg. (27.96) ein, erhält man
\[ \begin{align} \omega^4 - \omega^2\left[f_0^2 + \frac{8\Phi_0}{3d^2}\left(\frac{k_1^2d^2}{4} + \frac{k_2^2d^2}{4} + \frac{k_3^2d^2}{4}\right)\right] + \frac{8\Phi_0f_0^2}{9d^2}\underbrace{\left(\frac{k_1d}{2} + \frac{k_1d}{2} + \frac{k_1d}{2}\right)^2}_{\approx 0} &\approx 0\nonumber\\ \Leftrightarrow\omega^4 - \omega^2\left[f_0^2 + \frac{2\Phi_0}{3}\left(k_1^2 + k_2^2 + k_3^2\right)\right] &\approx 0.\nonumber\\ \Leftrightarrow\omega^4 - \omega^2\left[f_0^2 + \Phi\mathbf{k}^2\right] &\approx 0. \end{align} \]
Dabei wurde im letzten Schritt ausgenutzt, dass $k_1^2 + k_2^2 + k_3^2 = \frac{3}{2}\mathbf{k}^2$ gilt. Dies kann man sich klarmachen, indem man o. B. d. A. $\mathbf{k} \parallel \mathbf{i}_1$ annimmt, denn dann gelten $k_1 = k$, $k_2 = -k\sin\left(30^\circ\right) = -\frac{k}{2}$ sowie $k_3 = -\frac{k}{2}$. Dies führt auf die zwei Intertio-Schwere-Moden
\[ \begin{align} \omega^2 \approx f_0^2 + \Phi_0\left|\mathbf{k}\right|^2 \end{align} \]
sowie die geostrophische Mode
\[ \begin{align} \omega \approx 0. \end{align} \]
Bei weniger gut aufgelösten Wellen führt der Term $\propto\left(s_1c_1 + s_2c_2 + s_3c_3\right)^2$ jedoch dazu, dass die geostrophische Mode in zwei Teilmoden mit $\omega \not= 0$ aufspaltet. Dies ist problematisch und entspricht den Erkenntnissen von Abschn. 26.9.
Bei einer ebenen dreieckigen oder hexagonalen C-Gitter-Diskretisierung sind die Geschwindigkeitskomponenten um 120 bzw. 60$^\circ$ gegenüber einander rotiert, jedoch sind die Komponenten nicht am selben Ort definiert. Daher gilt die Bedingung Glg. (27.18) nicht ohne weiteres für die Vektorkomponenten, sondern erst nach einer Interpolation an Referenzorte, wofür hier die Zentren der Zellen verwendet werden. Diese Interpolation erfolgt mittels Operatoren
\[ \begin{align} \newtilde{u}^{(i)} \end{align} \]
mit $1 \leq i \leq 3$, wobei $u$ eine beliebige Vektorkomponente ist. Für jede Raumrichtung $i$ wird hier ein separater Mittelungsoperator zugelassen. Die Bedingungen Glg. (27.18) lassen sich also hier in der Form
\[ \begin{align} \newtilde{u}_1^{(1)} + \newtilde{u}_2^{(2)} + \newtilde{u}_3^{(3)} \hastobe 0, & {} & \newtilde{v}_1^{(1)} + \newtilde{v}_2^{(2)} + \newtilde{v}_3^{3} \hastobe 0\tag{27.102}\label{eq:linear_condition_trivariate_discrete} \end{align} \]
notieren. Diese Gleichungen sind als Forderungen an Vektorfelder $\mathbf{v}$ zu verstehen. Sollten Gradientefelder $\nabla\alpha$ die Bedingung Glg. (27.102) nicht, so hat man ein Problem, da die Berechnung von Gradientenfeldern auf dem C-Gitter festgelegt ist. Daher verwendet man als Forderung an den zu findenden Mittelungsoperator die Mittelung von Glg. (27.22):
\[ \begin{align} \newtilde{\delta_1\alpha}^{(1)} + \newtilde{\delta_2\alpha}^{(2)} + \newtilde{\delta_3\alpha}^{3} \stackrel{!}{=} 0\tag{27.103}\label{eq:linear_condition_trivariate_discrete_gradient} \end{align} \]
Für die Helmholtz-Aufteilung Glg.en (27.55) - (27.60) folgt durch Mittelung in den jeweils relevanten Raumrichtungen
\[ \begin{align} \newtilde{u}_1^{(1)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_3\psi}^{(1)} - \newtilde{\delta_2\psi}^{(1)}\right) + \newtilde{\delta_1\chi}^{(1)},\\ \newtilde{u}_2^{(2)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_1\psi}^{(2)} - \newtilde{\delta_3\psi}^{(2)}\right) + \newtilde{\delta_2\chi}^{(2)},\\ \newtilde{u}_3^{3} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_2\psi}^{3} - \newtilde{\delta_1\psi}^{3}\right) + \newtilde{\delta_3\chi}^{3},\\ \newtilde{v}_1^{(1)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_2\chi}^{(1)} - \newtilde{\delta_3\chi}^{(1)}\right) + \newtilde{\delta_1\psi}^{(1)},\\ \newtilde{v}_2^{(2)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_3\chi}^{(2)} - \newtilde{\delta_1\chi}^{(2)}\right) + \newtilde{\delta_2\psi}^{(2)},\\ \newtilde{v}_3^{3} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_1\chi}^{3} - \newtilde{\delta_2\chi}^{3}\right) + \newtilde{\delta_3\psi}^{3}. \end{align} \]
Hiermit sind die Glg.en (27.102) - (27.102) noch nicht erfüllt, vielmehr ist eine vorbereitende Mitttelung der Skalarfelder erforderlich:
\[ \begin{align} \newtilde{u}_1^{(1)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_3\newtilde{\psi}^{(2)}}^{(1)} - \newtilde{\delta_2\newtilde{\psi}^{3}}^{(1)}\right) + \newtilde{\delta_1\chi}^{(1)}\\ \newtilde{u}_2^{(2)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_1\newtilde{\psi}^{3}}^{(2)} - \newtilde{\delta_3\newtilde{\psi}^{(1)}}^{(2)}\right) + \newtilde{\delta_2\chi}^{(2)}\\ \newtilde{u}_3^{3} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_2\newtilde{\psi}^{(1)}}^{3} - \newtilde{\delta_1\newtilde{\psi}^{(2)}}^{3}\right) + \newtilde{\delta_3\chi}^{3}\\ \newtilde{v}_1^{(1)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_2\newtilde{\chi}^{3}}^{(1)} - \newtilde{\delta_3\newtilde{\chi}^{(2)}}^{(1)}\right) + \newtilde{\delta_1\psi}^{(1)}\\ \newtilde{v}_2^{(2)} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_3\newtilde{\chi}^{(1)}}^{(2)} - \newtilde{\delta_1\newtilde{\chi}^{3}}^{(2)}\right) + \newtilde{\delta_2\psi}^{(2)}\\ \newtilde{v}_3^{3} &= \frac{1}{\sqrt{3}}\left(\newtilde{\delta_1\newtilde{\chi}^{(2)}}^{3} - \newtilde{\delta_2\newtilde{\chi}^{(1)}}^{3}\right) + \newtilde{\delta_3\psi}^{3} \end{align} \]
Dies impliziert
\[ \begin{align} u_1 &= \frac{1}{\sqrt{3}}\left(\delta_3\newtilde{\psi}^{(2)} - \delta_2\newtilde{\psi}^{3}\right) + \delta_1\chi\tag{27.116}\label{eq:hex:helmholtz_0},\\ u_2 &= \frac{1}{\sqrt{3}}\left(\delta_1\newtilde{\psi}^{3} - \delta_3\newtilde{\psi}^{(1)}\right) + \delta_2\chi\tag{27.117}\label{eq:hex:helmholtz_1},\\ u_3 &= \frac{1}{\sqrt{3}}\left(\delta_2\newtilde{\psi}^{(1)} - \delta_1\newtilde{\psi}^{(2)}\right) + \delta_3\chi\tag{27.118}\label{eq:hex:helmholtz_2},\\ v_1 &= \frac{1}{\sqrt{3}}\left(\delta_2\newtilde{\chi}^{3} - \delta_3\newtilde{\chi}^{(2)}\right) + \delta_1\psi,\\ v_2 &= \frac{1}{\sqrt{3}}\left(\delta_3\newtilde{\chi}^{(1)} - \delta_1\newtilde{\chi}^{3}\right) + \delta_2\psi,\\ v_3 &= \frac{1}{\sqrt{3}}\left(\delta_1\newtilde{\chi}^{(2)} - \delta_2\newtilde{\chi}^{(1)}\right) + \delta_3\psi. \end{align} \]
Nun stellt man sich die Frage, wie der Mittelungsoperator $\newtilde{\psi}^{(j)}$ konkret aussehen muss. Für den in Glg. (27.73) eingeführten Operator gilt
\[ \begin{align} &\newoverline{\delta_1\psi}^{(1)} + \newoverline{\delta_1\psi}^{(2)} + \newoverline{\delta_1\psi}^{3}\nonumber\\ &= \psi\left(\mathbf{r}\right)\frac{2i}{d}\left[\sin\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_1d}{2}\right) + \sin\left(\frac{k_2d}{2}\right)\cos\left(\frac{k_2d}{2}\right) + \sin\left(\frac{k_3d}{2}\right)\cos\left(\frac{k_3d}{2}\right)\right]\nonumber\\ &= \psi\left(\mathbf{r}\right)\frac{i}{d}\left[2\sin\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_1d}{2}\right) + 2\sin\left(\frac{k_2d}{2}\right)\cos\left(\frac{k_2d}{2}\right) + 2\sin\left(\frac{k_3d}{2}\right)\cos\left(\frac{k_3d}{2}\right)\right]\nonumber\\ &= \psi\left(\mathbf{r}\right)\frac{i}{d}\left[\sin\left(k_1d\right) + \sin\left(k_2d\right) + \sin\left(k_3d\right)\right] \stackrel{\text{i. A.}}{\not=} 0. \end{align} \]
Dieser Operator erfüllt die Bedingung Glg. (27.103) also nicht. Daher definiert man einen zweiten Mittelungsoperator
\[ \begin{align} \newoverline{\psi}^{((1))} & \coloneqq \newoverline{\psi}^{(2, 3)} = \newoverline{\newoverline{\psi}^{3}}^{(2)} = \newoverline{\newoverline{\psi}^{(2)}}^{3} = \psi\left(\mathbf{r}\right)\cos\left(\frac{k_2d}{2}\right)\cos\left(\frac{k_3d}{2}\right).\text{ (und zyklisch)} \end{align} \]
Für diesen gilt
\[ \begin{align} \newoverline{\delta_1\alpha}^{((1))} + \newoverline{\delta_2\alpha}^{((2))} + \newoverline{\delta_3\alpha}^{(3)} &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\Big[\sin\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right)\cos\left(\frac{k_3d}{2}\right) + \cos\left(\frac{k_1d}{2}\right)\sin\left(\frac{k_2d}{2}\right)\cos\left(\frac{k_3d}{2}\right)\nonumber\\ & + \cos\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right)\sin\left(\frac{k_3d}{2}\right)\Big]\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\Big[\sin\left(\frac{k_1d + k_2d}{2}\right)\cos\left(\frac{k_3d}{2}\right) + \cos\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right)\sin\left(\frac{k_3d}{2}\right)\Big]\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\Big[\sin\left(\frac{k_1 + k_2 + k_3}{2}d\right) - \cos\left(\frac{k_1 + k_2}{2}d\right)\sin\left(\frac{k_3}{2}d\right)\nonumber\\ &+ \cos\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right)\sin\left(\frac{k_3d}{2}\right)\Big]. \end{align} \]
Mit $k_1 + k_2 + k_3 = 0$ folgt
\[ \begin{align} \newoverline{\delta_1\alpha}^{((1))} + \newoverline{\delta_2\alpha}^{((2))} + \newoverline{\delta_3\alpha}^{(3)} &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\Big[-\cos\left(\frac{k_1 + k_2}{2}d\right)\sin\left(\frac{k_3}{2}d\right) + \cos\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right)\sin\left(\frac{k_3d}{2}\right)\Big]\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_3}{2}d\right)\Big[\cos\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right) - \cos\left(\frac{k_1 + k_2}{2}d\right)\Big]\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_3}{2}d\right)\Big[\cos\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right) - \cos\left(\frac{k_1d}{2}\right)\cos\left(\frac{k_2d}{2}\right)\nonumber\\ &+ \sin\left(\frac{k_1d}{2}\right)\sin\left(\frac{k_2d}{2}\right)\Big]\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2d}{2}\right)\sin\left(\frac{k_3d}{2}\right) \stackrel{\text{i. A.}}{\not=} 0.\tag{27.125}\label{eq:thuburn_op_deriv_0} \end{align} \]
Auch dieser Operator erfüllt die Bedingung Glg. (27.103) also nicht. Man rechnet jedoch wiederum mit $k_1 + k_2 + k_3 = 0$ weiter:
\[ \begin{align} & \newoverline{\delta_1\alpha}^{((1))} + \newoverline{\delta_2\alpha}^{((2))} + \newoverline{\delta_3\alpha}^{(3)} = \alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2d}{2}\right)\sin\left(\frac{k_3d}{2}\right)\nonumber\\ &= \alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2d}{2}\right)\sin\left(-\frac{k_1 + k_2}{2}d\right) = -\alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2d}{2}\right)\sin\left(\frac{k_1 + k_2}{2}d\right)\nonumber\\ &= -\alpha\left(\mathbf{r}\right)\frac{2i}{d}\sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2d}{2}\right)\left[\sin\left(\frac{k_1}{2}d\right)\cos\left(\frac{k_2}{2}d\right) + \cos\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)\right]\nonumber\\ &= -\alpha\left(\mathbf{r}\right)\frac{2i}{d}\left[\sin\left(\frac{k_2}{2}d\right)\sin\left(\frac{k_1}{2}d\right)^2\cos\left(\frac{k_2}{2}d\right) + \sin\left(\frac{k_1}{2}d\right)\cos\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)^2\right]\nonumber\\ &= -\alpha\left(\mathbf{r}\right)\frac{i}{d}\left[\sin\left(k_2d\right)\sin\left(\frac{k_1}{2}d\right)^2 + \sin\left(k_1d\right)\sin\left(\frac{k_2}{2}d\right)^2\right]\nonumber\\ &= -\alpha\left(\mathbf{r}\right)\frac{i}{d}\left[\sin\left(k_2d\right) + \sin\left(k_1d\right) - \sin\left(k_2d\right)\cos\left(\frac{k_1}{2}d\right)^2 - \sin\left(k_1d\right)\cos\left(\frac{k_2}{2}d\right)^2\right]\nonumber \end{align} \] \[ \begin{align} &= -\alpha\left(\mathbf{r}\right)\frac{i}{d}\Big[\sin\left(k_2d\right) + \sin\left(k_1d\right) - \sin\left(k_2d\right)\cos\left(k_1\right) - \sin\left(k_1d\right)\cos\left(k_2d\right) + \sin\left(k_2d\right)\sin\left(\frac{k_1}{2}d\right)^2\nonumber\\ & + \sin\left(k_1d\right)\sin\left(\frac{k_2}{2}d\right)^2\Big]\nonumber\\ &= -\alpha\left(\mathbf{r}\right)\frac{i}{d}\Big[\sin\left(k_2d\right) + \sin\left(k_1d\right) - \sin\left[\left(k_1 + k_2\right)d\right] + \sin\left(k_2d\right)\sin\left(\frac{k_1}{2}d\right)^2 + \sin\left(k_1d\right)\sin\left(\frac{k_2}{2}d\right)^2\Big]\nonumber\\ &= -\alpha\left(\mathbf{r}\right)\frac{i}{d}\Big[\sin\left(k_2d\right) + \sin\left(k_1d\right) - \sin\left(-k_3d\right) + \sin\left(k_2d\right)\sin\left(\frac{k_1}{2}d\right)^2 + \sin\left(k_1d\right)\sin\left(\frac{k_2}{2}d\right)^2\Big]\nonumber\\ &= -\alpha\left(\mathbf{r}\right)\frac{i}{d}\Big[\sin\left(k_2d\right) + \sin\left(k_1d\right) + \sin\left(k_3d\right) + \sin\left(k_2d\right)\sin\left(\frac{k_1}{2}d\right)^2 + \sin\left(k_1d\right)\sin\left(\frac{k_2}{2}d\right)^2\Big]\nonumber\\ &= -\left(\newoverline{\delta_1\alpha}^{(1)} + \newoverline{\delta_1\alpha}^{(2)} + \newoverline{\delta_1\alpha}^{3}\right) - \alpha\left(\mathbf{r}\right)\frac{i}{d}\Big[\sin\left(k_2d\right)\sin\left(\frac{k_1}{2}d\right)^2 + \sin\left(k_1d\right)\sin\left(\frac{k_2}{2}d\right)^2\Big] \end{align} \]
Man rechnet nun
\[ \begin{align} &\sin\left(k_2d\right)\sin\left(\frac{k_1}{2}d\right)^2 + \sin\left(k_1d\right)\sin\left(\frac{k_2}{2}d\right)^2 = \sin\left(\frac{k_2}{2}d + \frac{k_2}{2}d\right)\sin\left(\frac{k_1}{2}d\right)^2 + \sin\left(\frac{k_1}{2}d + \frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)^2\nonumber\\ &= 2\sin\left(\frac{k_1}{2}d\right)^2\sin\left(\frac{k_2}{2}d\right)\cos\left(\frac{k_2}{2}d\right) + 2\sin\left(\frac{k_2}{2}d\right)^2\sin\left(\frac{k_1}{2}d\right)\cos\left(\frac{k_1}{2}d\right)\nonumber\\ &= \sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)\left[2\sin\left(\frac{k_1}{2}d\right)\cos\left(\frac{k_2}{2}d\right) + 2\cos\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)\right]\nonumber\\ &= \sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)\sin\left(\frac{k_1 + k_2}{2}d\right) = \sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)\sin\left(-\frac{k_3}{2}d\right)\nonumber\\ &= -\sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)\sin\left(\frac{k_3}{2}d\right). \end{align} \]
Somit folgt
\[ \begin{align} & \newoverline{\delta_1\alpha}^{((1))} + \newoverline{\delta_2\alpha}^{((2))} + \newoverline{\delta_3\alpha}^{(3)} = -\left(\newoverline{\delta_1\alpha}^{(1)} + \newoverline{\delta_1\alpha}^{(2)} + \newoverline{\delta_1\alpha}^{3}\right) + \alpha\left(\mathbf{r}\right)\frac{i}{d}\sin\left(\frac{k_1}{2}d\right)\sin\left(\frac{k_2}{2}d\right)\sin\left(\frac{k_3}{2}d\right)\nonumber\\ &\stackrel{\href{#eq:thuburn_op_deriv_0}{\text{Glg. (27.125)}}}{=} -\left(\newoverline{\delta_1\alpha}^{(1)} + \newoverline{\delta_1\alpha}^{(2)} + \newoverline{\delta_1\alpha}^{3}\right) + \frac{1}{2}\left(\newoverline{\delta_1\alpha}^{((1))} + \newoverline{\delta_2\alpha}^{((2))} + \newoverline{\delta_3\alpha}^{(3)}\right)\nonumber\\ &\Rightarrow\newoverline{\delta_1\alpha}^{((1))} + \newoverline{\delta_2\alpha}^{((2))} + \newoverline{\delta_3\alpha}^{(3)} = -\frac{1}{2}\left(\newoverline{\delta_1\alpha}^{(1)} + \newoverline{\delta_1\alpha}^{(2)} + \newoverline{\delta_1\alpha}^{3}\right).\tag{27.128}\label{eq:thuburn_op_deriv_1} \end{align} \]
Man macht nun den Ansatz
\[ \begin{align} \newtilde{\alpha}^{(j)} = \beta\newoverline{\alpha}^{(j)} + \left(1 - \beta\right)\newoverline{\alpha}^{((j))} \end{align} \]
mit $0 \leq \beta \leq 0$. Setzt man hier Glg. (27.128) ein, erhält man
\[ \begin{align} \newtilde{\delta_1\alpha}^{(1)} + \newtilde{\delta_2\alpha}^{(2)} + \newtilde{\delta_3\alpha}^{3} = \left(\beta - \frac{1}{2}\left(1 - \beta\right)\right)\left(\newoverline{\delta_1\alpha}^{(1)} + \newoverline{\delta_1\alpha}^{(2)} + \newoverline{\delta_1\alpha}^{3}\right) \stackrel{!}{=} 0. \end{align} \]
Dies ist erfüllt für
\[ \begin{align} \beta - \frac{1}{2}\left(1 - \beta\right) &= 0 \Leftrightarrow \beta - \frac{1}{2} + \frac{\beta}{2} = \frac{3\beta}{2} - \frac{1}{2} = 0\nonumber\\ \Leftrightarrow \beta = \frac{1}{3}. \end{align} \]
Somit gilt
\[ \begin{align} \newtilde{\alpha}^{(j)} = \frac{1}{3}\newoverline{\alpha}^{(j)} + \frac{2}{3}\newoverline{\alpha}^{((j))}.\tag{27.132}\label{eq:thuburn_op} \end{align} \]
Dieser Operator wird als Thuburn-Operator bezeichnet, da er von Thuburn in [37] erstmalig verwendet wurde.
Die Anwendung des Thuburn-Operators in den linearisierten shallow water equations auf dem hexagonalen Gitter Glg.en (27.81) - (27.84) führt auf
\[ \begin{align} \frac{\partial\Phi}{\partial t} + \frac{2}{3}\Phi_0\left(\delta_1u_1 + \delta_2u_2 + \delta_3u_3\right) &= 0,\\ \frac{\partial u_1}{\partial t} - \frac{f_0}{\sqrt{3}}\left(\newtilde{u}_2^{(3)} - \newtilde{u}_3^{(2)}\right) + \delta_1\Phi &= 0,\tag{27.134}\label{eq:x_1_momentum_thuburn_c-grid}\\ \frac{\partial u_2}{\partial t} - \frac{f_0}{\sqrt{3}}\left(\newtilde{u}_3^{(1)} - \newtilde{u}_1^{(3)}\right) + \delta_2\Phi &= 0,\\ \frac{\partial u_3}{\partial t} - \frac{f_0}{\sqrt{3}}\left(\newtilde{u}_1^{(2)} - \newtilde{u}_2^{(1)}\right) + \delta_3\Phi &= 0. \end{align} \]
Hieraus folgt unmittelbar
\[ \begin{align} \frac{\partial}{\partial t}\left(\newtilde{u}_1^{(1)} + \newtilde{u}_2^{(2)} + \newtilde{u}_3^{(3)}\right) = 0.\tag{27.137}\label{eq:thuburn_con_swe_c-grid} \end{align} \]
Die Dispersionsrelation dieses Gleichungssystems kann man analog zur Herleitung von Glg. (27.96) berechnen, indem man die Ersetzung
\[ \begin{align} c_1 \to a_1 \coloneqq\frac{c_1 + 2c_2c_3}{3}\text{ (und zyklisch)} \end{align} \]
vornimmt. Dies führt auf
\[ \begin{align} \omega^4 - \omega^2\left[\frac{f_0^2}{3}\left(a_1^2 + a_2^2 + a_3^2\right) + \frac{8\Phi_0}{3d^2}\left(s_1^2 + s_2^2 + s_3^2\right)\right] + \frac{8\Phi_0f_0^2}{9d^2}\left(s_1a_1 + s_2a_2 + s_3a_3\right)^2 = 0.\tag{27.139}\label{eq:disp_rel_hex_c_mod_0} \end{align} \]
Setzt man $\alpha = \exp\left(i\mathbf{k}\cdot\mathbf{r}\right)$ in Glg. (27.103) ein, erhält man
\[ \begin{align} \frac{2i}{d}\left(a_1s_1 + a_2s_2 + a_3s_3\right)\alpha = 0 \Leftrightarrow a_1s_1 + a_2s_2 + a_3s_3 = 0. \end{align} \]
Somit wird Glg. (27.139) zu
\[ \begin{align} \omega^4 - \omega^2\left[\frac{f_0^2}{3}\left(a_1^2 + a_2^2 + a_3^2\right) + \frac{8\Phi_0}{3d^2}\left(s_1^2 + s_2^2 + s_3^2\right)\right] = 0.\tag{27.141}\label{eq:thuburn_leads_to_geostrophic_mode} \end{align} \]
Verwendet man Glg. (27.132) zur Rekontruktion der tangentialen Geschwindigkeitskomponenten, so entsteht also eine stationäre geostrophische Mode.
Auf einem regelmäßigen Sechseckgitter mit der Gitterkonstante $d$ hat ein Dreieck die Fläche
\[ \begin{align} A_v = \frac{d}{2}l. \end{align} \]
Hierbei gilt
\[ \begin{align} l = \cos\left(30^\circ\right)d = \frac{\sqrt{3}}{2}d, \end{align} \]
woraus
\[ \begin{align} A_v = \frac{d^2\sqrt{3}}{4}\tag{27.144}\label{eq:hex_face_triangle} \end{align} \]
folgt. Somit kann man auf einem Dreieck (Superskript $t$), welches nach oben zeigt (Index $u$), die Divergenz laut Gauß'schem Satz in der Form
\[ \begin{align} D^{(t)}_u = -\frac{dv_1 + dv_2 + dv_3}{A_v} = -\frac{4}{\sqrt{3}d}\left(v_1 + v_2 + v_3\right) \end{align} \]
berechnen. Das Minuszeichen kommt daher, dass bei einem solchen Dreieck die Koordinatenachsen des j-Systems nach innen zeigen. Auf einem Dreieck, welches nach unten zeigt (Index $l$), ist dies andersherum, es gilt
\[ \begin{align} D^{(t)}_l = \frac{dv_1 + dv_2 + dv_3}{A_v} = \frac{4}{\sqrt{3}d}\left(v_1 + v_2 + v_3\right). \end{align} \]
Dies kann man in der Form
\[ \begin{align} D^{(t)}_{u, l} = \mp\frac{dv_1 + dv_2 + dv_3}{A_v} = \mp\frac{4}{\sqrt{3}d}\left(v_1 + v_2 + v_3\right) = \mp\frac{4}{\sqrt{3}d}\sum_{e \in t}v_e \end{align} \]
zusammenfassen, wobei mit $e$ die Kanten des Dreiecks bezeichnet werden. Analog gilt für die relative Vorticity laut Stokes'schem Satz
\[ \begin{align} \zeta^{(t)}_{u, l} = \pm\frac{du_1 + du_2 + du_3}{A_v} = \pm\frac{4}{\sqrt{3}d}\sum_{e \in t}u_e. \end{align} \]
$h$ bezeichne nun ein Hexagon. Dann bezeichnen $u \in h$ die drei nach oben gerichteten Dreiecke, welche einen Überlapp mit dem betrachteten Hexagon $h$ haben und analog $l \in h$ die drei nach unten gerichteten Dreiecke, welche sich mit $h$ überlappen. Mit diesen Bezeichnungen kann man
\[ \begin{align} \frac{\sqrt{3}d}{4}\left(\sum_{l \in h}D_l^{(t)} - \sum_{u \in h}D_u^{(t)}\right) &= \sum_{l \in h}\sum_{e \in l}v_e + \sum_{u \in h}\sum_{e \in u}v_e = \sum_{l \in h}\left(\sum_{e \in l, e \in h}v_e + \sum_{e \in l, e \not\in h}v_e\right) + \sum_{u \in h}\left(\sum_{e \in u, e \in h}v_e + \sum_{e \in u, e \not\in h}v_e\right) \end{align} \]
notieren. Beim letzten Umformungsschritt wurden die Kanten der Dreiecke in solche eingeteilt, welche auch Kanten des betrachteten Sechseckes sind ($e \in h$), und solche, bei denen dies nicht der Fall ist ($e \not\in h$). Ausklammern ergibt
\[ \begin{align} \frac{\sqrt{3}d}{4}\left(\sum_{l \in h}D_l^{(t)} - \sum_{u \in h}D_u^{(t)}\right) &= \sum_{l \in h}\sum_{e \in l, e \in h}v_e + \sum_{l \in h}\sum_{e \in l, e \not\in h}v_e + \sum_{u \in h}\sum_{e \in u, e \in h}v_e + \sum_{u \in h}\sum_{e \in u, e \not\in h}v_e. \end{align} \]
Wie man sieht, tragen alle Geschwindigkeitskomponenten mit positivem Vorzeichen zu $\sum_{l \in h}D_l^{(t)} - \sum_{u \in h}D_u^{(t)}$ bei. In den Summen $\sum_{l \in h}\sum_{e \in l, e \in h}$ und $\sum_{u \in h}\sum_{e \in u, e \in h}$ tritt jede Kante des Hexagons $h$ jeweils einmal auf. Somit gilt
\[ \begin{align} \frac{\sqrt{3}d}{4}\left(\sum_{l \in h}D_l^{(t)} - \sum_{u \in h}D_u^{(t)}\right) &= 2\sum_{e \in h}v_e + \sum_{l \in h}\sum_{e \in l, e \not\in h}v_e + \sum_{u \in h}\sum_{e \in u, e \not\in h}v_e. \end{align} \]
Diejenigen Kanten, die nicht an das betrachtete Hexagon grenzen, tragen jeweils einmal zu $\sum_{l \in h}D_l^{(t)} - \sum_{u \in h}D_u^{(t)}$ bei:
\[ \begin{align} \frac{\sqrt{3}d}{4}\left(\sum_{l \in h}D_l^{(t)} - \sum_{u \in h}D_u^{(t)}\right) &= 2\sum_{e \in h}v_e + \sum_{t \in h}\sum_{e \in t, e \not\in h}v_e \end{align} \]
Betrachte an dieser Stelle die Größe $\newtilde{v}_1^{(1)}$, welche im Zentrum des betrachteten Hexagons definiert ist. Diese ist eine Überlagerung aus vier Geschwindigkeitskomponenten, zwei davon liegen auf der $\mathbf{i}_1-$Achse, zwei senkrecht zu dieser. Erstere tragen mit dem Gewicht $1/3$ zu $\newtilde{v}_1^{(1)}$ bei, letztere mit dem Gewicht $1/6.$ Dies kann man zyklisch auch auf die anderen Raumrichtungen übertragen. Somit gilt
\[ \begin{align} \frac{\sqrt{3}d}{4}\left(\sum_{l \in h}D_l^{(t)} - \sum_{u \in h}D_u^{(t)}\right) &= 6\left(\newtilde{v}_1^{(1)} + \newtilde{v}_2^{(2)} + \newtilde{v}_3^{3}\right).\tag{27.153}\label{eq:checkerboard_triangular} \end{align} \]
Ist die Bedingung $\newtilde{v}_1^{(1)} + \newtilde{v}_2^{(2)} + \newtilde{v}_3^{3} = 0$ nicht erfüllt, so ist der Mittelwert der Divergenz der drei nach oben gerichteten Dreiecke i. A. ungleich dem Mittelwert der Divergenz der nach unten gerichteten Dreiecke. Dies führt zur Entstehung eines sogenannten Checkerboard-Patterns in der Divergenz. Dieser Lärm pflanzt sich von dort auch in andere Felder fort. Es ist also wichtig, die Bedingung $\newtilde{v}_1^{(1)} + \newtilde{v}_2^{(2)} + \newtilde{v}_3^{3} = 0$ immer einzuhalten.
Auf dem sechseckigen Gitter kann man die Gedankengänge vollständig analog auf die relative Vorticity übertragen, was auf
\[ \begin{align} \frac{\sqrt{3}d}{4}\left(\sum_{u \in h}\zeta_u^{(t)} - \sum_{l \in h}\zeta_l^{(t)}\right) &= 6\left(\newtilde{u}_1^{(1)} + \newtilde{u}_2^{(2)} + \newtilde{u}_3^{3}\right) \end{align} \]
führt. Nichteinhaltung der Bedingung $\newtilde{u}_1^{(1)} + \newtilde{u}_2^{(2)} + \newtilde{u}_3^{3} = 0$ führt also auch auf dem sechseckigen Gitter zu einem Checkerboard-Pattern (hier jedoch in der Vorticity) und ist daher auch hier zu beachten.
Im vorherigen Abschnitt wurde gezeigt, dass ein Checkerboard-Pattern auf dem sechseckigen Gitter in der auf Dreiecken berechneten Vorticity existiert, falls die Bedingung $\newtilde{u}_1^{(1)} + \newtilde{u}_2^{(2)} + \newtilde{u}_3^{3} = 0$ für das Windfeld $\mathbf{v}$ nicht eingehalten wird. Damit diese Bedingung an das Windfeld jedoch für alle Zeitschritte erfüllt ist, muss dies sowohl im Anfangszustand für das Windfeld gelten, als auch für alle Forcings $\mathbf{w}$. Simulationen werden üblicherweise jedoch mit Impulsdiffusion durchgeführt, um eine Aufstauung der kinetischen Energie auf der Gitterskala zu vermeiden. Diese diffusiven Terme bauen Gradienten auf der Gitterskala ab und können daher auch ein bestehendes Checkerboard-Pattern abbauen, solange dieses nicht zu stark ist. Dementsprechend kann auf die Forderung $\newtilde{u}_1^{(1)} + \newtilde{u}_2^{(2)} + \newtilde{u}_3^{3} = 0$ für das Anfangswindfeld $\mathbf{v}$ verzichtet werden. Wichtig ist jedoch folgendes:
Für alle Geschwindigkeitstendenzen $\mathbf{w}$ muss $\newtilde{w}_1^{(1)} + \newtilde{w}_2^{(2)} + \newtilde{w}_3^{3} = 0$ gelten.
Die analogen Tatsachen gelten auf dem dreieckigen Gitter.
Wendet man den Laplace-Operator auf Glg.en (27.116) - (27.118) an, erhält man
\[ \begin{align} \Delta u_1 &= \frac{1}{\sqrt{3}}\left(\delta_3\Delta\newtilde{\psi}^{(2)} - \delta_2\Delta \newtilde{\psi}^{3}\right) + \delta_1\Delta\chi,\tag{27.155}\label{eq:rhombus_curl_deriv_0}\\ \Delta u_2 &= \frac{1}{\sqrt{3}}\left(\delta_1\Delta\newtilde{\psi}^{3} - \delta_3\Delta \newtilde{\psi}^{(1)}\right) + \delta_2\Delta\chi,\tag{27.156}\label{eq:rhombus_curl_deriv_1}\\ \Delta u_3 &= \frac{1}{\sqrt{3}}\left(\delta_2\Delta\newtilde{\psi}^{(1)} - \delta_1\Delta \newtilde{\psi}^{(2)}\right) + \delta_3\Delta\chi.\tag{27.157}\label{eq:rhombus_curl_deriv_2} \end{align} \]
Da diskretisierte partielle Ableitungen und Mittelungsoperatoren vertauschen, folgt hieraus unmittelbar
\[ \begin{align} \newtilde{\Delta u_1}^{(1)} + \newtilde{\Delta u_2}^{(2)} + \newtilde{\Delta u_3}^{(3)} = 0. \end{align} \]
Diskretisiert man die Glg.en (27.62) - (27.64), erhält man
\[ \begin{align} \Delta u_1 &= \frac{1}{\sqrt{3}}\left(\delta_3\zeta - \delta_2\zeta\right) + \delta_1 D,\tag{27.159}\label{eq:rhombus_curl_deriv_3}\\ \Delta u_2 &= \frac{1}{\sqrt{3}}\left(\delta_1\zeta - \delta_3\zeta\right) + \delta_2 D,\tag{27.160}\label{eq:rhombus_curl_deriv_4}\\ \Delta u_3 &= \frac{1}{\sqrt{3}}\left(\delta_2\zeta - \delta_1\zeta\right) + \delta_3 D.\tag{27.161}\label{eq:rhombus_curl_deriv_5} \end{align} \]
Durch Vergleich stellt man fest, dass es im $\mathbf{i}-$System drei Arten von Vorticities gibt:
\[ \begin{align} \zeta_1 \coloneqq \Delta\newtilde{\psi}^{(1)}, & {} & \zeta_2 \coloneqq \Delta\newtilde{\psi}^{(2)}, & {} & \zeta_3 \coloneqq \Delta\newtilde{\psi}^{3} \end{align} \]
Es gilt somit
\[ \begin{align} \Delta u_1 &= \frac{1}{\sqrt{3}}\left(\delta_3\zeta_2 - \delta_2\zeta_3\right) + \delta_1 D,\tag{27.163}\label{eq:rhombus_curl_deriv_6}\\ \Delta u_2 &= \frac{1}{\sqrt{3}}\left(\delta_1\zeta_3 - \delta_3\zeta_1\right) + \delta_2 D,\tag{27.164}\label{eq:rhombus_curl_deriv_7}\\ \Delta u_3 &= \frac{1}{\sqrt{3}}\left(\delta_2\zeta_1 - \delta_1\zeta_2\right) + \delta_3 D.\tag{27.165}\label{eq:rhombus_curl_deriv_8} \end{align} \]
Für die Divergenz gilt nach Glg. (27.33)
\[ \begin{align} D &= \frac{2}{3}\left(\delta_1u_1 + \delta_2u_2 + \delta_3u_3\right). \end{align} \]
Setzt man dies in die Glg.en (27.163) - (27.165) ein, erhält man
\[ \begin{align} \Delta u_1 &= \frac{1}{\sqrt{3}}\left(\delta_3\zeta_2 - \delta_2\zeta_3\right) + \frac{2}{3}\left(\delta_1\delta_1u_1 + \delta_1\delta_2u_2 + \delta_1\delta_3u_3\right),\tag{27.167}\label{eq:hex_laplace_vec_0}\\ \Delta u_2 &= \frac{1}{\sqrt{3}}\left(\delta_1\zeta_3 - \delta_3\zeta_1\right) + \frac{2}{3}\left(\delta_2\delta_1u_1 + \delta_2\delta_2u_2 + \delta_2\delta_3u_3\right),\tag{27.168}\label{eq:hex_laplace_vec_1}\\ \Delta u_3 &= \frac{1}{\sqrt{3}}\left(\delta_2\zeta_1 - \delta_1\zeta_2\right) + \frac{2}{3}\left(\delta_3\delta_1u_1 + \delta_3\delta_2u_2 + \delta_3\delta_3u_3\right).\tag{27.169}\label{eq:hex_laplace_vec_2} \end{align} \]
Für die Vorticities $\zeta_i$ macht man, motiviert durch die Glg.en (27.49) - (27.51), die Ansätze
\[ \begin{align} \zeta_1 &= \frac{2}{\sqrt{3}}\left(\delta_2u_3 - \delta_3u_2\right),\tag{27.170}\label{eq:hex_curl_rhombus_0}\\ \zeta_2 &= \frac{2}{\sqrt{3}}\left(\delta_3u_1 - \delta_1u_3\right),\tag{27.171}\label{eq:hex_curl_rhombus_1}\\ \zeta_3 &= \frac{2}{\sqrt{3}}\left(\delta_1u_2 - \delta_2u_1\right).\tag{27.172}\label{eq:hex_curl_rhombus_2} \end{align} \]
Setzt man Glg. (27.170) in Glg. (27.167) ein, erhält man
\[ \begin{align} \Delta u_1 &= \frac{2}{3}\left(\delta_3\delta_3u_1 - \delta_3\delta_1u_3 - \delta_2\delta_1u_2 + \delta_2\delta_2u_1\right) + \frac{2}{3}\left(\delta_1\delta_1u_1 + \delta_1\delta_2u_2 + \delta_1\delta_3u_3\right) = \frac{2}{3}\left(\delta_1\delta_1u_1 + \delta_2\delta_2u_1 + \delta_3\delta_3u_1\right). \end{align} \]
Durch zyklische Vertauschungen erhält man die analoge Aussage auch für die anderen Komponenten. Somit sind die Ansätze der Glg.en (27.170) - (27.172) gerechtfertigt. Zwei benachbarte Dreiecke bilden zusammen einen Rhombus (ein Parallelogramm), ein solcher hat nach Glg. (27.144) die Fläche
\[ \begin{align} A_r = 2A_v = 2\frac{d^2\sqrt{3}}{4} = \frac{d^2\sqrt{3}}{2}. \end{align} \]
Somit entspricht $\zeta_i$ der mittels des Stokes'schen Satzes über einen Rhombus, dessen kurze Symmetrieachse parallel zur i-Achse liegt, berechneten Vorticity. Auf dem hexagonalen Gitter bilden die Rhomben also die natürlichen Mengen, auf denen der Stokes'sche Satz ausgewertet werden muss. Dies ist eine Abweichung von der eigentlich intuitiven Annahme, dass die Vorticity über das duale Gitter (das Dreicksgitter) berechnet werden muss.
Auf dem dreieckigen Gitter treten wie üblich Divergenz und Vorticity in vertauschten Rollen auf. Die Vorticity hat dort den eindeutigen Wert
\[ \begin{align} \zeta = \frac{2}{3}\left(\delta_1v_1 + \delta_2v_2 + \delta_3v_3\right) \end{align} \]
Dort treten jedoch drei Arten von Divergenzen auf, die jeweils mittels des Gauß'schen Satzes über die jeweiligen Rhomben berechnet werden:
\[ \begin{align} D_1 &= \frac{2}{\sqrt{3}}\left(\delta_3v_2 - \delta_2v_3\right)\\ D_2 &= \frac{2}{\sqrt{3}}\left(\delta_1v_3 - \delta_3v_1\right)\\ D_3 &= \frac{2}{\sqrt{3}}\left(\delta_2v_1 - \delta_1v_2\right) \end{align} \]
Dies wird auch durch die Diskretisierung des vektoriellen Laplace-Operators gemäß Glg. (B.54) unterstrichen, welche
\[ \begin{align} \Delta\mathbf{v} = \nabla\left(\nabla\cdot\mathbf{v}\right) - \nabla\times\left(\nabla\times\mathbf{v}\right) \end{align} \]
lautet. Wendet man dies auf die Vektorkomponente $u_1$ an, erhält man
\[ \begin{align} \Delta u_1 = \delta_1D_h - \delta_{i, \perp}\zeta_t. \end{align} \]
Hierbei ist $D_h$ die Divergenz auf Hexagonen und $\zeta_t$ die Vorticity auf Dreicken. $\delta_{i, \perp}$ bezeichnet den zentralen Differenzenquotienten in $\mathbf{j}_1-$Richtung. Für diesen gilt
\[ \begin{align} -\delta_{i, \perp}\zeta_t &= \frac{\sqrt{3}}{d}\left(\zeta_{t, l} - \zeta_{t, u}\right). \end{align} \]
Aus Glg. (27.163) folgt
\[ \begin{align} \Delta u_1 &= \delta_1D_h + \frac{\sqrt{3}}{d}\left(\zeta_{t, l} - \zeta_{t, u}\right) = \delta_1 D_h + \frac{1}{\sqrt{3}}\left(\delta_3\zeta_2 - \delta_2\zeta_3\right)\nonumber\\ & \Leftrightarrow \sqrt{3}\frac{\zeta_{t, l} - \zeta_{t, u}}{d} = \frac{\sqrt{3}}{3}\left(\delta_3\zeta_2 - \delta_2\zeta_3\right) \Leftrightarrow \sqrt{3}\frac{\zeta_{t, l} - \zeta_{t, u}}{d} = \frac{\sqrt{3}}{3}\frac{\zeta_{2, l} - \zeta_{2, u} - \zeta_{3, u} + \zeta_{3, l}}{d}\nonumber\\ & \Leftrightarrow \sqrt{3}\frac{\zeta_{t, l} - \zeta_{t, u}}{d} = \frac{\sqrt{3}}{3}\frac{\zeta_{2, l} + \zeta_{3, l} - \left(\zeta_{2, u} - \zeta_{3, u}\right)}{d}\nonumber\\ & \Leftrightarrow \sqrt{3}\frac{\zeta_{t, l} - \zeta_{t, u}}{d} = \frac{\sqrt{3}}{3}\frac{\zeta_{1, l} + \zeta_{2, l} + \zeta_{3, l} - \left(\zeta_{1, u} + \zeta_{2, u} - \zeta_{3, u}\right)}{d}. \end{align} \]
Die letzte Äquivalenzumformung folgt aus der Tatsache $\zeta_{1, l} = \zeta_{1, u}$. Dies ist mit den Definitionen
\[ \begin{align} \zeta_{t, u} &\coloneqq \frac{\zeta_{1, u} + \zeta_{2, u} + \zeta_{3, u}}{3},\\ \zeta_{t, l} &\coloneqq \frac{\zeta_{1, l} + \zeta_{2, l} + \zeta_{3, l}}{3} \end{align} \]
erfüllt. Die Rotation auf einem Dreieck ist also als Mittelwert der drei Rhomben zu berechnen, die sich mit diesem Dreieck überlappen.
Auf dem hexagonalen C-Gitter erfüllen Gradientenfelder $\mathbf{u} = \nabla\psi$ die Bedingung
\[ \begin{align} \newtilde{u}_1^{(1)} + \newtilde{u}_2^{(2)} + \newtilde{u}_3^{3} = 0. \end{align} \]
Auf dem dreieckigen C-Gitter hat ein Gradientenfeld $\mathbf{v} = \nabla\psi$ natürlicherweise Komponenten parallel zu den Basiselementen $\mathbf{j}_i$. Dies erfüllt
\[ \begin{align} \newtilde{v}_1^{(1)} + \newtilde{v}_2^{(2)} + \newtilde{v}_3^{3} = 0 \end{align} \]
nicht, wie man sich leicht klarmachen kann. Daher produzieren Gradientenfelder auf dem dreieckigen C-Gitter nach Glg. (27.153) ein Checkerboard-Pattern in der Divergenz, was unphysikalische numerische Stabilisatoren, wie zum Beispiel eine Divergenzdämpfung, erfordern würde. Daher wurde dieses Gitter in Abschn. 26.9.1 zu Recht ausgeschlossen. Diese Erkenntnis wurde erstmals von Gassmann in [24] vorgebracht.