Interplay of martensitic phase transformation and plastic slip in polycrystals

Inspired by key experimental and analytical results regarding Shape Memory Alloys (SMAs), we propose a modelling framework to explore the interplay between martensitic phase transformations and plastic slip in polycrystalline materials, with an eye towards computational efficiency. The resulting framework uses a convexified potential for the internal energy density to capture the stored energy associated with transformation at the meso-scale, and introduces kinetic potentials to govern the evolution of transformation and plastic slip. The framework is novel in the way it treats plasticity on par with transformation.

We implement the framework in the setting of anti-plane shear, using a staggered implicit/explict update: we first use a Fast-Fourier Transform (FFT) solver based on an Augmented Lagrangian formulation to implicitly solve for the full-field displacements of a simulated polycrystal, then explicitly update the volume fraction of martensite and plastic slip using their respective stick-slip type kinetic laws. We observe that, even in this simple setting with an idealized material comprising four martensitic variants and four slip systems, the model recovers a rich variety of SMA type behaviors. We use this model to gain insight into the isothermal behavior of stress-stabilized martensite, looking at the effects of the relative plastic yield strength, the memory of deformation history under non-proportional loading, and several others.

We extend the framework to the generalized 3-D setting, for which the convexified potential is a lower bound on the actual internal energy, and show that the fully implicit discrete time formulation of the framework is governed by a variational principle for mechanical equilibrium. We further propose an extension of the method to finite deformations via an exponential mapping. We implement the generalized framework using an existing Optimal Transport Mesh-free (OTM) solver. We then model the $\alpha$--$\gamma$ and $\alpha$--$\varepsilon$ transformations in pure iron, with an initial attempt in the latter to account for twinning in the parent phase. We demonstrate the scalability of the framework to large scale computing by simulating Taylor impact experiments, observing nearly linear (ideal) speed-up through 256 MPI tasks. Finally, we present preliminary results of a simulated Split-Hopkinson Pressure Bar (SHPB) experiment using the $\alpha$--$\varepsilon$ model.

Fabrication, Characterization, And Deformation of 3D Structural Meta-Materials

Current technological advances in fabrication methods have provided pathways to creating architected structural meta-materials similar to those found in natural organisms that are structurally robust and lightweight, such as diatoms. Structural meta-materials are materials with mechanical properties that are determined by material properties at various length scales, which range from the material microstructure (nm) to the macro-scale architecture (μm – mm). It is now possible to exploit material size effect, which emerge at the nanometer length scale, as well as structural effects to tune the material properties and failure mechanisms of small-scale cellular solids, such as nanolattices. This work demonstrates the fabrication and mechanical properties of 3-dimensional hollow nanolattices in both tension and compression. Hollow gold nanolattices loaded in uniaxial compression demonstrate that strength and stiffness vary as a function of geometry and tube wall thickness. Structural effects were explored by increasing the unit cell angle from 30° to 60° while keeping all other parameters constant; material size effects were probed by varying the tube wall thickness, t, from 200nm to 635nm, at a constant relative density and grain size. In-situ uniaxial compression experiments reveal an order-of-magnitude increase in yield stress and modulus in nanolattices with greater lattice angles, and a 150% increase in the yield strength without a concomitant change in modulus in thicker-walled nanolattices for fixed lattice angles. These results imply that independent control of structural and material size effects enables tunability of mechanical properties of 3-dimensional architected meta-materials and highlight the importance of material, geometric, and microstructural effects in small-scale mechanics. This work also explores the flaw tolerance of 3D hollow-tube alumina kagome nanolattices with and without pre-fabricated notches, both in experiment and simulation. Experiments demonstrate that the hollow kagome nanolattices in uniaxial tension always fail at the same load when the ratio of notch length (a) to sample width (w) is no greater than 1/3, with no correlation between failure occurring at or away from the notch. For notches with (a/w) > 1/3, the samples fail at lower peak loads and this is attributed to the increased compliance as fewer unit cells span the un-notched region. Finite element simulations of the kagome tension samples show that the failure is governed by tensile loading for (a/w) < 1/3 but as (a/w) increases, bending begins to play a significant role in the failure. This work explores the flaw sensitivity of hollow alumina kagome nanolattices in tension, using experiments and simulations, and demonstrates that the discrete-continuum duality of architected structural meta-materials gives rise to their flaw insensitivity even when made entirely of intrinsically brittle materials.

Techniques for strength measurement at high pressures and strain-rates using transverse waves

The study of the strength of a material is relevant to a variety of applications including automobile collisions, armor penetration and inertial confinement fusion. Although dynamic behavior of materials at high pressures and strain-rates has been studied extensively using plate impact experiments, the results provide measurements in one direction only. Material behavior that is dependent on strength is unaccounted for. The research in this study proposes two novel configurations to mitigate this problem.

The first configuration introduced is the oblique wedge experiment, which is comprised of a driver material, an angled target of interest and a backing material used to measure in-situ velocities. Upon impact, a shock wave is generated in the driver material. As the shock encounters the angled target, it is reflected back into the driver and transmitted into the target. Due to the angle of obliquity of the incident wave, a transverse wave is generated that allows the target to be subjected to shear while being compressed by the initial longitudinal shock such that the material does not slip. Using numerical simulations, this study shows that a variety of oblique wedge configurations can be used to study the shear response of materials and this can be extended to strength measurement as well. Experiments were performed on an oblique wedge setup with a copper impactor, polymethylmethacrylate driver, aluminum 6061-t6 target, and a lithium fluoride window. Particle velocities were measured using laser interferometry and results agree well with the simulations.

The second novel configuration is the y-cut quartz sandwich design, which uses the anisotropic properties of y-cut quartz to generate a shear wave that is transmitted into a thin sample. By using an anvil material to back the thin sample, particle velocities measured at the rear surface of the backing plate can be implemented to calculate the shear stress in the material and subsequently the strength. Numerical simulations were conducted to show that this configuration has the ability to measure the strength for a variety of materials.

Adsorption of Pb(II) and Cu(II) on α-quartz from aqueous solutions: influence of pH, ionic strength, and complexing ligands

Adsorption of aqueous Pb(II) and Cu(II) on α-quartz was studied as a function of time, system surface area, and chemical speciation. Experimental systems contained sodium as a major cation, hydroxide, carbonate, and chloride as major anions, and covered the pH range 4 to 8. In some cases citrate and EDTA were added as representative organic complexing agents. The adsorption equilibria were reached quickly, regardless of the system surface area. The positions of the adsorption equilibria were found to be strongly dependent on pH, ionic strength and concentration of citrate and EDTA. The addition of these non-adsorbing ligands resulted in a competition between chelation and adsorption. The experimental work also included the examination of the adsorption behavior of the doubly charged major cations Ca(II) and Mg(II) as a function of pH.

The theoretical description of the experimental systems was obtained by means of chemical equilibrium-plus-adsorption computations using two adsorption models: one mainly electrostatic (the James-Healy Model), and the other mainly chemical (the Ion Exchange-Surface Complex Formation Model). Comparisons were made between these two models.

The main difficulty in the theoretical predictions of the adsorption behavior of Cu(II) was the lack of the reliable data for the second hydrolysis constant(*β_2) The choice of the constant was made on the basis of potentiometric titratlons of Cu^(2+)

The experimental data obtained and the resulting theoretical observations were applied in models of the chemical behavior of trace metals in fresh oxic waters, with emphasis on Pb(II) and Cu(II).

Hydrogen-Atom Transfer Photochemistry of Tetrakis(µ-pyrophosphito)diplatinate(II)

In many senses, the hydrogen-atom transfer reactions observed with the triplet excited state of pyrophosphito-bridged platinum(II) dimers resemble the reactions of organic ketone nπ* states. The first two chapters describe our attempts to understand the reactivity differences between these two chromophores. Reactivity of the metal dimers is strongly regulated by the detailed nature of the ligands that ring the axial site, the hydrogen-abstraction center. A hydrogen-bonded network linking the ligands facilitates H-atom transfer quenching with alcohols through the formation of a hydrogen-bonded complex between the alcohol and a dimer. For substrates of equal C-H bond strength that lack a hydroxyl group (e.g., benzyl hydrocarbons), the quenching rate is several orders of magnitude slower.

The shape and size of the axial site, as determined by the ligands, also discriminate among quenchers by their steric characteristics. Very small quenchers quench slowly because of high entropies of activation, while very large ones have large enthalpic barriers. The two effects find a balance with quenchers of "just the right size."

The third chapter discusses the design of a mass spectrometer that uses positron annihilation to ionize neutral molecules. The mass spectrometer creates positron-molecule adducts whose annihilation produces fragmentation products that may yield information on the bonding of positrons in such complexes.

A. An investigation of the relation between the tensile strength and (Brinell) hardness of non-ferrous alloys. B. A design of a fatigue testing machine for a college laboratory

[no abstract]

Electron tunneling in proteins and water

The subject of this thesis is electronic coupling in donor-bridge-acceptor systems. In Chapter 2, ET properties of cyanide-bridged dinuclear ruthenium complexes were investigated. The strong interaction between the mixed-valent ruthenium centers leads to intense metal-to-metal charge transfer bands (MMCT). Hush analysis of the MMCT absorption bands yields the electronic-coupling strength between the metal centers (H_(AB)) and the total reorganization energy (λ). Comparison of ET kinetics to calculated rates shows that classical ET models fail to account for the observed kinetics and nuclear tunneling must be considered.

In Chapter 3, ET rates were measured in four ruthenium-modified highpotential iron-sulfur proteins (HiPIP), which were modified at position His50, His81, His42 and His18, respectively. ET kinetics for the His50 and His81 mutants are a factor of 300 different, while the donor-acceptor separation is nearly identical. PATHWAY calculations corroborate these measurements and highlight the importance of structural detail of the intervening protein matrix.

In Chapter 4, the distance dependence of ET through water bridges was measured. Photoinduced ET measurements in aqueous glasses at 77 K show that water is a poor medium for ET. Luminescence decay and quantum yield data were analyzed in the context of a quenching model that accounts for the exponential distance dependence of ET, the distance distribution of donors and acceptors embedded in the glass and the excluded volumes generated by the finite sizes of the donors and acceptors.

In Chapter 5, the pH-dependent excited state dynamics of ruthenium-modified amino acids were measured. The [Ru(bpy)_(3)] ^(2+) chromophore was linked to amino acids via an amide linkage. Protonation of the amide oxygen effectively quenches the excited state. In addition. time-resolved and steady-state luminescence data reveal that nonradiative rates are very sensitive to the protonation state and the structure of the amino acid moiety.

Strength of tantalum at high pressures through Richtmyer-Meshkov laser compression experiments and simulations

Strength at extreme pressures (>1 Mbar or 100 GPa) and high strain rates (106-108 s-1) of materials is not well characterized. The goal of the research outlined in this thesis is to study the strength of tantalum (Ta) at these conditions. The Omega Laser in the Laboratory for Laser Energetics in Rochester, New York is used to create such extreme conditions. Targets are designed with ripples or waves on the surface, and these samples are subjected to high pressures using Omega’s high energy laser beams. In these experiments, the observational parameter is the Richtmyer-Meshkov (RM) instability in the form of ripple growth on single-mode ripples. The experimental platform used for these experiments is the “ride-along” laser compression recovery experiments, which provide a way to recover the specimens having been subjected to high pressures. Six different experiments are performed on the Omega laser using single-mode tantalum targets at different laser energies. The energy indicates the amount of laser energy that impinges the target. For each target, values for growth factor are obtained by comparing the profile of ripples before and after the experiment. With increasing energy, the growth factor increased.

Engineering simulations are used to interpret and correlate the measurements of growth factor to a measure of strength. In order to validate the engineering constitutive model for tantalum, a series of simulations are performed using the code Eureka, based on the Optimal Transportation Meshfree (OTM) method. Two different configurations are studied in the simulations: RM instabilities in single and multimode ripples. Six different simulations are performed for the single ripple configuration of the RM instability experiment, with drives corresponding to laser energies used in the experiments. Each successive simulation is performed at higher drive energy, and it is observed that with increasing energy, the growth factor increases. Overall, there is favorable agreement between the data from the simulations and the experiments. The peak growth factors from the simulations and the experiments are within 10% agreement. For the multimode simulations, the goal is to assist in the design of the laser driven experiments using the Omega laser. A series of three-mode and four-mode patterns are simulated at various energies and the resulting growth of the RM instability is computed. Based on the results of the simulations, a configuration is selected for the multimode experiments. These simulations also serve as validation for the constitutive model and the material parameters for tantalum that are used in the simulations.

By designing samples with initial perturbations in the form of single-mode and multimode ripples and subjecting these samples to high pressures, the Richtmyer-Meshkov instability is investigated in both laser compression experiments and simulations. By correlating the growth of these ripples to measures of strength, a better understanding of the strength of tantalum at high pressures is achieved.

Strength, Deformation and Fracture in Metallic Nanostructures

An understanding of the mechanics of nanoscale metals and semiconductors is necessary for the safe and prolonged operation of nanostructured devices from transistors to nanowire- based solar cells to miniaturized electrodes. This is a fascinating but challenging pursuit because mechanical properties that are size-invariant in conventional materials, such as strength, ductility and fracture behavior, can depend critically on sample size when materials are reduced to sub- micron dimensions. In this thesis, the effect of nanoscale sample size, microstructure and structural geometry on mechanical strength, deformation and fracture are explored for several classes of solid materials. Nanocrystalline platinum nano-cylinders with diameters of 60 nm to 1 μm and 12 nm sized grains are fabricated and tested in compression. We find that nano-sized metals containing few grains weaken as sample diameter is reduced relative to grain size due to a change from deformation governed by internal grains to surface grain governed deformation. Fracture at the nanoscale is explored by performing in-situ SEM tension tests on nanocrystalline platinum and amorphous, metallic glass nano-cylinders containing purposely introduced structural flaws. It is found that failure location, mechanism and strength are determined by the stress concentration with the highest local stress whether this is at the structural flaw or a microstructural feature. Principles of nano-mechanics are used to design and test mechanically robust hierarchical nanostructures with structural and electrochemical applications. 2-photon lithography and electroplating are used to fabricate 3D solid Cu octet meso-lattices with micron- scale features that exhibit strength higher than that of bulk Cu. An in-situ SEM lithiation stage is developed and used to simultaneously examine morphological and electrochemical changes in Si-coated Cu meso-lattices that are of interest as high energy capacity electrodes for Li-ion batteries.

I. Thermochemistry and reaction kinetics of disolvated protons by ion cyclotron resonance spectroscopy. II. Thermochemical studies of small fluorocarbons by photoionization mass spectrometry

The disolvated proton, H(OH₂)₂⁺ is employed as a chemical reagent in low pressure (˂ 10^-5 torr) investigations by ion cyclotron resonance spectroscopy. Since termolecular reactions are absent at low pressure, disolvated protons are not generally observed. However H(OH₂)₂⁺ is produced in a sequence of bimolecular reactions in mixtures containing H₂O and one of a small number of organohalide precursors. Then a series of hydrated Lewis bases is produced by H₃O⁺ transfer from H(OH₂)₂⁺. In Chapter II, the relative stability of hydrated bases containing heteroatoms of both first and second row elements is determined from the preferred direction of H₃O⁺ transfer between BH(OH₂)⁺ complexes. S and P containing bases are shown to bind H₃O⁺ more weakly than O and N bases with comparable proton affinities. A simple model of hydrogen bonding is proposed to account for these observations.

H⁺ transfer from H(OH₂)₂⁺ to several Lewis bases also occurs at low pressure. In Chapter III the relative importance of H₃O⁺ transfer and H⁺ transfer from H(OH₂)₂⁺ to a series of bases is observed to be a function of base strength. Beginning with CH₃COOH, the weakest base for which H⁺ transfer is observed, the importance of H⁺ transfer increases with increasing proton affinity of the acceptor base. The nature of neutral products formed from H(OH₂)₂⁺ by loss of H⁺ is also considered.

Chapters IV and V deal with thermochemistry of small fluorocarbons determined by photoionization mass spectrometry. The enthalpy of formation of CF₂ is considered in Chapter IV. Photoionization of perfluoropropylene, perfluorocyclopropane, and trifluoromethyl benzene yield onsets for ions formed by loss of a CF₂ neutral fragment. Earlier determinations of ΔH^°_f298 (CF₂) are reinterpreted using updated thermochemical values and compared with results of this study. The heat of formation of neutral perfluorocyclopropane is also derived. Finally, the energetics of interconversion of perfluoropropylene and perfluorocyclopropane are considered for both the neutrals and their molecular ions.

In Chapter V the heats of formation of CF₃⁺ and CF₃I⁺are derived from photoionization of CF₃I. These are considered with respect to ion-molecule reactions observed in CF₃I monitored by the techniques of ion cyclotron resonance spectroscopy. Results obtained in previous experiments are also compared.

Strength of thin-walled elliptical cylinders supported at the minor axis

In this investigation it was found that the instability failure of curved sheet is nearly independent of the type of loading and is primarily a function of the maximum stress, radius-thickness ration and modulus of elasticity. A method of correlating the critical stress of thin sheet under several different types of loading is given. An explanation for the experimental critical stress of thin walled cylinders under bending being greater than that for pure compression is given. The strength of unstiffened thin walled circular nose sections under pure bending was found to be controlled by local instability of the section, rather than a large scale instability. The equation of local instability of curved sheet gives values which are in fair agreement with those found experimentally.

The strength of elliptical cylinders supported at the minor axis under bending plus shear loads is governed primarily by the bending strength, and is little effected by the sheer force unless the amount of shear is quite large with respect to the moment. The effect of increasing the amount of elliptically greatly reduces the bending and shear strength of nose sections. Under torsional loads the stress at buckling falls off as the ration of the major to minor axis increases but the failure stress decreases at a slower rate than the buckling stress. The length effect of semi-circular sections under torsion is similar to that of a circular tube, and can be obtained by Donnell's theoretical equation.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

Part I: Total yield measurements for the ^(21)Ne(α,n)^)(24)Mg reaction. Part II: A study of the ^(12)C(He,p)^(14)N reaction

PART I

The total cross-section for the reaction ²¹Ne(α, n)²⁴Mg has been measured in the energy range 1.49 Mev ≤ E_cm ≤ 2.6 Mev. The cross-section factor, S(O), for this reaction has been determined, by means of an optical model calculation, to be in the range 1.52 x 10¹² mb-Mev to 2.67 x 10¹² mb-Mev, for interaction radii in the range 5.0 fm to 6.6 fm. With S(O) ≈ 2 x 10¹² mb-Mev, the reaction ²¹Ne(α, n)²⁴Mg can produce a large enough neutron flux to be a significant astrophysical source of neutrons.

PART II

The reaction¹²C(³He, p)¹⁴N has been studied over the energy range 12 Mev ≤ E_lab ≤ 18 Mev. Angular distributions of the proton groups leading to the lowest seven levels in ¹⁴N were obtained.

Distorted wave calculations, based on two-nucleon transfer theory, were performed, and were found to be reliable for obtaining the value of the orbital angular momentum transferred. The present work shows that such calculations do not yield unambiguous values for the spectroscopic factors.