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.

EVALUATING INFRASTRUCTURE PERFORMANCE TO ASSIST FACILITY MANAGEMENT TOWARD SUSTAINABLE SYSTEMS

Civil infrastructure provides essential services for the development of both society and economy. It is very important to manage systems efficiently to ensure sound performance. However, there are challenges in information extraction from available data, which also necessitates the establishment of methodologies and frameworks to assist stakeholders in the decision making process. This research proposes methodologies to evaluate systems performance by maximizing the use of available information, in an effort to build and maintain sustainable systems. Under the guidance of problem formulation from a holistic view proposed by Mukherjee and Muga, this research specifically investigates problem solving methods that measure and analyze metrics to support decision making. Failures are inevitable in system management. A methodology is developed to describe arrival pattern of failures in order to assist engineers in failure rescues and budget prioritization especially when funding is limited. It reveals that blockage arrivals are not totally random. Smaller meaningful subsets show good random behavior. Additional overtime failure rate is analyzed by applying existing reliability models and non-parametric approaches. A scheme is further proposed to depict rates over the lifetime of a given facility system. Further analysis of sub-data sets is also performed with the discussion of context reduction. Infrastructure condition is another important indicator of systems performance. The challenges in predicting facility condition are the transition probability estimates and model sensitivity analysis. Methods are proposed to estimate transition probabilities by investigating long term behavior of the model and the relationship between transition rates and probabilities. To integrate heterogeneities, model sensitivity is performed for the application of non-homogeneous Markov chains model. Scenarios are investigated by assuming transition probabilities follow a Weibull regressed function and fall within an interval estimate. For each scenario, multiple cases are simulated using a Monte Carlo simulation. Results show that variations on the outputs are sensitive to the probability regression. While for the interval estimate, outputs have similar variations to the inputs. Life cycle cost analysis and life cycle assessment of a sewer system are performed comparing three different pipe types, which are reinforced concrete pipe (RCP) and non-reinforced concrete pipe (NRCP), and vitrified clay pipe (VCP). Life cycle cost analysis is performed for material extraction, construction and rehabilitation phases. In the rehabilitation phase, Markov chains model is applied in the support of rehabilitation strategy. In the life cycle assessment, the Economic Input-Output Life Cycle Assessment (EIO-LCA) tools are used in estimating environmental emissions for all three phases. Emissions are then compared quantitatively among alternatives to support decision making.

Effects of silicon content and cooling rate on mechanical properties of heavy section ductile cast iron

The effects of Si and cooling rate are investigated for their effect on the mechanical properties and microstructure. Three alloys were chosen with varying C and Si contents and an attempt to keep the remainder of the elements present constant. Within each heat, three test blocks were poured. Two blocks had chills – one with a fluid flowing through it to cool it (active chill) and one without the fluid (passive) – and the third block did not have a chill. Cooling curves were gathered and analyzed. The mechanical properties of the castings were correlated to the microstructure, cooling rate and Si content of each block. It was found that an increase in Si content increased the yield stress, tensile strength and hardness but decreased the impact toughness, elongation and Young’s modulus. The fast cooling rates produced by the chills caused a high nodule count in the castings along with a fine ferrite grain size and a high degree of nodularity. The fine microstructures, in turn, increased the strength and ductile to brittle transition temperature (DBTT) of the castings. The fast cooling rate was not adequate to overcome the dramatic increase in DBTT that is caused by the addition of Si.

Multimetallic complexes based on phosphine- and phosphine oxide- appended p-hydroquinones

Presented here, is the work done with a series of binucleating ligands based on phosphine and phosphine oxide appended p-hydroquinones and their reactions towards various metals sources. The long term goal of the project was to produce coordination polymers that would have novel electronic, magnetic, and optical properties which would be of use in the field of molecular electronics. Binucleating ligands contained a p-hydroquinone motif in which various phosphine- and phosphine oxide substituents have been placed in the ortho position relative to each of the hydroxy position were synthesized. A previously published synthetic method for such lugands utilized n-BuLi to form a phenyl lithium intermediate which was quenched with chlorodiphenylphosphine. This technique was also used to produce a ligand with diisopropylphosphine groups. Phosphine ligands, containing the same structural motif, were also generated using LDA as the lithiating agent. This technique was found to be higher yielding. Phosphine chalcogenide ligands were accessed by further oxidizing the low valent phosphorous centers with either hydrogen peroxide or with elemental sulfur. These ligands were characterized using multinuclear NMR, low and high resolution mass spectroscopy, FTIR, and single crystal X-ray diffraction. Their electrochemical properties were explored with cyclic voltammetry. The phosphine appended ligands were used in the synthesis of a several bimetallic complexes. It was found that the ligands readily reacted with NiCp2 and NiCp*2, displacing one of the cyclopentadiene (Cp) or pentamethylcyclopentadiene (Cp*) rings. A cyclopentadiene complexes, containing diisopropylphine, was readily oxidized by[FeCp2]PF6 to give a NMR silent mixed valence complex. Cyclic voltammetry of these complexes showed a number of reversible waves with a large potential separation. The mixed valence compounds also showed a large absorbance band in the NIR region which was assigned to be an intervalence charge transfer. The cyclic voltammetry and NIR spectroscopy suggest that these systems are very capable of efficient metal-to-metal charge transfer. These complexes were characterized by multinuclear NMR, single crystal X-ray diffraction, UV/VIS-NIR spectroscopy and elemental analysis. The phosphine oxide ligands were reacted with a variety of different metal sources but limited success was gained in obtaining single crystals, allowing structural characterization of these compounds. Single crystals were obtained from products generated by reacting the diphenylphosphine oxide ligand with (Bipy)Cu(NO3)2 and Cu(NO3)2. In all cases the ligand had been further oxidized to a 2,5-dihydroxy-1,4-benzoquinone motif. In the reaction between the diphenylphosphine oxide ligand and (Bipy)Cu(NO3)2 it was found that the phosphine oxide moiety was involved with intermolecular coordination leading to the formation of a one-dimensional polymer composed of a series of bimetallic complexes tethered together. When NaSbF6 was present in the reaction with (Bipy)Cu(NO3)2 a unique tetrametallic complex was formed. Here the phospine oxide moiety was oriented so that two bimetallic complexes were bound together. If only Cu(NO3)2 was present, a two-dimensional polymeric sheet was formed where the ligand was present in two different coordination modes. The electronic properties of these complexes remained to be assessed.

COMBUSTION INITIATION BY ELECTRICAL-DISCHARGE-INDUCED PLASMA IN LEAN AND DILUTE METHANE-AIR MIXTURE: EXPERIMENTAL AND MODELING INVESTIGATION

This dissertation represents experimental and numerical investigations of combustion initiation trigged by electrical-discharge-induced plasma within lean and dilute methane air mixture. This research topic is of interest due to its potential to further promote the understanding and prediction of spark ignition quality in high efficiency gasoline engines, which operate with lean and dilute fuel-air mixture. It is specified in this dissertation that the plasma to flame transition is the key process during the spark ignition event, yet it is also the most complicated and least understood procedure. Therefore the investigation is focused on the overlapped periods when plasma and flame both exists in the system. Experimental study is divided into two parts. Experiments in Part I focuses on the flame kernel resulting from the electrical discharge. A number of external factors are found to affect the growth of the flame kernel, resulting in complex correlations between discharge and flame kernel. Heat loss from the flame kernel to code ambient is found to be a dominant factor that quenches the flame kernel. Another experimental focus is on the plasma channel. Electrical discharges into gases induce intense and highly transient plasma. Detailed observation of the size and contents of the discharge-induced plasma channel is performed. Given the complex correlation and the multi-discipline physical/chemical processes involved in the plasma-flame transition, the modeling principle is taken to reproduce detailed transitions numerically with minimum analytical assumptions. Detailed measurement obtained from experimental work facilitates the more accurate description of initial reaction conditions. The novel and unique spark source considering both energy and species deposition is defined in a justified manner, which is the key feature of this Ignition by Plasma (IBP) model. The results of numerical simulation are intuitive and the potential of numerical simulation to better resolve the complex spark ignition mechanism is presented. Meanwhile, imperfections of the IBP model and numerical simulation have been specified and will address future attentions.

Spin-strain coupling in a 3-D transition metal oxides

ab-initio Hartree Fock (HF), density functional theory (DFT) and hybrid potentials were employed to compute the optimized lattice parameters and elastic properties of perovskite 3-d transition metal oxides. The optimized lattice parameters and elastic properties are interdependent in these materials. An interaction is observed between the electronic charge, spin and lattice degrees of freedom in 3-d transition metal oxides. The coupling between the electronic charge, spin and lattice structures originates due to localization of d-atomic orbitals. The coupling between the electronic charge, spin and crystalline lattice also contributes in the ferroelectric and ferromagnetic properties in perovskites. The cubic and tetragonal crystalline structures of perovskite transition metal oxides of ABO3 are studied. The electronic structure and the physics of 3-d perovskite materials is complex and less well considered. Moreover, the novelty of the electronic structure and properties of these perovskites transition metal oxides exceeds the challenge offered by their complex crystalline structures. To achieve the objective of understanding the structure and property relationship of these materials the first-principle computational method is employed. CRYSTAL09 code is employed for computing crystalline structure, elastic, ferromagnetic and other electronic properties. Second-order elastic constants (SOEC) and bulk moduli (B) are computed in an automated process by employing ELASTCON (elastic constants) and EOS (equation of state) programs in CRYSTAL09 code. ELASTCON, EOS and other computational algorithms are utilized to determine the elastic properties of tetragonal BaTiO3, rutile TiO2, cubic and tetragonal BaFeO3 and the ferromagentic properties of 3-d transition metal oxides. Multiple methods are employed to crosscheck the consistency of our computational results. Computational results have motivated us to explore the ferromagnetic properties of 3-d transition metal oxides. Billyscript and CRYSTAL09 code are employed to compute the optimized geometry of the cubic and tetragonal crystalline structure of transition metal oxides of Sc to Cu. Cubic crystalline structure is initially chosen to determine the effect of lattice strains on ferromagnetism due to the spin angular momentum of an electron. The 3-d transition metals and their oxides are challenging as the basis functions and potentials are not fully developed to address the complex physics of the transition metals. Moreover, perovskite crystalline structures are extremely challenging with respect to the quality of computations as the latter requires the well established methods. Ferroelectric and ferromagnetic properties of bulk, surfaces and interfaces are explored by employing CRYSTAL09 code. In our computations done on cubic TMOs of Sc-Fe it is observed that there is a coupling between the crystalline structure and FM/AFM spin polarization. Strained crystalline structures of 3-d transition metal oxides are subjected to changes in the electromagnetic and electronic properties. The electronic structure and properties of bulk, composites, surfaces of 3-d transition metal oxides are computed successfully.