Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-linear Simulation”

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.

514 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods

Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Robust optimization of flexible diafiltration systems for critical mineral separations

This paper provides major contributions in expanding the literature for membrane process design with critical mineral recovery applications and showcasing the importance of robust design techniques for reducing risks of underperformance in such systems. Here, a membrane process flowsheet featuring PrOMMiS membrane models for recovering lithium/cobalt from spent batteries is showcased, uncertainty in membrane sieving and localized fouling are considered, and robust designs are obtained using the PyROS toolset. This paper is intended for a general audience of researchers working in critical minerals, membranes, and optimization related areas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

X-ray linear dichroic tomography of crystallographic and topological defects

Abstract The functionality of materials is determined by their composition 1–4 and microstructure, that is, the distribution and orientation of crystalline grains, grain boundaries and the defects within them 5,6 . Until now, characterization techniques that map the distribution of grains, their orientation and the presence of defects have been limited to surface investigations, to spatial resolutions of a few hundred nanometres or to systems of thickness around 100 nm, thus requiring destructive sample preparation for measurements and preventing the study of system-representative volumes or the investigation of materials under operational conditions 7–15 . Here we present X-ray linear dichroic orientation tomography (XL-DOT), a quantitative, non-invasive technique that allows for an intragranular and intergranular characterization of extended polycrystalline and non-crystalline 16 materials in three dimensions. We present the detailed characterization of a polycrystalline sample of vanadium pentoxide (V 2 O 5 ), a key catalyst in the production of sulfuric acid 17 . We determine the nanoscale composition, microstructure and crystal orientation throughout the polycrystalline sample with 73 nm spatial resolution. We identify and characterize grains, as well as twist, tilt and twin grain boundaries. We further observe the creation and annihilation of topological defects promoted by the presence of volume crystallographic defects. The non-destructive and spectroscopic nature of our method opens the door to operando combined chemical and microstructural investigations 11,18 of functional materials, including energy, mechanical and quantum materials.

Science & Technology - Other Topics↗
DOE OSTI2024

Gyromagnetic nonlinear transmission line for radio frequency signal generation and pulse compression

Disclosed are non-linear transmission lines using ferromagnetic materials to generate ferromagnetic resonance oscillations. In one aspect, a non-linear transmission line apparatus is disclosed. The apparatus includes an outer conductor having a first side and a second internally facing side, and an inner conductor positioned internal to the non-linear transmission line apparatus. The apparatus further includes a ferromagnetic material surrounding the inner conductor, wherein the ferromagnetic material comprises nanoparticles of an ε-polymorph of iron oxide expressed as ε-Fe2O3. The apparatus also includes a first dielectric material positioned between the outer conductor and the inner conductor, the dielectric material in contact with both the ferromagnetic material and with the second internally facing side of the outer conductor, wherein the outer conductor, the inner conductor, the dielectric material and the ferromagnetic material form the nonlinear transmission line.

Schneider, Joseph Devin↗
DOE OSTI2026

Rapid Coal-Ash Characterization using Geophysical Methods & Machine Learning

Coal combustion products (CCP) are challenging to delineate in heterogeneous field settings. Conventional methods (test pits, coring, and laboratory analyses) are labor-intensive, slow, invasive, and provide sparse spatial coverage. This study evaluates whether rapid non-invasive geophysical screening methods—induced polarization (IP), magnetic susceptibility, and nuclear magnetic resonance (NMR) —combined with surface colorimetry (RGB_24), can discriminate CCP-soil mixtures and provide reliable estimates of CCP content. Laboratory measurements were collected on five CCP-soil mixtures (series) and modeled using (i) a linear baseline, (ii) a calibrated non-linear (power-mean) model, and (iii) a machine-learning (ML) Random Forest approach, with validation via leave-one-series-out and site-specific tests. Across the five series, individual signals—particularly IP and magnetic susceptibility—were strongly predictive of ash content but were consistently outperformed by combined models. The pooled calibrated non-linear and ML models captured the observed non-linearity and achieved high accuracy and precision, improving on linear fits. Colorimetry showed the weakest direct relationship with ash content for the tested samples but improved performance when included in multi-signal models. At pre-selected 3.5% decision threshold, calibrated and ML approaches yielded near-perfect classification (Matthews correlation coefficient ˜ 1), suggesting strong practical operability for field screening. Additionally, field-analog tests highlighted the role of endmembers—accuracy declined without access to end-member measurements but was largely recovered by collecting a minimal labeled pair for local recalibration. With end members, accuracy remained high. Globally trained models performed well on three operational unknowns; however, series-specific refits provided the most accurate predictions. Overall, these results highlight the potential of combining rapid geophysics and minimal local calibration for improved coal-ash delineation.

Peshtani, Klaudio↗
DOE OSTI2026

Modeling of a Four-Stage Linear Ionization Cooling Channel for a Muon Collider in G4Beamline

A previous study of an eight-stage rectilinear ionization cooling channel in the ICOOL software demonstrated a five-order-of-magnitude reduction in a muon beam’s 6D emittance. In this study, we look to compare the ways ICOOL and Muons, Inc.’s g4Beamline software model ionization cooling by comparing their modeling of the first four stages to this optimized cooling channel constructed in ICOOL. We begin by identifying the parameters used to construct the optimized ionization cooling channel in ICOOL. We then reconstruct this beam in g4Beamline with identical parameter specifications and simulate the cooling of an identical input beam. Finally, we compare the two simulations based on their beam transmission, longitudinal emittance, and transverse emittance along the channel length. Through this process, we demonstrate that G4Beamline accurately reproduces transverse cooling results but predicts systematically different longitudinal emittance evolution while maintaining similar overall cooling performance, reproducing a 97.9% reduction in 6D emittance over four stages.

Keeler, Dominic [Purdue U., West Lafayette] (ORCID↗
DOE OSTI2024

Surface science insight note: A linear algebraic approach to elucidate native films on Fe 3 O 4 surface

Standard materials are often used to obtain spectra that can be compared to those from unknown samples. Spectra measured from these known substances are also used as a means of computing sensitivity factors to allow quantification by X-ray photoelectron spectroscopy (XPS) of less well-defined materials. Spectra from known materials also provide line shapes suitable for inclusion in spectral models which, when fitted to spectra, permit the chemical state for a sample to be assessed. Both types of information depend on isolating photoemission signals from the inelastically scattered signal. In this Insight note, technical issues associated with the use of XPS of as received Fe 3 O 4 powder sample surface are discussed. The Insight note is designed to show how linear algebraic techniques applied to data collected from a sample marketed as pure Fe 3 O 4 powder are used to verify that XPS has been performed on chemistry representative of the sample. The methods described in this Insight note can further be utilized in elucidating complex XPS data obtained from thin films formed or evolved during cyclic/non-steady use of complex (electro)catalyst surfaces, especially in the presence of contaminants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Maximizing dynamic range and performance of anatase TiO 2 ECRAM through structure and programming

Here, in this study, we investigate the structure-dependent modulation characteristics of all-solid-state three-terminal electrochemical random-access memory (ECRAM) based on an anatase Li x TiO 2 channel. By directly comparing “asymmetric” and “symmetric” ECRAM device architectures, we reveal significant insight into the impact of a non-zero gate-drain open-circuit voltage and its influence on voltage vs. current-controlled gating. We also explore the impact of potentiation/depression write parameters on the symmetry, linearity, and dynamic range of the device response. Together, initial results from optimizing structure and programming approaches yielded unprecedented G max /G min ratios of >1,000 for ECRAM and hundreds of tunable memory states with excellent linearity and symmetry. Simulations based on these ECRAM devices further illustrate the promise of this analog memory technology, achieving near 2% classification error in the MNIST digit recognition benchmark for a range of training parameters compared to a theoretical best of 1.66% and outperforming other device models extracted from the literature.

AIHWKit↗
DOE OSTI2026

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2027

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗
DOE OSTI2025

Linear Dichroism of the Optical Properties of SnS and SnSe Van der Waals Crystals

Abstract Tin monochalcogenides SnS and SnSe, belonging to a family of Van der Waals crystals isoelectronic to black phosphorus, are known as environmentally friendly materials promising for thermoelectric conversion applications. However, they exhibit other desired functionalities, such as intrinsic linear dichroism of the optical and electronic properties originating from strongly anisotropic orthorhombic crystal structures. This property makes them perfect candidates for polarization‐sensitive photodetectors working in near‐infrared spectral range. A comprehensive study of the SnS and SnSe crystals is presented, performed by means of optical spectroscopy and photoemission spectroscopy, supported by ab initio calculations. The studies reveal the high sensitivity of the optical response of both materials to the incident light polarization, which is interpreted in terms of the electronic band dispersion and orbital composition of the electronic bands, dictating the selection rules. From the photoemission investigation the ionization potential, electron affinity and work function are determined, which are parameters crucial for the design of devices based on semiconductor heterostructures.

Chemistry↗
DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I↗
DOE OSTI2025

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗