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.

Buildings in the United States, account for over 68 percent of electricity consumed, 39 percent of total energy use, and 38 percent of the carbon dioxide emissions. By the year 2035, about 75% of the U.S. building sector will be either new or renovated. The energy efficiency requirements of current building codes would have a significant impact on future energy use, hence, one of the most widely accepted solutions to slowing the growth rate of GHG emissions and then reversing it involves a stringent adoption of building energy codes. A large number of building energy codes exist and a large number of studies which state the energy savings possible through code compliance. However, most codes are difficult to comprehend and require an extensive understanding of the code, the compliance paths, all mandatory and prescriptive requirements as well as the strategy to convert the same to energy model inputs. This paper provides a simplified solution for the entire process by providing an easy to use interface for code compliance and energy simulation through a spreadsheet based tool, the ECCO or the Energy Code COmpliance Tool. This tool provides a platform for a more detailed analysis of building codes as applicable to each and every individual building in each climate zone. It also facilitates quick building energy simulation to determine energy savings achieved through code compliance. This process is highly beneficial not only for code compliance, but also for identifying parameters which can be improved for energy efficiency. Code compliance is simplified through a series of parametric runs which generates the minimally compliant baseline building and 30% beyond code building. This tool is seen as an effective solution for architects and engineers for an initial level analysis as well as for jurisdictions as a front-end diagnostic check for code compliance.  

A method of determining nanoparticle temperature through fluorescence intensity levels is described. Intracellular processes are often tracked through the use of fluorescence tagging, and ideal temperatures for many of these processes are unknown. Through the use of fluorescence-based thermometry, cellular processes such as intracellular enzyme movement can be studied and their respective temperatures established simultaneously. Polystyrene and silica nanoparticles are synthesized with a variety of temperature-sensitive dyes such as BODIPY, rose Bengal, Rhodamine dyes 6G, 700, and 800, and Nile Blue A and Nile Red. Photographs are taken with a QImaging QM1 Questar EXi Retiga camera while particles are heated from 25 to 70 C and excited at 532 nm with a Coherent DPSS-532 laser. Photographs are converted to intensity images in MATLAB and analyzed for fluorescence intensity, and plots are generated in MATLAB to describe each dye's intensity vs temperature. Regression curves are created to describe change in fluorescence intensity over temperature. Dyes are compared as nanoparticle core material is varied. Large particles are also created to match the camera's optical resolution capabilities, and it is established that intensity values increase proportionally with nanoparticle size. Nile Red yielded the closest-fit model, with R2 values greater than 0.99 for a second-order polynomial fit. By contrast, Rhodamine 6G only yielded an R2 value of 0.88 for a third-order polynomial fit, making it the least reliable dye for temperature measurements using the polynomial model. Of particular interest in this work is Nile Blue A, whose fluorescence-temperature curve yielded a much different shape from the other dyes. It is recommended that future work describe a broader range of dyes and nanoparticle sizes, and use multiple excitation wavelengths to better quantify each dye's quantum efficiency. Further research into the effects of nanoparticle size on fluorescence intensity levels should be considered as the particles used here greatly exceed 2 ìm. In addition, Nile Blue A should be further investigated as to why its fluorescence-temperature curve did not take on a characteristic shape for a temperature-sensitive dye in these experiments.

The focus of this investigation is on the renewed assessment of nonlinear reduced order models (ROM) for the accurate prediction of the geometrically nonlinear response of a curved beam. In light of difficulties encountered in an earlier modeling effort, the various steps involved in the construction of the reduced order model are carefully reassessed. The selection of the basis functions is first addressed by comparison with the results of proper orthogonal decomposition (POD) analysis. The normal basis functions suggested earlier, i.e. the transverse linear modes of the corresponding flat beam, are shown in fact to be very close to the POD eigenvectors of the normal displacements and thus retained in the present effort. A strong connection is similarly established between the POD eigenvectors of the tangential displacements and the dual modes which are accordingly selected to complement the normal basis functions. The identification of the parameters of the reduced order model is revisited next and it is observed that the standard approach for their identification does not capture well the occurrence of snap-throughs. On this basis, a revised approach is proposed which is assessed first on the static, symmetric response of the beam to a uniform load. A very good to excellent matching between full finite element and ROM predicted responses validates the new identification procedure and motivates its application to the dynamic response of the beam which exhibits both symmetric and antisymmetric motions. While not quite as accurate as in the static case, the reduced order model predictions match well their full Nastran counterparts and support the reduced order model development strategy.

While much effort in Stirling engine development is placed on making the high-temperature region of the Stirling engine warmer, this research explores methods to lower the temperature of the cold region by improving heat transfer in the cooler. This paper presents heat transfer coefficients obtained for a Stirling engine heat exchanger with oscillatory flow. The effects of oscillating frequency and input heat rate on the heat transfer coefficients are evaluated and details on the design and development of the heat exchanger test apparatus are also explained. Featured results include the relationship between overall heat transfer coefficients and oscillation frequency which increase from 21.5 to 46.1 Wm-2K-1 as the oscillation frequency increases from 6.0 to 19.3 Hz. A correlation for the Nusselt number on the inside of the heat exchange tubes in oscillatory flow is presented in a concise, dimensionless form in terms of the kinetic Reynolds number as a result of a statistical analysis. The test apparatus design is proven to be successful throughout its implementation due to the usefulness of data and clear trends observed. The author is not aware of any other publicly-available research on a Stirling engine cooler to the extent presented in this paper. Therefore, the present results are analyzed on a part-by-part basis and compared to segments of other research; however, strong correlations with data from other studies are not expected. The data presented in this paper are part of a continuing effort to better understand heat transfer properties in Stirling engines as well as other oscillating flow applications.

Thermal interface materials (TIMs) are extensively used in thermal management applications especially in the microelectronics industry. With the advancement in microprocessors design and speed, the thermal management is becoming more complex. With these advancements in microelectronics, there have been parallel advancements in thermal interface materials. Given the vast number of available TIM types, selection of the material for each specific application is crucial. Most of the metrologies currently available on the market are designed to qualify TIMs between two perfectly flat surfaces, mimicking an ideal scenario. However, in realistic applications parallel surfaces may not be the case. In this study, a unique characterization method is proposed to address the need for TIMs characterization between non-parallel surfaces. Two different metrologies are custom-designed and built to measure the impact of tilt angle on the performance of TIMs. The first metrology, Angular TIM Tester, is based on the ASTM D5470 standard with flexibility to perform characterization of the sample under induced tilt angle of the rods. The second metrology, Bare Die Tilting Metrology, is designed to validate the performance of TIM under induced tilt angle between the bare die and the cooling solution in an "in-situ" package testing format. Several types of off-the-shelf thermal interface materials were tested and the results are outlined in the study. Data were collected using both metrologies for all selected materials. It was found that small tilt angles, up to 0.6°, have an impact on thermal resistance of all materials especially for in-situ testing. In addition, resistance change between 0° and the selected tilt angle was found to be in close agreement between the two metrologies for paste-based materials and phase-change material. However, a clear difference in the thermal performance of the tested materials was observed between the two metrologies for the gap filler materials.

The focus of this research is to investigate methods for material substitution for the purpose of re-engineering legacy systems that involves incomplete information about form, fit and function of replacement parts. The primary motive is to extract as much useful information about a failed legacy part as possible and use fuzzy logic rules for identifying the unknown parameter values. Machine elements can fail by any number of failure modes but the most probable failure modes based on the service condition are considered critical failure modes. Three main parameters are of key interest in identifying the critical failure mode of the part. Critical failure modes are then directly mapped to material properties. Target material property values are calculated from material property values obtained from the originally used material and from the design goals. The material database is searched for new candidate materials that satisfy the goals and constraints in manufacturing and raw stock availability. Uncertainty in the extracted data is modeled using fuzzy logic. Fuzzy member functions model the imprecise nature of data in each available parameter and rule sets characterize the imprecise dependencies between the parameters and makes decisions in identifying the unknown parameter value based on the incompleteness. A final confidence level for each material in a pool of candidate material is a direct indication of uncertainty. All the candidates satisfy the goals and constraints to varying degrees and the final selection is left to the designer's discretion. The process is automated by software that inputs incomplete data; uses fuzzy logic to extract more information and queries the material database with a constrained search for finding candidate alternatives.

Next generation gas turbines will be required to produce low concentrations of pollutants such as oxides of nitrogen (NOx), carbon monoxide (CO), and soot. In order to design gas turbines which produce lower emissions it is essential to have computational tools to help designers. Over the past few decades, computational fluid dynamics (CFD) has played a key role in the design of turbomachinary and will be heavily relied upon for the design of future components. In order to design components with the least amount of experimental rig testing, the ensemble of submodels used in simulations must be known to accurately predict the component's performance. The present work aims to validate a CFD model used for a reverse flow, rich-burn, quick quench, lean-burn combustor being developed at Honeywell. Initially, simulations are performed to establish a baseline which will help to assess impact to combustor performance made by changing CFD models. Rig test data from Honeywell is compared to these baseline simulation results. Reynolds averaged Navier-Stokes (RANS) and Large Eddy Simulation (LES) turbulence models are both used with the presumption that the LES turbulence model will better predict combustor performance. One specific model, the fuel spray model, is evaluated next. Experimental data of the fuel spray in an isolated environment is used to evaluate models for the fuel spray and a new, simpler approach for inputting the spray boundary conditions (BC) in the combustor is developed. The combustor is simulated once more to evaluate changes from the new fuel spray boundary conditions. This CFD model is then used in a predictive simulation of eight other combustor configurations. All computer simulations in this work were preformed with the commercial CFD software ANSYS FLUENT. NOx pollutant emissions are predicted reasonably well across the range of configurations tested using the RANS turbulence model. However, in LES, significant under predictions are seen. Causes of the under prediction in NOx concentrations are investigated. Temperature metrics at the exit of the combustor, however, are seen to be better predicted with LES.

In this work, we focused on the stability and reducibility of quasi-periodic systems. We examined the quasi-periodic linear Mathieu equation of the form x ̈+(ä+ϵ[cost+cosùt])x=0 The stability of solutions of Mathieu's equation as a function of parameter values (ä,ϵ) had been analyzed in this work. We used the Floquet type theory to generate stability diagrams which were used to determine the bounded regions of stability in the ä-ù plane for fixed ϵ. In the case of reducibility, we first applied the Lyapunov- Floquet (LF) transformation and modal transformation, which converted the linear part of the system into the Jordan form. Very importantly, quasi-periodic near-identity transformation was applied to reduce the system equations to a constant coefficient system by solving homological equations via harmonic balance. In this process we obtained the reducibility/resonance conditions that needed to be satisfied to convert a quasi-periodic system to a constant one.

It is possible in a properly controlled environment, such as industrial metrology, to make significant headway into the non-industrial constraints on image-based position measurement using the techniques of image registration and achieve repeatable feature measurements on the order of 0.3% of a pixel, or about an order of magnitude improvement on conventional real-world performance. These measurements are then used as inputs for a model optimal, model agnostic, smoothing for calibration of a laser scribe and online tracking of velocimeter using video input. Using appropriate smooth interpolation to increase effective sample density can reduce uncertainty and improve estimates. Use of the proper negative offset of the template function has the result of creating a convolution with higher local curvature than either template of target function which allows improved center-finding. Using the Akaike Information Criterion with a smoothing spline function it is possible to perform a model-optimal smooth on scalar measurements without knowing the underlying model and to determine the function describing the uncertainty in that optimal smooth. An example of empiric derivation of the parameters for a rudimentary Kalman Filter from this is then provided, and tested. Using the techniques of Exploratory Data Analysis and the "Formulize" genetic algorithm tool to convert the spline models into more accessible analytic forms resulted in stable, properly generalized, KF with performance and simplicity that exceeds "textbook" implementations thereof. Validation of the measurement includes that, in analytic case, it led to arbitrary precision in measurement of feature; in reasonable test case using the methods proposed, a reasonable and consistent maximum error of around 0.3% the length of a pixel was achieved and in practice using pixels that were 700nm in size feature position was located to within ± 2 nm. Robust applicability is demonstrated by the measurement of indicator position for a King model 2-32-G-042 rotameter.