Reservoir simulation is an important tool for decision making and field development management. It enables reservoir engineers to predict reservoir production performance, update an existing model to reproduce monitoring data, assess alternative field development scenarios and design robust production optimization strategies by taking into account the existing uncertainties. A big obstacle in automating model calibration and production optimization approaches is the massive computation required to predict the response of real reservoirs under proposed changes in the model inputs. To speed up reservoir response predictions without compromising accuracy, fast surrogate models have been proposed. These models are either derived by preserving the physics of the involved processes (e.g. mass balance equations) to provide reliable long-range predictions or are developed based solely on statistical relations, in which case they can only provide short-range predictions due to the absence of the physical processes that govern the long-term behavior of the reservoir. We present an alternative solution that combines the advantages of both statistics-based and physics-based methods by deriving the flow predictions in complex two-dimensional models from one-dimensional flow network models. The existing injection/production wells in the original model form the nodes or vertices of the flow network. Each pair of wells (nodes) in the flow network is connected using a one-dimensional numerical simulation model; hence, the entire reservoir is reduced to a connected network of one-dimensional simulation models where the coupling between the individual one-dimensional models is enforced at the nodes where network edges intersect. The proposed flow network model provides a useful and fast tool for characterizing inter-well connectivity, estimating drainage volume between each pair of wells, and predicting reservoir production over an extended period of time for optimization purposes. We estimate the parameters of the flow network model using a robust training approach to ensure that the flow network model reproduces the response of the original full model under a wide range of development strategies. This step helps preserve the flow network model's predictive power during the production optimization when development strategies can change at different iterations. The robust networks training and the subsequent production optimization iterations are computationally efficient as they are performed with the faster flow network model. We demonstrate the effectiveness and applicability of our proposed flow network modeling approach to rapid production optimization using two-phase waterflooding simulations in synthetic and benchmark models.

Data generated by complex, computational models can provide highly accurate predictions of hydrological and hydrodynamic data in multiple dimensions. Unfortunately, however, for large data sets, running these models is often timeconsuming and computationally expensive. Thus, finding a way to reduce the running time of these models, while still producing comparable results, is of notable interest. The Inference Engine is a proposed system for doing just this. It takes previously generated model data and uses them to predict additional data. Its performance, both accuracy and running time, has been compared to the performance of the actual models, in increasingly difficult data prediction tasks, and it is able, with sufficient accuracy, to quickly predict unknown model data.

abstract: This thesis focuses on the continued extension, validation, and application of combined thermal-structural reduced order models for nonlinear geometric problems. The first part of the thesis focuses on the determination of the temperature distribution and structural response induced by an oscillating flux on the top surface of a flat panel. This flux is introduced here as a simplified representation of the thermal effects of an oscillating shock on a panel of a supersonic/hypersonic vehicle. Accordingly, a random acoustic excitation is also considered to act on the panel and the level of the thermo-acoustic excitation is assumed to be large enough to induce a nonlinear geometric response of the panel. Both temperature distribution and structural response are determined using recently proposed reduced order models and a complete one way, thermal-structural, coupling is enforced. A steady-state analysis of the thermal problem is first carried out that is then utilized in the structural reduced order model governing equations with and without the acoustic excitation. A detailed validation of the reduced order models is carried out by comparison with a few full finite element (Nastran) computations. The computational expedience of the reduced order models allows a detailed parametric study of the response as a function of the frequency of the oscillating flux. The nature of the corresponding structural ROM equations is seen to be of a Mathieu-type with Duffing nonlinearity (originating from the nonlinear geometric effects) with external harmonic excitation (associated with the thermal moments terms on the panel). A dominant resonance is observed and explained. The second part of the thesis is focused on extending the formulation of the combined thermal-structural reduced order modeling method to include temperature dependent structural properties, more specifically of the elasticity tensor and the coefficient of thermal expansion. These properties were assumed to vary linearly with local temperature and it was found that the linear stiffness coefficients and the "thermal moment" terms then are cubic functions of the temperature generalized coordinates while the quadratic and cubic stiffness coefficients were only linear functions of these coordinates. A first validation of this reduced order modeling strategy was successfully carried out. / Dissertation/Thesis / M.S. Aerospace Engineering 2011

Composite materials hold numerous advantages over conventional aircraft grade metals. These include high stiffness/strength-to-weight ratios and beneficial stiffness coupling typically used for aeroelastic tailoring. Due to the complexity of modeling composites, designers often select safe, simple geometry and layup schedules for their wing/blade cross-sections. An example of this might be a box-beam made up of 4 laminates, all of which are quasi-isotropic. This results in neglecting more complex designs that might yield a more effective solution, but require a greater analysis effort. The present work aims to show that the incorporation of complex cross-sections are feasible in the early design process through the use of cross-sectional analysis in conjunction with Timoshenko beam theory. It is important to note that in general, these cross-sections can be inhomogeneous: made up of any number of various materials systems. In addition, these materials could all be anisotropic in nature. The geometry of the cross-sections can take on any shape. Through this reduced order modeling scheme, complex structures can be reduced to 1 dimensional beams. With this approach, the elastic behavior of the structure can be captured, while also allowing for accurate 3D stress and strain recovery. This efficient structural modeling would be ideal in the preliminary design optimization of a wing structure. Furthermore, in conjunction with an efficient unsteady aerodynamic model such as the doublet lattice method, the dynamic aeroelastic stability can also be efficiently captured. This work introduces a comprehensively verified, open source python API called AeroComBAT (Aeroelastic Composite Beam Analysis Tool). By leveraging cross-sectional analysis, Timoshenko beam theory, and unsteady doublet-lattice method, this package is capable of efficiently conducting linear static structural analysis, normal mode analysis, and dynamic aeroelastic analysis. AeroComBAT can have a significant impact on the design process of a composite structure, and would be ideally implemented as part of a design optimization. / Master of Science

No description available.

In this thesis, we develop a data-driven correction reduced order model (DDC-ROM) for numerical simulation of fluid flows. The general DDC-ROM involves two stages: (1) we apply ROM filtering (such as ROM projection) to the full order model (FOM) and construct the filtered ROM (F-ROM). (2) We use data-driven modeling to model the nonlinear interactions between resolved and unresolved modes, which solves the F-ROM's closure problem. In the DDC-ROM, a linear or quadratic ansatz is used in the data-driven modeling step. In this thesis, we propose a new cubic ansatz. To get the unknown coefficients in our ansatz, we solve an optimization problem that minimizes the difference between the FOM data and the ansatz. We test the new DDC-ROM in the numerical simulation of the one-dimensional Burgers equation with a small diffusion coefficient. Furthermore, we perform a cross-validation of the DDC-ROM to investigate whether it can be successful in computational settings that are different from the training regime. / M.S. / Practical engineering and scientific problems often require the repeated simulation of unsteady fluid flows. In these applications, the computational cost of high-fidelity full-order models can be prohibitively high. Reduced order models (ROMs) represent efficient alternatives to brute force computational approaches. In this thesis, we propose a data-driven correction ROM (DDC-ROM) in which available data and an optimization problem are used to model the nonlinear interactions between resolved and unresolved modes. In order to test the new DDC-ROM's predictability, we perform its cross-validation for the one-dimensional viscous Burgers equation and different training regimes.

No description available.

The need for low order models capable of performing damage identification has become apparent in many structural dynamics applications where structural health monitoring and damage prognosis programs are implemented. These programs require that damage identification routines have low computational requirements and be reliable with some quantifiable degree of accuracy. Response surface metamodels (RSMs) are proposed to fill this need. Popular in the fields of chemical and industrial engineering, RSMs have only recently been applied in the field of structural dynamics and to date there have been no studies which fully demonstrate the potential of these methods. In this thesis, several RSMs are developed in order to demonstrate the potential of the methodology. They are shown to be robust to noise (experimental variability) and have success in solving the damage identification problem, both locating and quantifying damage with some degree of accuracy, for both linear and nonlinear systems. A very important characteristic of the RSMs developed in this thesis is that they require very little information about the system in order to generate relationships between damage indicators and measureable system responses for both linear and nonlinear structures. As such, the potential of these methods for damage identification has been demonstrated and it is recommended that these methods be developed further. / Master of Science

Mathematical models describing real world processes are becoming increasingly complex to better match the dynamics of the true system. While this is a positive step towards more complete knowledge of our world, numerical evaluations of these models become increasingly computationally inefficient, requiring increased resources or time to evaluate. This has led to the need for simplified surrogates to these complex mathematical models. A growing surrogate modeling solution is with the usage of neural networks. Neural networks (NN) are known to generalize an approximation across a diverse dataset and minimize the solution along complex nonlinear boundaries. Additionally, these surrogate models can be found using only incomplete knowledge of the true dynamics. However, NN surrogates often suffer from a lack of interpretability, where the decisions made in the training process are not fully understood, and the roles of individual neurons are not well defined. We present two solutions towards this lack of interpretability. The first focuses on mimicking polynomial chaos (PC) modeling techniques, modifying the structure of a NN to produce polynomial approximations of the underlying dynamics. This methodology allows for an extractable meaning from the network and results in improvement in accuracy over traditional PC methods. Secondly, we examine the construction of a reduced order modeling scheme using NN autoencoders, guiding the decisions of the training process to better match the real dynamics. This guiding process is performed via a physics-informed (PI) penalty, resulting in a speed-up in training convergence, but still results in poor performance compared to traditional schemes. / Master of Science / The world is an elaborate system of relationships between diverse processes. To accurately represent these relationships, increasingly complex models are defined to better match what is physically seen. These complex models can lead to issues when trying to use them to predict a realistic outcome, either requiring immensely powerful computers to run the simulations or long amounts of time to present a solution. To fix this, surrogates or approximations to these complex models are used. These surrogate models aim to reduce the resources needed to calculate a solution while remaining as accurate to the more complex model as possible. One way to make these surrogate models is through neural networks. Neural networks try to simulate a brain, making connections between some input and output given to the network. In the case of surrogate modeling, the input is some current state of the true process, and the output is what is seen later from the same system. But much like the human brain, the reasoning behind why choices are made when connecting the input and outputs is often largely unknown. Within this paper, we seek to add meaning to neural network surrogate models in two different ways. In the first, we change what each piece in a neural network represents to build large polynomials (e.g., $x^5 + 4x^2 + 2$) to approximate the larger complex system. We show that the building of these polynomials via neural networks performs much better than traditional ways to construct them. For the second, we guide the choices made by the neural network by enforcing restrictions in what connections it can make. We do this by using additional information from the larger system to ensure the connections made focus on the most important information first before trying to match the less important patterns. This guiding process leads to more information being captured when the surrogate model is compressed into only a few dimensions compared to traditional methods. Additionally, it allows for a faster learning time compared to similar surrogate models without the information.

In this thesis we propose, analyze, and investigate numerically a novel two-level Galerkin reduced order model (2L-ROM) for the efficient and accurate numerical simulation of the steady Navier-Stokes equations. In the first step of the 2L-ROM, a relatively low-dimensional nonlinear system is solved. In the second step, the Navier-Stokes equations are linearized around the solution found in the first step, and a higher-dimensional system for the linearized problem is solved. We prove an error bound for the new 2L-ROM and compare it to the standard Galerkin ROM, or one-level ROM (1L-ROM), in the numerical simulation of the steady Burgers equation. The 2L-ROM significantly decreases (by a factor of 2 and even 3) the 1L-ROM computational cost, without compromising its numerical accuracy. / Master of Science / In this thesis we introduce a new method for efficiently and accurately simulating fluid flow, the Navier-Stokes equations, called the two-level Galerkin reduced order model (2L-ROM). The 2L-ROM involves solving a relatively low-dimensional nonlinear system in the first step, followed by a higher-dimensional linearized system in the second step. We show that this method produces highly accurate results while significantly reducing computational costs compared to previous methods. We provide a comparison between the 2L-ROM and the standard Galerkin ROM, or one-level ROM (1L-ROM), by modeling the steady Burgers equation, as an example. Our results demonstrate that the 2L-ROM reduces the computational cost of the 1L-ROM by a factor of 2 to 3 without sacrificing accuracy.