Continuous orthonormalization for boundary value problems

From MaRDI portal





Continuous orthonormalization describes an initial value method for linear two-point boundary value problems \(y'=A(t)y+r(t)\), \(0<t<1\), \(By(0)=\alpha_ 0\), \(g(y(1))=0\), which provides an orthogonal basis for the solution space at all points of the interval. The equations of continuous orthonormalization are derived with elementary projection arguments to provide geometric insight and motivate some modification of the earlier algorithm by \textit{A. Davey} [J. Comput. Phys. 51, 343-356 (1983; Zbl 0516.65065)]. The method is then applied to some oscillatory and stiff boundary value problems to demonstrate that it is simple to use, problem independent, and as adaptive as the initial value code which is used to integrate the equations of continuous orthonormalization. Three numerical examples, including the Orr-Sommerfeld eigenvalue problem, are given to illustrate the method.





Describes a project that uses

Uses Software






This page was built for publication: Continuous orthonormalization for boundary value problems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1082049)