Implementation of the conjugate heat transfer code in KIVA-4

KIVA is an open Computational Fluid Dynamics (CFD) source code that is capable to compute the transient two and three-dimensional chemically reactive fluid flows with spray. The latest version in the family of KIVA codes is the KIVA-4 which is capable of handling the unstructured mesh. This project focuses on the implementation of the Conjugate Heat Transfer code (CHT) in KIVA-4. The previous version of KIVA code with conjugate heat transfer code has been developed at Michigan Technological University by Egel Urip and is be used in this project. During the first phase of the project, the difference in the code structure between the previous version of KIVA and the KIVA-4 has been studied, which is the most challenging part of the project. The second phase involves the reverse engineering where the CHT code in previous version is extracted and implemented in KIVA-4 according to the new code structure. The validation of the implemented code is performed using a 4-valve Pentroof engine case. A solid cylinder wall has been developed using GRIDGEN which surrounds 3/4th of the engine cylinder and heat transfer to the solid wall during one engine cycle (0-720 Crank Angle Degree) is compared with that of the reference result. The reference results are nothing but the same engine case run in the previous version with the original code developed by Egel. The results of current code are very much comparable to that of the reference results which verifies that successful implementation of the CHT code in KIVA-4.

Sulfoxides as an intramolecular sulfenylating agent for indoles and diverse applications of the sulfide-sulfoxide redox cycle in organic chemistry

This dissertation involves study of various aspects of sulfoxide chemistry. Specifically designed t-butyl and propanenitrile sulfoxides tethered to indole-2-carboxamide were used as a source of intramolecular sulfenylating agents to synthesize novel indolo[3,2-b]-1-5-benzothiazepinones which are structurally analogous to the other biologically active benzothiazepinones. This study reveals that the intramolecular cyclization of sulfoxide follows an electrophilic sulfenylation (Sulfoxide Electrophilic Sulfenylation, SES) reaction pathway. Evidence of the absence of sulfenic acid as a transient reactive intermediate in such intramolecular cyclization is also provided. In another study, sulfoxide was used as a “protecting group” of thioether to synthesize 8-membered, indole substituted, thiazocine-2-acetic acid derivative via Ring Closing Metathesis (RCM). Protection (oxidation) of inert (to RCM) sulfide to sulfoxide followed by RCM produced cyclized product in good yields. Deprotection (reduction) of sulfoxide was achieved using Lawessons Reagent (L.R.). Application of the sulfide-sulfoxide redox cycle to solve the existing difficulties in using RCM methodology to thioethers is illustrated. A new design of a “molecular brake”, based on the sulfide-sulfoxide redox cycle is described. N-Ar rotation in simple isoindolines is controlled by the oxidation state of the proximate sulfur atom. Sulfide [S(II)] shows “free” [brake OFF] N-Ar rotation whereas sulfoxide displayed hindered [brake ON] N-Ar rotation. The semi-empirical molecular orbital (PM3) calculations revealed concerted pyramidalization of amidic nitrogen with N-Ar rotation.

Numerical modeling of the failure mechanisms in Si thin film anode for Li-ion batteries

In recent times, the demand for the storage of electrical energy has grown rapidly for both static applications and the portable electronics enforcing the substantial improvement in battery systems, and Li-ion batteries have been proven to have maximum energy storage density in all rechargeable batteries. However, major breakthroughs are required to consummate the requirement of higher energy density with lower cost to penetrate new markets. Graphite anode having limited capacity has become a bottle neck in the process of developing next generation batteries and can be replaced by higher capacity metals such as Silicon. In the present study we are focusing on the mechanical behavior of the Si-thin film anode under various operating conditions. A numerical model is developed to simulate the intercalation induced stress and the failure mechanism of the complex anode structure. Effect of the various physical phenomena such as diffusion induced stress, plasticity and the crack propagation are investigated to predict better performance parameters for improved design.

Public preferences of the Great Lakes Environment : a Lake Michigan pilot study

The Great Lakes watershed is home to over 40 million people, and the health of the Great Lakes ecosystem is vital to the overall economic, societal, and environmental health of the U.S. and Canada. However, environmental issues related to them are sometimes overlooked. Policymakers and the public face the challenges of balancing economic benefits with the need to conserve and/or replenish regional natural resources to ensure long term prosperity. From the literature review, nine critical stressors of ecological services were delineated, which include pollution and contamination, agricultural erosion, non-native species, degraded recreational resources, loss of wetlands habitat, climate change, risk of clean water shortage, vanishing sand dunes, and population overcrowding; this list was validated through a series of stakeholder discussions and focus groups in Grand Rapids. Focus groups were conducted in Grand Rapids to examine the awareness of, concern with, and willingness to expend resources on these stressors. Stressors that the respondents have direct contact with tend to be the most important. The focus group results show that concern related to pollution and contamination is much higher than for any of the other stressors. Low responses to climate change result in recommendations for outreach programs.

Exploration of the role of serum factors in maintaining bone mass during hibernation in black bears

Disuse osteoporosis is a condition in which reduced mechanical loading (e.g. bed-rest, immobilization, or paralysis) results in unbalanced bone turnover. The American black bear is a unique, naturally occurring model for the prevention of disuse osteoporosis. Bears remain mostly inactive for up to half a year of hibernation annually, yet they do not lose bone mechanical strength or structural properties throughout hibernation. The long-term goal of this study is to determine the biological mechanism through which bears maintain bone during hibernation. This mechanism could pinpoint new signaling pathway targets for the development of drugs for osteoporosis prevention. In this study, bone specific alkaline phosphatase (BSALP), a marker of osteoblast activity, and tartrate resistant acid phosphatase (TRACP), a marker of osteoclast number, were quantified in the serum of hibernating and active black bears. BSALP and TRACP decreased during hibernation, suggesting a balanced reduction in bone turnover. This decrease in BSALP and TRACP were correlated positively to serum adiponectin and inversely to serum neuropeptide Y, suggesting a possible role of these hormones in suppressing bone turnover during hibernation. Osteocalcin (OCN) and undercarboxylated OCN increased dramatically in the serum of hibernating bears. These increases were inversely correlated with adiponectin, glucose, and serotonin, suggesting that OCN may have a unique role in energy homeostasis during hibernation. Finally, MC3T3-E1 osteoblasts were cultured in the serum from active and hibernating bears, and seasonal cell responses were quantified. Cells cultured in serum from hibernating bears had a reduced caspase-3/7 response, and more living cells, after apoptotic threat. The caspase-3/7 response was positively correlated to serum adiponectin and to gene expression of OCN and Runx2, suggesting that reduced caspase-3/7 activity may be related to the reduced differentiation potential of osteoblasts in hibernation serum, and that adiponectin is a potential effector hormone. In summary, the activities of osteoblasts and osteoclasts are reduced during hibernation in bears. This reduced turnover is due, in part, to hormonal control. Further study of potential effectors adiponectin and neuropeptide Y may provide insight into the biological mechanism through which bears maintain bone throughout hibernation.

Cholesterol granuloma of the petrous apex: benefit of computer-aided surgery

The following is an analysis of the role of computer aided surgery by infralabyrinthine-subcochlear approach to the petrous apex for cholesterol granulomas with hearing preservation. In a retrospective case review from 1996 to 2008 six patients were analysed in our tertiary referral centre, otorhinolaryngology outpatient clinic. Excellent intraoperative localisation of the carotid artery, facial nerve and the entrance into the cholesterol cyst of the bone by means of the navigation system was seen. Additionally, the operation time decreased from an initial 4 h down to 2 h. The application of computer-aided surgery allows intraoperative monitoring of the position of the tip of the microsurgical instruments in case of a rare disease and in the delicate area of the petrous apex giving a high security level.

Statistics of the received power for free space optical channels

Free space optical (FSO) communication links can experience extreme signal degradation due to atmospheric turbulence induced spatial and temporal irradiance fuctuations (scintillation) in the laser wavefront. In addition, turbulence can cause the laser beam centroid to wander resulting in power fading, and sometimes complete loss of the signal. Spreading of the laser beam and jitter are also artifacts of atmospheric turbulence. To accurately predict the signal fading that occurs in a laser communication system and to get a true picture of how this affects crucial performance parameters like bit error rate (BER) it is important to analyze the probability density function (PDF) of the integrated irradiance fuctuations at the receiver. In addition, it is desirable to find a theoretical distribution that accurately models these ?uctuations under all propagation conditions. The PDF of integrated irradiance fuctuations is calculated from numerical wave-optic simulations of a laser after propagating through atmospheric turbulence to investigate the evolution of the distribution as the aperture diameter is increased. The simulation data distribution is compared to theoretical gamma-gamma and lognormal PDF models under a variety of scintillation regimes from weak to very strong. Our results show that the gamma-gamma PDF provides a good fit to the simulated data distribution for all aperture sizes studied from weak through moderate scintillation. In strong scintillation, the gamma-gamma PDF is a better fit to the distribution for point-like apertures and the lognormal PDF is a better fit for apertures the size of the atmospheric spatial coherence radius ρ0 or larger. In addition, the PDF of received power from a Gaussian laser beam, which has been adaptively compensated at the transmitter before propagation to the receiver of a FSO link in the moderate scintillation regime is investigated. The complexity of the adaptive optics (AO) system is increased in order to investigate the changes in the distribution of the received power and how this affects the BER. For the 10 km link, due to the non-reciprocal nature of the propagation path the optimal beam to transmit is unknown. These results show that a low-order level of complexity in the AO provides a better estimate for the optimal beam to transmit than a higher order for non-reciprocal paths. For the 20 km link distance it was found that, although minimal, all AO complexity levels provided an equivalent improvement in BER and that no AO complexity provided the correction needed for the optimal beam to transmit. Finally, the temporal power spectral density of received power from a FSO communication link is investigated. Simulated and experimental results for the coherence time calculated from the temporal correlation function are presented. Results for both simulation and experimental data show that the coherence time increases as the receiving aperture diameter increases. For finite apertures the coherence time increases as the communication link distance is increased. We conjecture that this is due to the increasing speckle size within the pupil plane of the receiving aperture for an increasing link distance.

Habitat Associations of the Golden-Winged Warbler in Honduras

The global population of the Neotropical migrant Golden-winged Warbler (Vermivora chrysoptera) has declined steadily over the past fifty years. While factors influencing this decline have been well researched on the breeding grounds, little is known about the distribution and habitat requirements of this warbler on its stationary non-breeding range. Recent efforts to quantify the non-breeding habitat requirements of this warbler have focused on Colombia and Costa Rica, though the species ranges as far north as the Yucatan Peninsula, Mexico. To address the gap in knowledge from the northern portion of the non-breeding range, I conducted 80 serial point-count surveys targeting Goldenwinged Warblers at eight field sites in Honduras, Central America. I found that Goldenwinged Warblers occupy a greater variety of habitats than previously recognized, including pine-oak forest and semi-deciduous broadleaf forest. I also documented habitat associations that have not been observed in other parts of the non-breeding range with respect to elevation, rainfall, and spatial segregation by sex. These results demonstrate the need to consider the entire non-breeding range in conservation planning, as Goldenwinged Warbler habitat associations appear to vary regionally.

Evaluation of the performance and cost-effectiveness of pavement sections containing open-graded base courses

Moisture induced distresses have been the prevalent distress type affecting the deterioration of both asphalt and concrete pavement sections. While various surface techniques have been employed over the years to minimize the ingress of moisture into the pavement structural sections, subsurface drainage components like open-graded base courses remain the best alternative in minimizing the time the pavement structural sections are exposed to saturated conditions. This research therefore focuses on assessing the performance and cost-effectiveness of pavement sections containing both treated and untreated open-graded aggregate base materials. Three common roadway aggregates comprising of two virgin aggregates and one recycled aggregate were investigated using four open-ended gradations and two binder types. Laboratory tests were conducted to determine the hydraulic, mechanical and durability characteristics of treated and untreated open-graded mixes made from these three aggregate types. Results of the experimental program show that for the same gradation and mix design types, limestone samples have the greatest drainage capacity, stability to traffic loads and resistance to degradation from environmental conditions like freeze-thaw. However, depending on the gradation and mix design used, all three aggregate types namely limestone, natural gravel and recycled concrete can meet the minimum coefficient of hydraulic conductivity required for good drainage in most pavements. Tests results for both asphalt and cement treated open-graded samples indicate that a percent air void content within the range of 15-25 will produce a treated open-graded base course with sufficient drainage capacity and also long term stability under both traffic and environmental loads. Using the new Mechanistic and Empirical Design Guide software, computer simulations of pavement performance were conducted on pavement sections containing these open-graded base aggregate base materials to determine how the MEPDG predicted pavement performance is sensitive to drainage. Using three truck traffic levels and four climatic regions, results of the computer simulations indicate that the predicted performance was not sensitive to the drainage characteristics of the open-graded base course. Based on the result of the MEPDG predicted pavement performance, the cost-effectiveness of the pavement sections with open-graded base was computed on the assumption that the increase service life experienced by these sections was attributed to the positive effects of subsurface drainage. The two cost analyses used gave two contrasting results with the one indicating that the inclusion of open-graded base courses can lead to substantial savings.

Extracellular matrix of the charophycean green algae

A comprehensive knowledge of cell wallstructure and function throughout the plant kingdom is essential to understanding cell wall evolution. The fundamental understanding of the charophycean green algal cell wall is broadening. The similarities and differences that exist between land plant and algal cell walls provide opportunities to understand plant evolution. A variety of polymers previously associated with higher plants were discovered in the charophycean green algae (CGA), including homogalacturonans, cross-linking glycans, arabinogalactan protein, β-glucans, and cellulose. The cellulose content of CGA cell walls ranged from 6% to 43%, with the higher valuescomparable to that found in the primary cell wall of land plants (20-30%). (1,3)β-glucans were found in the unicellular Chlorokybus atmophyticus, Penium margaritaceum, and Cosmarium turpini, the unbranched filamentous Klebsormidium flaccidum, and the multicellular Chara corallina. The discovery of homogalacturonan in Penium margaritaceum representsthe first confirmation of land plant-type pectinsin desmids and the second rigorous characterization of a pectin polymer from the charophycean algae. Homogalacturonan was also indicated from the basal species Chlorokybus atmophyticus and Klebsormidium flaccidum. There is evidence of branched pectins in Cosmarium turpini and linkage analysis suggests the presence of type I rhamnogalacturonan (RGI). Cross-linking β-glucans are associated with cellulose microfibrils during land plant cell growth, and were found in the cell wall of CGA. The evidence of mixed-linkage glucan (MLG) in the 11 charophytesis both suprising and significant given that MLG was once thought to be specific to some grasses. The organization and structure of Cosmarium turpini and Chara corallina MLG was found to be similar to that of Equisetumspp., whereas the basal species of the CGA, Chlorokybus atmophyticus and Klebsormidium flaccidum, have unique organization of alternating of 3- and 4-linkages. The significance of this result on the evolution of the MLG synthetic pathway has yet to be determined. The extracellular matrix (ECM) of Chlorokybus atmophyticus, Klebsormidium flaccidum, and Spirogyra spp. exhibits significant biochemical diversity, ranging from distinct “land plant” polymers to polysaccharides unique to these algae. The neutral sugar composition of Chlorokybus atmophyticus hot water extract and Spirogyra extracellular polymeric substance (EPS), combined with antibody labeling results, revealed the distinct possibility of an arabinogalactan protein in these organisms. Polysaccharide analysis of Zygnematales (desmid) EPS, indicated a probable range of different EPS backbones and substitution patterns upon the core portions of the molecules. Desmid EPS is predominately composed of a complex matrix of branched, uronic acid containing polysaccharides with ester sulfate substitutions and, as such, has an almost infinite capacity for various hydrogen bonding, hydrophobic interaction and ionic cross-bridging motifs, which characterize their unique function in biofilms. My observations support the hypothesis that members of the CGA represent the phylogenetic line that gave rise to vascular plants and that the primary cell wall of vascular plants many have evolved directly from structures typical of the cell wall of filamentous green algae found in the charophycean green algae.

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.

Length of the annular regime for condensing flows inside a horizontal channel - the experimental determination of its values and its trends

The experiments observe and measure the length of the annular regime in fully condensing quasi-steady (steady-in-the-mean) flows of pure FC-72 vapor in a horizontal condenser (rectangular cross-section of 2 mm height, 15 mm width, and 1 m length). The sides and top of the duct are made of clear plastic that allows flow visualization. The experimental system in which this condenser is used is able to control and achieve different quasi-steady mass flow rates, inlet pressures, and wall cooling conditions (by adjustment of the temperature and flow rate of the cooling water flowing underneath the condensing-plate). The reported correlations and measurements for the annular length are also vital information for determining the length of the annular regime and proposing extended correlation (covering many vapors and a larger parameter set than the experimentally reported version here) by ongoing independent modeling and computational simulation approach.

Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.

Muffler characterization with implementation of the finite element method and experimental techniques

Determining how an exhaust system will perform acoustically before a prototype muffler is built can save the designer both a substantial amount of time and resources. In order to effectively use the simulation tools available it is important to understand what is the most effective tool for the intended purpose of analysis as well as how typical elements in an exhaust system affect muffler performance. An in-depth look at the available tools and their most beneficial uses are presented in this thesis. A full parametric study was conducted using the FEM method for typical muffler elements which was also correlated to experimental results. This thesis lays out the overall ground work on how to accurately predict sound pressure levels in the free field for an exhaust system with the engine properties included. The accuracy of the model is heavily dependent on the correct temperature profile of the model in addition to the accuracy of the source properties. These factors will be discussed in detail and methods for determining them will be presented. The secondary effects of mean flow, which affects both the acoustical wave propagation and the flow noise generation, will be discussed. Effective ways for predicting these secondary effects will be described. Experimental models will be tested on a flow rig that showcases these phenomena.