Matching Items (11)
Filtering by
- All Subjects: Molecular Dynamics
- All Subjects: Aerospace Engineering
- Creators: Jiao, Yang
Description
The study of deflagration to detonation transition (DDT) in explosives is of prime importance with regards to insensitive munitions (IM). Critical damage owing to thermal or shock stimuli could translate to significant loss of life and material. The present study models detonation and deflagration of a commonly used granular explosive: cyclotetramethylene-tetranitramine, HMX. A robust literature review is followed by computational modeling of gas gun and DDT tube test data using the Sandia National Lab three-dimensional multi-material Eulerian hydrocode CTH. This dissertation proposes new computational practices and models that aid in predicting shock stimulus IM response. CTH was first used to model experimental data sets of DDT tubes from both Naval Surface Weapons Center and Los Alamos National Laboratory which were initiated by pyrogenic material and a piston, respectively. Analytical verification was performed, where possible, for detonation via empirical based equations at the Chapman Jouguet state with errors below 2.1%, and deflagration via pressure dependent burn rate equations. CTH simulations include inert, history variable reactive burn and Arrhenius models. The results are in excellent agreement with published HMX detonation velocities. Novel additions include accurate simulation of the pyrogenic material BKNO3 and the inclusion of porosity in energetic materials. The treatment of compaction is especially important in modeling precursory hotspots, caused by hydrodynamic collapse of void regions or grain interactions, prior to DDT of granular explosives. The CTH compaction model of HMX was verified within 11% error via a five pronged validation approach using gas gun data and employed use of a newly generated set of P-α parameters for granular HMX in a Mie-Gruneisen Equation of State. Next, the additions of compaction were extended to a volumetric surface burning model of HMX and compare well to a set of empirical burn rates. Lastly, the compendium of detonation and deflagration models was applied to the aforementioned DDT tubes and demonstrate working functionalities of all models, albeit at the expense of significant computational resources. A robust hydrocode methodology is proposed to make use of the deflagration, compaction and detonation models as a means to predict IM response to shock stimulus of granular explosive materials.
ContributorsMahon, Kelly Susan (Author) / Lee, Taewoo (Thesis advisor) / Herrmann, Marcus (Committee member) / Chen, Kangping (Committee member) / Jiao, Yang (Committee member) / Huang, Huei-Ping (Committee member) / Arizona State University (Publisher)
Created2015
Description
Materials with unprecedented properties are necessary to make dramatic changes in current and future aerospace platforms. Hybrid materials and composites are increasingly being used in aircraft and spacecraft frames; however, future platforms will require an optimal design of novel materials that enable operation in a variety of environments and produce known/predicted damage mechanisms. Nanocomposites and nanoengineered composites with CNTs have the potential to make significant improvements in strength, stiffness, fracture toughness, flame retardancy and resistance to corrosion. Therefore, these materials have generated tremendous scientific and technical interest over the past decade and various architectures are being explored for applications to light-weight airframe structures. However, the success of such materials with significantly improved performance metrics requires careful control of the parameters during synthesis and processing. Their implementation is also limited due to the lack of complete understanding of the effects the nanoparticles impart to the bulk properties of composites. It is common for computational methods to be applied to explain phenomena measured or observed experimentally. Frequently, a given phenomenon or material property is only considered to be fully understood when the associated physics has been identified through accompanying calculations or simulations.
The computationally and experimentally integrated research presented in this dissertation provides improved understanding of the mechanical behavior and response including damage and failure in CNT nanocomposites, enhancing confidence in their applications. The computations at the atomistic level helps to understand the underlying mechanochemistry and allow a systematic investigation of the complex CNT architectures and the material performance across a wide range of parameters. Simulation of the bond breakage phenomena and development of the interface to continuum scale damage captures the effects of applied loading and damage precursor and provides insight into the safety of nanoengineered composites under service loads. The validated modeling methodology is expected to be a step in the direction of computationally-assisted design and certification of novel materials, thus liberating the pace of their implementation in future applications.
The computationally and experimentally integrated research presented in this dissertation provides improved understanding of the mechanical behavior and response including damage and failure in CNT nanocomposites, enhancing confidence in their applications. The computations at the atomistic level helps to understand the underlying mechanochemistry and allow a systematic investigation of the complex CNT architectures and the material performance across a wide range of parameters. Simulation of the bond breakage phenomena and development of the interface to continuum scale damage captures the effects of applied loading and damage precursor and provides insight into the safety of nanoengineered composites under service loads. The validated modeling methodology is expected to be a step in the direction of computationally-assisted design and certification of novel materials, thus liberating the pace of their implementation in future applications.
ContributorsSubramanian, Nithya (Author) / Chattopadhyay, Aditi (Thesis advisor) / Dai, Lenore (Committee member) / Jiao, Yang (Committee member) / Liu, Yongming (Committee member) / Rajadas, John (Committee member) / Arizona State University (Publisher)
Created2018
Description
Improved knowledge connecting the chemistry, structure, and properties of polymers is necessary to develop advanced materials in a materials-by-design approach. Molecular dynamics (MD) simulations can provide tremendous insight into how the fine details of chemistry, molecular architecture, and microstructure affect many physical properties; however, they face well-known restrictions in their applicable temporal and spatial scales. These limitations have motivated the development of computationally-efficient, coarse-grained methods to investigate how microstructural details affect thermophysical properties. In this dissertation, I summarize my research work in structure-based coarse-graining methods to establish the link between molecular-scale structure and macroscopic properties of two different polymers. Systematically coarse-grained models were developed to study the viscoelastic stress response of polyurea, a copolymer that segregates into rigid and viscous phases, at time scales characteristic of blast and impact loading. With the application of appropriate scaling parameters, the coarse-grained models can predict viscoelastic properties with a speed up of 5-6 orders of magnitude relative to the atomistic MD models. Coarse-grained models of polyethylene were also created to investigate the thermomechanical material response under shock loading. As structure-based coarse-grained methods are generally not transferable to states different from which they were calibrated at, their applicability for modeling non-equilibrium processes such as shock and impact is highly limited. To address this problem, a new model is developed that incorporates many-body interactions and is calibrated across a range of different thermodynamic states using a least square minimization scheme. The new model is validated by comparing shock Hugoniot properties with atomistic and experimental data for polyethylene. Lastly, a high fidelity coarse-grained model of polyethylene was constructed that reproduces the joint-probability distributions of structural variables such as the distributions of bond lengths and bond angles between sequential coarse-grained sites along polymer chains. This new model accurately represents the structure of both the amorphous and crystal phases of polyethylene and enabling investigation of how polymer processing such as cold-drawing and bulk crystallization affect material structure at significantly larger time and length scales than traditional molecular simulations.
ContributorsAgrawal, Vipin (Author) / Oswald, Jay (Thesis advisor) / Peralta, Pedro (Committee member) / Chamberlin, Ralph (Committee member) / Solanki, Kiran (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2016
Description
A hybrid molecular dynamics (MD) simulation framework is developed to emulate mechanochemical reaction of mechanophores in epoxy-based nanocomposites. Two different force fields, a classical force field and a bond order based force field are hybridized to mimic the experimental processes from specimen preparation to mechanical loading test. Ultra-violet photodimerization for mechanophore synthesis and epoxy curing for thermoset polymer generation are successfully simulated by developing a numerical covalent bond generation method using the classical force field within the framework. Mechanical loading tests to activate mechanophores are also virtually conducted by deforming the volume of a simulation unit cell. The unit cell deformation leads to covalent bond elongation and subsequent bond breakage, which is captured using the bond order based force field. The outcome of the virtual loading test is used for local work analysis, which enables a quantitative study of mechanophore activation. Through the local work analysis, the onset and evolution of mechanophore activation indicating damage initiation and propagation are estimated; ultimately, the mechanophore sensitivity to external stress is evaluated. The virtual loading tests also provide accurate estimations of mechanical properties such as elastic, shear, bulk modulus, yield strain/strength, and Poisson’s ratio of the system. Experimental studies are performed in conjunction with the simulation work to validate the hybrid MD simulation framework. Less than 2% error in estimations of glass transition temperature (Tg) is observed with experimentally measured Tgs by use of differential scanning calorimetry. Virtual loading tests successfully reproduce the stress-strain curve capturing the effect of mechanophore inclusion on mechanical properties of epoxy polymer; comparable changes in Young’s modulus and yield strength are observed in experiments and simulations. Early damage signal detection, which is identified in experiments by observing increased intensity before the yield strain, is captured in simulations by showing that the critical strain representing the onset of the mechanophore activation occurs before the estimated yield strain. It is anticipated that the experimentally validated hybrid MD framework presented in this dissertation will provide a low-cost alternative to additional experiments that are required for optimizing material design parameters to improve damage sensing capability and mechanical properties.
In addition to the study of mechanochemical reaction analysis, an atomistic model of interphase in carbon fiber reinforced composites is developed. Physical entanglement between semi-crystalline carbon fiber surface and polymer matrix is captured by introducing voids in multiple graphene layers, which allow polymer matrix to intertwine with graphene layers. The hybrid MD framework is used with some modifications to estimate interphase properties that include the effect of the physical entanglement. The results are compared with existing carbon fiber surface models that assume that carbon fiber has a crystalline structure and hence are unable to capture the physical entanglement. Results indicate that the current model shows larger stress gradients across the material interphase. These large stress gradients increase the viscoplasticity and damage effects at the interphase. The results are important for improved prediction of the nonlinear response and damage evolution in composite materials.
In addition to the study of mechanochemical reaction analysis, an atomistic model of interphase in carbon fiber reinforced composites is developed. Physical entanglement between semi-crystalline carbon fiber surface and polymer matrix is captured by introducing voids in multiple graphene layers, which allow polymer matrix to intertwine with graphene layers. The hybrid MD framework is used with some modifications to estimate interphase properties that include the effect of the physical entanglement. The results are compared with existing carbon fiber surface models that assume that carbon fiber has a crystalline structure and hence are unable to capture the physical entanglement. Results indicate that the current model shows larger stress gradients across the material interphase. These large stress gradients increase the viscoplasticity and damage effects at the interphase. The results are important for improved prediction of the nonlinear response and damage evolution in composite materials.
ContributorsKoo, Bonsung (Author) / Chattopadhyay, Aditi (Thesis advisor) / Dai, Lenore (Committee member) / Jiang, Hanqing (Committee member) / Jiao, Yang (Committee member) / Rajadas, John (Committee member) / Arizona State University (Publisher)
Created2017
Description
Applications such as heat exchangers, surface-based cellular structures, rotating blades, and waveguides rely on thin metal walls as crucial constituent elements of the structure. The design freedom enabled by laser powder bed fusion has led to an interest in exploiting this technology to further the performance of these components, many of which retain their as-built surface morphologies on account of their design complexity. However, there is limited understanding of how and why mechanical properties vary by wall thickness for specimens that are additively manufactured and maintain an as-printed surface finish. Critically, the contributions of microstructure and morphology to the mechanical behavior of thin wall laser powder bed fusion structures have yet to be systematically identified and decoupled. This work focuses on elucidating the room temperature quasi-static tensile and high cycle fatigue properties of as-printed, thin-wall Inconel 718 fabricated using laser powder bed fusion, with the aim of addressing this critical gap in the literature. Wall thicknesses studied range from 0.3 - 2.0 mm, and the effects of Hot Isostatic Pressing are also examined, with sheet metal specimens used as a baseline for comparison. Statistical analyses are conducted to identify the significance of the dependence of properties on wall thickness and Hot Isostatic Pressing, as well as to examine correlations of these properties to section area, porosity, and surface roughness. A thorough microstructural study is complemented with a first-of-its-kind study of surface morphology to decouple their contributions and identify underlying causes for observed changes in mechanical properties. This thesis finds that mechanical properties in the quasi-static and fatigue framework do not see appreciable declines until specimen thickness is under 0.75 mm in thickness. The added Hot Isostatic Pressing heat treatment effectively closed pores, recrystallized the grain structure, and provided a more homogenous microstructure that benefits the modulus, tensile strength, elongation, and fatigue performance at higher stresses. Stress heterogeneities, primarily caused by surface defects, negatively affected the thinner specimens disproportionately. Without the use of the Hot Isostatic Pressing, the grain structure remained much more refined and benefitted the yield strength and fatigue endurance limit.
ContributorsParadise, Paul David (Author) / Bhate, Dhruv (Thesis advisor) / Chawla, Nikhilesh (Committee member) / Azeredo, Bruno (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2022
Description
Intelligent engineering designs require an accurate understanding of material behavior, since any uncertainties or gaps in knowledge must be counterbalanced with heightened factors of safety, leading to overdesign. Therefore, building better structures and pushing the performance of new components requires an improved understanding of the thermomechanical response of advanced materials under service conditions. This dissertation provides fundamental investigations of several advanced materials: thermoset polymers, a common matrix material for fiber-reinforced composites and nanocomposites; aluminum alloy 7075-T6 (AA7075-T6), a high-performance aerospace material; and ceramic matrix composites (CMCs), an advanced composite for extreme-temperature applications. To understand matrix interactions with various interfaces and nanoinclusions at their fundamental scale, the properties of thermoset polymers are studied at the atomistic scale. An improved proximity-based molecular dynamics (MD) technique for modeling the crosslinking of thermoset polymers is carefully established, enabling realistic curing simulations through its ability to dynamically and probabilistically perform complex topology transformations. The proximity-based MD curing methodology is then used to explore damage initiation and the local anisotropic evolution of mechanical properties in thermoset polymers under uniaxial tension with an emphasis on changes in stiffness through a series of tensile loading, unloading, and reloading experiments. Aluminum alloys in aerospace applications often require a fatigue life of over 109 cycles, which is well over the number of cycles that can be practically tested using conventional fatigue testing equipment. In order to study these high-life regimes, a detailed ultrasonic cycle fatigue study is presented for AA7075-T6 under fully reversed tension-compression loading. The geometric sensitivity, frequency effects, size effects, surface roughness effects, and the corresponding failure mechanisms for ultrasonic fatigue across different fatigue regimes are investigated. Finally, because CMCs are utilized in extreme environments, oxidation plays an important role in their degradation. A multiphysics modeling methodology is thus developed to address the complex coupling between oxidation, mechanical stress, and oxygen diffusion in heterogeneous carbon fiber-reinforced CMC microstructures.
ContributorsSchichtel, Jacob (Author) / Chattopadhyay, Aditi (Thesis advisor) / Dai, Lenore (Committee member) / Ghoshal, Anindya (Committee member) / Huang, Huei-Ping (Committee member) / Jiao, Yang (Committee member) / Oswald, Jay (Committee member) / Arizona State University (Publisher)
Created2022
Description
High-entropy alloys (HEAs) is a new class of materials which have been studied heavily due to their special mechanical properties. HEAs refers to alloys with multiple equimolar or nearly equimolar elements. HEAs show exceptional and attractive properties currently absent from conventional alloys, which make them the center of intense investigation. HEAs obtain their properties from four core effects that they exhibit and most of the work on them have been dedicated to study their mechanical properties. In contrast, little or no research have gone into studying the functional or even thermal properties of HEAs. Some HEAs have also shown exceptional or very high melting points. According to the definition of HEAs, Si-Ge-Sn alloys with equal or comparable concentrations of the three group IV elements belong to the category of HEAs. Thus, the equimolar components of Si-Ge-Sn alloys probably allow their atomic structures to display the same fundamental effects of metallic HEAs. The experimental fabrication of such alloys has been proven to be very difficult, which is mainly due to differences between the properties of their constituent elements, as indicated from their binary phase diagrams. However, previous computational studies have shown that SiGeSn HEAs have some very interesting properties, such as high electrical conductivity, low thermal conductivity and semiconducting properties. In this work, going for a complete characterization of the SiGeSn HEA properties, the melting point of this alloy is studied using classical molecular dynamics (MD) simulations and density functional theory (DFT) calculations. The aim is to investigate the effects of high Sn content in this alloy on the melting point compared with the traditional SiGe alloys. Classical MD simulations results strongly indicates that none of the available empirical potentials is able to predict accurate or reasonable melting points for SiGeSn HEAs and most of its subsystems. DFT calculations results show that SiGeSn HEA have a melting point which represent the mean value of its constituent elements and that no special deviations are found. This work contributes to the study of SiGeSn HEA properties, which can serve as guidance before the successful experimental fabrication of this alloy.
ContributorsAlqaisi, Ahmad Madhat Odeh (Author) / Hong, Qi-Jun (Thesis advisor) / Zhuang, Houlong (Thesis advisor) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2023
Description
Dissimilar metal joints such as aluminum-steel joints are extensively used in automobile, naval and aerospace applications and these are subjected to corrosive environmental and mechanical loading resulting in eventual failure of the structural joints. In the case of aluminum alloys under aggressive environment, the damage accumulation is predominantly due to corrosion and is accelerated in presence of other metals. During recent years several approaches have been employed to develop models to assess the metal removal rate in the case of galvanic corrosion. Some of these models are based on empirical methods such as regression analysis while others are based on quantification of the ongoing electrochemical processes. Here, a numerical model for solving the Nernst- Planck equation, which captures the electrochemical process, is implemented to predict the galvanic current distribution and, hence, the corrosion rate of a galvanic couple. An experimentally validated numerical model for an AE44 (Magnesium alloy) and mild steel galvanic couple, available in the literature, is extended to simulate the mechano- electrochemical process in order to study the effect of mechanical loading on the galvanic current density distribution and corrosion rate in AE44-mild steel galvanic couple through a multiphysics field coupling technique in COMSOL Multiphysics®. The model is capable of tracking moving boundariesy of the corroding constituent of the couple by employing Arbitrary Langrangian Eulerian (ALE) method.Results show that, when an anode is under a purely elastic deformation, there is no apparent effect of mechanical loading on the electrochemical galvanic process. However, when the applied tensile load is sufficient to cause a plastic deformation, the local galvanic corrosion activity at the vicinity of the interface is increased remarkably. The effect of other factors, such as electrode area ratios, electrical conductivity of the electrolyte and depth of the electrolyte, are studied. It is observed that the conductivity of the electrolyte significantly influences the surface profile of the anode, especially near the junction. Although variations in electrolyte depth for a given galvanic couple noticeably affect the overall corrosion, the change in the localized corrosion rate at the interface is minimal. Finally, we use the model to predict the current density distribution, rate of corrosion and depth profile of aluminum alloy 7075-stainless steel 316 galvanic joints, which are extensively used in maritime structures.
ContributorsMuthegowda, Nitin Chandra (Author) / Solanki, Kiran N (Thesis advisor) / Rykaczewski, Konrad (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2015
Description
Polyurea is a highly versatile material used in coatings and armor systems to protect against extreme conditions such as ballistic impact, cavitation erosion, and blast loading. However, the relationships between microstructurally-dependent deformation mechanisms and the mechanical properties of polyurea are not yet fully understood, especially under extreme conditions. In this work, multi-scale coarse-grained models are developed to probe molecular dynamics across the wide range of time and length scales that these fundamental deformation mechanisms operate. In the first of these models, a high-resolution coarse-grained model of polyurea is developed, where similar to united-atom models, hydrogen atoms are modeled implicitly. This model was trained using a modified iterative Boltzmann inversion method that dramatically reduces the number of iterations required. Coarse-grained simulations using this model demonstrate that multiblock systems evolve to form a more interconnected hard phase, compared to the more interrupted hard phase composed of distinct ribbon-shaped domains found in diblock systems. Next, a reactive coarse-grained model is developed to simulate the influence of the difference in time scales for step-growth polymerization and phase segregation in polyurea. Analysis of the simulated cured polyurea systems reveals that more rapid reaction rates produce a smaller diameter ligaments in the gyroidal hard phase as well as increased covalent bonding connecting the hard domain ligaments as evidenced by a larger fraction of bridging segments and larger mean radius of gyration of the copolymer chains. The effect that these processing-induced structural variations have on the mechanical properties of the polymer was tested by simulating uniaxial compression, which revealed that the higher degree of hard domain connectivity leads to a 20% increase in the flow stress. A hierarchical multiresolution framework is proposed to fully link coarse-grained molecular simulations across a broader range of time scales, in which a family of coarse-grained models are developed. The models are connected using an incremental reverse–mapping scheme allowing for long time scale dynamics simulated at a highly coarsened resolution to be passed all the way to an atomistic representation.
ContributorsLiu, Minghao (Author) / Oswald, Jay (Thesis advisor) / Muhich, Christopher (Committee member) / Jiang, Hanqing (Committee member) / Peralta, Pedro (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2020
Description
Accurate knowledge and understanding of thermal conductivity is very important in awide variety of applications both at microscopic and macroscopic scales. Estimation,however varies widely with respect to scale and application. At a lattice level, calcu-lation of thermal conductivity of any particular alloy require very heavy computationeven for a relatively small number of atoms. This thesis aims to run conventionalmolecular dynamic simulations for a particular supercell and then employ a machinelearning based approach and compare the two in hopes of developing a method togreatly reduce computational costs as well as increase the scale and time frame ofthese systems. Conventional simulations were run using interatomic potentials basedon density function theory-basedab initiocalculations. Then deep learning neuralnetwork based interatomic potentials were used run similar simulations to comparethe two approaches.
ContributorsDabir, Anirudh (Author) / Zhuang, Houlong (Thesis advisor) / Nian, Qiong (Committee member) / Jiao, Yang (Committee member) / Arizona State University (Publisher)
Created2021