Computational prediction of the influence of crosslink distribution on the thermo-mechanical properties of epoxies

Experimental studies on epoxies report that the microstructure consists of highly crosslinked localized regions connected with a dispersed phase of low crosslink density. The various thermo-mechanical properties of epoxies might be affected by the crosslink distribution. But as experiments cannot report the exact number of crosslinked covalent bonds present in the structure, molecular dynamics is thus being used in this work to determine the influence of crosslink distribution on thermo-mechanical properties. Molecular dynamics and molecular mechanics simulations are used to establish wellequilibrated molecular models of EPON 862-DETDA epoxy system with a range of crosslink densities and various crosslink distributions. Crosslink distributions are being varied by forming differently crosslinked localized clusters and then by forming different number of crosslinks interconnecting the clusters. Simulations are subsequently used to predict the volume shrinkage, thermal expansion coefficients, and elastic properties of each of the crosslinked systems. The results indicate that elastic properties increase with increasing levels of overall crosslink density and the thermal expansion coefficient decreases with overall crosslink density, both above and below the glass transition temperature. Elastic moduli and coefficients of linear thermal expansion values were found to be different for systems with same overall crosslink density but having different crosslink distributions, thus indicating an effect of the epoxy nanostructure on physical properties. The values of thermo-mechanical properties for all the crosslinked systems are within the range of values reported in literature.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

Molecular modeling of EPON 862-DETDA polymer

EPON 862 is an epoxy resin which is cured with the hardening agent DETDA to form a crosslinked epoxy polymer and is used as a component in modern aircraft structures. These crosslinked polymers are often exposed to prolonged periods of temperatures below glass transition range which cause physical aging to occur. Because physical aging can compromise the performance of epoxies and their composites and because experimental techniques cannot provide all of the necessary physical insight that is needed to fully understand physical aging, efficient computational approaches to predict the effects of physical aging on thermo-mechanical properties are needed. In this study, Molecular Dynamics and Molecular Minimization simulations are being used to establish well-equilibrated, validated molecular models of the EPON 862-DETDA epoxy system with a range of crosslink densities using a united-atom force field. These simulations are subsequently used to predict the glass transition temperature, thermal expansion coefficients, and elastic properties of each of the crosslinked systems for validation of the modeling techniques. The results indicate that glass transition temperature and elastic properties increase with increasing levels of crosslink density and the thermal expansion coefficient decreases with crosslink density, both above and below the glass transition temperature. The results also indicate that there may be an upper limit to crosslink density that can be realistically achieved in epoxy systems. After evaluation of the thermo-mechanical properties, a method is developed to efficiently establish molecular models of epoxy resins that represent the corresponding real molecular structure at specific aging times. Although this approach does not model the physical aging process, it is useful in establishing a molecular model that resembles the physically-aged state for further use in predicting thermo-mechanical properties as a function of aging time. An equation has been predicted based on the results which directly correlate aging time to aged volume of the molecular model. This equation can be helpful for modelers who want to study properties of epoxy resins at different levels of aging but have little information about volume shrinkage occurring during physical aging.

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.

AMPLITUDE-VERSUS-ANGLE ANALYSIS AND WIDE-ANGLE-INVERSION OF CROSSWELL SEISMIC DATA IN A CARBONATE RESERVOIR

Crosswell data set contains a range of angles limited only by the geometry of the source and receiver configuration, the separation of the boreholes and the depth to the target. However, the wide angles reflections present in crosswell imaging result in amplitude-versus-angle (AVA) features not usually observed in surface data. These features include reflections from angles that are near critical and beyond critical for many of the interfaces; some of these reflections are visible only for a small range of angles, presumably near their critical angle. High-resolution crosswell seismic surveys were conducted over a Silurian (Niagaran) reef at two fields in northern Michigan, Springdale and Coldspring. The Springdale wells extended to much greater depths than the reef, and imaging was conducted from above and from beneath the reef. Combining the results from images obtained from above with those from beneath provides additional information, by exhibiting ranges of angles that are different for the two images, especially for reflectors at shallow depths, and second, by providing additional constraints on the solutions for Zoeppritz equations. Inversion of seismic data for impedance has become a standard part of the workflow for quantitative reservoir characterization. Inversion of crosswell data using either deterministic or geostatistical methods can lead to poor results with phase change beyond the critical angle, however, the simultaneous pre-stack inversion of partial angle stacks may be best conducted with restrictions to angles less than critical. Deterministic inversion is designed to yield only a single model of elastic properties (best-fit), while the geostatistical inversion produces multiple models (realizations) of elastic properties, lithology and reservoir properties. Geostatistical inversion produces results with far more detail than deterministic inversion. The magnitude of difference in details between both types of inversion becomes increasingly pronounced for thinner reservoirs, particularly those beyond the vertical resolution of the seismic. For any interface imaged from above and from beneath, the results AVA characters must result from identical contrasts in elastic properties in the two sets of images, albeit in reverse order. An inversion approach to handle both datasets simultaneously, at pre-critical angles, is demonstrated in this work. The main exploration problem for carbonate reefs is determining the porosity distribution. Images of elastic properties, obtained from deterministic and geostatistical simultaneous inversion of a high-resolution crosswell seismic survey were used to obtain the internal structure and reservoir properties (porosity) of Niagaran Michigan reef. The images obtained are the best of any Niagaran pinnacle reef to date.

Local degradation of structural, mechanical, electrical and chemical properties of membrane electrode assembly in polymer electrolyte fuel cell

Polymer electrolyte fuel cell (PEMFC) is promising source of clean power in many applications ranging from portable electronics to automotive and land-based power generation. However, widespread commercialization of PEMFC is primarily challenged by degradation. The mechanisms of fuel cell degradation are not well understood. Even though the numbers of installed units around the world continue to increase and dominate the pre-markets, the present lifetime requirements for fuel cells cannot be guarantee, creating the need for a more comprehensive knowledge of material’s ageing mechanism. The objective of this project is to conduct experiments on membrane electrode assembly (MEA) components of PEMFC to study structural, mechanical, electrical and chemical changes during ageing and understanding failure/degradation mechanism. The first part of this project was devoted to surface roughness analysis on catalyst layer (CL) and gas diffusion layer (GDL) using surface mapping microscopy. This study was motivated by the need to have a quantitative understanding of the GDL and CL surface morphology at the submicron level to predict interfacial contact resistance. Nanoindentation studies using atomic force microscope (AFM) were introduced to investigate the effect of degradation on mechanical properties of CL. The elastic modulus was decreased by 45 % in end of life (EOL) CL as compare to beginning of life (BOL) CL. In another set of experiment, conductive AFM (cAFM) was used to probe the local electric current in CL. The conductivity drops by 62 % in EOL CL. The future task will include characterization of MEA degradation using Raman and Fourier transform infrared (FTIR) spectroscopy. Raman spectroscopy will help to detect degree of structural disorder in CL during degradation. FTIR will help to study the effect of CO in CL. XRD will be used to determine Pt particle size and its crystallinity. In-situ conductive AFM studies using electrochemical cell on CL to correlate its structure with oxygen reduction reaction (ORR) reactivity

Measuring the elastic modulus of polymers using the atomic force microscope

Testing a new method of nanoindentation using the atomic force microscope (AFM) was the purpose of this research. Nanoindentation is a useful technique to study the properties of materials on the sub-micron scale. The AFM has been used as a nanoindenter previously; however several parameters needed to obtain accurate results, including tip radius and cantilever sensitivity, can be difficult to determine. To solve this problem, a new method to determine the elastic modulus of a material using the atomic force microscope (AFM) has been proposed by Tang et al. This method models the cantilever and the sample as two springs in a series. The ratio of the cantilever spring constant (k) to diameter of the tip (2a) is treated in the model as one parameter (α=k/2a). The value of a, along with the cantilever sensitivity, are determined on two reference samples with known mechanical properties and then used to find the elastic modulus of an unknown sample. To determine the reliability and accuracy of this technique, it was tested on several polymers. Traditional depth-sensing nanoindentation was preformed for comparison. The elastic modulus values from the AFM were shown to be statistically similar to the nanoindenter results for three of the five samples tested.

MICROWAVE ABSORPTION PROPERTIES OF TIRES

The waste tire is belonging to insoluble high polymer elastic materials. It takes hundreds of years to resolve the macromolecules of waste tire into the standard which does not pollute the environment. More and more waste tires are air stored which causes space occupation and mosquito-breeding in the places that will spread diseases. The disposal methods include landfill, stockpiles, dumping and incising into particles. However, all these methods are not technically and economically efficient. The trend for the development of waste tire treatment processes is low cost, on-site, and high product recovery at high energy efficiency. In this project, microwave energy has been applied for treatment of the waste tire in laboratory scale. Experimental conditions were varied in order to find the optimum processing parameters such as temperature and atmosphere. The microwave absorption capability of waste tire rubber was investigated by measuring its dielectric properties from room temperature to 800°C in stagnant air and pure nitrogen atmospheres, respectively, at both 915 and 2466MHz.The dielectric parameters data increase steadily at temperatures below 400°C. At temperatures above 400°C, the relative dielectric loss factor and relative dielectric constant begin to decrease. This is due to the solid phase of tire rubber begins to transform to gas phase and the release of volatiles. The calculations of microwave half-power depth and penetration depth of waste tire rubber show that the pyrolysis process significantly improves the microwave absorption capability of the waste tire rubber at low temperatures. The calculated reflection loss of the waste tire rubber suggests that its maximum microwave absorption can be obtained when the rubber has a thickness of 25mm at 915MHz. The sample dimension has a significant effect on the overall performance of microwave absorption in waste tire during pyrolysis and thus on the efficiency of microwave waste tire rubber pyrolysis.