appearing in Eqn. 58 and Eqn. 63 is computed by consideration of the Courant number (
) and Von Neumann stability conditions:
The local time-step appearing in Eqn. 58 and Eqn. 63 is computed by consideration of the Courant number () and Von Neumann stability conditions:
| (66) |
where 
is the maximum eigenvalue of the system, 
is the Courant number, 
is the cell volume, 
is the Von Neumann number (
) and 
is a characteristic cell length scale.