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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Material point generation from an STL file

In this short note we describe a method for generating material points from a file in the stere olithography (STL) format commonly used in computer aided design (CAD), 3-D printing, and computer aided manufacturing applications. In this file format the surface of an object is described by triangles in three-dimensional space. The coordinates of the triangle vertexes are provided in the file. In this note we assume that the object of interest is inside of the enclosing set of triangles. Similar methods can also be used to generate material points outside of the enclosing triangles.

42 ENGINEERING↗
DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗
NASA NTRS2005

Laser And Nonlinear Optical Materials For Laser Remote Sensing

NASA remote sensing missions involving laser systems and their economic impact are outlined. Potential remote sensing missions include: green house gasses, tropospheric winds, ozone, water vapor, and ice cap thickness. Systems to perform these measurements use lanthanide series lasers and nonlinear devices including second harmonic generators and parametric oscillators. Demands these missions place on the laser and nonlinear optical materials are discussed from a materials point of view. Methods of designing new laser and nonlinear optical materials to meet these demands are presented.

Barnes, Norman P.↗
NASA NTRS2000

Materials Compatibility Testing in Concentrated Hydrogen Peroxide

Materials test methods from the 1960's have been used as a starting point in evaluating materials for today's space launch vehicles. These established test methods have been modified to incorporate today's analytical laboratory equipment. The Orbital test objective was to test a wide range of materials to incorporate the revolution in polymer and composite materials that has occurred since the 1960's. Testing is accomplished in 3 stages from rough screening to detailed analytical tests. Several interesting test observations have been made during this testing and are included in the paper. A summary of the set-up, test and evaluation of long-term storage sub-scale tanks is also included. This sub-scale tank test lasted for a 7-month duration prior to being stopped due to a polar boss material breakdown. Chemical evaluations of the hydrogen peroxide and residue left on the polar boss surface identify the material breakdown quite clearly. The paper concludes with recommendations for future testing and a specific effort underway within the industry to standardize the test methods used in evaluating materials.

Boxwell, R.↗
NASA NTRS2009

Heat rejection sublimator

A sublimator includes a sublimation plate having a thermal element disposed adjacent to a feed water channel and a control point disposed between at least a portion of the thermal element and a large pore substrate. The control point includes a sintered metal material. A method of dissipating heat using a sublimator includes a sublimation plate having a thermal element and a control point. The thermal element is disposed adjacent to a feed water channel and the control point is disposed between at least a portion of the thermal element and a large pore substrate. The method includes controlling a flow rate of feed water to the large pore substrate at the control point and supplying heated coolant to the thermal element. Sublimation occurs in the large pore substrate and the controlling of the flow rate of feed water is independent of time. A sublimator includes a sublimation plate having a thermal element disposed adjacent to a feed water channel and a control point disposed between at least a portion of the thermal element and a large pore substrate. The control point restricts a flow rate of feed water from the feed water channel to the large pore substrate independent of time.

Dingell, Charles W.↗
DOE OSTI2021

Point Adversarial Self-Mining: A Simple Method for Facial Expression Recognition

In this article, we propose a simple yet effective approach, called point adversarial self mining (PASM), to improve the recognition accuracy in facial expression recognition (FER). Unlike previous works focusing on designing specific architectures or loss functions to solve this problem, PASM boosts the network capability by simulating human learning processes: providing updated learning materials and guidance from more capable teachers. Specifically, to generate new learning materials, PASM leverages a point adversarial attack method and a trained teacher network to locate the most informative position related to the target task, generating harder learning samples to refine the network. The searched position is highly adaptive since it considers both the statistical information of each sample and the teacher network capability. Other than being provided new learning materials, the student network also receives guidance from the teacher network. After the student network finishes training, the student network changes its role and acts as a teacher, generating new learning materials and providing stronger guidance to train a better student network. Here, the adaptive learning materials generation and teacher/student update can be conducted more than one time, improving the network capability iteratively. Extensive experimental results validate the efficacy of our method over the existing state of the arts for FER.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fuel Properties of Oxymethylene Ethers with Terminating Groups from Methyl to Butyl

Oxymethylene ethers (OMEs) have been studied as possible additives or replacements for diesel fuels. Typically, studies have considered only methyl-terminated OMEs. Recent structure-property relationship models suggest that extended-alkyl OMEs may provide improvements to many of the properties of methyl-terminated OMEs that make them less suitable as diesel fuel blendstocks. In this work, we describe the synthesis and characterization of 16 different OMEs with methyl, ethyl, propyl, butyl, isopropyl, and isobutyl terminating alkyl groups with varying oxymethylene chain length. Indicated Cetane Number, Lower Heating Value, Flash Point, Density, Viscosity, Vapor Pressure, and Oxidative Stability are tested via ASTM standard methods. Additionally, Water Solubility, Boiling Point, seal material compatibility, and sooting propensity (via the Yield Sooting Index) are measured for these fuels. For diesel compatibility, all tested OMEs except smaller methyl and ethyl OMEs, and the branched isopropyl OME, meet cetane number requirements. Further, extending the alkyl end group increases the heating value, but all OMEs, due to their oxygen content, have heating values less than diesel; despite this, all OMEs show significant reductions in soot production per unit heating value. Only the heaviest OMEs meet diesel viscosity requirements, and most are higher density than diesel. OMEs with larger alkyl groups show the highest stability under accelerated auto-oxidation conditions. Increases in alkyl group length cause order of magnitude reduction in water solubility, from hundreds of g/L for methyl terminated OMEs to hundreds of mg/L for butyl terminated OMEs. Limited seal material testing indicates that PEEK polymers are unaffected by OMEs; while extended alkyl groups may improve compatibility with FKM (Viton), other common elastomers (NBR, silicone) remain incompatible with all tested OMEs. Overall, it is found that methyl-terminated OMEs exhibit the most potential for soot reduction, but OMEs with larger propyl and butyl terminating alkyl groups show improved compatibility with existing diesel systems.

09 BIOMASS FUELS↗
DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗
NASA NTRS1992

Advanced composites structural concepts and materials technologies for primary aircraft structures: Structural response and failure analysis

Non-linear analysis methods were adapted and incorporated in a finite element based DIAL code. These methods are necessary to evaluate the global response of a stiffened structure under combined in-plane and out-of-plane loading. These methods include the Arc Length method and target point analysis procedure. A new interface material model was implemented that can model elastic-plastic behavior of the bond adhesive. Direct application of this method is in skin/stiffener interface failure assessment. Addition of the AML (angle minus longitudinal or load) failure procedure and Hasin's failure criteria provides added capability in the failure predictions. Interactive Stiffened Panel Analysis modules were developed as interactive pre-and post-processors. Each module provides the means of performing self-initiated finite elements based analysis of primary structures such as a flat or curved stiffened panel; a corrugated flat sandwich panel; and a curved geodesic fuselage panel. This module brings finite element analysis into the design of composite structures without the requirement for the user to know much about the techniques and procedures needed to actually perform a finite element analysis from scratch. An interactive finite element code was developed to predict bolted joint strength considering material and geometrical non-linearity. The developed method conducts an ultimate strength failure analysis using a set of material degradation models.

Dorris, William J.↗
NASA NTRS1995

Measurement of Density, Sound Velocity, Surface Tension, and Viscosity of Freely Suspended Supercooled Liquids

Non-contact methods have been implemented in conjunction with levitation techniques to carry out the measurement of the macroscopic properties of liquids significantly cooled below their nominal melting point. Free suspension of the sample and remote methods allow the deep excursion into the metastable liquid state and the determination of its thermophysical properties. We used this approach to investigate common substances such as water, o-terphenyl, succinonitrile, as well as higher temperature melts such as molten indium, aluminum and other metals. Although these techniques have thus far involved ultrasonic, electromagnetic, and more recently electrostatic levitation, we restrict our attention to ultrasonic methods in this paper. The resulting magnitude of maximum thermal supercooling achieved have ranged between 10 and 15% of the absolute temperature of the melting point for the materials mentioned above. The physical properties measurement methods have been mostly novel approaches, and the typical accuracy achieved have not yet matched their standard equivalent techniques involving contained samples and invasive probing. They are currently being refined, however, as the levitation techniques become more widespread, and as we gain a better understanding of the physics of levitated liquid samples.

Supercooling Liquids↗
NASA NTRS1979

Accurate determination of work in three-point bend tests

Article presents procedure where correction curve accounts for coincidental displacement and simplifies data analysis in three point bend test in field of materials testing. Method is applicable to any test in above field regardless of load displacement.

Buzzard, R. J.↗
DOE OSTI2022

Fully self-consistent finite-temperature $GW$ in Gaussian Bloch orbitals for solids

In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc ⁢$GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc⁢ $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Starting-point-independent quantum Monte Carlo calculations of iron oxide

Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. In this study, we examine the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

Neutron and high speed photogrammetric arcjet diagnosis

Two methods for real time internal diagnostics of arcjet engines are described. One method uses cold, thermal, or epithermal neutrons. Cold neutrons are used to detect the presence and location of hydrogenous propellants. Thermal neutrons are used to delineate the edge contours of anode and cathode surfaces and to measure stress/strain. Epithermal neutrons are used to measure temperatures on arcjet surfaces, bulk material temperatures, and point temperatures in bulk materials. It is found that this method, with an exposure time of 10 min, produces at temperature accuracy for W or Re of + or - 2.5 C. The other method uses visible-light high-speed photogrammetry to obtain images of the transient behavior of the arc during start-up and to relate this behavior to electrial supply characteristics such as voltage, current, and ripple.

Stewart, P. A. E.↗
DOE OSTI2022

A Novel Rotating Flexure-Test Technique for Brittle Materials with Circular Geometries

The four-point bend test is one of the simplest and often the preferred flexural strength evaluation method for brittle materials. In this loading mode, fracture often initiates from a critical surface (or subsurface) flaw when subjected to a tensile stress state. However, if the critical flaw exists on the compression side of the test specimen, it may not activate to grow a crack and hence the resulting flexural strength will be higher than the true value. The goal of this study is to measure the true flexural strength of a solid or hollow cylindrical brittle specimen by ensuring that failure occurs at its weakest point by rotating along its longitudinal axis, thereby exposing and activating its critical surface flaw during a four-point bend test. A novel test fixture has been designed and fabricated, and the true flexural strengths of cylindrical brittle and quasi-brittle tubular specimens have been measured and compared to existing experimental data obtained through traditional four-point bend tests. In this work, experimental results showcase the orientation dependance on flexural strength for various materials. Additionally, similarities between experimental findings and those available in literature, including observations of fracture surfaces and relationships between surface roughness and material strength, are discussed.

36 MATERIALS SCIENCE↗