Search NASA⌕ Search

SEARCH · Search NASA

Results for “material point method”

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 361 records · Page 20

Fire safety experiments on MIR Orbital Station

The process of heterogeneous combustion of most materials under zero-g without forced motion of air is practically impossible. However, ventilation is required to support astronauts' life and cool equipment. The presence of ventilation flows in station compartments at accidental ignition can cause a fire. An additional, but exceedingly important parameter of the fire risk of solid materials under zero-g is the minimum air gas velocity at which the extinction of materials occurs. Therefore, the conception of fire safety can be based on temporarily lowering the intensity of ventilation and even turning it off. The information on the limiting conditions of combustion under natural conditions is needed from both scientific and practical points of view. It will enable us to judge the reliability of results of ground-based investigations and develop a conception of fire safety of inhabited sealed compartments of space stations to by provided be means of nontraditional and highly-effective methods without both employing large quantities of fire-extinguishing compounds and hard restrictions on use of polymers. In this connection, an experimental installation was created to study the process of heterogeneous combustion of solid non-metals and to determine the conditions of its extinction under microgravity. This installation was delivered to the orbital station 'Mir' and the cosmonauts Viktorenko and Kondakova performed initial experiments on it in late 1994. The experimental installation consists of a combustion chamber with an electrical systems for ignition of samples, a device for cleaning air from combustion products, an air suction unit, air pipes and a control panel. The whole experiment is controlled by telemetry and recorded with two video cameras located at two different places. Besides the picture, parameters are recorded to determine the velocity of the air flow incoming to the samples, the time points of switching on/off the devices, etc. The combustion chamber temperature is also controlled. The main objectives of experiments of this series were as follows: (1) verification of the reliability of the installation in orbital flight; (2) verification of the experimental procedure; and (3) investigation of combustion of two types of materials under microgravity at various velocities of the incoming air flow.

Egorov, S. D.↗

Enhancing Cluster Identification in Atom Probe Tomography Data Using Transfer Learning

Atom Probe Tomography (APT) is a powerful technique for visualizing the atomic-scale distribution of solutes in materials, but quantitative cluster analysis of APT datasets remains a challenge due to the need for subjective parameter selection in clustering algorithms. While distance-based and density-based methods such as HDBSCAN are widely used, their performance is highly sensitive to user-defined parameters, which undermines reproducibility and accuracy. This study proposes an image-based, deep learning-aided workflow for automating parameter selection and cluster detection in APT data analysis. By projecting 3D APT point clouds onto 2D planes, we leverage pretrained convolutional neural networks (ConvNeXt-Tiny and ResNet-50) through transfer learning to predict the number of clusters present in synthetic datasets. The output is used to guide K-means clustering and estimate HDBSCAN parameters, specifically minimum cluster size and minimum sample points. This approach reduces reliance on manual parameter tuning, improving consistency and scalability. The methodology demonstrates the feasibility of using image-based deep learning for interpreting complex spatial patterns in APT data, enabling faster and more objective analysis. The complete workflow and code are made publicly available to support reproducibility and future research.

Density-based clustering↗

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗

Uncovering multiscale structure-property correlations via active learning in scanning tunneling microscopy

Atomic arrangements and local sub-structures fundamentally influence emergent material functionalities. These structures are conventionally probed using spatially resolved studies and the property correlations are deciphered by a researcher based on sequential explorations, thereby limiting the efficiency and scope. Here we demonstrate a multi-scale Bayesian deep-learning based framework that automatically correlates material structure with its electronic properties using scanning tunneling microscopy (STM) measurements in real-time. Its predictions are used to autonomously direct exploration toward regions of the sample that optimize a given material property. This method is deployed on a low-temperature ultra-high vacuum STM to understand the structure-property relationship in a europium-based semimetal, EuZn 2 As 2 , a promising candidate relevant to magnetism-driven topological phenomena. The framework employs a sparse-sampling approach to efficiently construct the scalar-property space using minimal measurements, about 1–10% of the data required in standard hyperspectral methods. Moreover, we formulate the problem hierarchically across length scales, implementing autonomous workflow to locate mesoscopic and atomic structures that correspond to a target material property. This framework offers the choice to design scalar-property from the spectroscopic data to steer sample exploration. Our findings reveal correlations of the electronic properties unique to surface terminations, local defect density, and point defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Automated optimization techniques for aircraft synthesis

Application of numerical optimization techniques to automated conceptual aircraft design is examined. These methods are shown to be a general and efficient way to obtain quantitative information for evaluating alternative new vehicle projects. Fully automated design is compared with traditional point design methods and time and resource requirements for automated design are given. The NASA Ames Research Center aircraft synthesis program (ACSYNT) is described with special attention to calculation of the weight of a vehicle to fly a specified mission. The ACSYNT procedures for automatically obtaining sensitivity of the design (aircraft weight, performance and cost) to various vehicle, mission, and material technology parameters are presented. Examples are used to demonstrate the efficient application of these techniques.

Vanderplaats, G. N.↗

Radial basis function network learns ceramic processing and predicts related strength and density

Radial basis function (RBF) neural networks were trained using the data from 273 Si3N4 modulus of rupture (MOR) bars which were tested at room temperature and 135 MOR bars which were tested at 1370 C. Milling time, sintering time, and sintering gas pressure were the processing parameters used as the input features. Flexural strength and density were the outputs by which the RBF networks were assessed. The 'nodes-at-data-points' method was used to set the hidden layer centers and output layer training used the gradient descent method. The RBF network predicted strength with an average error of less than 12 percent and density with an average error of less than 2 percent. Further, the RBF network demonstrated a potential for optimizing and accelerating the development and processing of ceramic materials.

Cios, Krzysztof J.↗

Failure of Harmonic Gears During Verification of a Two-Axis Gimbal for the Mars Reconnaissance Orbiter Spacecraft

The Mars Reconnaissance Orbiter (MRO) spacecraft has three two-axis gimbal assemblies that support and move the High Gain Antenna and two solar array wings. The gimbal assemblies are required to move almost continuously throughout the mission's seven-year lifetime, requiring a large number of output revolutions for each actuator in the gimbal assemblies. The actuator for each of the six axes consists of a two-phase brushless dc motor with a direct drive to the wave generator of a size-32 cup-type harmonic gear. During life testing of an actuator assembly, the harmonic gear teeth failed completely, leaving the size-32 harmonic gear with a maximum output torque capability less than 10% of its design capability. The investigation that followed the failure revealed limitations of the heritage material choices that were made for the harmonic gear components that had passed similar life requirements on several previous programs. Additionally, the methods used to increase the stiffness of a standard harmonic gear component set, while accepted practice for harmonic gears, is limited in its range. The stiffness of harmonic gear assemblies can be increased up to a maximum stiffness point that, if exceeded, compromises the reliability of the gear components for long life applications.

Johnson, Michael R.↗

Method of making macrocrystalline or single crystal semiconductor material

A macrocrystalline or single crystal semiconductive material is formed from a primary substrate including a single crystal or several very large crystals of a relatively low melting material. This primary substrate is deposited on a base such as steel or ceramic, and it may be formed from such metals as zinc, cadmium, germanium, aluminum, tin, lead, copper, brass, magnesium silicide, or magnesium stannide. These materials generally have a melting point below about 1000 C and form on the base crystals the size of fingernails or greater. The primary substrate has an epitaxial relationship with a subsequently applied layer of material, and because of this epitaxial relationship, the material deposited on the primary substrate will have essentially the same crystal size as the crystals in the primary substrate. If required, successive layers are formed, each of a material which has an epitaxial relationship with the previously deposited layer, until a layer is formed which has an epitaxial relationship with the semiconductive material. This layer is referred to as the epitaxial substrate, and its crystals serve as sites for the growth of large crystals of semiconductive material. The primary substrate is passivated to remove or otherwise convert it into a stable or nonreactive state prior to deposition of the seconductive material.

Shlichta, P. J.↗

Sudbury project (University of Muenster-Ontario Geological Survey): Sr-Nd in heterolithic breccias and gabbroic dikes

One major objective of our Sudbury project was to define origin and age of the huge breccia units below and above the Sudbury Igneous Complex (SIC). The heterolithic Footwall Breccia (FB) represents a part of the uplifted crater floor. It contains subrounded fragments up to several meters in size and lithic fragments with shock features (greater than 10 GPa) embedded into a fine- to medium-grained matrix. Epsilon(sub Nd)-epsilon(sub Sr) relationships point to almost exclusively parautochthonous precursor lithologies. The different textures of the matrix reflect the metamorphic history of the breccia layer; thermal annealing by the overlying hot impact melt sheet (SIC) at temperatures greater than 1000 C resulted in melting of the fine crushed material, followed by an episode of metasomatic K-feldspar growth and, finally, formation of low-grade minerals such as actinolite and chlorite. Isotope relationships in the Onaping breccias (Gray and Green Member) are much more complex. All attempts to date the breccia formation failed: Zircons are entirely derived from country rocks and lack the pronounced Pb loss caused by the heat of the slowly cooling impact melt sheet (SIC). Rb-Sr techniques using either lithic fragments of different shock stages or the thin slab method, set time limits for the apparently pervasive alkali mobility in these suevitic breccias. The data array and the intercept in the plots point to a major Rb-Sr fractionation around 1.54 Ga ago. This model age is in the same range as the age obtained for the metasomatic matrix of the FB. Rb-Sr dating of a shock event in impact-related breccias seems to be possible only if their matrix had suffered total melting by the hot melt sheet (FB) or if they contain a high fraction of impact melt (suevitic Onaping breccias), whereas the degree of shock metamorphism in rock or lithic fragments plays a minor role. In the Sudbury case, however, the impact melt in the seuvitic breccias is devitrified and recrystallized, which changed Rb/Sr ratios quite drastically. Therefore, the Onaping breccias give only age limits for alteration and low-grade metamorphism. The Sm-Nd system was not reset during the Sudbury event; clasts as well as the matrix in the FB and in the Onaping breccias show preimpact 'Archean' Nd isotope signatures.

Buhl, D.↗

Unzipping hBN with ultrashort mid-infrared pulses

Manipulating the nanostructure of materials is critical for numerous applications in electronics, magnetics, and photonics. However, conventional methods such as lithography and laser writing require cleanroom facilities or leave residue. We describe an approach to creating atomically sharp line defects in hexagonal boron nitride (hBN) at room temperature by direct optical phonon excitation with a mid-infrared pulsed laser from free space. We term this phenomenon “unzipping” to describe the rapid formation and growth of a crack tens of nanometers wide from a point within the laser-driven region. Formation of these features is attributed to the large atomic displacement and high local bond strain produced by strongly driving the crystal at a natural resonance. This process occurs only via coherent phonon excitation and is highly sensitive to the relative orientation of the crystal axes and the laser polarization. Its cleanliness, directionality, and sharpness enable applications such as polariton cavities, phonon-wave coupling, and in situ flake cleaving.

36 MATERIALS SCIENCE↗

Fracture mechanics; Proceedings of the Seventeenth National Symposium, Albany, NY, August 7-9, 1984

The present conference gives attention to topics in the application of fracture mechanics, subcritical crack growth phenomena, fracture testing methods, ductile fracture behavior, and fracture mechanisms and their analysis. Specific papers treat the resistance curve approach to composite materials characterization, fracture toughness in ductile iron and cast steel, hold-time effects in elevated temperature fatigue crack propagation, creep crack growth under nonsteady conditions, viscoplastic fatigue in a superalloy at elevated temperatures, fracture testing with arc bend specimens, one-point bend impact test application, and a compact mode II fracture specimen. Also discussed are the computation of stable crack growth using the J-integral, the use of plastic energy dissipation to characterize crack growth, the extension of surface cracks under cyclic loading, the minimum time criterion for crack instability in structural materials, dynamic crack propagation and branching under biaxial loading, and boundary layer effects in cracked bodies.

Underwood, J. M.↗

Improved silicon carbide for advanced heat engines

This is the third annual technical report for the program entitled, Improved Silicon Carbide for Advanced Heat Engines, for the period February 16, 1987 to February 15, 1988. The objective of the original program was the development of high strength, high reliability silicon carbide parts with complex shapes suitable for use in advanced heat engines. Injection molding is the forming method selected for the program because it is capable of forming complex parts adaptable for mass production on an economically sound basis. The goals of the revised program are to reach a Weibull characteristic strength of 550 MPa (80 ksi) and a Weibull modulus of 16 for bars tested in 4-point loading. Two tasks are discussed: Task 1 which involves materials and process improvements, and Task 2 which is a MOR bar matrix to improve strength and reliability. Many statistically designed experiments were completed under task 1 which improved the composition of the batches, the mixing of the powders, the sinter and anneal cycles. The best results were obtained by an attritor mixing process which yielded strengths in excess of 550 MPa (80 ksi) and an individual Weibull modulus of 16.8 for a 9-sample group. Strengths measured at 1200 and 1400 C were equal to the room temperature strength. Annealing of machined test bars significantly improved the strength. Molding yields were measured and flaw distributions were observed to follow a Poisson process. The second iteration of the Task 2 matrix experiment is described.

Whalen, Thomas J.↗

Improved boundary layer heat transfer calculations near a stagnation point

A thermal design of a solar receiver has been developed for the solutions of problems involving phase-change thermal energy storage and natural convection loss. Two dimensional axisymmetrical solidification and melting of materials contained between two concentric cylinders of finite length has been studied for thermal energy storage analysis. For calculation of free convection loss inside receiver cavity, two dimensional axisymmetrical, laminar, transient free convection including radiation effects has been studied using integral/finite difference method. Finite difference equations are derived for the above analysis subject to constant or variable material properties, initial conditions, and boundary conditions. The validity of the analyses has been substantiated by comparing results of the present general method with available analytic solutions or numerical results reported in the literature. Both explicit and implicit schemes are tested in phase change analysis with different number of nodes ranging from 4 to 18. The above numerical methods have been applied to the existing solar receiver analyzing computer code as additional subroutines. The results were computed for one of the proposed Brayton cycle solar receiver models running under the actual environmental conditions. Effect of thermal energy storage on the thermal behavior of the receiver has been estimated. Due to the thermal energy storage, about 65% reduction on working gas outlet temperature fluctuation has been obtained; however, maximum temperature of thermal energy storage containment has been increased about 18%. Also, effect of natural convection inside a receiver cavity on the receiver heat transfer has been analyzed. The finding indicated that thermal stratification occurs during the sun time resulting in higher receiver temperatures at the outlet section of the gas tube, and lower temperatures at the inlet section of the gas tube when compared with the results with no natural convection. Due to heat supply from the air during the shade time, minimum temperature has been increased, while maximum temperature has been reduced due to convection loss to air. Consequently, cyclic temperature fluctuation has been reduced 29% for working gas and 16% for thermal energy storage containment. On the other hand, despite the presence of the natural convection the time-averaged temperatures for receiver components were found to be similar for two cases with/without natural convection (maximum difference was 1.8%).

Ahn, Kyung Hwan↗

Investigation into Spectroscopic Techniques for Thermal Barrier Coating Spall Detection

Spectroscopic methods are proposed for detection of thermal barrier coating (TBC) spallation from engine hot zone components. These methods include absorption and emission of airborne marker species originally embedded in the TBC bond coat. In this study, candidate marker materials for this application were evaluated. Thermochemical analysis of candidate marker materials combined with additional constraints such as toxicity and uniqueness to engine environment, provided a short list of four potential species: platinum, copper oxide, zinc oxide. and indium. The melting point of indium was considered to be too low for serious consideration. The other three candidate marker materials, platinum, copper oxide, and zinc oxide were placed in a high temperature furnace and emission and absorption properties were measured over a temperature range from 800-1400 C and a spectral range from 250 to 18000 nm. Platinum did not provide the desired response, likely due to the low vapor Pressure of the metallic species and the low absorption of the oxide species. It was also found, however. that platinum caused a broadening of the carbon dioxide absorption at 4300 nm. The nature of this effect is not known. Absorption and emission caused by sodium and potassium impurities in the platinum were found in the platinum tests. Zinc oxide did not provide the desired response, again, most likely due to the low vapor pressure of the metallic species and the low absorption of the oxide species. Copper oxide generated two strongly temperature dependent absorption peaks at 324.8 and 327.4 nm. The melting point of copper oxide was determined to be too low for serious consideration as marker material.

deGroot, Wim↗

Real-Time Optical Monitoring and Simulations of Gas Phase Kinetics in InN Vapor Phase Epitaxy at High Pressure

Understanding the kinetics of nucleation and coalescence of heteroepitaxial thin films is a crucial step in controlling a chemical vapor deposition process, since it defines the perfection of the heteroepitaxial film both in terms of extended defect formation and chemical integrity of the interface. The initial nucleation process also defines the film quality during the later stages of film growth. The growth of emerging new materials heterostructures such as InN or In-rich Ga(x)In(1-x)N require deposition methods operating at higher vapor densities due to the high thermal decomposition pressure in these materials. High nitrogen pressure has been demonstrated to suppress thermal decomposition of InN, but has not been applied yet in chemical vapor deposition or etching experiments. Because of the difficulty with maintaining stochiometry at elevated temperature, current knowledge regarding thermodynamic data for InN, e.g., its melting point, temperature-dependent heat capacity, heat and entropy of formation are known with far less accuracy than for InP, InAs and InSb. Also, no information exists regarding the partial pressures of nitrogen and phosphorus along the liquidus surfaces of mixed-anion alloys of InN, of which the InN(x)P(1-x) system is the most interesting option. A miscibility gap is expected for InN(x)P(1-x) pseudobinary solidus compositions, but its extent is not established at this point by experimental studies under near equilibrium conditions. The extension of chemical vapor deposition to elevated pressure is also necessary for retaining stoichiometric single phase surface composition for materials that are characterized by large thermal decomposition pressures at optimum processing temperatures.

Dietz, Nikolaus↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

An Approximate Ablative Thermal Protection System Sizing Tool for Entry System Design

A computer tool to perform entry vehicle ablative thermal protection systems sizing has been developed. Two options for calculating the thermal response are incorporated into the tool. One, an industry-standard, high-fidelity ablation and thermal response program was integrated into the tool, making use of simulated trajectory data to calculate its boundary conditions at the ablating surface. Second, an approximate method that uses heat of ablation data to estimate heat shield recession during entry has been coupled to a one-dimensional finite-difference calculation that calculates the in-depth thermal response. The in-depth solution accounts for material decomposition, but does not account for pyrolysis gas energy absorption through the material. Engineering correlations are used to estimate stagnation point convective and radiative heating as a function of time. The sizing tool calculates recovery enthalpy, wall enthalpy, surface pressure, and heat transfer coefficient. Verification of this tool is performed by comparison to past thermal protection system sizings for the Mars Pathfinder and Stardust entry systems and calculations are performed for an Apollo capsule entering the atmosphere at lunar and Mars return speeds.

Dec, John A.↗

Development and Application of Computational/In Vitro Toxicological Methods for Chemical Hazard Risk Reduction of New Materials for Advanced Weapon Systems

The development of quantitative structure-activity relationship (QSAR) is essential for reducing the chemical hazards of new weapon systems. The current collaboration between HEST (toxicology research and testing), MLPJ (computational chemistry) and PRS (computational chemistry, new propellant synthesis) is focusing R&D efforts on basic research goals that will rapidly transition to useful products for propellant development. Computational methods are being investigated that will assist in forecasting cellular toxicological end-points. Models developed from these chemical structure-toxicity relationships are useful for the prediction of the toxicological endpoints of new related compounds. Research is focusing on the evaluation tools to be used for the discovery of such relationships and the development of models of the mechanisms of action. Combinations of computational chemistry techniques, in vitro toxicity methods, and statistical correlations, will be employed to develop and explore potential predictive relationships; results for series of molecular systems that demonstrate the viability of this approach are reported. A number of hydrazine salts have been synthesized for evaluation. Computational chemistry methods are being used to elucidate the mechanism of action of these salts. Toxicity endpoints such as viability (LDH) and changes in enzyme activity (glutahoione peroxidase and catalase) are being experimentally measured as indicators of cellular damage. Extrapolation from computational/in vitro studies to human toxicity, is the ultimate goal. The product of this program will be a predictive tool to assist in the development of new, less toxic propellants.

Frazier, John M.↗