Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomistic computations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations

We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.

Good, Brian S.↗
NASA NTRS1997

Computational Modeling Method for Superalloys

Computer modeling based on theoretical quantum techniques has been largely inefficient due to limitations on the methods or the computer needs associated with such calculations, thus perpetuating the notion that little help can be expected from computer simulations for the atomistic design of new materials. In a major effort to overcome these limitations and to provide a tool for efficiently assisting in the development of new alloys, we developed the BFS method for alloys, which together with the experimental results from previous and current research that validate its use for large-scale simulations, provide the ideal grounds for developing a computationally economical and physically sound procedure for supplementing the experimental work at great cost and time savings.

Bozzolo, Guillermo↗
NASA NTRS2009

A Model for Predicting Thermoelectric Properties of Bi2Te3

A parameterized orthogonal tight-binding mathematical model of the quantum electronic structure of the bismuth telluride molecule has been devised for use in conjunction with a semiclassical transport model in predicting the thermoelectric properties of doped bismuth telluride. This model is expected to be useful in designing and analyzing Bi2Te3 thermoelectric devices, including ones that contain such nano - structures as quantum wells and wires. In addition, the understanding gained in the use of this model can be expected to lead to the development of better models that could be useful for developing other thermoelectric materials and devices having enhanced thermoelectric properties. Bi2Te3 is one of the best bulk thermoelectric materials and is widely used in commercial thermoelectric devices. Most prior theoretical studies of the thermoelectric properties of Bi2Te3 have involved either continuum models or ab-initio models. Continuum models are computationally very efficient, but do not account for atomic-level effects. Ab-initio models are atomistic by definition, but do not scale well in that computation times increase excessively with increasing numbers of atoms. The present tight-binding model bridges the gap between the well-scalable but non-atomistic continuum models and the atomistic but poorly scalable ab-initio models: The present tight-binding model is atomistic, yet also computationally efficient because of the reduced (relative to an ab-initio model) number of basis orbitals and flexible parameterization of the Hamiltonian.

Lee, Seungwon↗
NASA NTRS2003

Heats of Segregation of BCC Binaries from Ab Initio and Quantum Approximate Calculations

We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.

Good, Brian S.↗
NASA NTRS2004

Heats of Segregation of BCC Binaries from ab Initio and Quantum Approximate Calculations

We compare dilute-limit heats of segregation for selected BCC transition metal binaries computed using ab initio and quantum approximate energy methods. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent LMTO-based parameters. Quantum approximate segregation energies are computed with and without atomistic relaxation, while the ab initio calculations are performed without relaxation. Results are discussed within the context of a segregation model driven by strain and bond-breaking effects. We compare our results with full-potential quantum calculations and with available experimental results.

Good, Brian S.↗
NASA NTRSDate not supplied

Molten Salt Electrolytes: Computational Analysis of Transport, Electrochemical Properties and Designing New Mixtures

Electric aircraft propulsion has gained a traction over the last decade due to possible high-energy density electrochemistries with reasonable cycle life and identification of non-flammable chemistries that can guarantee much better safety. Commercial Li-ion batteries cannot reach such high capacities because of weight limitations that arise from the widely used intercalation electrodes. Also, use of organic electrolytes make them highly susceptible to fire on exposure to air and humidity. Currently, use of Lithium metal anode along with conversion electrochemistries such as Li-O2 and Li-S are being pursued to achieve such high energy densities. Safer inorganic solid-state electrolytes, recently discovered Water-in-Salt electrolytes, molten salt electrolytes are some of the alternatives being pursued for a safer/non-flammable battery. Of these safer alternatives, molten-salt electrolytes offer some attractive properties: liquid at operating temperatures resulting in better electrode-wetting, stable interface with Li-metal anodes and good ionic conductivity; their only drawback being high operating temperatures. In this talk, we will discuss some of the computational studies, driven via atomistic simulations, of the properties of molten-salt electrolytes including transport mechanism and interface stability. We will also discuss predicting electrochemical properties of these high-temperature electrolytes. Further, we will discuss a thermodynamic approach to predicting new molten-salt eutectics with lower melting points.

Molten Salt Electrolytes↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
NASA NTRSDate not supplied

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗
NASA NTRS2008

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗
NASA NTRS2008

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.↗
NASA NTRS1984

Computer simulation of surface and film processes

All the investigations which were performed employed in one way or another a computer simulation technique based on atomistic level considerations. In general, three types of simulation methods were used for modeling systems with discrete particles that interact via well defined potential functions: molecular dynamics (a general method for solving the classical equations of motion of a model system); Monte Carlo (the use of Markov chain ensemble averaging technique to model equilibrium properties of a system); and molecular statics (provides properties of a system at T = 0 K). The effects of three-body forces on the vibrational frequencies of triatomic cluster were investigated. The multilayer relaxation phenomena for low index planes of an fcc crystal was analyzed also as a function of the three-body interactions. Various surface properties for Si and SiC system were calculated. Results obtained from static simulation calculations for slip formation were presented. The more elaborate molecular dynamics calculations on the propagation of cracks in two-dimensional systems were outlined.

Tiller, W. A.↗
NASA NTRS2019

Computational Materials Techniques for Thermal Protection Solutions: Materials and Process Design

Integrated computational materials techniques that span the atomistic and continuum scales have the potential to aid the design and manufacturing of thermal protection materials. Two cases demonstrating the practical application of these methods are discussed. Case one examines the selection of a high temperature coating for carbon/carbon, with the target application being a solar thermal propulsion heat exchanger. The performance of various refractory metal and metal-carbide coatings is characterized considering extreme thermal (3500 K) and chemical (hydrogen flows) conditions. The recession rate, hydrogen leakage, and likelihood of mechanical failure are characterized and provide directions for further experimental investigation. Case two examines the process optimization of a heat shield material composed of a woven silica fiber preform and cyanate ester resin. Frequently, internal voids were found to be present in this composite after the resin infusion and curing stages of manufacturing. Using the manufacturing conditions, computations indicate that both water adsorption and resin cure shrinkage are contributing factors to void formation. The results suggest an alternative process configuration for curing that would mitigate voids.

thermal protection↗
NASA NTRS2018

Computational Materials Techniques for Thermal Protection Solutions: Materials and Process Design

Integrated computational materials techniques that span the atomistic and continuum scales have the potential to aid the design and manufacturing of thermal protection materials. Two cases demonstrating the practical application of these methods are discussed. Case one examines the selection of a high temperature coating for carbon/carbon, with the target application being a solar thermal propulsion heat exchanger. The performance of various refractory metal and metal-carbide coatings is characterized considering extreme thermal (3500 K) and chemical (hydrogen flows) conditions. The recession rate, hydrogen leakage, and likelihood of mechanical failure are characterized and provide directions for further experimental investigation. Case two examines the process optimization of a heat shield material composed of a woven silica fiber preform and cyanate ester resin. Frequently, internal voids were found to be present in this composite after the resin infusion and curing stages of manufacturing. Using the manufacturing conditions, computations indicate that both water adsorption and resin cure shrinkage are contributing factors to void formation. The results suggest an alternative process configuration for curing that would mitigate voids.

Haskins, Justin↗
NASA NTRS2019

Computational Materials Techniques for Thermal Protection Solutions: Materials and Process Design

Integrated computational materials techniques that span the atomistic and continuum scales have the potential to aid the design and manufacturing of thermal protection materials. Two cases demonstrating the practical application of these methods are discussed. Case one examines the selection of a high temperature coating for carbon/carbon, with the target application being a solar thermal propulsion heat exchanger. The performance of various refractory metal and metal-carbide coatings is characterized considering extreme thermal (3500 degrees Kelvin) and chemical (hydrogen flows) conditions. The recession rate, hydrogen leakage, and likelihood of mechanical failure are characterized and provide directions for further experimental investigation. Case two examines the process optimization of a heat shield material composed of a woven silica fiber preform and cyanate ester resin. Frequently, internal voids were found to be present in this composite after the resin infusion and curing stages of manufacturing. Using the manufacturing conditions, computations indicate that both water adsorption and resin cure shrinkage are contributing factors to void formation. The results suggest an alternative process configuration for curing that would mitigate voids.

Thermal Protection↗
NASA NTRS2002

Atomistic Modeling of Semiconductors: Si, C, and 3C-SiC

An ongoing task of the Computational Materials Group (CMG) at the NASA Glenn Research Center is to enhance the role of atomistic simulations based on quantum-approximate methods in the study of new materials and their properties. One of the main goals of the activity continues to be breaching limitations that arise from the natural balance between accuracy, range of application, and computational simplicity. Whether that balance can be maintained while breaking new ground depends on the methods available with a minimum of constraints and limitations for the study of the energetics of arbitrary systems. The main tool used in CMG research, the Bozzolo- Ferrante-Smith (BFS) method for alloys, has no inherent constraint in its formulation, a feature that has allowed for successful research on various topics. In this article, we report on the latest development of the CMG program, namely, the extension and application of the BFS method to compound semiconductors, a departure from our previous research based primarily on metallic alloys.

Bozzolo, Guillermo↗
NASA NTRS2003

Atomistic Modeling of Nanostructures via the BFS Quantum Approximate Method

Ideally, computational modeling techniques for nanoscopic physics would be able to perform free of limitations on the type and number of elements, while providing comparable accuracy when dealing with bulk or surface problems. Computational efficiency is also desirable, if not mandatory, for properly dealing with the complexity of typical nano-strucured systems. A quantum approximate technique, the BFS method for alloys, which attempts to meet these demands, is introduced for the calculation of the energetics of nanostructures. The versatility of the technique is demonstrated through analysis of diverse systems, including multi-phase precipitation in a five element Ni-Al-Ti-Cr-Cu alloy and the formation of mixed composition Co-Cu islands on a metallic Cu(III) substrate.

Bozzolo, Guillermo↗
NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗
NASA NTRS2004

adwTools Developed: New Bulk Alloy and Surface Analysis Software for the Alloy Design Workbench

A suite of atomistic modeling software, called the Alloy Design Workbench, has been developed by the Computational Materials Group at the NASA Glenn Research Center and the Ohio Aerospace Institute (OAI). The main goal of this software is to guide and augment experimental materials research and development efforts by creating powerful, yet intuitive, software that combines a graphical user interface with an operating code suitable for real-time atomistic simulations of multicomponent alloy systems. Targeted for experimentalists, the interface is straightforward and requires minimum knowledge of the underlying theory, allowing researchers to focus on the scientific aspects of the work. The centerpiece of the Alloy Design Workbench suite is the adwTools module, which concentrates on the atomistic analysis of surfaces and bulk alloys containing an arbitrary number of elements. An additional module, adwParams, handles ab initio input for the parameterization used in adwTools. Future modules planned for the suite include adwSeg, which will provide numerical predictions for segregation profiles to alloy surfaces and interfaces, and adwReport, which will serve as a window into the database, providing public access to the parameterization data and a repository where users can submit their own findings from the rest of the suite. The entire suite is designed to run on desktop-scale computers. The adwTools module incorporates a custom OAI/Glenn-developed Fortran code based on the BFS (Bozzolo- Ferrante-Smith) method for alloys, ref. 1). The heart of the suite, this code is used to calculate the energetics of different compositions and configurations of atoms.

Bozzolo, Guillermo↗