Designing and implementing coupled Earth System Models (ESMs) is a challenge for climate scientists and software engineers alike. Coupled models incorporate two or more independent numerical models into a single application, allowing for the simulation of complex feedback effects. As ESMs increase in sophistication, incorporating higher fidelity models of geophysical processes, developers are faced with the issue of managing increasing software complexity. Recently, reusable coupling software has emerged to aid developers in building coupled models. Effective reuse of coupling infrastructure means increasing the number of coupling functions reused, minimizing code duplication, reducing the development time required to couple models, and enabling flexible composition of coupling infrastructure with existing constituent model implementations. Despite the widespread availability of software packages that provide coupling infrastructure, effective reuse of coupling technologies remains an elusive goal: coupling models is effort-intensive, often requiring weeks or months of developer time to work through implementation details, even when starting from a set of existing software components. Coupling technologies are never used in isolation: they must be integrated with multiple existing constituent models to provide their primary services, such as model-to-model data communication and transformation. Unfortunately, the high level of interdependence between coupling concerns and scientific concerns has resulted in high interdependence between the infrastructure code and the scientific code within a model’s implementation. These dependencies are a source of complexity which tends to reduce reusability of coupling infrastructure. This dissertation presents mechanisms for increasing modeler productivity based on improving reuse of coupling infrastructure and raising the level of abstraction at which modelers work. This dissertation argues that effective reuse of coupling technologies can be achieved by decomposing existing coupling technologies into a salient set of implementation-independent features required for coupling high-performance models, increasing abstraction levels at which model developers work, and facilitating integration of coupling infrastructure with constituent models via component-based modularization of coupling features. The contributions of this research include: (1) a comprehensive feature model that identifies the multi-dimensional design space of coupling technologies used in high-performance Earth System Models, (2) Cupid, a domain-specific language and compiler for specifying coupling configurations declaratively and generating their implementations automatically, and (3) Component-based Coupling Operators (CC-Ops), a modular approach to code reuse of coupling infrastructure based on component technologies for high-performance scientific settings. The Cupid domain-specific language is evaluated by specifying a coupling configuration for an example fluid dynamics model and measuring the amount of code generated by the Cupid compiler compared to a hand coded version. The CC-Op approach is evaluated by implementing several CC-Ops using an existing high-performance component framework and measuring performance in terms of scalability and overhead.

[Truncated abstract] Rocks, as the host to natural chains of coupled thermal, hydraulic and mechanical processes, are heterogeneous at a variety of length scales, and in their mechanical properties, as well as in the hydraulic and thermal transport properties. Rock heterogeneity affects the ultimate hydro-carbon recovery or geothermal energy production. This heterogeneity has been considered one important and difficult problem that needs to be taken into account for its effect on the coupled processes. The aim of this thesis is to investigate the effect of rock heterogeneity on multi-physical processes. A fully coupled finite element model, hereinafter referred to as a porosity-based model (PBM) was developed to characterise the thermal-hydraulic-mechanical (THM) coupling processes. The development of the PBM consists of a two-staged workflow. First, based on poromechanics, porosity, one of the inherent rock properties, was derived as a variant function of the thermal, hydraulic and mechanical effects. Then, empirical relations or experimental results, correlating porosity with the mechanical, hydraulic and thermal properties, were incorporated as the coupling effects. In the PBM, the bulk volume of the model is assumed to be changeable. The rate of the volumetric strain was derived as the difference of two parts: the first part is the change in volume per unit of volume and per unit of time (this part was traditionally considered the rate of volumetric strain); and the second is the product of the first part and the volumetric strain. The second part makes the PBM a significant advancement of the models reported in the literature. ... impact of the rock heterogeneity on the hydro-mechanical responses because of the requirement of large memory and long central processing unit (CPU) time for the 3D applications. In the 2D PBM applications, as the thermal boundary condition applied to the rock samples containing some fractures, the pore pressure is generated by the thermal gradient. Some pore pressure islands can be generated as the statistical model and the digital image model are applied to characterise the initial porosity distribution. However, by using the homogeneous model, this phenomenon cannot be produced. In the 3D PBM applications, the existing fractures become the preferential paths for the fluid flowing inside the numerical model. The numerical results show that the PBM is sufficiently reliable to account for the rock mineral distribution in the hydro-mechanical coupling processes. The applications of the statistical method and the digital image processing technique make it possible to visualise the rock heterogeneity effect on the pore pressure distribution and the heat dissipation inside the rock model. Monitoring the fluid flux demonstrates the impact of the rock heterogeneity on the fluid product, which concerns petroleum engineering. The overall fluid flux (OFF) is mostly overestimated when the rock and fluid properties are assumed to be homogeneous. The 3D PBM application is an example. As the rock is heterogeneous, the OFF by the digital core is almost the same as that by the homogeneous model (this is due to that some fractures running through the digital core become the preferential path for the fluid flow), and around 1.5 times of that by the statistical model.

Die vorliegende Arbeit stellt eine mehrskalige Modellierungs-Strategie für die Beschreibung magnetorheologischer Elastomere (MRE) vor. Diese ermöglicht die Betrachtung von MRE sowohl auf der Mikroskala, wo die heterogene Mikrostruktur bestehend aus Partikeln und Matrix explizit aufgelöst ist, als auch auf der Makroskala, in welcher das MRE als homogener magnetisch aktiver Körper aufzufassen ist. Auf beiden Skalen kommt dabei eine Kontinuumsformulierung des gekoppelten magneto-mechanischen Feldproblems mit Gültigkeit für finite Deformationen zum Einsatz, wobei die Lösung des Systems partieller Differentialgleichungen mittels der Finite-Elemente-Methode erfolgt. Ausgehend von einer experimentellen Charakterisierung der Konstituenten werden Materialmodelle für die elastomere Matrix sowie Carbonyleisen- und Neodym-Eisen-Bor-Partikel formuliert und mittels dieser Daten kalibriert. Im nächsten Schritt erfolgt die Analyse des effektiven Verhaltens hart- und weichmagnetischer MRE auf Basis von numerischen Homogenisierungen verschiedener mikroskopischer Partikelverteilungen und den Materialmodellen für die Konstituenten. Um weiterhin die effiziente Simulation makroskopischer MRE-Proben und -Bauteile zu ermöglichen, ist daran anschließend die Entwicklung und Parametrisierung eines Makromodells ausgehend von mikroskopisch generierten Datensätzen beschrieben. Mit diesem für isotrope, weichmagnetische und elastische MRE gültigen Modell werden abschließend Simulationen des magnetostriktiven sowie des magnetorheologischen Effektes verschiedener Proben durchgeführt. / In this contribution, a strategy for the multiscale modeling of magnetorheological elastomers (MREs) is presented. It allows to consider these materials on the microscopic scale, where the heterogeneous microstructure consisting of an elastomer matrix and embedded magnetizable particles is explicitly resolved, as well as the macroscopic scale, where the MRE is considered to be a homogeneous magneto-active body. On both scales, a continuum formulation of the coupled magneto-mechanical boundary value problem valid for finite strains is applied. The solution of the system of partial differential equations is calculated by using the finite element method. Starting with an experimental characterization of the individual constituents, constitutive models for the elastomer matrix as well as carbonyl iron and neodymium-iron-boron particles are formulated and adjusted to experimental data. In a next step, basic effective properties of magnetically soft and hard MREs are analyzed by using a computational homogenization scheme, where different geometrical arrangements of the particles on the microscale are considered. In order to enable the efficient simulation of macroscopic MRE samples and components, the developement and parametrization of a macroscopic model based on a microscopically generated data basis is described. With this model which is applicable for isotropic, magnetically soft and elastic MREs, simulations of the magnetostrictive and magnetorheological effects of several sample geometries are performed.

Indiana University-Purdue University Indianapolis (IUPUI) / A wave rotor constant volume combustor (WRCVC) was designed and built as a collaborative work of Rolls Royce LibertyWorks, Indiana University-Purdue University at Indianapolis (IUPUI), and Purdue University, and ran experimental tests at Purdue's Zucrow Laboratories in 2009. Instrumentation of the WRCVC rig inlet flow included temperature and pressure transducers upstream of the venturi and at the fuel delivery plane. Other instrumentation included exhaust pressures and temperatures. In addition, ion sensors, dynamic pressure sensors, and accelerometers were used to instrument the rotating hardware. The rig hardware included inlet guide vanes directly in front of the rotating hardware, which together with concern for damage potential, prevented use of any pressure transducers at the entrance to the rotor. For this reason, a complete understanding of the conditions at the WRCVC inlet is unavailable, requiring simulations of the WRCVC to estimate the inlet pressure at a specific operating condition based on airflow. The operation of a WRCVC rig test is a sequence of events over a short time span. These events include introduction of the main air flow followed by time-sequenced delivery of fuel, lighting of the ignition source, and the combustion sequence. The fast changing conditions in the rig inlet hardware make necessary a time-dependent computation of the rig inlet section in order to simulate the overall rig operation. The chosen method for computing inlet section temperature and pressure was a time-dependent lumped volume model of the inlet section hardware, using a finite difference modified Euler predictor-corrector method for computing the continuity and energy equations. This is coupled with perfect gas prediction of venturi air and fuel flow rates, pressure drag losses at the fuel nozzles, pressure losses by mass addition of the fuel or nitrogen purge, friction losses at the inlet guide vanes, and a correlation of the non-dimensional flow characteristics of the WRCVC. The flow characteristics of the WRCVC are computed by varying the non-dimensional inlet stagnation pressure and the WRCVC's operational conditions, assuming constant rotational speed and inlet stagnation temperature. This thesis documents the creation of a computer simulation of the entire WRCVC rig, to understand the pressure losses in the inlet system and the dynamic coupling of the inlet section and the WRCVC, so that an accurate prediction of the WRCVC rotor inlet conditions can be computed. This includes the computational development of the WRCVC upstream rig dynamic model, the background behind supporting computations, and results for one test sequence. The computations provide a clear explanation of why the pressures at the rotor inlet differ so much from the upstream measured values. The pressure losses correlate very well with the computer predictions and the dynamic response tracks well with the estimation of measured airflow. A simple Fortran language computer program listing is included, which students can use to simulate charging or discharging of a container.

Indiana University-Purdue University Indianapolis (IUPUI) / Some problems of complex systems are internal to the system whereas other problems exist peripherally; two such problems will be explored in this thesis. First, is the issue of excessive jerk from instantaneous velocity demand changes produced by an adaptive cruise control system. Calculations will be demonstrated and an example control solution will be proposed in Chapter 3. Second, is the issue of a non-perfect plant, called an uncertain or corrupted plant. In initial control analysis, the adaptive cruise control systems are assumed to have a perfect plant; that is to say, the plant always behaves as commanded. In reality, this is seldom the case. Plant corruption may come from a variation in performance through use or misuse, or from noise or imperfections in the sensor signal data. A model for plant corruption is introduced and methods for analysis and compensation are explored in Chapter 4. To facilitate analysis, Chapter 2 discusses the concept of system identification, an order reduction tool which is employed herein. Adaptive cruise control systems are also discussed with special emphasis on the situations most likely to employ jerk limitation.