Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete Element 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.

700 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Multilevel Conditional Disturbance Rejection Control for Satellite Attitude Tracking

Recently, the conditional disturbance rejection controller (CDRC) was proposed to improve control performance by leveraging disturbances that have beneficial effects on the system. However, it only considers disturbances acting on the state variable directly influenced by them. Although fast convergence of this state can be achieved with the CDRC, it may unintentionally affect the convergence of the output (i.e., the primary state). Here, in this article, a multilevel CDRC is proposed to enhance satellite attitude control performance by accounting for the effect of disturbances on both attitude (output) and angular velocity. The extended state observer is employed to estimate the lumped disturbance, including the modeling errors and external disturbances. Then, a backstepping-based controller with the multilevel disturbance rejection law (ML-DRL) is designed for attitude tracking. The ML-DRL is developed to improve the control performance by using a disturbance with a damping effect on both attitude and velocity. Faster convergence of attitude and velocity can be achieved by conditionally compensating for the disturbance. The stability of the proposed control method is analyzed by demonstrating that the errors are bounded as time tends to infinity. The attitude control performance of the proposed method is evaluated through numerical examples conducted using the MATLAB/Simulink Multibody tool.

Active disturbance rejection control (ADRC)↗
NASA NTRS2021

Mathematical Characterization of Battery Models

The purpose of this document is to demonstrate the use of the Extended Kalman Filter as a tool for battery state estimation and the estimation of battery state of charge. The mathematical details based on the equivalent circuit model are presented followed by an electrochemical engineering model. A simplified first-order model is used to demonstrate the procedure followed by second and third-order models. Next a simplified electrochemistry model is presented along with observer development. State observability is calculated for the simpler equivalent circuit models and the simplified electrochemistry model. An outline of the battery model parameter identification method is presented, and model performance based on experimental and flight data is demonstrated.

Battery↗
DOE OSTI2024

Uniaxial stress effect on the electronic structure of quantum materials

Uniaxial stress has proven to be a powerful experimental tuning parameter for effectively controlling lattice, charge, orbital, and spin degrees of freedom in quantum materials. In addition, its ability to manipulate the symmetry of materials has garnered significant attention. Recent technical progress to combine uniaxial stress cells with quantum oscillation and angle-resolved photoemission techniques allowed to study the electronic structure as function of uniaxial stress. This review provides an overview on experimental advancements in methods and examines studies on diverse quantum materials, encompassing the semimetal WTe2, the unconventional superconductor Sr2RuO4, Fe-based superconductors, and topological materials.

FOS: Physical sciences↗
NASA NTRS2026

Dynamic Regulation of Sub-Atmospheric Pressure for Constant and Cyclic Gas Loads During Testing of Spacesuit Components

New iterations of the various subsystems within the spacesuit will benefit from a new sorbent technology. For example, sorbents used in the Trace Contaminant Control (TCC) and the Rapid Cycle Amine (RCA) systems within the Exploration Portable Life Support System (xPLSS), part of the Extra-Vehicular Mobility Unit (xEMU). For proper validation, the sub-atmospheric pressure needs to be maintained within a simulated 2 ft 3 spacesuit volume. The traditional approach, which involves placing the entire vent loop within a hypobaric chamber, is not practical for the widespread testing of components and prototypes, as these specialized chambers are costly and not widely available. A two-stage regulator based sub-atmospheric pressure system is presented in this work. The method regulates the pressure based on pressure differentials between the pressure regulation pump, the test system, and the ambient environment. The performance of this system is demonstrated using data collected during a 150+ hour non-regenerative TCC sorbent evaluation, several 24-hour regenerative TCC sorbent evaluations with different regeneration cycles, and an 8-hour RCA sorbent evaluation, all under xEMU operating conditions. For regulation, the inlet regulator was set to propagate a small leak to increase the stability of the system as the gas load from the testing system changes. It was also found that controlling the input flow rate that replaces the material lost during regeneration is critical for ensuring the stability of the system. This system provided excellent pressure regulation, without adjustments, for constant loads during shorter time evaluations (hours), while longer time evaluations (days) are easily obtainable with periodic regulator adjustments. A second method is under development to address longer-term stability and to automate pressure regulation. It will incorporate flow and pressure measurements to dynamically control a set of electronically controlled proportional valves for adjusting the pumping speed and air bleed supply.

Nicholas F Materer↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
NASA NTRSDate not supplied

Dry Pressed, High Areal Loading Electrode Architectures Enabled by Holey Graphene

For future electric aviation, advanced battery cell chemistry beyond lithium ion batteries are required to meet mission requirements. High energy density battery concepts such as lithium-sulfur (Li-S) and lithiumoxygen (Li-O2) chemistries are being intensively investigated to realize their extraordinary theoretical promise in terms of energy density. Most fabrication methods of cathodes for these novel battery chemistries followed a conventional approach. In this approach, the active material is mixed with a polymer binder and a conductive carbon in a high-boiling organic solvent to form a slurry, followed by casting onto a current collector and solvent evaporation. The process is usually lengthy and poses environmental hazards due to the use of organic solvents.

Lin Yi↗
DOE OSTI2025

SUDA: A SUrface Dust Analyser for Compositional Mapping of the Galilean Moon Europa

Abstract The Surface Dust Analyser (SUDA) is a mass spectrometer onboard the Europa Clipper mission for investigating the surface composition of the Galilean moon Europa. Atmosphereless planetary moons such as the Galilean satellites are wrapped into a ballistic dust exosphere populated by tiny samples from the moon’s surface produced by impacts of fast micrometeoroids. SUDA will measure the composition of such surface ejecta during close flybys of Europa to obtain key chemical signatures for revealing the satellite’s composition such as organic molecules and salts, history, and geological evolution. Because of their ballistic orbits, detected ejecta can be traced back to the surface with a spatial resolution roughly equal to the instantaneous altitude of the spacecraft. SUDA is a Time-Of-Flight (TOF), reflectron-type impact mass spectrometer, optimized for a high mass resolution which only weakly depends on the impact location. The instrument will measure the mass, speed, charge, elemental, molecular, and isotopic composition of impacting grains. The instrument’s small size of$268 ~\mathrm {mm} \times 250 ~\mathrm {mm} \times 171$ 268 mm × 250 mm × 171 $~\mathrm {mm}$ mm , radiation-hard design, and rather large sensitive area of 220 cm 2 matches well the challenging demands of the Clipper mission.

Astronomy & Astrophysics↗
NASA NTRS2026

Pulsed Electroacoustic Analysis Capabilities at NASA Glenn Research Center

The High Voltage Testing Laboratory at NASA Glenn Research Center is equipped with a pulsed electroacoustic (PEA) analyzer to study space charge phenomena in dielectric materials. Dielectric materials such as those used in high voltage insulation systems are of great interest due to their crucial role in withstanding high voltages and avoiding electrical breakdown in electrified aircraft power systems. Current research is underway for the development of novel lightweight electrical insulation and the thorough evaluation of these materials is required. The PEA instrumentation at NASA GRC is a nondestructive method to evaluate the space charge phenomena which can give a unique insight on the feasibility and lifetime of material systems. The instrumentation capability and initial results are reported here.

pulsed electroacoustic↗
DOE OSTI2026

Reduced salen in a monomeric Fe complex

We report the first monomeric formally M I transition metal salen complex without ion pairing. A combination of spectroscopy, crystallography, and computations indicates that reduction of an Fe II salen complex gives substantial reduced ligand character. We also present a computational method for accurately predicting 57 Fe Mössbauer parameters of Fe salen complexes.

Feldman, Dana M. [Yale University, New Haven, CT (↗
DOE OSTI2024

Utilities Perspective on Protection Challenges with High IBR Penetration

Utilities have seen rapid increase in Solar, Wind, and Battery Energy Storage Resources that interconnect to their electric system through Inverters. Inverter-Based Resources (IBRs) have fault current characteristics that are unlike the fault current response of traditional rotating-machine-based generators, which is well known and repeatable. IBR’s non-traditional fault current behavior is due to the IBR control scheme, which is configured to provide a clean AC output but also protect the inverter’s sensitive power electronics devices from damage, one source of which is overcurrent. This results in low fault current magnitude, low or no negative sequence current injection, the variability of sequence component currents, the variability of voltage with respect to current angles, and the lack of inertia. The control scheme also results in a fault current response that can vary between manufactures and between models of the same manufacturer. High penetration of IBRs can adversely affect the protection schemes applied in areas with high penetration of IBRs. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems. This is mainly because the existing protection elements and practices have been designed based on the fault current response of conventional rotating machines. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This presentation identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. Instead of facing on one issue, we are looking broadly on all the challenges that system protection has experienced with high penetration of IBRs. Based on the IBR response from various utilities during real fault events and gathering perspective from different utility SMEs via questionnaire, the presentation summarizes on gathered data and internal experiences.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

The first high-redshift cavity power measurements of cool-core galaxy clusters with the International LOFAR Telescope

Radio-mode feedback associated with the active galactic nuclei (AGNs) at the cores of galaxy clusters injects a large amount of energy into the intracluster medium (ICM), offsetting radiative losses through X-ray emission. This mechanism prevents the ICM from rapidly cooling down and fueling extreme starburst activity as it accretes onto the central galaxies, and it is therefore a key ingredient in the evolution of galaxy clusters. However, the influence and mode of feedback at high redshifts (z ∼ 1) remains largely unknown. Low-frequency sub-arcsecond-resolution radio observations taken with the International LOFAR Telescope have demonstrated their ability to assist X-ray observations with constraining the energy output from the AGNs (or “cavity power”) in galaxy clusters, thereby enabling research at higher redshifts than before. In this pilot project, we tested this hybrid method on a high-redshift (0.6 < z < 1.3) sample of 13 galaxy clusters for the first time with the aim of verifying the performance of this method at these redshifts and providing the first estimates of the cavity power associated with the central AGN for a sample of distant clusters. We were able to detect clear radio lobes in three out of 13 galaxy clusters at redshifts of 0.7 < z < 0.9, and we used these detections in combination with ICM pressures surrounding the radio lobes obtained from standard profiles to calculate the corresponding cavity powers of the AGNs. Combining our results with the literature, the current data appear to suggest that the average cavity power peaked at a redshift ofz ∼ 0.4 and slowly decreases toward higher redshifts. However, we require more and tighter constraints on the cavity volume and a better understanding of our observational systematics to confirm any deviation of the cavity power trend from a constant level.

Astronomy & Astrophysics↗
DOE OSTI2026

EvoDiffMol: evolutionary diffusion framework for 3D molecular design with optimized properties

Designing molecules with specific target properties remains a fundamental challenge in computational chemistry. While existing approaches show promise, most rely on simplified representations like SMILES strings or 2D graphs that lack essential three-dimensional geometric information. We present EvoDiffMol, a computational framework that integrates evolutionary algorithms with three-dimensional diffusion models for property-driven molecular generation. The method operates through adaptive evolutionary optimization, where population-based selection guides the generation process toward desired property landscapes. EvoDiffMol supports both unconstrained molecular design and scaffold-constrained generation that preserves fixed substructures while optimizing complementary regions. Comprehensive evaluation demonstrates exceptional performance, achieving the highest drug-likeness score (0.94) among all compared state-of-the-art methods while maintaining excellent validity, uniqueness, and novelty. Beyond single property optimization, the framework demonstrates flexible multi-property optimization capabilities, simultaneously controlling multiple molecular descriptors including synthetic accessibility, lipophilicity, topological polar surface area, and clinically relevant ADMET properties such as cardiotoxicity (hERG) and intestinal permeability (Caco-2). This adaptability spans from simple descriptors to practical pharmaceutical endpoints without requiring complete model retraining. The framework achieves precise control over target property values, generating molecules with properties closely matching specified targets for both single and multiple descriptors. Scaffold-constrained experiments preserve fixed molecular cores while maintaining effective property optimization. The three-dimensional representation offers advantages in maintaining structural validity during iterative optimization, with potential for geometry-aware applications in materials science and drug discovery.

3D molecular generation↗
NASA NTRS1966

Spacecraft Antenna Pointing

In examining the characteristic inaccuracies of a pre-programmed high-gain antenna pointing system, the effects of spacecraft attitude errors, as embodied in pitch, yaw and roll limit cycle activity, are of major concern. If the pointing accuracy requirements for a particular mission allow the use of a single degree of freedom pointing system, it is doubtful whether a method of correcting antenna position for spacecraft attitude errors would be worthwhile. This is because the pointing error due to fitting a plane curve to a three-dimensional trajectory would probably be at least as large as the additional pointing errors due to attitude drift within the deadbands. Further, trying to correct pointing for pitch, yaw, and roll drifts by adjusting a single hinge angle is relatively ineffective.

ERROR SIGNAL↗
NASA NTRS2025

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path↗
NASA NTRS2004

Energetics of Single Substitutional Impurities in NiTi

Shape-memory alloys are of considerable current interest, with applications ranging from stents to Mars rover components. In this work, we present results on the energetics of single substitutional impurities in B2 NiTi. Specifically, energies of Pd, Pt, Zr and Hf impurities at both Ni and Ti sites are computed. All energies are computed using the CASTEP ab initio code, and, for comparison, using the quantum approximate energy method of Bozzolo, Ferrante and Smith. Atomistic relaxation in the vicinity of the impurities is investigated via quantum approximate Monte Carlo simulation, and in cases where the relaxation is found to be important, the resulting relaxations are applied to the ab initio calculations. We compare our results with available experimental work.

Binary Alloys↗
NASA NTRS1998

High Temperature - Thin Film Strain Gages Based on Alloys of Indium Tin Oxide

A stable, high temperature strain gage based on reactively sputtered indium tin oxide (ITO) was demonstrated at temperatures up to 1050 C. These strain sensors exhibited relatively large, negative gage factors at room temperature and their piezoresistive response was both linear and reproducible when strained up to 700 micro-in/in. When cycled between compression and tension, these sensors also showed very little hysteresis, indicating excellent mechanical stability. Thin film strain gages based on selected ITO alloys withstood more than 50,000 strain cycles of +/- 500 micro-in/in during 180 hours of testing in air at 1000 C, with minimal drift at temperature. Drift rates as low as 0.0009%/hr at 1000 C were observed for ITO films that were annealed in nitrogen at 700 C prior to strain testing. These results compare favorably with state of the art 10 micro-m thick PdCr films deposited by NASA, where drift rates of 0.047%/hr at 1050 C were observed. Nitrogen annealing not only produced the lowest drift rates to date, but also produce the largest dynamic gage factors (G = 23.5). These wide bandgap, semiconductor strain sensors also exhibited moderately low temperature coefficients of resistance (TCR) at temperatures up to 1100 C, when tested in a nitrogen ambient. A TCR of +230 ppm/C over the temperature range 200 C < T < 500 C and a TCR of -469 ppm/C over the temperature range 600 C < T < 1100 C was observed for the films tested in nitrogen. However, the resistivity behavior changed considerably when the same films were tested in oxygen ambients. A TCR of -1560 ppm/C was obtained over the temperature range of 200 C < T < 1100 C. When similar films were protected with an overcoat or when ITO films were prepared with higher oxygen contents in the plasma, two distinct TCR's were observed. At T < 800 C, a linear TCR of -210 ppm/C was observed and at T > 800 C, a linear TCR of -2170 DDm/C was observed. The combination of a moderately low TCR and a relatively large gage factor make these semiconducting oxide films promising candidates for the active strain elements in high temperature thin film strain gages, particularly in applications where static strain measurement is desired.

Otto J Gregory↗
DOE OSTI2026

Vapor-phase pillarization of MXenes for engineering hierarchical interlayer porosity

MXenes, a family of two-dimensional (2D) multilamellar materials, possess excellent thermal and electronic properties for a range of applications. Their use in heterogeneous catalysis, however, is limited by the low surface area resulting from stacked layers. Pillarization with inorganic oxides can create more open, mesoporous MXene structures, improving accessibility for guest species to diffuse, reside or react in the space between 2D layers. A previous liquid-phase pillarization method, however, involves excessive use of solvent-based precursors and multiple processing steps. Here, we report a vapor-phase pillarization (VPP) strategy to introduce pillars, exemplified by silica pillars, with high pillar precursor usage efficiency and a simplified processing workflow. The resulting silica-pillared mesoporous MXene exhibits significantly increased surface area and porosity. These textural properties can be easily tuned by the VPP synthesis conditions. When applied as a ruthenium (Ru) catalyst support for the hydrogenolysis of low-density polyethylene (LDPE), the silica-pillared MXene enabled high Ru dispersion and catalytic activity. This study highlights the potential of the VPP method for engineering mesoporous, 2D MXene materials and demonstrates the effectiveness of mesoporous MXene as a catalyst support in overcoming mass transport and active-site accessibility challenges in heterogeneous catalysis involving bulky substances, such as plastics upcycling.

Luo, Song [University of Delaware, Newark, DE (Uni↗
DOE OSTI2024

Improvements in the utilization of calcium carbonate in promoting sustainability and environmental health

Calcium carbonate (CaCO 3 ) is an incredibly abundant mineral on Earth, with over 90% of it being found in the lithosphere. To address the CO 2 crisis and combat ocean acidification, it is essential to produce more CaCO 3 using various synthetic methods. Additionally, this approach can serve as a substitute for energy-intensive processes like cement production. By doing so, we have the potential to not only reverse the damage caused by climate change but also protect biological ecosystems and the overall environment. The key lies in maximizing the utilization of CaCO 3 in various human activities, paving the way for a more sustainable future for our planet.

Chemistry↗