Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

Shear band evolution in Zr/Hf-based bulk metallic glasses under static and dynamic indentations

Bulk metallic glasses (BMGs) exhibit superior mechanical properties as compared with other conventional materials and have been proposed for numerous engineering and technological applications. Zr/Hf-based BMGs or tungsten reinforced BMG composites are considered as a potential replacement for depleted uranium armor-piercing projectiles because of their ability to form localized shear bands during impact, which has been known to be the dominant plastic deformation mechanism in BMGs. However, in conventional tensile, compressive and bending tests, limited ductility has been observed because of fracture initiation immediately following the shear band formation. To fully investigate shear band characteristics, indentation tests that can confine the deformation in a limited region have been pursued. In this thesis, a detailed investigation of thermal stability and mechanical deformation behavior of Zr/Hf-based BMGs is conducted. First, systematic studies had been implemented to understand the influence of relative compositions of Zr and Hf on thermal stability and mechanical property evolution. Second, shear band evolution under indentations were investigated experimentally and theoretically. Three kinds of indentation studies were conducted on BMGs in the current study. (a) Nano-indentation to determine the mechanical properties as a function of Hf/Zr content. (b) Static Vickers indentation on bonded split specimens to investigate the shear band evolution characteristics beneath the indention. (c) Dynamic Vickers indentation on bonded split specimens to investigate the influence of strain rate. It was found in the present work that gradually replacing Zr by Hf remarkably increases the density and improves the mechanical properties. However, a slight decrease in glass forming ability with increasing Hf content has also been identified through thermodynamic analysis although all the materials in the current study were still found to be amorphous. Many indentation studies have revealed only a few shear bands surrounding the indent on the top surface of the specimen. This small number of shear bands cannot account for the large plastic deformation beneath the indentations. Therefore, a bonded interface technique has been used to observe the slip-steps due to shear band evolution. Vickers indentations were performed along the interface of the bonded split specimen at increasing loads. At small indentation loads, the plastic deformation was primarily accommodated by semi-circular primary shear bands surrounding the indentation. At higher loads, secondary and tertiary shear bands were formed inside this plastic zone. A modified expanding cavity model was then used to predict the plastic zone size characterized by the shear bands and to identify the stress components responsible for the evolution of the various types of shear bands. The applicability of various hardness—yield-strength ( H −σγ ) relationships currently available in the literature for bulk metallic glasses (BMGs) is also investigated. Experimental data generated on ZrHf-based BMGs in the current study and those available elsewhere on other BMG compositions were used to validate the models. A modified expanding-cavity model, employed in earlier work, was extended to propose a new H −σγ relationship. Unlike previous models, the proposed model takes into account not only the indenter geometry and the material properties, but also the pressure sensitivity index of the BMGs. The influence of various model parameters is systematically analyzed. It is shown that there is a good correlation between the model predictions and the experimental data for a wide range of BMG compositions. Under dynamic Vickers indentation, a decrease in indentation hardness at high loading rate was observed compared to static indentation hardness. It was observed that at equivalent loads, dynamic indentations produced more severe deformation features on the loading surface than static indentations. Different from static indentation, two sets of widely spaced semi-circular shear bands with two different curvatures were observed. The observed shear band pattern and the strain rate softening in indentation hardness were rationalized based on the variations in the normal stress on the slip plane, the strain rate of shear and the temperature rise associated with the indentation deformation. Finally, a coupled thermo-mechanical model is proposed that utilizes a momentum diffusion mechanism for the growth and evolution of the final spacing of shear bands. The influence of strain rate, confinement pressure and critical shear displacement on the shear band spacing, temperature rise within the shear band, and the associated variation in flow stress have been captured and analyzed. Consistent with the known pressure sensitive behavior of BMGs, the current model clearly captures the influence of the normal stress in the formation of shear bands. The normal stress not only reduces the time to reach critical shear displacement but also causes a significant temperature rise during the shear band formation. Based on this observation, the variation of shear band spacing in a typical dynamic indentation test has been rationalized. The temperature rise within a shear band can be in excess of 2000K at high strain rate and high confinement pressure conditions. The associated drop in viscosity and flow stress may explain the observed decrease in fracture strength and indentation hardness. The above investigations provide valuable insight into the deformation behavior of BMGs under static and dynamic loading conditions. The shear band patterns observed in the above indentation studies can be helpful to understand and model the deformation features under complex loading scenarios such as the interaction of a penetrator with armor. Future work encompasses (1) extending and modifying the coupled thermo-mechanical model to account for the temperature rise in quasistatic deformation; and (2) expanding this model to account for the microstructural variation-crystallization and free volume migration associated with the deformation.

Dynamic modeling, simulation and control design of a parafoil-payload system for ship launched aerial delivery system (SLADS)

The objective of this research was to develop a high-fidelity dynamic model of a parafoilpayload system with respect to its application for the Ship Launched Aerial Delivery System (SLADS). SLADS is a concept in which cargo can be transfered from ship to shore using a parafoil-payload system. It is accomplished in two phases: An initial towing phase when the glider follows the towing vessel in a passive lift mode and an autonomous gliding phase when the system is guided to the desired point. While many previous researchers have analyzed the parafoil-payload system when it is released from another airborne vehicle, limited work has been done in the area of towing up the system from ground or sea. One of the main contributions of this research was the development of a nonlinear dynamic model of a towed parafoil-payload system. After performing an extensive literature review of the existing methods of modeling a parafoil-payload system, a five degree-of-freedom model was developed. The inertial and geometric properties of the system were investigated to predict accurate results in the simulation environment. Since extensive research has been done in determining the aerodynamic characteristics of a paraglider, an existing aerodynamic model was chosen to incorporate the effects of air flow around the flexible paraglider wing. During the towing phase, it is essential that the parafoil-payload system follow the line of the towing vessel path to prevent an unstable flight condition called ‘lockout’. A detailed study of the causes of lockout, its mathematical representation and the flight conditions and the parameters related to lockout, constitute another contribution of this work. A linearized model of the parafoil-payload system was developed and used to analyze the stability of the system about equilibrium conditions. The relationship between the control surface inputs and the stability was investigated. In addition to stability of flight, one more important objective of SLADS is to tow up the parafoil-payload system as fast as possible. The tension in the tow cable is directly proportional to the rate of ascent of the parafoil-payload system. Lockout instability is more favorable when tow tensions are large. Thus there is a tradeoff between susceptibility to lockout and rapid deployment. Control strategies were also developed for optimal tow up and to maintain stability in the event of disturbances.

Time-domain models for power system stability and unbalance

It is an important and difficult challenge to protect modern interconnected power system from blackouts. Applying advanced power system protection techniques and increasing power system stability are ways to improve the reliability and security of power systems. Phasor-domain software packages such as Power System Simulator for Engineers (PSS/E) can be used to study large power systems but cannot be used for transient analysis. In order to observe both power system stability and transient behavior of the system during disturbances, modeling has to be done in the time-domain. This work focuses on modeling of power systems and various control systems in the Alternative Transients Program (ATP). ATP is a time-domain power system modeling software in which all the power system components can be modeled in detail. Models are implemented with attention to component representation and parameters. The synchronous machine model includes the saturation characteristics and control interface. Transient Analysis Control System is used to model the excitation control system, power system stabilizer and the turbine governor system of the synchronous machine. Several base cases of a single machine system are modeled and benchmarked against PSS/E. A two area system is modeled and inter-area and intra-area oscillations are observed. The two area system is reduced to a two machine system using reduced dynamic equivalencing. The original and the reduced systems are benchmarked against PSS/E. This work also includes the simulation of single-pole tripping using one of the base case models. Advantages of single-pole tripping and comparison of system behavior against three-pole tripping are studied. Results indicate that the built-in control system models in PSS/E can be effectively reproduced in ATP. The benchmarked models correctly simulate the power system dynamics. The successful implementation of a dynamically reduced system in ATP shows promise for studying a small sub-system of a large system without losing the dynamic behaviors. Other aspects such as relaying can be investigated using the benchmarked models. It is expected that this work will provide guidance in modeling different control systems for the synchronous machine and in representing dynamic equivalents of large power systems.

Thickness measurement of dynamic thin liquid films generated by plug-annular flow in non-wetting microchannels

Surface tension forces are significant at millimeter length-scales, causing profoundly different flow morphologies in microchannels than in macroscale flows. The existence and morphology of thin liquid films is particularly relevant for predicting performance and operational stability of devices containing microscale two phase flows. Analytical, computational, and experimental methods previously employed in the study of thin liquid films are discussed. Thicknesses before and after a novel film morphology, referred to as a `shock,' are measured with a novel film thickness measurement technique that uses confocal microscopy. Film thicknesses predicted by previous work are compared to experimental results. Methods for increasing the accuracy of the confocal film thickness measurement technique are discussed.