Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 469 records · Page 26

Towards Efficient Uncertainty estimation in deep learning for robust energy prediction in crystal materials

In material science, recent studies have started to explore the potential of using deep learning to improve property prediction from high-fidelity simulations, e.g, density functional theory (DFT). However, the design spaces are sometimes too large and intractable to sample completely. This results in a critical question that is how to evaluate the confidence and robustness of the prediction. In this paper, we propose an efficient approach to estimate uncertainty in deep learning using a single forward pass and then apply it for robust prediction of the total energy in crystal lattice structures. Our approach is built upon the deep kernel learning (DKL) that originally introduces to leverage the expressiveness of deep neural networks as input with a probabilistic prediction of Gaussian processes (GPs) as output. Existing DKL methods have difficulties in the accuracy of predictive uncertainty, training stability, and scaling to large datasets, which lead to significant barriers in real-world applications. We propose to address these challenges by using an inducing point approximate GP in feature space combined with spectral normalization as a regularization. We finally demonstrate our robust performance on an artificial example and a real-world application from materials chemistry.

Bi, Sirui↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd [1] and associated criteria (KnGLL > 0.05) can be accurately used to predict zones of breakdown in the flow, and associated errors for surface quantities. Furthermore, when studying reacting flow, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown in minimal (see Fig. 1). Finally, we studied the impact of the grid refinement on the surface quantities in DSMC.

DSMC↗

Comparison Between DSMC and CFD for Hypersonic Planetary Entry Simulations

Hypersonic planetary entry flows span a wide range of Knudsen numbers between rarefied and continuum flows. While computational fluid dynamics (CFD) techniques cannot provide an accurate solution for flows in the rarefied regime, the direct simulation Monte Carlo (DSMC) method is capable of providing accurate solutions for flows in both in the rarefied and continuum regimes but becomes prohibitively expensive as the Knudsen number decreases. For the purpose of thermal protection systems (TPS) design and post-flight reconstruction, various selected points along an entry trajectory are often solved using hypersonic solvers. The quantities of interest that are obtained from that exercise are generally surface quantities, such as pressure, heat flux and enthalpy. Then, material response solvers are used to either design the heat shield to an optimal thickness based on a choice of material, or to provide in-depth heating profiles through the material at various select locations, and compare with flight instrumentation such as the ones that flew on NASA’s two most recent Mars missions, MSL and Mars2020. While most of the heating is generally experienced during the continuum part of the entry, the heating within the rarefied regime is significant for some atmospheres, and hence the flow solutions need to be computed using the DSMC method. Ensuring consistency between hypersonic CFD and the DSMC is crucial so that reliable surface quantities can be passed to material response solvers. Studies were performed to compare the two methods at various select locations, for both non-reacting argon flows as well reacting CO2/N2 flows. Preliminary conclusions show that, for non-reacting flows, the agreement between the two methods for surface heating is excellent (within expected uncertainties) for a freestream Knudsen number of 0.0006, and gets progressively worse as the Knudsen number increases to 0.06. Continuum breakdown analyses were performed and showed that, in general, the Gradient Length Local Knudsen number (KnGLL) from Boyd and associated criterion (KnGLL > 0.05) seems conservative in predicting zones of breakdown in the flow, and associated errors for surface quantities. Updated criteria of KnGLL = 2.0 and 0.5 appear to be more appropriate for surface and flow quantities, respectively. Furthermore, when studying reacting flows, our studies showed that while flow quantities are highly dependent on relaxation parameters and chemistry rates, it is possible to obtain a good agreement for surface heating, as long as the continuum breakdown is minimal.

DSMC↗

Electron Microscopy Approaches to Unraveling the Structure of Amorphous Materials

Determining atomic structures in crystalline materials—where atoms are arranged in rigid, periodic lattices—has been highly successful using probes such as electrons, X-rays, and neutrons. In contrast, amorphous materials, despite their ubiquity and technological importance, remain far more challenging to characterize with comparable accuracy and precision. This review highlights existing, emerging, and potential (scanning) transmission electron microscopy ((S)TEM) techniques for probing short- and medium-range order in amorphous materials. Approaches ranging from high-resolution (S)TEM imaging and selected electron diffraction pattern to four-dimensional STEM (4D-STEM) based pair distribution function, fluctuation electron microscopy, tomography, ptychography, and spectroscopic methods are discussed, emphasizing their ability to provide complementary insights across multiple length scales—from sub-angstrom local environments to nanometer-scale correlations. Here, we further explore the promise of multimodal and correlative strategies, as well as the growing role of machine learning and physics-informed AI in enabling real-time, quantitative interpretation of complex structural signatures. Together, these advances point toward a future where electron microscopy not only reveals the hidden order in amorphous systems but also establishes robust structure–property relationships, paving the way for materials innovation in disordered matter.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Raman structural studies of the nickel electrode

Raman spectroscopy is sensitive to empirically controlled nickel electrode structural variations, and has unique potential for structural characterization of these materials. How the structure relates to electrochemical properties is examined so that the latter can be more completely understood, controlled, and optimized. Electrodes were impregnated and cycled, and cyclic voltammetry is being used for electrochemical characterization. Structural variation was observed which has escaped detection using other methods. Structural changes are induced by: (1) cobalt doping, (2) the state of change or discharge, (3) the preparation conditions and type of buffer used, and (4) the formation process. Charged active mass has an NiOOH-type structure, agreeing with X-ray diffraction results. Discharged active mass, however, is not isostructural with beta-Ni(OH)2. Chemically prepared alpha phases are not isostructural either. A disordered structural model, containing point defects, is proposed for the cycled materials. This model explains K(+) incorporation. Band assignments were made and spectra interpreted for beta-Ni(OH)2, electrochemical NiOOH and chemically precipitated NiOOH.

Cornilsen, B. C.↗

Simultaneous elastic and photoelastic calibration of birefringent orthotropic model materials

Transmission photoelastic analysis of composite models has been shown in recent years to be a viable stress analysis method for composite structures. However, difficulties in the fabrication of the model materials require the development of efficient mechanical and optical calibration procedures. In this paper, the methods for measuring elastic and photoelastic constants of composites are briefly reviewed. Then a new method, utilizing a relatively small calibration specimen, is described. In this method, electrical resistance strain gauges are mounted at selected points of a half-plane model and photoelastic measurements and strain readings are simultaneously taken. The data are processed by a least-squares procedure to yield the elastic and photoelastic constants. Comparison of results obtained by the proposed method with the results from conventional tests shows good agreement.

Prabhakaran, R.↗

Fracture of composites using singular finite elements

Six noded quarter point natural isoparametric triangular elements are employed to obtain displacement and stress distributions around the crack tip in unidirectional graphite epoxy. The performance of the quarter point element for predicting crack growth direction in unidirectional composite materials is studied for different loading conditions and fiber orientations. The finite element results are compared to analytical solutions of anisotropic elasticity where available. Notched Iosipescu shear specimens are analyzed to study fracture behavior under shear loading. Two methods for determining stress intensity factors are presented.

Choksi, G. N.↗

Variable Depth Bragg Peak Method for Single Event Effects Testing

Traditionally, accelerator SEE testing is accomplished by removing the tops of packages so that the IC chips are accessible to heavy ions. However, ICs in some advanced packages cannot be de-lidded so a different approach is used that involves grinding and/or chemically etching away part of the package and the chip from the back side. The parts are then tested from the back side with ions having sufficient range to reach the sensitive volume. More recently, the entire silicon substrate in an SOI/SRAM was removed, making it possible to use low-energy ions with shorter ranges. Where removal of part of the package is not possible, facilities at Michigan State, NASA Space Radiation Laboratory, GANIL (France) and GSI (Germany) offer high-energy heavy ions with long ranges so that the ions can reach the devices' sensitive volumes without much change in the LET. Unfortunately, a run will typically involve only one ion species having a single energy and LET due to the long time it takes to tune a new energy. The Variable Depth Bragg Peak (VDBP) method is similar to the above method in that it involves the use of high-energy heavy ions that are able to pass through the packaging material and reach the device, obviating the need to remove the package. However, the method provides a broad range of LETs from a single ion by inserting degraders in the beam that modify the ion energy and, therefore, the LET. The crux of the method involves establishing a fiduciary point for degrader thickness, i.e., where the Bragg peak is located precisely at the sensitive volume in the device, for which the measured SEU cross-section and the ion LET are both also maxima and can be calculated using a Monte-Carlo program, TRIM. Once the fiduciary point has been established, calibrated high density polyethylene (HDPE) degraders are inserted into or removed from the beam to vary the ion LET at the device in a known manner. After each change of degrader thickness, the SEU cross-section is measured and the corresponding LET calculated from the change in degrader thickness. That information is used to generate a plot of cross-section as a function of ion LET. The advantages of this approach are that the part does not have to be de-lidded and a broad range of LETs is available from a single heavy ion without having to go to non-normal angles of incidence to change the "effective" LET. As we will show, it is possible to obtain an entire curve of cross-section versus LET using just two or three ions. Fig. 1 shows curves of cross-section vs LET for a Freescale 4 Mbit SOI/SRAM measured at the 88" Cyclotron at Berkeley and at NSRL. The open symbols are the data obtained from Berkeley for top-side and back-side irradiation. The solid data points are for the data obtained at NSRL using a device for which the package was intact. The data are for Iron and Gold and cover a range of LETs from 4 MeV.cm2/mg to 84 MeV.cm2/mg. The agreement between the data obtained from Berkeley and from NSRL is excellent, demonstrating that the VDBP method is capable of providing accurate values of cross-section versus LET, at least for the 4 Mbit SRAM. Details of the technique will be included in the final presentation.

Buchner, S.↗

Abstract for CRADA between NETL and the University of Massachusetts Lowell (AGMT-1275)

The demand for printed electronics over conventional electronics is growing each year due to their low costs, low material waste during manufacturing, and compatibility with flexible substrates. The development of inks for printed electronics is a challenge because each printing technology requires different viscosities and curing behaviors to create devices with controlled structures at the micron scale, or smaller. Carbon is a great candidate for fabricating printed electronics devices because it is inexpensive and earth abundant, easily functionalized to impart miscibility with solvent systems, and has electrical properties that can be tuned from insulating to conducting. Despite these advantages, carbons are rarely used for printed electronics because there is a scarcity of inks available for non-contact printing technologies. To address these challenges, the National Energy Technology Laboratory (NETL) will collaborate with the University of Massachusetts Lowell (Participant) through its Printed Electronics Research Collaborative (PERC) to develop engineered carbons and ink formulations optimized specifically for printing radio- and microwave-frequency electronic devices using ink jet, aerosol jet, syringe dispensing, and other additive manufacturing techniques. NETL will process coal, petroleum, and other carbon feedstocks to make engineered carbons with the appropriate composition, microstructure, surface functionalization, and particle morphologies needed to stabilize ink suspensions and impart useful electronic properties to printed devices. The University of Massachusetts Lowell will utilize these engineered carbons to determine which solvent systems and ink additives are needed to achieve the appropriate viscosity, stability, and printability to make useful ink formulations. The University of Massachusetts Lowell will then utilize these inks to manufacture printed resistors and conductors using inkjet, aerosol jet, and syringe dispensing printing methods and will characterize device performance using four-point probe resistivity measurements, contact angle measurements to assess wettability, adhesion evaluation to various substrates, and focused ion beam and scanning electron microscopy for surface and density investigations of printed devices.

36 MATERIALS SCIENCE↗

A methodology for small scale rural land use mapping in semi-arid developing countries using orbital imagery. Part 5: Experimental and operational techniques of mapping land use

The author has identified the following significant results. Scope of the preprocessing techniques was restricted to standard material from the EROS Data Center accompanied by some enlarging procedures and the use of the diazo process. Investigation has shown that the most appropriate sampling strategy for this study is the stratified random technique. A viable sampling procedure, together with a method for determining minimum number of sample points in order to test results of any interpretation are presented.

Vangenderen, J. L.↗

Sandwich beam compressive test method

Tests presented show that the sandwich beam in four-point bending could be used to obtain reliable compressive elastic constants for graphite/polyimide laminates although some difficulties were encountered measuring composite compressive strengths. Data were obtained for ultimate stress, ultimate strain, Young's modulus, and Poisson's ratio values in 117K, room temperature, and 589K test environments for the HTS1/PMR-15 material system. A total of 36 compressive tests on various laminates were performed. Also, 24 tensile tests on the same laminate orientations were performed to obtain input to the analytical portion of this program. A portion of the beam test section was analyzed using a linear elastic finite element computer program to predict the stress state in the graphite/polyimide composite. The influence of the honeycomb core of this stress state was analytically determined.

Shuart, M. J.↗

FEM/SINDA: Combining the strengths of NASTRAN, SINDA, I-DEAS, and PATRAN for thermal and structural analysis

This paper describes the interface/integration between FEM/SINDA, a general purpose geometry driven thermal analysis code, and the FEM software: I-DEAS, PATRAN, and NASTRAN. FEM/SINDA brings together the advantages of the finite element method to model arbitrary geometry and anisotropic materials and SINDA's finite difference capability to model thermal properties, loads, and boundary conditions that vary with time or temperature. I-DEAS and PATRAN thermal entities are directly supported since FEM/SINDA uses the nodes of the FEM model as the point at which the temperature is determined. Output from FEM/SINDA (as well as the FEM/SINDA input deck) can be used directly by NASTRAN for structural analysis.

Zarda, P. Richard↗

A Dividend in Food Safety

When NASA faced the problem of how and what to feed an astronaut in a sealed capsule under weightless conditions while planning for manned space mission, they enlisted the aid of The Pillsbury Company. There were two principal concerns: barring crumbs of food that might contaminate the spacecraft's atmosphere or float their way into sensitive instruments; and assuring absolute freedom from potentially catastrophic disease-producing bacteria and toxins. Pillsbury quickly solved the first concern, but the other part of the problem was not as easy. They found that with using standard methods, there was no way to be assured there would not be any bacteria. It was concluded that the only way to succeed was to establish control over the entire process, the raw materials, the processing environment and the people involved. Pillsbury developed the Hazard Analysis and Critical Control Point (HACCP) concept. The HACCP is designed to prevent food safety problems rather than to catch them after they have occurred. Three other government agencies are taking preliminary steps toward extending HACCP to meat/poultry and seafood inspection operations. Today, Pillsbury plants are still operating under HACCP.

Source record↗

Health Monitoring of Composite Material Structures Using a Vibrometry Technique

Non-destructive evaluation (NDE) methods for quantifying and locating damage are essential for inspecting structures to ensure safety and reliability. Transmittance function monitoring is a potentially new NDE technique being tested as a tool to detect, quantify, and locate damage on flexible structures. The technique has a large spatial range that is practical for detecting damage on large composite material structures such as a reusable launch vehicle. The Transmittance Function (TF) theory is based on structural dynamics principles that define how vibration at one point in a structure is related to a force at another point. This relationship is called the Frequency Response Function (FRF). A Transmittance Function (TF) is derived as the ratio of FRFs, and can detect damage because the FRFs change due to damage. If one excitation is used for the testing, the force does not need to be measured to compute the TF. In the damage detection procedure, the structure is subjected to wide-band vibration and TFs are computed between different accelerometers to detect changes in the structure, presumably due to damage. In the first year of the project the TF method was tested on a bolted panel, a curved panel, and beams, all made of fiberglass. It was shown that damage could be detected using low frequency vibration, 250 to 1,250 Hz. The technique is sensitive to damage, but it requires storage of historical or pre-damage TFs for the healthy structure. This would become a large data storage requirement for large structures. Thus one objective for the second year of the project was to eliminate the need to store historical data. The second year report gives details of how storage of historical data was eliminated. Further results of testing panel structures are also given.

Schulz, Mark J.↗

Machine learning the metastable phase diagram of covalently bonded carbon

Abstract Conventional phase diagram generation involves experimentation to provide an initial estimate of the set of thermodynamically accessible phases and their boundaries, followed by use of phenomenological models to interpolate between the available experimental data points and extrapolate to experimentally inaccessible regions. Such an approach, combined with high throughput first-principles calculations and data-mining techniques, has led to exhaustive thermodynamic databases (e.g. compatible with the CALPHAD method), albeit focused on the reduced set of phases observed at distinct thermodynamic equilibria. In contrast, materials during their synthesis, operation, or processing, may not reach their thermodynamic equilibrium state but, instead, remain trapped in a local (metastable) free energy minimum, which may exhibit desirable properties. Here, we introduce an automated workflow that integrates first-principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the metastable phases to construct “metastable” phase diagrams for materials far-from-equilibrium. Using carbon as a prototypical system, we demonstrate automated metastable phase diagram construction to map hundreds of metastable states ranging from near equilibrium to far-from-equilibrium (400 meV/atom). We incorporate the free energy calculations into a neural-network-based learning of the equations of state that allows for efficient construction of metastable phase diagrams. We use the metastable phase diagram and identify domains of relative stability and synthesizability of metastable materials. High temperature high pressure experiments using a diamond anvil cell on graphite sample coupled with high-resolution transmission electron microscopy (HRTEM) confirm our metastable phase predictions. In particular, we identify the previously ambiguous structure of n -diamond as a cubic-analog of diaphite-like lonsdaelite phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-modal Free-moving Data Fusion (MFDF) v1.0

Multi-modal Free-moving Data Fusion (MFDF) is a python package that implements low and high-level functionality for performing qualitative and quantitative gamma-ray imaging analyses. Such analyses enable the use of gamma-ray spectrometers and/or imagers to estimate the distribution and quantity of radiological materials in an environment. Free-moving 3D imaging requires additional information about the 3D trajectory and orientation of the system, however all of these methods can be applied to 2D static gamma-ray imaging as well. The package includes methods (MLEM and MAP) for distributed source reconstruction, methods (PSL) for point source reconstruction, the ability to ingest and utilize quantitative detector response functions (for absolute analyses), and the ability to perform 3D estimation of dose-rates from quantitative reconstructions.

Joshi, Tenzing↗

Comparison of the surface charge behavior of commercial silicon nitride and silicon carbide powders

The adsorption and desorption of protons from aqueous solution onto the surfaces of a variety of commercial silicon carbide and silicon nitride powders has been examined using a surface titration methodology. This method provides information on some colloidal characteristics, such as the point of zero charge (pzc) and the variation of proton adsorption with dispersion pH, useful for the prediction of optimal ceramic-processing conditions. Qualitatively, the magnitude of the proton adsorption from solution reveals small differences among all of the materials studied. However, the results show that the pzc for the various silicon nitride powders is affected by the powder synthesis route. Complementary investigations have shown that milling can also act to shift the pzc exhibited by silicon nitride powder. Also, studies of the role of the electrolyte in the development of surface charge have indicated no evidence of specific adsorption of ammonium ion on either silicon nitride or silicon carbide powders.

Whitman, Pamela K.↗

Determination of the structural damping coefficients of six full-scale helicopter rotor blades of different materials and methods of construction

Results of an experimental investigation of the structural damping of six full-scale helicopter rotor blades, made to determine the variation of structural damping with materials and methods of construction, are presented. The damping of the blades was determined for the first three flapwise bending modes, first chordwise bending mode, and first torsion mode. The contribution of structural damping to the total damping of the blades is discussed for several aerodynamic conditions in order to point out situations where structural damping is significant.

WINGS, ROTATING - EXPERIMENTAL STUDIES↗