next up previous
Next: The (3+1)D Interpolated Rectilinear Up: Finite Difference Schemes for Previous: The Cubic Rectilinear Scheme

The Octahedral Scheme

The grid for an octahedral scheme is constructed from two superimposed rectilinear grids; if the points of the first grid are located at cube corners, then the points of the second will occur at the centers of the cubes defined by the first. The relevant difference scheme on an octahedral grid can be written as

\begin{displaymath}\begin{split}U_{i,j,k}(n+1)+U_{i,j,k}(n-1) &= \frac{3}{4}\lam...
...g)\\ &\quad+\left(2-8\lambda^{2}\right)U_{i,j,k}(n) \end{split}\end{displaymath} (A.30)

for $ i$, $ j$ and $ k$ which are either all even or all odd integers. Now, we have taken the spacing between nearest neighbors to be $ \Delta$, so the indices $ i$, $ j$ and $ k$ refer to a point with coordinates $ x=i\Delta/\sqrt{3}$, $ y=j\Delta/\sqrt{3}$ and $ z=k\Delta/\sqrt{3}$. The amplification polynomial equation is again of the form (A.5), with

$\displaystyle B_{\mbox{{\scriptsize\boldmath$\beta$}}} = -2\left(1+3\lambda^{2}...


$\displaystyle F_{\mbox{{\scriptsize\boldmath$\beta$}}} = 3\left(\cos(\frac{\bet...

and it is again easy to determine that

$\displaystyle \max_{\mbox{{\scriptsize\boldmath$\beta$}}}F_{\mbox{{\scriptsize\...
...ox{{\scriptsize\boldmath$\beta$}}}F_{\mbox{{\scriptsize\boldmath$\beta$}}} = -6$    

which are the same as the bounds in the cubic rectilinear case. We again have that

$\displaystyle \lambda\leq\frac{1}{\!\!\sqrt{3}}$   (for Von Neumann stability)    

Thus the stability bound coincides with the passivity bound for the mesh implementation. For $ \lambda = 1/\sqrt{3}$, instabilities may appear at any spatial frequency triplets $ \beta$$ = [\beta_{x}, \beta_{y}, \beta_{z}]^{T}$ where each component is either 0 or $ \pm\sqrt{3}\pi/\Delta$.

The computational and add densities are given by

$\displaystyle \rho_{oct}$ $\displaystyle = \frac{3\sqrt{3}v_{0}}{4\Delta^{4}}$   $\displaystyle \sigma_{oct} = \frac{27\sqrt{3}v_{0}}{4\Delta^{4}}$   for$\displaystyle \hspace{0.3in} v_{0}>\sqrt{3}\gamma$    
$\displaystyle \rho^{s}_{oct}$ $\displaystyle = \frac{9\gamma}{8\Delta^{4}}$   $\displaystyle \sigma^{s}_{oct} = \frac{9\gamma}{\Delta^{4}}$   for$\displaystyle \hspace{0.3in} v_{0}=\sqrt{3}\gamma$    

At the stability limit, the scheme can be divided into two mutually exclusive subschemes; plots of numerical dispersion are shown in Figure A.9(b) and (c). It is interesting to note that there is no dispersion error along the six axial directions; this should be compared with the cubic rectilinear scheme, for which wave propagation is dispersionless along the diagonal directions (there are eight such directions).

Figure A.9: The octahedral scheme (A.29)-- (a) numerical grid and connections, where grey/white coloring of points indicates a division into mutually exclusive subschemes at the stability bound; (b) $ v_{\mbox{{\scriptsize\boldmath$\beta$}}, phase}/\gamma$ for the scheme at the stability bound $ \lambda = 1/\sqrt{3}$, for a spherical surface with $ \Vert$$ \beta$$ \Vert _{2} = \pi/(2\Delta)$--the shading is normalized over the surface so that white corresponds to no dispersion error, and black to the maximum error over the surface (which is 5 per cent in this case). (c) Contour plots of the $ v_{\mbox{{\scriptsize\boldmath$\beta$}}, phase}/\gamma$ for various cross-sections of the space of spatial frequencies $ \beta$; contours indicate successive deviations of 2 per cent from the ideal value of 1 which is obtained at spatial DC.
...{4}{\tiny {$\beta_{x}$}}
\end{picture} \end{center} \vspace{0.3in}

next up previous
Next: The (3+1)D Interpolated Rectilinear Up: Finite Difference Schemes for Previous: The Cubic Rectilinear Scheme
Stefan Bilbao 2002-01-22