GLOBAL DISSIPATIVITY FOR A-STABLE METHODS

摘要

This paper concerns the discretization of the initial value problem u(t) = f(u) under the three structural conditions: (i) f:C-N --> C-N, Re(f(u), u) less than or equal to a - bu(2), a greater than or equal to 0, b > 0 for all u epsilon CN; (ii) f:C-N --> C-N, Re(f(u), u) < 0 for u epsilon C-NB(0, R) and R > 0; (iii) f:W --> H, Re(f(w)(H) less than or equal to a - bw(2)(H), a greater than or equal to 0, b > 0 for all w epsilon W for complex Hilbert spaces W subset of or equal to H. Dahlquist's G-stability theory is used to show that linear multistep and one-leg methods yield dissipative discretizations for all f satisfying (i) if and only if the method (rho,sigma) is A-stable. Extensions of G-stability theory are made to find necessary and sufficient conditions on (rho,sigma) for similar properties to hold in cases (ii) and (iii). In every case, conditions are found for the strict contractivity of solutions for large initial data, and bounds for the asymptotic rate of decay are calculated in cases (i) and (iii). [References: 17]