next up previous contents
Next: 3.2.3 Stability Up: 3.2 Solving the Level Previous: 3.2.1 Upwind Scheme

3.2.2 Lax-Friedrichs Scheme

To solve the LS equation for non-convex Hamiltonians the Lax-Friedrichs scheme can be applied [87,110]. Some artificial dissipation is introduced to obtain stability. The approximation to the Hamiltonian for this scheme is written as

$\displaystyle {\hat{H}}({\vec{p}},{\Phi},{V}({\vec{p}})) := {V}({\vec{p}}) \nab...
...rac{\partial^{+}_{i}{\Phi}({{\vec{p}}})-\partial^{-}_{i}{\Phi}({{\vec{p}}})}{2}$ (3.12)

using the central gradient approximation

$\displaystyle \nabla^0{\Phi}({{\vec{p}}}):= \sqrt{\sum_{i=1}^{D}\left(\frac{\pa...
...l^{-}_{i}{\Phi}({{\vec{p}}})+\partial^{+}_{i}{\Phi}({{\vec{p}}})}{2}\right)^2}.$ (3.13)

$ \partial^{\pm}_{i}{\Phi}$ are again the one-sided approximations defined in (3.10). Again, depending on whether the second order term in (3.10) is included or not, the scheme is second or first order accurate in space, respectively [110].

The constants $ {\alpha}_i$ are dissipation coefficients which should fulfill

$\displaystyle {\alpha}_i \geq \max_{{\vec{x}}} \left\vert \frac{\partial{H}}{\partial q_i}\right\vert$   with$\displaystyle \ q_i:=\frac{\partial {\Phi}}{\partial{x}_i}$ (3.14)

in order to guarantee stability [85,110].

The choice of these dissipation coefficients is the main problem for the applicability of this scheme. Larger values lead to unnecessary smoothing of the surface, while smaller values cause numerical instabilities. The analytical evaluation of the right-hand side in (3.14) to find the optimal values for the constants $ {\alpha}_i$ is usually not possible for general surface velocities. However, often the optimal values cannot be straightforwardly evaluated, because the velocity field may depend on the particle transport which further depends on the surface geometry. Hence, the likely best way to find the optimal values is to test different values and analyze the results [98].


next up previous contents
Next: 3.2.3 Stability Up: 3.2 Solving the Level Previous: 3.2.1 Upwind Scheme

Otmar Ertl: Numerical Methods for Topography Simulation