A major concern in the operation of present-day gas turbine engines is the ingestion of hot mainstream gas into rotor-stator disk cavities of the high-pressure turbine stages. Although the engines require high gas temperature at turbine entry for good performance efficiency, the ingested gas shortens the lives of the cavity internals, particularly that of the rotor disks. Steps such as installing seals at the disk rims and injecting purge (secondary) air bled from the compressor discharge into the cavities are implemented to reduce the gas ingestion. Although there are advantages to the above-mentioned steps, the performance of a gas turbine engine is diminished by the purge air bleed-off. This then requires that the cavity sealing function be achieved with as low a purge air supply rate as possible. This, in turn, renders imperative an in-depth understanding of the pressure and velocity fields in the main gas path and within the disk cavities. In this work, experiments were carried out in a model 1.5-stage (stator-rotor-stator) axial air turbine to study the ingestion of main air into the aft, rotor-stator, disk cavity. The cavity featured rotor and stator rim seals with radial clearance and axial overlap and an inner labyrinth seal. First, time-average static pressure distribution was measured in the main gas path upstream and downstream of the rotor as well as in the cavity to ensure that a nominally steady run condition had been achieved. Main gas ingestion was determined by measuring the concentration distribution of tracer gas (CO2) in the cavity. To map the cavity fluid velocity field, particle image velocimetry was employed. Results are reported for two main air flow rates, two rotor speeds, and four purge air flow rates.

Hydrogen embrittlement (HE) is a phenomenon that affects both the physical and chemical properties of several intrinsically ductile metals. Consequently, understanding the mechanisms behind HE has been of particular interest in both experimental and modeling research. Discrepancies between experimental observations and modeling results have led to various proposals for HE mechanisms. Therefore, to gain insights into HE mechanisms in iron, this dissertation aims to investigate several key issues involving HE such as: a) the incipient crack tip events; b) the cohesive strength of grain boundaries (GBs); c) the dislocation-GB interactions and d) the dislocation mobility.

The crack tip, which presents a preferential trap site for hydrogen segregation, was examined using atomistic methods and the continuum based Rice-Thompson criterion as sufficient concentration of hydrogen can alter the crack tip deformation mechanism. Results suggest that there is a plausible co-existence of the adsorption induced dislocation emission and hydrogen enhanced decohesion mechanisms. In the case of GB-hydrogen interaction, we observed that the segregation of hydrogen along the interface leads to a reduction in cohesive strength resulting in intergranular failure. A methodology was further developed to quantify the role of the GB structure on this behavior.

GBs play a fundamental role in determining the strengthening mechanisms acting as an impediment to the dislocation motion; however, the presence of an unsurmountable barrier for a dislocation can generate slip localization that could further lead to intergranular crack initiation. It was found that the presence of hydrogen increases the strain energy stored within the GB which could lead to a transition in failure mode. Finally, in the case of body centered cubic metals, understanding the complex screw dislocation motion is critical to the development of an accurate continuum description of the plastic behavior. Further, the presence of hydrogen has been shown to drastically alter the plastic deformation, but the precise role of hydrogen is still unclear. Thus, the role of hydrogen on the dislocation mobility was examined using density functional theory and atomistic simulations. Overall, this dissertation provides a novel atomic-scale understanding of the HE mechanism and development of multiscale tools for future endeavors.

Interest in Micro Aerial Vehicle (MAV) research has surged over the past decade. MAVs offer new capabilities for intelligence gathering, reconnaissance, site mapping, communications, search and rescue, etc. This thesis discusses key modeling and control aspects of flapping wing MAVs in hover. A three degree of freedom nonlinear model is used to describe the flapping wing vehicle. Averaging theory is used to obtain a nonlinear average model. The equilibrium of this model is then analyzed. A linear model is then obtained to describe the vehicle near hover. LQR is used to as the main control system design methodology. It is used, together with a nonlinear parameter optimization algorithm, to design a family multivariable control system for the MAV. Critical performance trade-offs are illuminated. Properties at both the plant output and input are examined. Very specific rules of thumb are given for control system design. The conservatism of the rules are also discussed. Issues addressed include

What should the control system bandwidth be vis--vis the flapping frequency (so that averaging the nonlinear system is valid)?

When is first order averaging sufficient? When is higher order averaging necessary?

When can wing mass be neglected and when does wing mass become critical to model?

This includes how and when the rules given can be tightened; i.e. made less conservative.

As robots are increasingly migrating out of factories and research laboratories and into our everyday lives, they should move and act in environments designed for humans. For this reason, the need of anthropomorphic movements is of utmost importance. The objective of this thesis is to solve the inverse kinematics problem of redundant robot arms that results to anthropomorphic configurations. The swivel angle of the elbow was used as a human arm motion parameter for the robot arm to mimic. The swivel angle is defined as the rotation angle of the plane defined by the upper and lower arm around a virtual axis that connects the shoulder and wrist joints. Using kinematic data recorded from human subjects during every-day life tasks, the linear sensorimotor transformation model was validated and used to estimate the swivel angle, given the desired end-effector position. Defining the desired swivel angle simplifies the kinematic redundancy of the robot arm. The proposed method was tested with an anthropomorphic redundant robot arm and the computed motion profiles were compared to the ones of the human subjects. This thesis shows that the method computes anthropomorphic configurations for the robot arm, even if the robot arm has different link lengths than the human arm and starts its motion at random configurations.

In order to achieve higher gas turbine efficiency, the main gas temperature at turbine inlet has been steadily increased from approximately 900°C to about 1500°C over the last few decades. This temperature is higher than the maximum acceptable temperature for turbine internals. The hot main gas may get ingested into the space between rotor and stator, the rotor-stator disk cavity in a stage because of the pressure differential between main gas annulus and the disk cavity. To reduce this ingestion, the disk cavity is equipped with a rim seal; additionally, secondary (purge) air is supplied to the cavity. Since the purge air is typically bled off the compressor discharge, this reducing the overall gas turbine efficiency, much research has been carried out to estimate the minimum purge flow necessary (c_w,min) for complete sealing of disk cavities.

In this work, experiments have been performed in a subscale single-stage axial turbine featuring vanes, blades and an axially-overlapping radial-clearance seal at the disk cavity rim. The turbine stage is also equipped with a labyrinth seal radially inboard. The stage geometry and the experimental conditions were such that the ingestion into the disk cavity was driven by the pressure asymmetry in the main gas annulus. In the experiments, time-averaged static pressure was measured at several locations in the main annulus and in the disk cavity; the pressure differential between a location on the vane platform close to lip (this being the rim seal part on the stator) and a location in the 'seal region' in the cavity is considered to be the driving potential for both ingestion and egress. Time-averaged volumetric concentration of the tracer gas (CO₂) in the purge air supplied was measured at multiple radial locations on the stator surface. The pressure and ingestion data were then used to calculate the ingestion and egress discharge coefficients for a range of purge flow rates, employing a simple orifice model of the rim seal. For the experiments performed, the egress discharge coefficient increased and the ingestion discharge coefficient decreased with the purge air flow rate. A method for estimation of c_w,min is also proposed.

Origami and kirigami, the technique of generating three-dimensional (3D) structures from two-dimensional (2D) flat sheets, are now more and more involved in scientific and engineering fields. Therefore, the development of tools for their theoretical analysis becomes more and more important. Since much effort was paid on calculations based on pure mathematical consideration and only limited effort has been paid to include mechanical properties, the goal of my research is developing a method to analyze the mechanical behavior of origami and kirigami based structures. Mechanical characteristics, including nonlocal effect and fracture of the structures, as well as elasticity and plasticity of materials are studied. For calculation of relative simple structures and building of structures’ constitutive relations, analytical approaches were used. For more complex structures, finite element analysis (FEA), which is commonly applied as a numerical method for the analysis of solid structures, was utilized. The general study approach is not necessarily related to characteristic size of model. I believe the scale-independent method described here will pave a new way to understand the mechanical response of a variety of origami and kirigami based structures under given mechanical loading.

The stability of nanocrystalline microstructural features allows structural materials to be synthesized and tested in ways that have heretofore been pursued only on a limited basis, especially under dynamic loading combined with temperature effects. Thus, a recently developed, stable nanocrystalline alloy is analyzed here for quasi-static (<100 s-1) and dynamic loading (103 to 104 s-1) under uniaxial compression and tension at multiple temperatures ranging from 298-1073 K. After mechanical tests, microstructures are analyzed and possible deformation mechanisms are proposed. Following this, strain and strain rate history effects on mechanical behavior are analyzed using a combination of quasi-static and dynamic strain rate Bauschinger testing. The stable nanocrystalline material is found to exhibit limited flow stress increase with increasing strain rate as compared to that of both pure, coarse grained and nanocrystalline Cu. Further, the material microstructural features, which includes Ta nano-dispersions, is seen to pin dislocation at quasi-static strain rates, but the deformation becomes dominated by twin nucleation at high strain rates. These twins are pinned from further growth past nucleation by the Ta nano-dispersions. Testing of thermal and load history effects on the mechanical behavior reveals that when thermal energy is increased beyond 200 °C, an upturn in flow stress is present at strain rates below 104 s-1. However, in this study, this simple assumption, established 50-years ago, is shown to break-down when the average grain size and microstructural length-scale is decreased and stabilized below 100nm. This divergent strain-rate behavior is attributed to a unique microstructure that alters slip-processes and their interactions with phonons; thus enabling materials response with a constant flow-stress even at extreme conditions. Hence, the present study provides a pathway for designing and synthesizing a new-level of tough and high-energy absorbing materials.

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.

High Pressure Superheater 1 (HPSH1) is the first heat exchange tube bank inside the Heat Recovery Steam Generator (HRSG) to encounter exhaust flue gas from the gas turbine of a Combined Cycle Power Plant. Steam flowing through the HPSH1 gains heat from the flue gas prior to entering the steam turbine. During cold start-ups, rapid temperature changes in operating condition give rise to significant temperature gradients in the thick-walled components of HPSH1 (manifolds, links, and headers). These temperature gradients produce thermal-structural stresses in the components. The resulting high cycle fatigue is a major concern as this can lead to premature failure of the components. The main objective of this project was to address the thermal-structural stress field induced in HPSH1 during a typical cold start-up transient. To this end, computational fluid dynamics (CFD) was used to carry out the thermal-fluid analysis of HPSH1. The calculated temperature distributions in the component walls were the primary inputs for the finite element (FEA) model that performed structural analysis. Thermal-structural analysis was initially carried out at full-load steady state condition in order to gain confidence in the CFD and FEA methodologies. Results of the full-load steady state thermal-fluid analysis were found in agreement with the temperature values measured at specific locations on the outer surfaces of the inlet links and outlet manifold. It was found from the subsequent structural analysis that peak effective stresses were located at the connecting regions of the components and were well below the allowed stress values. Higher temperature differences were observed between the thick-walled HPSH1 components during the cold start-up transient as compared to the full-load steady state operating condition. This was because of the rapid temperature changes that occurred, especially in the steam temperature at the HPSH1 entry, and the different rates of heating or cooling for components with different wall thicknesses. Results of the transient thermal-fluid analysis will be used in future to perform structural analysis of the HPSH1. The developed CFD and FEA models are capable of analyzing various other transients (e.g., hot start-up and shut-down) and determine their influence on the durability of plant components.

Advanced composites are being widely used in aerospace applications due to their high stiffness, strength and energy absorption capabilities. However, the assurance of structural reliability is a critical issue because a damage event will compromise the integrity of composite structures and lead to ultimate failure. In this dissertation a novel homogenization based multiscale modeling framework using semi-analytical micromechanics is presented to simulate the response of textile composites. The novelty of this approach lies in the three scale homogenization/localization framework bridging between the constituent (micro), the fiber tow scale (meso), weave scale (macro), and the global response. The multiscale framework, named Multiscale Generalized Method of Cells (MSGMC), continuously bridges between the micro to the global scale as opposed to approaches that are top-down and bottom-up. This framework is fully generalized and capable of modeling several different weave and braids without reformulation. Particular emphasis in this dissertation is placed on modeling the nonlinearity and failure of both polymer matrix and ceramic matrix composites.