Topological changes in 2D simplicial meshes for the simulation of fractures

A method for the introduction of strong discontinuities into a mesh will be developed. This method, applicable to a number of eXtended Finite Element Methods (XFEM) with intra-element strong discontinuities will be demonstrated with one specific method: the Generalized Cohesive Element (GCE) method. The algorithm utilizes a subgraph mesh representation which may insert the GCE either adaptively during the course of the analysis or a priori. Using this subgraphing algorithm, the insertion time is O(n) to the number of insertions. Numerical examples are presented demonstrating the advantages of the subgraph insertion method.

Charge induced enhancement of adsorption for hydrogen storage materials

The rising concerns about environmental pollution and global warming have facilitated research interest in hydrogen energy as an alternative energy source. To apply hydrogen for transportations, several issues have to be solved, within which hydrogen storage is the most critical problem. Lots of materials and devices have been developed; however, none is able to meet the DOE storage target. The primary issue for hydrogen physisorption is a weak interaction between hydrogen and the surface of solid materials, resulting negligible adsorption at room temperature. To solve this issue, there is a need to increase the interaction between the hydrogen molecules and adsorbent surface. In this study, intrinsic electric dipole is investigated to enhance the adsorption energy. The results from the computer simulation of single ionic compounds with hydrogen molecules to form hydrogen clusters showed that electrical charge of substances plays an important role in generation of attractive interaction with hydrogen molecules. In order to further examine the effects of static interaction on hydrogen adsorption, activated carbon with a large surface area was impregnated with various ionic salts including LiCl, NaCl, KCl, KBr, and NiCl and their performance for hydrogen storage was evaluated by using a volumetric method. Corresponding computer simulations have been carried out by using DFT (Density Functional Theory) method combined with point charge arrays. Both experimental and computational results prove that the adsorption capacity of hydrogen and its interaction with the solid materials increased with electrical dipole moment. Besides the intrinsic dipole, an externally applied electric field could be another means to enhance hydrogen adsorption. Hydrogen adsorption under an applied electric field was examined by using porous nickel foil as electrodes. Electrical signals showed that adsorption capacity increased with the increasing of gas pressure and external electric voltage. Direct measurement of the amount of hydrogen adsorption was also carried out with porous nickel oxides and magnesium oxides using the piezoelectric material PMN-PT as the charge supplier due to the pressure. The adsorption enhancement from the PMN-PT generated charges is obvious at hydrogen pressure between 0 and 60 bars, where the hydrogen uptake is increased at about 35% for nickel oxide and 25% for magnesium oxide. Computer simulation reveals that under the external electric field, the electron cloud of hydrogen molecules is pulled over to the adsorbent site and can overlap with the adsorbent electrons, which in turn enhances the adsorption energy Experiments were also carried out to examine the effects of hydrogen spillover with charge induced enhancement. The results show that the overall storage capacity in nickel oxide increased remarkably by a factor of 4.

MULTISCALE EXAMINATION AND MODELING OF ELECTRON TRANSPORT IN NANOSCALE MATERIALS AND DEVICES

For half a century the integrated circuits (ICs) that make up the heart of electronic devices have been steadily improving by shrinking at an exponential rate. However, as the current crop of ICs get smaller and the insulating layers involved become thinner, electrons leak through due to quantum mechanical tunneling. This is one of several issues which will bring an end to this incredible streak of exponential improvement of this type of transistor device, after which future improvements will have to come from employing fundamentally different transistor architecture rather than fine tuning and miniaturizing the metal-oxide-semiconductor field effect transistors (MOSFETs) in use today. Several new transistor designs, some designed and built here at Michigan Tech, involve electrons tunneling their way through arrays of nanoparticles. We use a multi-scale approach to model these devices and study their behavior. For investigating the tunneling characteristics of the individual junctions, we use a first-principles approach to model conduction between sub-nanometer gold particles. To estimate the change in energy due to the movement of individual electrons, we use the finite element method to calculate electrostatic capacitances. The kinetic Monte Carlo method allows us to use our knowledge of these details to simulate the dynamics of an entire device— sometimes consisting of hundreds of individual particles—and watch as a device ‘turns on’ and starts conducting an electric current. Scanning tunneling microscopy (STM) and the closely related scanning tunneling spectroscopy (STS) are a family of powerful experimental techniques that allow for the probing and imaging of surfaces and molecules at atomic resolution. However, interpretation of the results often requires comparison with theoretical and computational models. We have developed a new method for calculating STM topographs and STS spectra. This method combines an established method for approximating the geometric variation of the electronic density of states, with a modern method for calculating spin-dependent tunneling currents, offering a unique balance between accuracy and accessibility.

ADVANCED COMPUTATIONAL MODELING OF THE INTERNAL STRUCTURE OF SMART WIND-TURBINE BLADES

Implementation of stable aeroelastic models with the ability to capture the complex features of Multi concept smartblades is a prime step in reducing the uncertainties that come along with blade dynamics. The numerical simulations of fluid structure interaction can thus be used to test a realistic scenarios comprising of full-scale blades at a reasonably low computational cost. A code which was a combination of two advanced numerical models was designed and was run with the help of paralell HPC supercomputer platform. The first model was based on a variation of dimensional reduction technique proposed by Hodges and Yu. This model was the one to record the structural response of heterogenous composite blades. This technique reduces the geometrical complexities of the heterogenous blade section into a stiffness matrix for an equivalent beam. This derived equivalent 1-D strain energy matrix is similar to the actual 3-D strain energy matrix in an asymptotic sense. As this 1-D matrix helps in accurately modeling the blade structure as a 1-D finite element problem, this substantially redues the computational effort and subsequently the computational cost that are required to model the structural dynamics at each step. Second model comprises of implementation of the Blade Element Momentum Theory. In this approach we map all the velocities and the forces with the help of orthogonal matrices that help in capturing the large deformations and the effects of rotations in calculating the aerodynamic forces. This ultimately helps us to take into account the complex flexo torsional deformations. In this thesis we have succesfully tested these computayinal tools developed by MTU’s research team lead by for the aero elastic analysis of wind-turbine blades. The validation in this thesis is majorly based on several experiments done on NREL-5MW blade, as this is widely accepted as a benchmark blade in the wind industry. Along with the use of this innovative model the internal blade structure was also changed to add up to the existing benefits of the already advanced numerical models.

CHARACTERIZATION OF ASPHALT MATERIALS CONTAINING BIO OIL FROM MICHIGAN WOOD

The objective of this research is to develop sustainable wood-blend bioasphalt and characterize the atomic, molecular and bulk-scale behavior necessary to produce advanced asphalt paving mixtures. Bioasphalt was manufactured from Aspen, Basswood, Red Maple, Balsam, Maple, Pine, Beech and Magnolia wood via a 25 KWt fast-pyrolysis plant at 500 °C and refined into two distinct end forms - non-treated (5.54% moisture) and treated bioasphalt (1% moisture). Michigan petroleum-based asphalt, Performance Grade (PG) 58-28 was modified with 2, 5 and 10% of the bioasphalt by weight of base asphalt and characterized with the gas chromatography-mass spectroscopy (GC-MS), Fourier Transform Infra-red (FTIR) spectroscopy and the automated flocculation titrimetry techniques. The GC-MS method was used to characterize the Carbon-Hydrogen-Nitrogen (CHN) elemental ratio whiles the FTIR and the AFT were used to characterize the oxidative aging performance and the solubility parameters, respectively. For rheological characterization, the rotational viscosity, dynamic shear modulus and flexural bending methods are used in evaluating the low, intermediate and high temperature performance of the bio-modified asphalt materials. 54 5E3 (maximum of 3 million expected equivalent standard axle traffic loads) asphalt paving mixes were then prepared and characterized to investigate their laboratory permanent deformation, dynamic mix stiffness, moisture susceptibility, workability and constructability performance. From the research investigations, it was concluded that: 1) levo, 2, 6 dimethoxyphenol, 2 methoxy 4 vinylphenol, 2 methyl 1-2 cyclopentandione and 4-allyl-2, 6 dimetoxyphenol are the dominant chemical functional groups; 2) bioasphalt increases the viscosity and dynamic shear modulus of traditional asphalt binders; 3) Bio-modified petroleum asphalt can provide low-temperature cracking resistance benefits at -18 °C but is susceptible to cracking at -24 °C; 3) Carbonyl and sulphoxide oxidation in petroleum-based asphalt increases with increasing bioasphalt modifiers; 4) bioasphalt causes the asphaltene fractions in petroleum-based asphalt to precipitate out of the solvent maltene fractions; 5) there is no definite improvement or decline in the dynamic mix behavior of bio-modified mixes at low temperatures; 6) bio-modified asphalt mixes exhibit better rutting performance than traditional asphalt mixes; 7) bio-modified asphalt mixes have lower susceptibility to moisture damage; 8) more field compaction energy is needed to compact bio-modified mixes.

LIFE CYCLE ASSESSMENT OF BIOFUEL PRODUCED FROM ALGAE

Algae are considered a promising source of biofuels in the future. However, the environmental impact of algae-based fuel has high variability in previous LCA studies due to lack of accurate data from researchers and industry. The National Alliance for Advanced Biofuels and Bioproducts (NAABB) project was designed to produce and evaluate new technologies that can be implemented by the algal biofuel industry and establish the overall process sustainability. The MTU research group within NAABB worked on the environmental sustainability part of the consortium with UOP-Honeywell and with the University of Arizona (Dr. Paul Blowers). Several life cycle analysis (LCA) models were developed within the GREET Model and SimaPro 7.3 software to quantitatively assess the environment viability and sustainability of algal fuel processes. The baseline GREET Harmonized algae life cycle was expanded and replicated in SimaPro software, important differences in emission factors between GREET/E-Grid database and SimaPro/Ecoinvent database were compared, and adjustments were made to the SimaPro analyses. The results indicated that in most cases SimaPro has a higher emission penalty for inputs of electricity, chemicals, and other materials to the algae biofuels life cycle. A system-wide model of algae life cycle was made starting with preliminary data from the literature, and then progressed to detailed analyses based on inputs from all NAABB research areas, and finally several important scenarios in the algae life cycle were investigated as variations to the baseline scenario. Scenarios include conversion to jet fuel instead of biodiesel or renewable diesel, impacts of infrastructure for algae cultivation, co-product allocation methodology, and different usage of lipid-extracted algae (LEA). The infrastructure impact of algae cultivation is minimal compared to the overall life cycle. However, in the scenarios investigating LEA usage for animal feed instead of internal recycling for energy use and nutrient recovery the results reflect the high potential variability in LCA results. Calculated life cycle GHG values for biofuel production scenarios where LEA is used as animal feed ranged from a 55% reduction to 127% increase compared to the GREET baseline scenario depending on the choice of feed meal. Different allocation methods also affect LCA results significantly. Four novel harvesting technologies and two extraction technologies provided by the NAABB internal report have been analysis using SimaPro LCA software. The results indicated that a combination of acoustic extraction and acoustic harvesting technologies show the most promising result of all combinations to optimize the extraction of algae oil from algae. These scenario evaluations provide important insights for consideration when planning for the future of an algae-based biofuel industry.

Processing and characterization of PbSnTe-based thermoelectric materials made by mechanical alloying

The research reported in this dissertation investigates the processes required to mechanically alloy Pb1-xSnxTe and AgSbTe2 and a method of combining these two end compounds to result in (y)(AgSbTe2)–(1 - y)(Pb1-xSnxTe) thermoelectric materials for power generation applications. In general, traditional melt processing of these alloys has employed high purity materials that are subjected to time and energy intensive processes that result in highly functional material that is not easily reproducible. This research reports the development of mechanical alloying processes using commercially available 99.9% pure elemental powders in order to provide a basis for the economical production of highly functional thermoelectric materials. Though there have been reports of high and low ZT materials fabricated by both melt alloying and mechanical alloying, the processing-structure-properties-performance relationship connecting how the material is made to its resulting functionality is poorly understood. This is particularly true for mechanically alloyed material, motivating an effort to investigate bulk material within the (y)(AgSbTe2)–(1 - y)(Pb1-xSnx- Te) system using the mechanical alloying method. This research adds to the body of knowledge concerning the way in which mechanical alloying can be used to efficiently produce high ZT thermoelectric materials. The processes required to mechanically alloy elemental powders to form Pb1-xSnxTe and AgSbTe2 and to subsequently consolidate the alloyed powder is described. The composition, phases present in the alloy, volume percent, size and spacing of the phases are reported. The room temperature electronic transport properties of electrical conductivity, carrier concentration and carrier mobility are reported for each alloy and the effect of the presence of any secondary phase on the electronic transport properties is described. An mechanical mixing approach for incorporating the end compounds to result in (y)(AgSbTe2)–(1-y)(Pb1-xSnxTe) is described and when 5 vol.% AgSbTe2 was incorporated was found to form a solid solution with the Pb1-xSnxTe phase. An initial attempt to change the carrier concentration of the Pb1-xSnxTe phase was made by adding excess Te and found that the carrier density of the alloys in this work are not sensitive to excess Te. It has been demonstrated using the processing techniques reported in this research that this material system, when appropriately doped, has the potential to perform as highly functional thermoelectric material.

Development and characterization of fluorescent pH sensors based on porous silica and hydrogel support matrices

The hydrogen ion activity (pH) is a very important parameter in environment monitoring, biomedical research and other applications. Optical pH sensors have several advantages over traditional potentiometric pH measurement, such as high sensitivity, no need of constant calibration, easy for miniaturization and possibility for remote sensing. Several pH indicators has been successfully immobilized in three different solid porous materials to use as pH sensing probes. The fluorescent pH indicator fluorescein-5-isothiocyanate (FITC) was covalently bound onto the internal surface of porous silica (pore size ~10 nm) and retained its pH sensitivity. The excited state pK* a of FITC in porous silica (5.58) was slightly smaller than in solution (5.68) due to the free silanol groups (Si-OH) on the silica surface. The pH sensitive range for this probe is pH 4.5 - 7.0 with an error less than 0.1 pH units. The probe response was reproducible and stable for at least four month, stored in DI water, but exhibit a long equilibrium of up to 100 minutes. Sol-gel based pH sensors were developed with immobilization of two fluorescent pH indicators fluorescein-5-(and-6)-sulfonic acid, trisodium salt (FS) and 8-hydroxypyrene- 1,3,6-trisulfonic acid (HPTS) through physical entrapment. Prior to immobilization, the indicators were ion-paired with a common surfactant hexadecyltrimethylammonium bromide (CTAB) in order to prevent leaching. The sol-gel films were synthesized through the hydrolysis of two different precursors, ethyltriethoxysilane (ETEOS) and 3- glycidoxypropyltrimethoxysilane (GPTMS) and deposited on a quartz slide through spin coating. The pK a of the indicators immobilized in sol-gel films was much smaller than in solutions due to silanol groups on the inner surface of the sol-gel films and ammonium groups from the surrounding surfactants. Unlike in solution, the apparent pK a of the indicators in sol-gel films increased with increasing ionic strength. The equilibrium time for these sensors was within 5 minutes (with film thickness of ~470 nm). Polyethylene glycol (PEG) hydrogel was of interest for optical pH sensor development because it is highly proton permeable, transparent and easy to synthesize. pH indicators can be immobilized in hydrogel through physical entrapment and copolymerization. FS and HPTS ion-pairs were physically entrapped in hydrogel matrix synthesized via free radical initiation. For covalent immobilization, three indicators, 6,8-dihydroxypyrene-1,3- disulfonic acid (DHPDS), 2,7-dihydroxynaphthalene-3,6-disulfonic acid (DHNDS) and cresol red were first reacted with methacrylic anhydride (MA) to form methacryloylanalogs for copolymerization. These hydrogels were synthesized in aqueous solution with a redox initiation system. The thickness of the hydrogel film is controlled as ~ 0.5 cm and the porosity can be adjusted with the percentage of polyethylene glycol in the precursor solutions. The pK a of the indicators immobilized in the hydrogel both physically and covalently were higher than in solution due to the medium effect. The sensors are stable and reproducible with a short equilibrium time (less than 4 minutes). In addition, the color change of cresol red immobilized hydrogel is vivid from yellow (acidic condition) to purple (basic condition). Due to covalently binding, cresol red was not leaching out from the hydrogel, making it a good candidate of reusable "pH paper".

A study of acoustic emission technique for concrete damage detection

The Acoustic emission (AE) technique, as one of non-intrusive and nondestructive evaluation techniques, acquires and analyzes the signals emitting from deformation or fracture of materials/structures under service loading. The AE technique has been successfully applied in damage detection in various materials such as metal, alloy, concrete, polymers and other composite materials. In this study, the AE technique was used for detecting crack behavior within concrete specimens under mechanical and environmental frost loadings. The instrumentations of the AE system used in this study include a low-frequency AE sensor, a computer-based data acquisition device and a preamplifier linking the AE sensor and the data acquisition device. The AE system purchased from Mistras Group was used in this study. The AE technique was applied to detect damage with the following laboratory tests: the pencil lead test, the mechanical three-point single-edge notched beam bending (SEB) test, and the freeze-thaw damage test. Firstly, the pencil lead test was conducted to verify the attenuation phenomenon of AE signals through concrete materials. The value of attenuation was also quantified. Also, the obtained signals indicated that this AE system was properly setup to detect damage in concrete. Secondly, the SEB test with lab-prepared concrete beam was conducted by employing Mechanical Testing System (MTS) and AE system. The cumulative AE events and the measured loading curves, which both used the crack-tip open displacement (CTOD) as the horizontal coordinate, were plotted. It was found that the detected AE events were qualitatively correlated with the global force-displacement behavior of the specimen. The Weibull distribution was vii proposed to quantitatively describe the rupture probability density function. The linear regression analysis was conducted to calibrate the Weibull distribution parameters with detected AE signals and to predict the rupture probability as a function of CTOD for the specimen. Finally, the controlled concrete freeze-thaw cyclic tests were designed and the AE technique was planned to investigate the internal frost damage process of concrete specimens.