On the stability of a finite difference scheme for the nonlinear silicon problem
摘要
This paper introduces a finite difference numeric scheme designed for the computational solution of the nonlinear heat conduction problem in a homogeneous silicon rod. It is shown that the numerical method is consistent with the conservation law, is conditionally stable and convergent. Focusing on the implicit schemes, an efficient Newton multigrid method with Gauss-Seidel red-black smoother is developed. Computational experiments confirm the stability theory as well as the robustness of the multigrid method.