TO Sebastien

Zhu zhiwen

Dear Sebastien,

continue to the last question, On SHEN's paper(J.of Comp.Phys.1998) the scheme which using the ( 3*u^(n+1)-4*n^n+u^(n-1))/2dt for the time derivative ,he said that the linear parabolic operator(Laplace operator)is approximated by a second-order backward scheme which appears to be more stable than the crank-nicholson scheme ,I donot know how he deal with the viscocity term,I didnot mean that the viscosity is varable,in generally,the central difference is used for the viscosity term.but on his paper,it seems not,

Thanks