Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization.

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.

Synthesis, structure, and stability of metal-organic frameworks

Metal-organic frameworks (MOFs) obtained much attention because of their unusual structures and properties as well as their potential applications. This dissertation research was focused on (1) the effects of synthesis conditions on the structures of MOFs, (2) the thermal stability of MOFs, (3) pressure-induced amorphization, and (4) the effect of high-valent ions on the structure of a MOF. This research demonstrated that the crystal structure of MOF-5 could be controlled by drying solvents. If the vacuum solvent is dimethylformamide (DMF), the crystal structure of MOF-5 is tetragonal. In contrast, if the DMF is displaced by CH2Cl2 before the vacuum, the obtained MOF-5 occupies a cubic structure. Furthermore, it was found that the tetragonal MOF-5 exhibited a mediate surface area (300-1000 m2/g). The surface area of tetragonal MOF-5 is also dependent on Zn(NO3)2/H2BDC (H2BDC: terephthalic acid) molar ratios used for its synthesis. The optimum ratio is 1.38, at which synthesized tetragonal MOF-5 exhibits the highest crystallinity and surface area (1297 m2/g). The thermal stability and decomposition of MOF-5 were systematically investigated. The thermal decomposition of cubic and tetragonal MOF-5s resulted in the same products: CO2, benzene, amorphous carbon, and crystal ZnO. The thermal decomposition is due to breaking carboxylic bridges between benzene rings and Zn4O clusters. Identifying structural relationships between crystalline and noncrystalline states is of fundamental interest in materials research. Currently, amorphization of solid materials at ambient temperature requires an ultra-high pressure (several GPa). However, this research demonstrated that MOF-5 and IRMOF-8 can be irreversibly amorphized at ambient temperature by employing a low compressing pressure of 3.5 MPa, which is 100 times lower than that required for amorphization of other solids. Furthermore, the pressure-induced amorphization (PIA) of MOFs is strongly dependent on the changeability of bond angles. If the geometric structure of a MOF can allow bond angles to be changed without breaking bonds, it can easily be amorphized by compression. This can explain why MOF-5 and IRMOF-8 can easily be amorphized via compression than Cu-BTC. It is generally recognized that zeolitic imidazolate frameworks (ZIFs) occupy much higher stability than other types of MOFs. The representative of ZIFs is Zn(2-methylimidazole)2 (ZIF-8) exhibiting high-decomposition temperature and high chemical resistance to various solvents. However, so far, it is still unknown whether the high stability of ZIF-8 can be challenged by ions, which is important for its modification by doping ions. In this research, we performed aqueous salt solution treatment on ZIF-8, and the results showed that anions (Cl¯ and NO3¯) in a solution exhibited no effect on the crystal structure of ZIF-8. However, the effect of cations (in a solution) on structure of ZIF-8 strongly depends on the cation valences. The univalent metal cations showed no effect on the structure of ZIF-8, whereas the bivalent or higher-valent metal cations caused the collapse of ZIF-8 crystal structure. Therefore, structure stability of ZIF-8 is considered when it is subjected to the application, in which high-valent metal cations are involved.

STUDY OF SPARK IGNITION ENGINE COMBUSTION MODEL FOR THE ANALYSIS OF CYCLIC VARIATION AND COMBUSTION STABILITY AT LEAN OPERATING CONDITIONS

A fundamental combustion model for spark-ignition engine is studied in this report. The model is implemented in SIMULINK to simulate engine outputs (mass fraction burn and in-cylinder pressure) under various engine operation conditions. The combustion model includes a turbulent propagation and eddy burning processes based on literature [1]. The turbulence propagation and eddy burning processes are simulated by zero-dimensional method and the flame is assumed as sphere. To predict pressure, temperature and other in-cylinder variables, a two-zone thermodynamic model is used. The predicted results of this model match well with the engine test data under various engine speeds, loads, spark ignition timings and air fuel mass ratios. The developed model is used to study cyclic variation and combustion stability at lean (or diluted) combustion conditions. Several variation sources are introduced into the combustion model to simulate engine performance observed in experimental data. The relations between combustion stability and the introduced variation amount are analyzed at various lean combustion levels.

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.

Patch selection by wintering white-tailed deer : consequences for eastern hemlock regeneration at local and landscape scales

Tsuga canadensis (eastern hemlock) is a highly shade-tolerant, late-successional, and long-lived conifer species found throughout eastern North America. It is most often found in pure or nearly pure stands, because highly acidic and nutrient poor forest floor conditions are thought to favor T. canadensis regeneration while simultaneously limiting the establishment of some hardwood species with greater nutrient requirements. Once a common species, T. canadensis is currently experiencing widescale declines across its range. The hemlock woolly adelgid (Adelges tsugae) is decimating the population across its eastern distribution. Across the Upper Great Lakes region, where the adelgid is currently being held at bay by cold winter temperatures, T. canadensis has been experiencing failures in regeneration attributed, in part, to herbivory by white-tailed deer (Odocoileus virginianus). Deer utilize T. canadensis stands as winter habitat in areas of high snow depth. Tsuga canadensis, once a major component of these forests, currently exists at just a fraction of its pre-settlement abundance due to historic logging and contemporary forest management practices, and what remains is found in small remnant patches surrounded by second- and third-growth deciduous forests. The deer population across the region, however, is likely double that of pre-European settlement times. In this dissertation I explore the relationship between white-tailed deer use of T. canadensis as winter habitat and the effect this use is having on regeneration and forest succession. For this research I quantified stand composition and structure and abiotic variables of elevation and snow depth in 39 randomly selected T. canadensis stands from across the western Upper Peninsula of Michigan. I also quantified composition and the configuration of the landscapes surrounding these stands. I measured relative deer use of T. canadensis stands as pellet group piles deposited in each stand during each of three consecutive winters, 2005-06, 2006-07, and 2007-08. The results of this research suggest that deer use of T. canadensis stands as winter habitat is influenced primarily by snow depth, elevation, and the composition and configuration of the greater landscapes surrounding these stands. Specifically, stands with more heterogeneous landscapes surrounding them (i.e., a patchy mosaic of conifer, deciduous, and open cover) had higher relative deer use than stands surrounded by homogenous deciduous forest cover. Additionally, the intensity of use and the number of stands used was greater in years with higher average snow depth. Tsuga canadensis regeneration in these stands was negatively associated with deer use and Acer saccharum (sugar maple) basal area. Of the 39 stands, 17 and 22 stands had no T. canadensis regeneration in small and large sapling categories, respectively. Acer saccharum was the most common understory tree species, and the importance of A. saccharum in the understory (stems < 10 cm dbh) of the stands was positively associated with overstory A. saccharum dominance. Tsuga canadensis establishment was associated with high-decay coarse woody debris and moss, and deciduous leaf litter inputs in these stands may be limiting access to these important microsites. Furthermore, A. saccharum is more tolerant to the effects of deer herbivory than T. canadensis, giving A. saccharum a competitive advantage in stands being utilized as winter habitat by deer. My research suggests that limited microsite availability, in conjunction with deer herbivory, may be leading to an erosion in T. canadensis patch stability and an altered successional trajectory toward one of A. saccharum dominance, an alternately stable climax species.

PMMA-CNT matrices for vacuum electronic, biosensing and energy applications

Carbon nanotubes (CNTs) are interesting materials with extraordinary properties for various applications. Here, vertically-aligned multiwalled CNTs (VA-MWCNTs) are grown by our dual radio frequency plasma enhanced chemical vapor deposition (PECVD). After optimizing the synthesis processes, these VA-MWCNTs were fabricated in to a series of devices for applications in vacuum electronics, glucose biosensors, glucose biofuel cells, and supercapacitors In particular, we have created the so-called PMMA-CNT matrices (opened-tip CNTs embedded in poly-methyl methacrylate) that are promising components in a novel energy sensing, generation and storage (SGS) system that integrate glucose biosensors, biofuel cells, and supercapacitors. The content of this thesis work is described as follows: 1. We have first optimized the synthesis of VA-MWCNTs by our PECVD technique. The effects of CH4 flow rate and growth duration on the lengths of these CNTs were studied. 2. We have characterized these VA-MWCNTs for electron field emission. We noticed that as grown CNTs suffers from high emission threshold, poor emission density and poor long-term stability. We attempted a series of experiments to understand ways to overcome these problems. First, we decrease the screening effects on VA-MWCNTs by creating arrays of self-assembled CNT bundles that are catalyst-free and opened tips. These bundles are found to enhance the field emission stability and emission density. Subsequently, we have created PMMA-CNT matrices that are excellent electron field emitters with an emission threshold field of more than two-fold lower than that of the as-grown sample. Furthermore, no significant emission degradation was observed after a continuous emission test of 40 hours (versus much shorter tests in reported literatures). Based on the new understanding we learnt from the PMMA-CNT matrices, we further created PMMA-STO-CNT matrices by embedding opened-tip VA-MWCNTs that are coated with strontium titanate (SrTiO3) with PMMA. We found that the PMMA-STO-CNT matrices have all the desired properties of the PMMA-CNT matrices. Furthermore, PMMA-STO-CNT matrices offer much lower emission threshold field, about five-fold lower than that of as grown VA-MWCNTs. The new understandings we obtained are important for practical application of VA-MWCNTs in field emission devices. 3. Subsequently, we have functionalized PMMA-CNT matrices for glucose biosensing. Our biosensor was developed by immobilized glucose oxidase (GOχ) on the opened-tip CNTs exposed on the matrices. The durability, stability and sensitivity of the biosensor were studied. In order to understand the performance of miniaturized glucose biosensors, we have then investigated the effect of working electrode area on the sensitivity and current level of our biosensors. 4. Next, functionalized PMMA-CNT matrices were utilized for energy generation and storage. We found that PMMA-CNT matrices are promising component in glucose/O2 biofuel cells (BFCs) for energy generation. The construction of these BFCs and the effect of the electrode area on the power density of these BFCs were investigated. Then, we have attempted to use PMMA-CNT matrices as supercapacitors for energy storage devices. The performance of these supercapacitors and ways to enhance their performance are discussed. 5. Finally, we further evaluated the concept of energy SGS system that integrated glucose biosensors, biofuel cells, and supercapacitors. This SGS system may be implantable to monitor and control the blood glucose level in our body.

Influence of Mechanical and Thermal Boundary Conditions on Stabilizing/Destabilizing Mechanisms in Evaporating Liquid Films

Liquid films, evaporating or non-evaporating, are ubiquitous in nature and technology. The dynamics of evaporating liquid films is a study applicable in several industries such as water recovery, heat exchangers, crystal growth, drug design etc. The theory describing the dynamics of liquid films crosses several fields such as engineering, mathematics, material science, biophysics and volcanology to name a few. Interfacial instabilities typically manifest by the undulation of an interface from a presumed flat state or by the onset of a secondary flow state from a primary quiescent state or both. To study the instabilities affecting liquid films, an evaporating/non-evaporating Newtonian liquid film is subject to a perturbation. Numerical analysis is conducted on configurations of such liquid films being heated on solid surfaces in order to examine the various stabilizing and destabilizing mechanisms that can cause the formation of different convective structures. These convective structures have implications towards heat transfer that occurs via this process. Certain aspects of this research topic have not received attention, as will be obvious from the literature review. Static, horizontal liquid films on solid surfaces are examined for their resistance to long wave type instabilities via linear stability analysis, method of normal modes and finite difference methods. The spatiotemporal evolution equation, available in literature, describing the time evolution of a liquid film heated on a solid surface, is utilized to analyze various stabilizing/destabilizing mechanisms affecting evaporating and non-evaporating liquid films. The impact of these mechanisms on the film stability and structure for both buoyant and non-buoyant films will be examined by the variation of mechanical and thermal boundary conditions. Films evaporating in zero gravity are studied using the evolution equation. It is found that films that are stable to long wave type instabilities in terrestrial gravity are prone to destabilization via long wave instabilities in zero gravity.