Discrete-time orthogonal spline collocation methods for Schrodinger-type problems

Time-dependent Schrodinger and vibration problems are of much physical and engineering significance. In this dissertation, we study the fully discrete high-order approximate solution of initial-boundary value problems for linear and nonlinear Schrodinger equations and vibration equation in two space variables on a unit square. We employ orthogonal spline collocation (OSC) with C¹ piecewise polynomials of arbitrary degree r/ge3 in each space variable for the spatial discretization, then use Crank-Nicolson and alternating direction implicit (ADI) techniques for the time discretization. For linear Schrodinger problems, we formulate Crank- Nicolson and ADI OSC schemes, and verify the existence and uniqueness of the approximations. We prove that the schemes are stable and are of second-order accuracy in time and (r+1-m)-order in space in the H^m-norm, m=0,1, which are optimal order estimate. We consider vibrating plate problems with three kinds of boundary conditions, in which the plate is (a) clamped on all sides; (b) hinged on all sides; and (c) hinged on two vertical sides and clamped on two horizontal sides. We reformulate the vibration equation as a system of Schrodinger equations in two new variables, and consider Crank-Nicolson schemes for this system. Stability and optimal order H^m-error estimates, m=1,2, are derived. ADI schemes are formulated for problems with boundary conditions (b) and (c), and stability results and optimal order H^m-error estimates, m=1,2, are derived for a typical case. Also an optimal order L2- error estimate is obtained for the ADI approximation of problem with the boundary condition (b). For nonlinear Schrodinger problems, we present extrapolated Crank-Nicolson and ADI OSC schemes, and examine existence and uniqueness of the approximations. The schemes involve the solution of linear algebraic systems at each time step. We derive the stability and an optimal order L²-error estimate for the Crank- Nicolson scheme. Matrix forms of the ADI schemes are given, and for their implementation, algorithms are described which cost O(N²) operations per time step, where N is the number of partition intervals in each space direction. For the vibration problem with boundary condition (a), for which no ADI scheme has been formulated, we employ the capacitance matrix method to determine the Crank- Nicolson approximation for a special case, with cost O(N²log₂N) per time step. Results of numerical experiments confirm the theoretical analyses.