| (62) |
An Euler implicit discretization in time of the governing equations is combined with a Newton-type linearization of all fluxes to produce the following linearized system in delta form:
|
| (62) |
where 
. The center and off-diagonal coefficient matrices, 
and 
, and the residual vector 
are given by:
| (63) |
| (64) |
| (65) |