3. 2. 4 Laplace Equation

One method typically used in simulation is the dual graph method. As an introductory example, the numerical solution of the Laplace equation is shown. In this case the functional $ \phi$ is the gradient and $ g = 0$ .

For linear shape functions, the function $ f(\mathbf{e})$ along one edge can be written as

$\displaystyle f(\mathbf{e}) = \sum_{EV}^{\underline{e}} q f(\bullet, \underline{e}) \; ,$ (3.34)

The function $ f(\mathbf{e}, \mathbf{v})$ denotes a shape function defined on the edge-skeleton of the cell complex. These functions have local support on one edge $ \mathbf{e}$ and yield unity in the vertex $ \mathbf{v}$ . Consequently, the gradient is constant along the edge and, therefore, the evaluation point of the function as described in () does not play any role.

$\displaystyle \Phi(\mathbf{e}, \mathbf{v}) := \frac{1}{l_ {\underline{2}}} \sum_{EV}^{\underline{e}} \mathcal{O}(\underline{e}, \bullet) q \;$ (3.35)

The complete expression for the discretized Laplace equation yields

$\displaystyle R(\mathbf{v}) := m(\mathbf{v}) \Delta(f(\mathbf{v})) = 0 \; ,$    
$\displaystyle R = [m\mathcal{L}f] = \sum_{VE}^{\underline{v}} \mathcal{O}(\bull...
...frac{A}{l} \sum_{EV}^{\underline{e}} \mathcal{O}(\underline{e}, \bullet) q \; .$ (3.36)

Figure 3.8: Topological neighborhood required for finite volume schemes.
\includegraphics[width=10cm]{DRAWINGS/fvm_topo.eps}

Michael 2008-01-16