Show simple item record

dc.contributor.authorWyatt, Sarah Aliceen_US
dc.date.accessioned2014-03-14T20:37:46Z
dc.date.available2014-03-14T20:37:46Z
dc.date.issued2009-05-04en_US
dc.identifier.otheretd-05202009-144059en_US
dc.identifier.urihttp://hdl.handle.net/10919/33042
dc.description.abstractDynamical systems are mathematical models characterized by a set of differential or difference equations. Due to the increasing demand for more accuracy, the number of equations involved may reach the order of thousands and even millions. With so many equations, it often becomes computationally cumbersome to work with these large-scale dynamical systems. Model reduction aims to replace the original system with a reduced system of significantly smaller dimension which will still describe the important dynamics of the large-scale model. Interpolation is one method used to obtain the reduced order model. This requires that the reduced order model interpolates the full order model at selected interpolation points. Reduced order models are obtained through the Krylov reduction process, which involves solving a sequence of linear systems. The Iterative Rational Krylov Algorithm (IRKA) iterates this Krylov reduction process to obtain an optimal $\mathcal{H}_2$ reduced model. Especially in the large-scale setting, these linear systems often require employing inexact solves. The aim of this thesis is to investigate the impact of inexact solves on interpolatory model reduction.

We considered preconditioning the linear systems, varying the stopping tolerances, employing GMRES and BiCG as the inexact solvers, and using different initial shift selections. For just one step of Krylov reduction, we verified theoretical properties of the interpolation error. Also, we found a linear improvement in the subspace angles between the inexact and exact subspaces provided that a good shift selection was used. For a poor shift selection, these angles often remained of the same order regardless of how accurately the linear systems were solved. These patterns were reflected in $\mathcal{H}_2$ and $\mathcal{H}_{\infty}$ errors between the inexact and exact subspaces, since these errors improved linearly with a good shift selection and were typically of the same order with a poor shift. We found that the shift selection also influenced the overall model reduction error between the full model and inexact model as these error norms were often several orders larger when a poor shift selection was used. For a given shift selection, the overall model reduction error typically remained of the same order for tolerances smaller than $1 \times 10^{-3}$, which suggests that larger tolerances for the inexact solver may be used without necessarily augmenting the model reduction error. With preconditioned linear systems as well as BiCG, we found smaller errors between the inexact and exact models while the order of the overall model reduction error remained the same. With IRKA, we observed similar patterns as with just one step of Krylov reduction. However, we also found additional benefits associated with using an initial guess in the inexact solve and by varying the tolerance of the inexact solve.

en_US
dc.publisherVirginia Techen_US
dc.relation.haspartetd_updated.pdfen_US
dc.rightsI hereby certify that, if appropriate, I have obtained and attached hereto a written permission statement from the owner(s) of each third party copyrighted matter to be included in my thesis, dissertation, or project report, allowing distribution as specified below. I certify that the version I submitted is the same as that approved by my advisory committee. I hereby grant to Virginia Tech or its agents the non-exclusive license to archive and make accessible, under the conditions specified below, my thesis, dissertation, or project report in whole or in part in all forms of media, now or hereafter known. I retain all other ownership rights to the copyright of the thesis, dissertation or project report. I also retain the right to use in future works (such as articles or books) all or part of this thesis, dissertation, or project report.en_US
dc.subjectH₂ Approximationen_US
dc.subjectRational Kryloven_US
dc.subjectinterpolationen_US
dc.subjectmodel reductionen_US
dc.titleInexact Solves in Interpolatory Model Reductionen_US
dc.typeThesisen_US
dc.contributor.departmentMathematicsen_US
dc.description.degreeMaster of Scienceen_US
thesis.degree.nameMaster of Scienceen_US
thesis.degree.levelmastersen_US
thesis.degree.grantorVirginia Polytechnic Institute and State Universityen_US
thesis.degree.disciplineMathematicsen_US
dc.contributor.committeechairGugercin, Serkanen_US
dc.contributor.committeememberBeattie, Christopher A.en_US
dc.contributor.committeememberde Sturler, Ericen_US
dc.identifier.sourceurlhttp://scholar.lib.vt.edu/theses/available/etd-05202009-144059/en_US
dc.date.sdate2009-05-20en_US
dc.date.rdate2009-05-27
dc.date.adate2009-05-27en_US


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record