In this thesis we conduct a numerical study of the 1D viscous Burgers' equation and several Reduced Order Models (ROMs) over a range of parameter values. This study is motivated by the need for robust reduced order models that can be used both for design and control. Thus the model should first, allow for selection of optimal parameter values in a trade space and second, identify impacts from changes of parameter values that occur during development, production and sustainment of the designs. To facilitate this study we apply a Finite Element Method (FEM) and where applicable, the Group Finite Element Method (GFE) due its demonstrated stability and reduced complexity over the standard FEM. We also utilize Proper Orthogonal Decomposition (POD) as a model reduction technique and modifications of POD that include Global POD, and the sensitivity based modifications Extrapolated POD and Expanded POD. We then use a single baseline parameter in the parameter range to develop a ROM basis for each method above and investigate the error of each ROM method against a full order "truth" solution for the full parameter range. / Master of Science

This thesis presents a study of reduced-order models based on proper orthogonal decomposition applied to non-isothermal transport phenomena in °uidized beds. A numerical °ow solver called Multiphase Flow with Interphase eXchanges (MFIX) was used to generate a database of solution snapshots for proper orthogonal decomposi- tion (POD). Using POD, time independent basis functions were extracted from the data and the governing equations of the numerical solver were projected onto the basis functions to generate reduced-order models. A reduced-order model was constructed that simulates multi-phase isothermal and non-isothermal °ow. In the reduced-order models (ROMs) the large number of partial di®erential equations were replaced by a much smaller number of ordinary di®erential equations. These reduced-order models were applied to two reference cases, a time extrapolation case and a time-dependent period boundary condition case. Three additional acceleration techniques were devel- oped to further improve computational e±ciency of the POD based ROM: 1) Database splitting, 2) Freezing the matrix of the linear system and 3) Time step adjustment. Detailed numerical analysis of both the full-order model, MFIX and the POD-based ROM, including estimating the number of operations and the CPU time per iteration, was performed as part of this study. The results of this investigation show that the reduced-order models are capable of producing qualitatively accurate results with less than 5% error with a two-order of magnitude reduction of computational costs.

An investigation into the responses of microbeams to electric actuations is presented. Attention is focused mainly on the use of microbeams in two important MEMS-based devices: capacitive microswitches and resonant microsensors. Nonlinear models are developed to simulate the behavior of the microbeams in each device. The models account for mid-plane stretching, an applied axial load, a DC electrostatic force, and, for the case of resonant sensors, an AC harmonic force. Further, a novel method that uses a reduced-order model is introduced for simulating the behavior of microbeams under a DC electrostatic force. The presented method shows attractive features, like for example, a high stability near the pull-in and a low computational cost. Thus, it can be of significant benefit to the development of MEMS design software. The static behavior of microbeams under electrostatic forces is studied using two methods. One method employs a shooting technique for solving the boundary-value problem that governs the static behavior. The second method is based on solving an algebraic system of equations obtained from the reduced-order model. Further, the eigenvalue problem describing the vibrations of a microbeam around its statically deflected position is solved using a shooting method to obtain the microbeam mode shapes and natural frequencies. The dynamic behavior of resonant microbeams is also investigated. A perturbation method, the method of multiple scales, is used to obtain two first-order nonlinear ordinary-differential equations that describe the amplitude and phase of the response and its stability. The results show that an inaccurate representation of the system nonlinearities may lead to an erroneous prediction of the nonlinear resonance frequency of a microbeam. The case of three-to-one internal resonance between the lowest two modes is treated. Finally, the reduced-order model is used to study the dynamic behavior of the electrostatically actuated microbeams. The proposed models are validated by comparing their results with experimental results available in the literature. / Master of Science

Lean-premixed combustion has the advantage of low emissions for modern gas turbines, but it is susceptible to thermoacoustic instabilities, which can result in large amplitude pressure oscillations in the combustion chamber. The thermoacoustic limit cycle is generated by the unsteady heat release dynamics coupled to the combustor acoustics. In this dissertation, we focused on reduced-order modeling of the dynamics of a laminar premixed flame. From first principles of combustion dynamics, a physically-based, reduced-order, nonlinear model was developed based on the proper orthogonal decomposition technique and generalized Galerkin method. In addition, the describing function for the flame was measured experimentally and used to identify an empirical nonlinear flame model. Furthermore, a linear acoustic model was developed and identified for the Rijke tube experiment. Closed-loop thermoacoustic modeling using the first principles flame model coupled to the linear acoustics successfully reproduced the linear instability and predicted the thermoacoustic limit cycle amplitude. With the measured experimental flame data and the modeled linear acoustics, the describing function technique was applied for limit cycle analysis. The thermoacoustic limit cycle amplitude was predicted with reasonable accuracy, and the closed-loop model also predicted the performance for a phase shift controller. Some problems found in the predictions for high heat release cases were documented. / Ph. D.

This dissertation proposes a modular fast direct (MFD) analysis method for a class of problems involving a large fixed platform region and a smaller, variable design region. A modular solution algorithm is obtained by first decomposing the problem geometry into platform and design regions. The two regions are effectively detached from one another using basic equivalence concepts. Equivalence principles allow the total system model to be constructed in terms of independent interaction modules associated with the platform and design regions. These modules include interactions with the equivalent surface that bounds the design region. This dissertation discusses how to analyze (fill and factor) each of these modules separately and how to subsequently compose the solution to the original system using the separately analyzed modules. The focus of this effort is on surface integral equation formulations of electromagnetic scattering from conductors and dielectrics. In order to treat large problems, it is necessary to work with sparse representations of the underlying system matrix and other, related matrices. Fortunately, a number of such representations are available. In the following, we will primarily use the adaptive cross approximation (ACA) to fill the multilevel simply sparse method (MLSSM) representation of the system matrix. The MLSSM provides a sparse representation that is similar to the multilevel fast multipole method. Solutions to the linear systems obtained using the modular analysis strategies described above are obtained using direct methods based on the local-global solution (LOGOS) method. In particular, the LOGOS factorization provides a data sparse factorization of the MLSSM representation of the system matrix. In addition, the LOGOS solver also provides an approximate sparse factorization of the inverse of the system matrix. The availability of the inverse eases the development of the MFD method. Because the behavior of the LOGOS factorization is critical to the development of the proposed MFD method, a significant part of this dissertation is devoted to providing additional analyses, improvements, and characterizations of LOGOS-based direct solution methods. These further developments of the LOGOS factorization algorithms and their application to the development of the MFD method comprise the most significant contributions of this dissertation.

Les microsystèmes, qui sont réalisés à partir de technologies micro-électroniques, connaissent un essor scientifique et technologique important grâce à leurs applications qui sont de plus en plus présentes dans la vie courante. Un des microsystèmes très en vogue est le transducteur ultrasonore capacitif micro-usiné, couramment appelé CMUT. Il est utilisé pour transmettre ou réceptionner des ondes ultrasonores et son domaine d’application est très vaste puisqu’on le trouve dans des sondes d’imagerie médicale, dans des hauts parleurs ultra directifs, pour le contrôle non destructif de matériaux… Dans la plupart des applications la puissance acoustique émise par le CMUT doit être très élevée ce qui implique que le CMUT va être utilisé en régime non-linéaire. En outre, même en utilisant des procédés de fabrication avancés, la microplaque mobile constituant le CMUT possède une déformation géométrique dans son état de repos. Il faudra par conséquent tenir compte des non-linéarités et des imperfections géométriques lors de l’analyse statique et dynamique du CMUT.Dans ce travail le modèle multiphysique d’un CMUT est développé en tenant compte des non-linéarités géométriques et électrostatiques ainsi que de la déflexion initiale de la microplaque. Les équations différentielles du mouvement de la microplaque, issues de la théorie des plaques de von Kármán, sont discrétisées spatialement en utilisant la méthode différentielle quadratique. La réponse statique d’un CMUT a été analysée à partir de simulations numériques et d’essais expérimentaux, en considérant des plaques planes et des plaques courbes et on montre qu’une déflexion initiale de la plaque conduit à une augmentation de la tension de pull-in. Le comportement dynamique non-linéaire du CMUT est analysé en discrétisant la variable temporelle et en utilisant la méthode des différences finies. En utilisant la technique de continuation arclength, nous déterminons la réponse en fréquence non-linéaire du CMUT. Suivant la valeur de la tension DC, le CMUT aura un comportement raidissant ou assouplissant. Une validation expérimentale du modèle numérique est réalisée en utilisant des microplaques planes et des microplaques courbes. En particulier nous montrons que l’utilisation de microplaques courbes, dues aux imperfections géométriques, change la réponse en fréquence du CMUT, passant d’un comportement raidissant à un comportement assouplissant, augmente le domaine de bi-stabilité et modifie la topologie de bifurcation.Le modèle numérique est par la suite étendu afin d’analyser les effets du film d’air sur le comportement dynamique de la microplaque en couplant les équations mécaniques du CMUT avec les équations de Reynolds du fluide. Les fréquences de résonance du problème multiphysique sont obtenues par résolution d'un système linéaire amorti. La validation expérimentale et numérique du modèle est effectuée en déterminant les fréquences de résonance du CMUT à des pressions différentes. Nous montrons que l’air comprimé change la réponse dynamique du CMUT par l’ajout d’une raideur et d’un amortissement. La diminution de la pression conduit à une diminution de la fréquence de résonance du système couplé et tend vers la fréquence de résonance de la microplaque. D'autre part la réponse en fréquence du système devient non-linéaire due à la diminution du coefficient d'amortissement. A la pression atmosphérique, on montre que le CMUT a un comportement non-linéaire de type assouplissant lorsque les excitations sont élevées. Le modèle numérique développé est un outil efficace pour analyser les CMUTs et augmenter leurs performanaces. / Micro Electro Mechanical Systems (MEMS) have attracted the interest of scientists and engineers thanks to the variety of their applications and their significant roles in our real life. One of the most important microsystems is the capacitive micromachined ultrasonic transducer (CMUT), which is used for transmitting ultrasonic waves, for instance in medical imaging and therapy. In such applications, a high-transmitted acoustic power is needed which implies driving the CMUT in the nonlinear regime. Moreover, from a manufacturing point of view, the fabrication of a CMUT with a flat surface is extremely difficult even with the recent advances in the fabrication process. Modeling this type of microsystem while including the main sources of nonlinearities and geometric imperfections is a challenging step in understanding its static and dynamic behavior.In this thesis, a multiphysics model of imperfect CMUTs is developed taking into account the geometric and electrostatic nonlinearities. The governing equations of motions are derived from the von Kármán plate theory and spatially discretized using the Differential Quadrature Method (DQM). For the static response, numerical simulations and experimental characterizations have been conducted on flat and curved CMUTs, showing that a positive initial deflection leads to an increase in the pull-in voltage. The nonlinear dynamic behavior of a CMUT is studied by discretizing the time variable using the Finite Difference Method (FDM). The nonlinear frequency and force responses have been determined by combining FDM with the arclength continuation technique. It is shown that the CMUT can exhibit a hardening or softening behavior depending on the DC voltage. An experimental validation of the numerical model is performed for the case of flat and curved microplates. We demonstrate that the geometric imperfection modifies the nonlinear frequency response of a CMUT from hardening to softening, increases its bistability domain and permits the tuning of its bifurcation topology.The numerical model is extended to investigate the effect of an air film on the dynamic behavior of the microplate by coupling the nonlinear mechanical equations with the Reynolds equation. The complex resonance frequencies of the multi-physical problem are determined by solving the damped linear system. An experimental and numerical validation of the model is performed by determining the resonance frequencies at several static pressures. We demonstrate that the air film is able to modify the dynamic response of the CMUT by adding stiffness and damping. By decreasing the static pressure, the resonance frequency of the coupled problem decreases and becomes closer to the natural resonance frequency of the microplate. Moreover, the frequency response of the system becomes nonlinear due the decrease in the damping coefficient. At atmospheric pressure, the softening type behavior of the CMUT is obtained by applying high excitation levels. The presented numerical model is a very efficient tool to understand the nonlinear dynamic behavior CMUTs and to enhance their performances.

abstract: This investigation focuses on the development of uncertainty modeling methods applicable to both the structural and thermal models of heated structures as part of an effort to enable the design under uncertainty of hypersonic vehicles. The maximum entropy-based nonparametric stochastic modeling approach is used within the context of coupled structural-thermal Reduced Order Models (ROMs). Not only does this strategy allow for a computationally efficient generation of samples of the structural and thermal responses but the maximum entropy approach allows to introduce both aleatoric and some epistemic uncertainty into the system. While the nonparametric approach has a long history of applications to structural models, the present investigation was the first one to consider it for the heat conduction problem. In this process, it was recognized that the nonparametric approach had to be modified to maintain the localization of the temperature near the heat source, which was successfully achieved. The introduction of uncertainty in coupled structural-thermal ROMs of heated structures was addressed next. It was first recognized that the structural stiffness coefficients (linear, quadratic, and cubic) and the parameters quantifying the effects of the temperature distribution on the structural response can be regrouped into a matrix that is symmetric and positive definite. The nonparametric approach was then applied to this matrix allowing the assessment of the effects of uncertainty on the resulting temperature distributions and structural response. The third part of this document focuses on introducing uncertainty using the Maximum Entropy Method at the level of finite element by randomizing elemental matrices, for instance, elemental stiffness, mass and conductance matrices. This approach brings some epistemic uncertainty not present in the parametric approach (e.g., by randomizing the elasticity tensor) while retaining more local character than the operation in ROM level. The last part of this document focuses on the development of “reduced ROMs” (RROMs) which are reduced order models with small bases constructed in a data-driven process from a “full” ROM with a much larger basis. The development of the RROM methodology is motivated by the desire to optimally reduce the computational cost especially in multi-physics situations where a lack of prior understanding/knowledge of the solution typically leads to the selection of ROM bases that are excessively broad to ensure the necessary accuracy in representing the response. It is additionally emphasized that the ROM reduction process can be carried out adaptively, i.e., differently over different ranges of loading conditions. / Dissertation/Thesis / Doctoral Dissertation Mechanical Engineering 2019

abstract: Recent studies of the occurrence of post-flutter limit cycle oscillations (LCO) of the F-16 have provided good support to the long-standing hypothesis that this phenomenon involves a nonlinear structural damping. A potential mechanism for the appearance of nonlinearity in the damping are the nonlinear geometric effects that arise when the deformations become large enough to exceed the linear regime. In this light, the focus of this investigation is first on extending nonlinear reduced order modeling (ROM) methods to include viscoelasticity which is introduced here through a linear Kelvin-Voigt model in the undeformed configuration. Proceeding with a Galerkin approach, the ROM governing equations of motion are obtained and are found to be of a generalized van der Pol-Duffing form with parameters depending on the structure and the chosen basis functions. An identification approach of the nonlinear damping parameters is next proposed which is applicable to structures modeled within commercial finite element software. The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses. / Dissertation/Thesis / Masters Thesis Mechanical Engineering 2015

abstract: This investigation develops small-size reduced order models (ROMs) that provide an accurate prediction of the response of only part of a structure, referred to as component-centric ROMs. Four strategies to construct such ROMs are presented, the first two of which are based on the Craig-Bampton Method and start with a set of modes for the component of interest (the β component). The response in the rest of the structure (the α component) induced by these modes is then determined and optimally represented by applying a Proper Orthogonal Decomposition strategy using Singular Value Decomposition. These first two methods are effectively basis reductions techniques of the CB basis. An approach based on the “Global - Local” Method generates the “global” modes by “averaging” the mass property over α and β comp., respectively (to extract a “coarse” model of α and β) and the “local” modes orthogonal to the “global” modes to add back necessary “information” for β. The last approach adopts as basis for the entire structure its linear modes which are dominant in the β component response. Then, the contributions of other modes in this part of the structure are approximated in terms of those of the dominant modes with close natural frequencies and similar mode shapes in the β component. In this manner, the non-dominant modal contributions are “lumped” onto the dominant ones, to reduce the number of modes for a prescribed accuracy. The four approaches are critically assessed on the structural finite element model of a 9-bay panel with the modal lumping-based method leading to the smallest sized ROMs. Therefore, it is extended to the nonlinear geometric situation and first recast as a rotation of the modal basis to achieve unobservable modes. In the linear case, these modes completely disappear from the formulation owing to orthogonality. In the nonlinear case, however, the generalized coordinates of these modes are still present in the nonlinear terms of the observable modes. A closure-type algorithm is then proposed to eliminate the unobserved generalized coordinates. This approach, its accuracy and computational savings, was demonstrated on a simple beam model and the 9-bay panel model. / Dissertation/Thesis / Doctoral Dissertation Mechanical Engineering 2017

abstract: This investigation is focused on the consideration of structural uncertainties in nearly-straight pipes conveying fluid and on the effects of these uncertainties on the dynamic response and stability of those pipes. Of interest more specifically are the structural uncertainties which affect directly the fluid flow and its feedback on the structural response, e.g., uncertainties on/variations of the inner cross-section and curvature of the pipe. Owing to the complexity of introducing such uncertainties directly in finite element models, it is desired to proceed directly at the level of modal models by randomizing simultaneously the appropriate mass, stiffness, and damping matrices. The maximum entropy framework is adopted to carry out the stochastic modeling of these matrices with appropriate symmetry constraints guaranteeing that the nature, e.g., divergence or flutter, of the bifurcation is preserved when introducing uncertainty. To support the formulation of this stochastic ROM, a series of finite element computations are first carried out for pipes with straight centerline but inner radius varying randomly along the pipe. The results of this numerical discovery effort demonstrate that the dominant effects originate from the variations of the exit flow speed, induced by the change in inner cross-section at the pipe end, with the uncertainty on the cross-section at other locations playing a secondary role. Relying on these observations, the stochastic reduced order model is constructed to model separately the uncertainty in inner cross-section at the pipe end and at other locations. Then, the fluid related mass, damping, and stiffness matrices of this stochastic reduced order model (ROM) are all determined from a single random matrix and a random variable. The predictions from this stochastic ROM are found to closely match the corresponding results obtained with the randomized finite element model. It is finally demonstrated that this stochastic ROM can easily be extended to account for the small effects due to uncertainty in pipe curvature. / Dissertation/Thesis / Masters Thesis Mechanical Engineering 2017