SEARCH · Search NASA

Results for “points”

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.

SourceResearch connectionFrom yearThrough year

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

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗

DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗

DOE OSTI2025

Impact of chemical ordering on thermodynamic properties of point defects and Xe substitutional in U-10Mo

The accurate knowledge of defect energetics is critical to understanding the aging and irradiation behavior of U-10Mo nuclear fuel, which is selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). Furthermore, using hybrid molecular dynamics and Monte Carlo (MDMC) simulation, we studied the impact of chemical ordering on the formation energies of vacancies, interstitials, and the solution energy of the Xe substitutionals. Instead of forming a random solid solution (RSS), substantial short-range-order (SRO) develops in U-10Mo, particularly at low temperatures. Mo atoms are found to repel each other and prefer U-rich local atomic environments within the 1st nearest neighbor (1NN) cutoff. Compared to the case of a RSS, the state with equilibrated Mo ordering shifts the distributions of vacancy and interstitial formation energies due to the dependence of defect energies on the local atomic environment, without a clear effect on Xe solution energy. In the operation temperature range (100–250 °C) of U-Mo fuels, neglecting SRO can lead to an inaccurate estimate of thermal equilibrium point defect concentrations by over an order of magnitude and incorrectly predict the preference among different types of dumbbells, highlighting the critical importance of accounting for the impact of chemical ordering for accurate atomistic calculations of defect properties.

Molecular dynamics and Monte Carlo (MDMC)↗

DOE OSTI2025

Leveraging explainable AI to characterize floating-point exceptions in linear solvers

Linear solver packages are central to many scientific, engineering, and machine learning applications. When floating-point exceptions occur in these solvers, e.g., division by zero or overflow, numerical results are compromised and become unreliable. Existing static and dynamic analysis tools can detect such exceptions, but they do not explain why the exceptions occur in terms of the solver inputs. Here, we present a study to characterize the inputs that cause numerical exceptions in linear solver packages. Our approach uses explainable AI (XAI) to find the most relevant characteristics of input matrices that explain the occurrence of exceptions in the solvers. Since training data in this domain is scarce, we perform extensive data gathering and data augmentation to obtain exception-inducing inputs. Our approach uses a repair strategy on the features blamed by XAI to validate that such features indeed explain the exceptions. We compare the LIME and SHAP XAI techniques using a dozen matrix features with three classifiers. We evaluate the approach on three widely used linear solver packages and find that some input characteristics can explain the occurrence of exceptions 100% of the time, in specific solvers and preconditioners.

Explainable AI↗

DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗

DOE OSTI2026

Insights from Initial Engineering Designs of Point Source Capture at Industrial Facilities

Initial engineering design studies examining the application of state-of-the-art carbon capture technology at industrial plants contain generally overlooked real-world design considerations for near-term deployment of point source capture (PSC). The implementation of PSC across a wide range of industrial applications presents unique challenges associated with fluctuating CO<sub>2</sub> concentrations, flue gas composition, and utility and land availability. In this article, seven recent industrial retrofit PSC projects are reviewed to investigate the impact of site-specific factors on project design and cost. Across the seven projects, three capture technology classes and four industrial applications are considered, allowing insight into industry-specific opportunities for PSC technology synergy. Common challenges across projects and proposed design solutions are highlighted to propagate ideas and solutions to close the technology gaps and accelerate learning rates.

FEED Studies↗

DOE OSTI2025

Bi-stability and period-doubling cascade of frequency combs in exceptional-point lasers

Abstract Recent studies have demonstrated that a laser can self-generate frequency combs when tuned near an exceptional point (EP), where two cavity modes coalesce. These EP combs induce periodic modulation of the population inversion in the gain medium, and their repetition rate is independent of the laser cavity’s free spectral range. In this work, we perform a stability analysis that reveals two notable properties of EP combs, bi-stability and a period-doubling cascade. The period-doubling cascade enables halving of the repetition rate while maintaining the comb’s total bandwidth, presenting opportunities for the design of highly compact frequency comb generators.

exceptional points↗

DOE OSTI2025

Remote sensing images, DEM, and point clouds associated with “Accuracy evaluation of cost-effective 3D reconstruction approaches for hydrobiogeochemical processes in non-perennial stream riverbeds”

This data package is associated with the publication “Accuracy evaluation of cost-effective 3D reconstruction approaches for hydrobiogeochemical processes in non-perennial stream riverbeds” published in Frontiers in Environmental Science, Environmental Informatics and Remote Sensing (Bao et al., 2026; doi: 10.3389/fenvs.2026.1725258). This data package includes the drone photos for a section of Umtanum Creek in Washington, Unted States. The photos were used to reconstruct the 3-dimensional (3D) digital elevation model (DEM) of the riverbed for the investigated stream section. The reconstruction results from four approaches are provided: (1) unoccupied aerial vehicle (UAV, colloquially known as drone) imagery-based Structure-from-Motion (SfM), (2) a machine learning-based 3D reconstruction model, Visual Geometry Grounded Deep Structure from Motion (VGGSfM), (3) Visual Geometry Grounded Transformer for long sequence of images (VGGT-Long), and (4) handheld smartphone LiDAR scanning. The ground truth measurements by tripod-mounted optical level kit and ground control points GPS locations for evaluating the accuracy of the four reconstruction approaches are also provided in this data package. A preliminary version of this data package was published in October 2025 at the time of manuscript submission. It was updated in March 2026, at the time of manuscript acceptance, to include additional metadata (this readme, data dictionary, and file level metadata). The data did not change. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) 8 folders; (2) the detailed flight configuration html files; (3) field metadata; (4) a readme; (5) a data dictionary; and (6) file-level metadata. The folders “2024_10_18_d01” and “2024_10_18_d02” contain the original drone photos for the two drone flights (d01 and d02) on October 18, 2024. The reconstruction results from each of the approaches are in the folders called “ODM_SfM”, “VGGSfM”, “VGGTLong”, and “LiDAR”. The ground truth measurements are in the folder called “optical_level_kit”. Lastly, results comparing the different approaches are in the folder called “comparisons”. All files are .csv, .html, .jpg, .obj, .txt, and .npy. For information on using the .obj and .npy files, see the readme files within the same folder as the files.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

DOE OSTI2024

(Doublon) Benchmarking of Different Inverse Point Kinetics Implementations for an Autocorrected Reactimeter Algorithm

In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed, such as the Transient Heatsink Overpower Response capsule (THOR), the Transient Water Irradiation System for TREAT (TWIST), and Sirius. Each has provided valuable data for materials performance and reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments can cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used for code verification, validation and uncertainty quantification. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the verification, validation and uncertainty quantification approach. To overcome this problem, the autocorrected reactimeter algorithm (ACRA) has been developed. This approach infers a time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e., when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for developing an improved version of the ACRA that does not require frozen configurations. To develop this new algorithm, we need a precise and unbiased implementation of the inverse point kinetic equations (IPKEs) as any error in the reactivity evaluation will be propagated into the choice of the optimal spatial correction. Indeed, the previous reactimeter algorithm would use approximations, such as a negligible flux amplitude derivative, to focus on rapidity. For the numerical validation of ACRA, we aim at absolute error under for reactivity derived from signals similar to the one of this study. In this summary, we test eight different IPKE implementations. Each will process a mockup signal built for this study, similar to those that the future ACRA will process. Each reactivity output will be compared to the reference reactivity that has been used to generate the mockup signal. The implementation minimizing the difference with the reference reactivity will be used in the development of a new ACRA formulation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Atomically Revealing Bulk Point Defect Dynamics in Hydrogen‐Driven γ‐Fe 2 O 3 → Fe 3 O 4 → FeO Transformation

Understanding how point defects in the bulk govern redox transformations is essential for advancing hydrogen-based metal production and designing high-performance oxide materials. This study reveals the atomic-scale mechanisms driving hydrogen-induced reduction of γ-Fe 2 O 3 to Fe 3 O 4 , focusing on how bulk vacancy dynamics dictate structural evolution and reaction kinetics. A key finding is the pronounced contrast in defect behavior between the two oxides: in γ-Fe 2 O 3 , intrinsic Fe vacancies promote oxygen vacancy clustering, destabilizing the local lattice and driving nanopore formation. In contrast, Fe 3 O 4 exhibits a higher oxygen vacancy formation energy and lacks intrinsic Fe vacancies, suppressing vacancy aggregation and maintaining a dense, pore-free structure. This divergence governs distinct reduction pathways—γ-Fe 2 O 3 undergoes an interface-reaction-limited transformation confined to the γ-Fe 2 O 3 /Fe 3 O 4 boundary, while Fe 3 O 4 supports a uniform increase in oxygen vacancy concentration, enabling bulk-phase reduction to lower-oxide FeO. Integrated in situ electron microscopy and density functional theory modeling uncover a vacancy-mediated mechanism, where synergistic cation-anion vacancy dynamics steer microstructure evolution and phase progression. These insights highlight the critical role of vacancy dynamics in controlling oxide reactivity and offer a pathway toward vacancy engineering to enhance reduction kinetics in hydrogen metallurgy and to tailor porosity, reactivity, and structural resilience in oxide-based catalysts and energy materials.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Converting Second‐Order Saddle Points to Transition States: New Principles for the Design of 4π Photoswitches

Abstract Molecular photoswitches have demonstrated potential for storing solar energy at the molecular level, with power densities comparable to commercial batteries and hydroelectric energy storage. However, development of efficient photoswitches is hindered by limitations in cyclability and optical properties of existing materials. We here demonstrate that certain limitations in photoswitches based on electrocyclizations stem from the issue of controlling competition between Woodward‐Hoffmann allowed and forbidden pathways. Our approach moves beyond the traditional view of activation barriers and reveals that second‐order saddle points are crucial in dictating the competition between disrotatory and conrotatory pathways. These insights suggest new opportunities to manipulate the competition between these pathways through geometric constraints, fundamentally altering the connectivity of the potential energy surface. Our study also emphasizes the necessity of multi‐reference methods and the need to conduct higher‐dimensional explorations for competing pathways beyond photoswitch design.

Chemistry↗

DOE OSTI2025

A two-loop four-point form factor at function level

Recently, the maximally-helicity-violating four-point form factor for the chiral stress-energy tensor in planar $\mathcal{N}$ = 4 super Yang-Mills was computed to three loops at the level of the symbol associated with multiple polylogarithms. It exhibits antipodal self-duality, or invariance under the combined action of a kinematic map and reversing the ordering of letters in the symbol. Here we lift the two-loop form factor from symbol level to function level. We provide an iterated representation of the function’s derivatives (coproducts). In order to do so, we find a three-parameter limit of the five-parameter phase space where the symbol’s letters are all rational. We also use function-level information about dihedral symmetries and the soft, collinear, and factorization limits, as well as limits governed by the form-factor operator product expansion (FFOPE). We provide plots of the remainder function on several kinematic slices, and show that the result is compatible with the FFOPE data. We further verify that antipodal self-duality is valid at two loops beyond the level of the symbol.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

The three-point form factor of Tr ϕ 3 to six loops

We study the three-point form factor of the length-three half-BPS operator (Tr ϕ 3 ) in planar $\mathcal{N}$ = 4 Super-Yang-Mills theory, using analyticity and integrability methods. We find that the functions describing the form factor in perturbation theory live in the same restrictive space of multiple polylogarithms as the one describing the form factor of the stress-tensor operator (Tr ϕ 2 ). Furthermore, we find that the leading-order data in the collinear limit provided by the form factor operator product expansion (FFOPE) is enough to fix the form factor uniquely, at least through six loops. We perform various tests of our results using the subleading FFOPE corrections. We also analyze the form factor in the Regge limit where two Mandelstam invariants are large; we obtain a compact representation for the form factor in this limit which is valid to all orders in the coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

DK/Dπ scattering and an exotic virtual bound state at the SU(3) flavour symmetric point from lattice QCD

Abstract ElasticS-wave scattering of a charm meson with a light pseudoscalar meson inJ P = 0 + is investigated in the flavour$$ \overline{\textbf{3}} $$ 3 ¯ ,6and$$ \overline{\textbf{15}} $$ 15 ¯ sectors at the SU(3) f flavour point using lattice QCD, working on three volumes withm π ≈ 700 MeV. Large bases of interpolating operators are employed to extract finite-volume spectra, which are subsequently used with the Lüscher method to provide constraints on infinite-volume scattering amplitudes. Examining the singularities of the amplitudes, theS-wave amplitude in the flavour$$ \overline{\textbf{3}} $$ 3 ¯ sector is found to contain a deeply bound state, strongly coupled to elastic threshold, corresponding to the$$ {J}^P={0}^{+}\ {D}_{s0}^{\ast }(2317) $$ J P = 0 + D s 0 ∗ 2317 . In the exotic flavour6sector a virtual bound state is found at$$ \sqrt{s_{\textrm{pole}}} $$ s pole = 2510 – 2610 MeV, roughly 40 – 140 MeV below threshold, whereas the$$ \overline{\textbf{15}} $$ 15 ¯ channel shows weak repulsion.

Physics↗

DOE OSTI2024

Development of 52Mn Labeled Trastuzumab for Extended Time Point PET Imaging of HER2

Abstract Purpose Due to their long circulation time in the blood, monoclonal antibodies (mAbs) such as trastuzumab, are usually radiolabeled with long-lived positron emitters for the development of agents for Positron Emission Tomography (PET) imaging. Manganese-52 ( 52 Mn, t 1/2 = 5.6 d, β + = 29.6%, E(β ave ) = 242 keV) is suitable for imaging at longer time points providing a complementary technique to Zirconium-89 ( 89 Zr, t 1/2 = 3.3 d, β + = 22.7%, E(β ave ) = 396 keV)) because of its long half-life and low positron energy. To exploit these properties, we aimed to investigate suitable bifunctional chelators that could be readily conjugated to antibodies and labeled with 52 Mn under mild conditions using trastuzumab as a proof-of-concept. Procedures Trastuzumab was incubated with S-2-(4-isothiocyanatobenzyl)-1,4,7,10-tetraazacyclododecane tetraacetic acid (p-SCN-Bn-DOTA), 1-Oxa-4,7,10-tetraazacyclododecane-5-S-(4-isothiocyantobenzyl)-4,7,10-triacetic acid (p-SCN-Bn-Oxo-DO3A), and 3,6,9,15-tetraazabicyclo[9.3.1] pentadeca-1(15),11,13-triene-4-S-(4-isothiocyanatobenzyl)-3,6,9-triacetic acid (p-SCN-Bn-PCTA) at a tenfold molar excess. The immunoconjugates were purified, combined with [ 52 Mn]MnCl 2 at different ratios, and the labeling efficiency was assessed by iTLC. The immunoreactive fraction of the radiocomplex was determined through a Lindmo assay. Cell studies were conducted in HER2 + (BT474) and HER2- (MDA-MB-468) cell lines followed by in vivo studies. Results Trastuzumab-Oxo-DO3A was labeled within 30 min at 37 °C with a radiochemical yield (RCY) of 90 ± 1.5% and with the highest specific activity of the chelators investigated of 16.64 MBq/nmol. The labeled compound was purified with a resulting radiochemical purity of > 98% and retained a 67 ± 1.2% immunoreactivity. DOTA and PCTA immunoconjugates resulted in < 50 ± 2.5% (RCY) with similar specific activity. Mouse serum stability studies of [ 52 Mn]Mn-Oxo-DO3A-trastuzumab showed 95% intact complex for over 5 days. Cell uptake studies showed higher uptake in HER2 + (12.51 ± 0.83% /mg) cells compared to HER2- (0.85 ± 0.10%/mg) cells. PET images of mice bearing BT474 tumors showed high tumor uptake that was consistent with the biodistribution (42.02 ± 2.16%ID/g, 14 d) compared to MDA-MB-468 tumors (2.20 ± 0.80%ID/g, 14 d). Additionally, both models exhibited low bone uptake of < 1% ID/g. Conclusion The bifunctional chelator p-SCN-Bn-Oxo-DO3A is promising for the development of 52 Mn radiopharmaceuticals as it was easily conjugated, radiolabeled at mild conditions, and illustrated stability for a prolonged duration both in vitro and in vivo . High-quality PET/CT images of [ 52 Mn]Mn-Oxo-DO3A-trastuzumab were obtained 14 d post-injection. This study illustrates the potential of [ 52 Mn]Mn-Oxo-DO3A for the evaluation of antibodies using PET imaging.

Omweri, James M.↗

DOE OSTI2024

Three-dimensional continuum point cloud method for large deformation and its verification

This study presents a strong form based meshfree collocation method, which is named Continuum Point Cloud Method, to solve nonlinear field equations derived from classical mechanics for deformed bodies in three-dimensional Euclidean space. The method and its implementation are benchmarked against a nonlinear vector field using manufactured solutions. The analysis of mechanical fields firstly focuses on the study of St. Venant Kirchhoff and compressible neo-Hookean materials. Results for various initial boundary value problems are presented, including benchmark cases involving unidirectional tension and simple shear. Subsequently, the study concludes with an analysis of a displacement-controlled simulation of a compressible neo-Hookean material, specifically a bar that is pulled to 50% of its original length and rotated 90°. The pure tension case yields a 1.5% error in displacement between computed and expected values and a combined tension and torsion loading case provides further insight into material behavior under complex loading conditions. The resulting normal axial and transverse stress-strain curves are also presented. Lastly, the consistency and robustness of the proposed nonlinear numerical schemes are successfully demonstrated through various numerical experiments.

Compressible neo-Hookean materials↗

DOE OSTI2025

A direct-adjoint approach for material point model calibration with application to plasticity

Here, this paper proposes a new approach for the calibration of material parameters in local elastoplastic constitutive models. The calibration is posed as a constrained optimization problem, where the constitutive model evolution equations for a single material point serve as constraints. The objective function quantifies the mismatch between the stress predicted by the model and corresponding experimental measurements. To improve calibration efficiency, a novel direct-adjoint approach is presented to compute the Hessian of the objective function, which enables the use of second-order optimization algorithms. Automatic differentiation is used for gradient and Hessian computations. Two numerical examples are employed to validate the Hessian matrices and to demonstrate that the Newton–Raphson algorithm consistently outperforms gradient-based algorithms such as L-BFGS-B.

36 MATERIALS SCIENCE↗