This paper considers questions regarding conditioning of, and numerical methods, for certain differential algebraic equations subject to initial and boundary conditions. The approach taken is that of separating “differential” and “algebraic” solution components, at least theoretically. This yields conditioning results for differential algebraic boundary value problems in terms of “pure” differential problems, for which existing theory is well developed. The process is carried out for problems with (global) index 1 or 2. For semi-explicit boundary value problems of index 1 (where solution components are separated) a convergence theorem is given for a special class of collocation methods. For general index 1 problems advantages and disadvantages of certain symmetric difference schemes are discussed. For initial value problems with index 2 the use of BDF schemes is examined, with a summary of conditions for their successful and stable utilization. Finally, the present considerations and analysis are applied to two problems involving differential algebraic equations arising in semiconductor device simulation.
Read full abstract