Issues in Interpolatory Model Reduction: Inexact Solves, Second-order Systems and DAEs

dc.contributor.authorWyatt, Sarah Aliceen
dc.contributor.committeechairGugercin, Serkanen
dc.contributor.committeememberde Sturler, Ericen
dc.contributor.committeememberBorggaard, Jeffrey T.en
dc.contributor.committeememberBeattie, Christopher A.en
dc.contributor.departmentMathematicsen
dc.date.accessioned2014-03-14T20:11:54Zen
dc.date.adate2012-05-25en
dc.date.available2014-03-14T20:11:54Zen
dc.date.issued2012-05-01en
dc.date.rdate2012-05-25en
dc.date.sdate2012-05-11en
dc.description.abstractDynamical systems are mathematical models characterized by a set of differential or difference equations. Model reduction aims to replace the original system with a reduced system of significantly smaller dimension that still describes the important dynamics of the large-scale model. Interpolatory model reduction methods define a reduced model that interpolates the full model at selected interpolation points. The reduced model may be obtained through a Krylov reduction process or by using the Iterative Rational Krylov Algorithm (IRKA), which iterates this Krylov reduction process to obtain an optimal ℋ₂ reduced model. This dissertation studies interpolatory model reduction for first-order descriptor systems, second-order systems, and DAEs. The main computational cost of interpolatory model reduction is the associated linear systems. Especially in the large-scale setting, inexact solves become desirable if not necessary. With the introduction of inexact solutions, however, exact interpolation no longer holds. While the effect of this loss of interpolation has previously been studied, we extend the discussion to the preconditioned case. Then we utilize IRKA's convergence behavior to develop preconditioner updates. We also consider the interpolatory framework for DAEs and second-order systems. While interpolation results still hold, the singularity associated with the DAE often results in unbounded model reduction errors. Therefore, we present a theorem that guarantees interpolation and a bounded model reduction error. Since this theorem relies on expensive projectors, we demonstrate how interpolation can be achieved without explicitly computing the projectors for index-1 and Hessenberg index-2 DAEs. Finally, we study reduction techniques for second-order systems. Many of the existing methods for second-order systems rely on the model's associated first-order system, which results in computations of a 2𝑛 system. As a result, we present an IRKA framework for the reduction of second-order systems that does not involve the associated 2𝑛 system. The resulting algorithm is shown to be effective for several dynamical systems.en
dc.description.degreePh. D.en
dc.identifier.otheretd-05112012-103400en
dc.identifier.sourceurlhttp://scholar.lib.vt.edu/theses/available/etd-05112012-103400/en
dc.identifier.urihttp://hdl.handle.net/10919/27668en
dc.publisherVirginia Techen
dc.relation.haspartWyatt_SA_D_2012.pdfen
dc.rightsIn Copyrighten
dc.rights.urihttp://rightsstatements.org/vocab/InC/1.0/en
dc.subjectSecond-order Systemsen
dc.subjectInexact Solvesen
dc.subjectKrylov reductionen
dc.subjectDAEsen
dc.titleIssues in Interpolatory Model Reduction: Inexact Solves, Second-order Systems and DAEsen
dc.typeDissertationen
thesis.degree.disciplineMathematicsen
thesis.degree.grantorVirginia Polytechnic Institute and State Universityen
thesis.degree.leveldoctoralen
thesis.degree.namePh. D.en

Files

Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
Wyatt_SA_D_2012.pdf
Size:
801.59 KB
Format:
Adobe Portable Document Format