Ingestion of high temperature mainstream gas into the rotor-stator cavities of a gas turbine is one of the major problems faced by the turbine designers. The ingested gas heats up rotor disks and induces higher thermal stresses on them, giving rise to durability concern. Ingestion is usually reduced by installing seals on the rotor and stator rims and by purging the disk cavity by secondary air bled from the compressor discharge. The geometry of the rim seals and the secondary air flow rate, together, influence the amount of gas that gets ingested into the cavities. Since the amount of secondary air bled off has a negative effect on the gas turbine thermal efficiency, one goal is to use the least possible amount of secondary air. This requires a good understanding of the flow and ingestion fields within a disk cavity. In the present study, the mainstream gas ingestion phenomenon has been experimentally studied in a model single-stage axial flow gas turbine. The turbine stage featured vanes and blades, and rim seals on both the rotor and stator. Additionally, the disk cavity contained a labyrinth seal radially inboard which effectively divided the cavity into a rim cavity and an inner cavity. Time-average static pressure measurements were obtained at various radial positions within the disk cavity, and in the mainstream gas path at three axial locations at the outer shroud spread circumferentially over two vane pitches. The time-average static pressure in the main gas path exhibited a periodic asymmetry following the vane pitch whose amplitude diminished with increasing distance from the vane trailing edge. The static pressure distribution increased with the secondary air flow rate within the inner cavity but was found to be almost independent of it in the rim cavity. Tracer gas (CO2) concentration measurements were conducted to determine the sealing effectiveness of the rim seals against main gas ingestion. For the rim cavity, the sealing effectiveness increased with the secondary air flow rate. Within the inner cavity however, this trend reversed -this may have been due to the presence of rotating low-pressure flow structures inboard of the labyrinth seal.

Shock loading is a complex phenomenon that can lead to failure mechanisms such as strain localization, void nucleation and growth, and eventually spall fracture. Studying incipient stages of spall damage is of paramount importance to accurately determine initiation sites in the material microstructure where damage will nucleate and grow and to formulate continuum models that account for the variability of the damage process due to microstructural heterogeneity. The length scale of damage with respect to that of the surrounding microstructure has proven to be a key aspect in determining sites of failure initiation. Correlations have been found between the damage sites and the surrounding microstructure to determine the preferred sites of spall damage, since it tends to localize at and around the regions of intrinsic defects such as grain boundaries and triple points. However, considerable amount of work still has to be done in this regard to determine the physics driving the damage at these intrinsic weak sites in the microstructure. The main focus of this research work is to understand the physical mechanisms behind the damage localization at these preferred sites. A crystal plasticity constitutive model is implemented with different damage criteria to study the effects of stress concentration and strain localization at the grain boundaries. A cohesive zone modeling technique is used to include the intrinsic strength of the grain boundaries in the simulations. The constitutive model is verified using single elements tests, calibrated using single crystal impact experiments and validated using bicrystal and multicrystal impact experiments. The results indicate that strain localization is the predominant driving force for damage initiation and evolution. The microstructural effects on theses damage sites are studied to attribute the extent of damage to microstructural features such as grain orientation, misorientation, Taylor factor and the grain boundary planes. The finite element simulations show good correlation with the experimental results and can be used as the preliminary step in developing accurate probabilistic models for damage nucleation.

The focus of this investigation includes three aspects. First, the development of nonlinear reduced order modeling techniques for the prediction of the response of complex structures exhibiting "large" deformations, i.e. a geometrically nonlinear behavior, and modeled within a commercial finite element code. The present investigation builds on a general methodology, successfully validated in recent years on simpler panel structures, by developing a novel identification strategy of the reduced order model parameters, that enables the consideration of the large number of modes needed for complex structures, and by extending an automatic strategy for the selection of the basis functions used to represent accurately the displacement field. These novel developments are successfully validated on the nonlinear static and dynamic responses of a 9-bay panel structure modeled within Nastran. In addition, a multi-scale approach based on Component Mode Synthesis methods is explored. Second, an assessment of the predictive capabilities of nonlinear reduced order models for the prediction of the large displacement and stress fields of panels that have a geometric discontinuity; a flat panel with a notch was used for this assessment. It is demonstrated that the reduced order models of both virgin and notched panels provide a close match of the displacement field obtained from full finite element analyses of the notched panel for moderately large static and dynamic responses. In regards to stresses, it is found that the notched panel reduced order model leads to a close prediction of the stress distribution obtained on the notched panel as computed by the finite element model. Two enrichment techniques, based on superposition of the notch effects on the virgin panel stress field, are proposed to permit a close prediction of the stress distribution of the notched panel from the reduced order model of the virgin one. A very good prediction of the full finite element results is achieved with both enrichments for static and dynamic responses. Finally, computational challenges associated with the solution of the reduced order model equations are discussed. Two alternatives to reduce the computational time for the solution of these problems are explored.

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.

This dissertation presents methods for addressing research problems that currently can only adequately be solved using Quality Reliability Engineering (QRE) approaches especially accelerated life testing (ALT) of electronic printed wiring boards with applications to avionics circuit boards. The methods presented in this research are generally applicable to circuit boards, but the data generated and their analysis is for high performance avionics. Avionics equipment typically requires 20 years expected life by aircraft equipment manufacturers and therefore ALT is the only practical way of performing life test estimates. Both thermal and vibration ALT induced failure are performed and analyzed to resolve industry questions relating to the introduction of lead-free solder product and processes into high reliability avionics. In chapter 2, thermal ALT using an industry standard failure machine implementing Interconnect Stress Test (IST) that simulates circuit board life data is compared to real production failure data by likelihood ratio tests to arrive at a mechanical theory. This mechanical theory results in a statistically equivalent energy bound such that failure distributions below a specific energy level are considered to be from the same distribution thus allowing testers to quantify parameter setting in IST prior to life testing. In chapter 3, vibration ALT comparing tin-lead and lead-free circuit board solder designs involves the use of the likelihood ratio (LR) test to assess both complete failure data and S-N curves to present methods for analyzing data. Failure data is analyzed using Regression and two-way analysis of variance (ANOVA) and reconciled with the LR test results that indicating that a costly aging pre-process may be eliminated in certain cases. In chapter 4, vibration ALT for side-by-side tin-lead and lead-free solder black box designs are life tested. Commercial models from strain data do not exist at the low levels associated with life testing and need to be developed because testing performed and presented here indicate that both tin-lead and lead-free solders are similar. In addition, earlier failures due to vibration like connector failure modes will occur before solder interconnect failures.

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.

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.

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.