Reduced-Order Modeling of Complex Engineering and Geophysical Flows: Analysis and Computations

Simple item page
dc.contributor.authorWang, Zhuen
dc.contributor.committeechairIliescu, Traianen
dc.contributor.committeememberBurns, John A.en
dc.contributor.committeememberBorggaard, Jeffrey T.en
dc.contributor.committeememberLin, Taoen
dc.contributor.committeememberZietsman, Lizetteen
dc.contributor.departmentMathematicsen
dc.date.accessioned2014-03-14T20:11:27Zen
dc.date.adate2012-05-14en
dc.date.available2014-03-14T20:11:27Zen
dc.date.issued2012-04-17en
dc.date.rdate2012-05-14en
dc.date.sdate2012-05-02en
dc.description.abstractReduced-order models are frequently used in the simulation of complex flows to overcome the high computational cost of direct numerical simulations, especially for three-dimensional nonlinear problems. Proper orthogonal decomposition, as one of the most commonly used tools to generate reduced-order models, has been utilized in many engineering and scientific applications. Its original promise of computationally efficient, yet accurate approximation of coherent structures in high Reynolds number turbulent flows, however, still remains to be fulfilled. To balance the low computational cost required by reduced-order modeling and the complexity of the targeted flows, appropriate closure modeling strategies need to be employed. In this dissertation, we put forth two new closure models for the proper orthogonal decomposition reduced-order modeling of structurally dominated turbulent flows: the dynamic subgrid-scale model and the variational multiscale model. These models, which are considered state-of-the-art in large eddy simulation, are carefully derived and numerically investigated. Since modern closure models for turbulent flows generally have non-polynomial nonlinearities, their efficient numerical discretization within a proper orthogonal decomposition framework is challenging. This dissertation proposes a two-level method for an efficient and accurate numerical discretization of general nonlinear proper orthogonal decomposition closure models. This method computes the nonlinear terms of the reduced-order model on a coarse mesh. Compared with a brute force computational approach in which the nonlinear terms are evaluated on the fine mesh at each time step, the two-level method attains the same level of accuracy while dramatically reducing the computational cost. We numerically illustrate these improvements in the two-level method by using it in three settings: the one-dimensional Burgers equation with a small diffusion parameter, a two-dimensional flow past a cylinder at Reynolds number Re = 200, and a three-dimensional flow past a cylinder at Reynolds number Re = 1000. With the help of the two-level algorithm, the new nonlinear proper orthogonal decomposition closure models (i.e., the dynamic subgrid-scale model and the variational multiscale model), together with the mixing length and the Smagorinsky closure models, are tested in the numerical simulation of a three-dimensional turbulent flow past a cylinder at Re = 1000. Five criteria are used to judge the performance of the proper orthogonal decomposition reduced-order models: the kinetic energy spectrum, the mean velocity, the Reynolds stresses, the root mean square values of the velocity fluctuations, and the time evolution of the proper orthogonal decomposition basis coefficients. All the numerical results are benchmarked against a direct numerical simulation. Based on these numerical results, we conclude that the dynamic subgrid-scale and the variational multiscale models are the most accurate. We present a rigorous numerical analysis for the discretization of the new models. As a first step, we derive an error estimate for the time discretization of the Smagorinsky proper orthogonal decomposition reduced-order model for the Burgers equation with a small diffusion parameter. The theoretical analysis is numerically verified by two tests on problems displaying shock-like phenomena. We then present a thorough numerical analysis for the finite element discretization of the variational multiscale proper orthogonal decomposition reduced-order model for convection-dominated convection-diffusion-reaction equations. Numerical tests show the increased numerical accuracy over the standard reduced-order model and illustrate the theoretical convergence rates. We also discuss the use of the new reduced-order models in realistic applications such as airflow simulation in energy efficient building design and control problems as well as numerical simulation of large-scale ocean motions in climate modeling. Several research directions that we plan to pursue in the future are outlined.en
dc.description.degreePh. D.en
dc.identifier.otheretd-05022012-012117en
dc.identifier.sourceurlhttp://scholar.lib.vt.edu/theses/available/etd-05022012-012117/en
dc.identifier.urihttp://hdl.handle.net/10919/27504en
dc.publisherVirginia Techen
dc.relation.haspartWang_Z_D_2012.pdfen
dc.rightsIn Copyrighten
dc.rights.urihttp://rightsstatements.org/vocab/InC/1.0/en
dc.subjectvariational multiscaleen
dc.subjectdynamic subgrid-scale modelen
dc.subjecttwo-level algorithmen
dc.subjectapproximate deconvolutionen
dc.subjectfinite elementsen
dc.subjectnumerical analysisen
dc.subjectProper orthogonal decompositionen
dc.subjectreduced-order modelingen
dc.subjectlarge eddy simulationen
dc.subjecteddy viscosityen
dc.subjectTurbulenceen
dc.titleReduced-Order Modeling of Complex Engineering and Geophysical Flows: Analysis and Computationsen
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
Thumbnail Image
Name:
Wang_Z_D_2012.pdf
Size:
8.42 MB
Format:
Adobe Portable Document Format
Download

Collections

Doctoral Dissertations