SEARCH · Engineering Papers

Results for “linear algebra”

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.

92 records · Page 2

All sources NASA NTRS DOE OSTI

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 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 OSTI2025

Strong Kitaev Interaction in BaCo 2⁢ (AsO 4 ) 2

The inelastic neutron scattering results and their analysis unequivocally point to a dominant Kitaev interaction in the honeycomb-lattice cobaltate BaCo 2 ⁢(AsO 4 ) 2 . Our anisotropic-exchange model closely describes all available neutron scattering data in the material’s field-polarized phase. Furthermore, the density-matrix renormalization group results for our model are in close accord with the unusual double-zigzag magnetic order and the low in-plane saturation field of BaCo 2 ⁢(AsO 4 ) 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1993

A SIMULINK Environment for Flight Dynamics and Control Analysis - Application to the DHC-2 'Beaver', Part 1. Implementation of a Model Library in SIMULINK

The design of advanced Automatic Aircraft Control Systems (AACS's) can be improved upon considerably if the designer can access all models and tools required for control system design and analysis through a graphical user-interface, from within one software environment. This MSc-thesis presents the first step in the development of such an environment, which is currently being done at the Section for Stability and Control of Delft University of Technology, Faculty of Aerospace Engineering. The environment is implemented within the commercially available software package MATLAB/SIMULINK. The report consists of two parts. Part I gives a detailed description of the AACS design environment. The heart of this environment is formed by the SIMULINK implementation of a nonlinear aircraft model in block-diagram format. The model has been worked out for the old laboratory aircraft of the Faculty, the De Havilland DHC-2 'Beaver', but due to its modular structure, it can easily be adapted for other aircraft. Part I also describes MATLAB programs which can be applied for finding steady-state trimmed-flight conditions and for linearization of the aircraft model, and it shows how the built-in simulation routines of SIMULINK have been used for open-loop analysis of the aircraft dynamics. Apart from the implementation of the models and tools, a thorough treatment of the theoretical backgrounds is presented. Part II of this report presents a part of an autopilot design process for the 'Beaver' aircraft, which clearly demonstrates the power and flexibility of the AACS design environment from part I. Evaluations of all longitudinal and lateral control laws by means of nonlinear simulations are treated in detail. The AACS design environment from part I proved to be a very useful tool for designing the control laws of the 'Beaver' autopilot within a very tight time-schedule. The autopilot design process itself will be used as a guideline for future AACS research at the Faculty of Aerospace Engineering. Flight tests of the 'Beaver' autopilot, done after evaluating the control laws in the SIMULINK package, proved to be quite successful. In the future, the AACS design package will evolve into a standardized, integrated design environment which can be applied to virtually any type of aircraft. The AACS design cycle will be shortened further by developing tools for automatically porting control laws from the MATLAB/SIMULINK environment to a piloted real-time flight simulator and the Flight Control Computers of the aircraft.

Aircraft Modelling↗

NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗

DOE OSTI2025

Erosion of high-Z refractory coatings for helicon plasma source window

Proto-MPEX (Materials Plasma Exposure eXperiment), a linear plasma device (LPD), using a high power radio frequency (RF) (⩾100 kW, 13.56 MHz) helicon plasma source, has suffered RF sheath-induced window erosion. The sputtered impurities from the window surface transport toward the downstream target affecting plasma-material interaction studies. The rectified sheath voltage formed was high enough to cause window erosion by light ions due to its low-Z components (e.g. Si 3 N 4 and AlN). Hence, we are proposing impurity mitigation strategies, which involve the application of a Faraday screen to lower the rectified sheath voltage and the application of high-Z refractory coatings on the existing window plasma-facing surface. In this work, we report the erosion of two different coatings, i.e., tantalum oxide (Ta 2 O 5 ) and hafnium oxide (HfO 2 ) on a silicon nitride (Si 3 N 4 ) window material under conditions of low-energy deuterium (D) ion impact, well below the Ta 2 O 5 sputtering energy threshold (i.e. 250 eV). The test samples were RF-biased and exposed to a high D-ion fluence (∼10 26 m −2 ) in Plasma Interaction with Surface Component Experimental Station (PISCES-A) LPD. The experimentally measured sputtering yields of the high-Z refractory coatings below the sputtering energy threshold of Ta 2 O 5 indicate an increased erosion due to the plasma impurities, especially due to oxygen. Improving the vacuum level could reduce the oxygen impurities and increase the lifespan of the coated helicon window for plasma operation. Post-surface analysis indicates no preferential erosion of oxygen or enrichment of the high-Z surface component. This is likely due to an effective energy transfer from the impurity ion for sputtering of the high-Z component in the coating. With the reduced rectified sheath voltage at the window surface and improved vacuum conditions, the proposed high-Z refractory coatings offer a promising solution for reducing window erosion in future MPEX plasma operations at ORNL.

RF bias↗

DOE OSTI2025

Geospatial modeling of near subsurface temperatures of the contiguous United States for assessment of materials degradation

Abstract Understanding subsurface temperature variations is crucial for assessing material degradation in underground structures. This study maps subsurface temperatures across the contiguous United States for depths from 50 to 3500 m, comparing linear interpolation, gradient boosting (LightGBM), neural networks, and a novel hybrid approach combining linear interpolation with LightGBM. Results reveal heterogeneous temperature patterns both horizontally and vertically. The hybrid model performed best achieving a root mean square error of 2.61 °C at shallow depths (50–350 m). Model performance generally decreased with depth, highlighting challenges in deep temperature prediction. State-level analyses emphasized the importance of considering local geological factors. This study provides valuable insights for designing efficient underground facilities and infrastructure, underscoring the need for depth-specific and region-specific modeling approaches in subsurface temperature assessment.

Science & Technology - Other Topics↗

NASA NTRS1985

Battery-Charge-State Model

Charge-state model for lead/acid batteries proposed as part of effort to make equivalent of fuel gage for battery-powered vehicles. Models based on equations that approximate observable characteristics of battery electrochemistry. Uses linear equations, easier to simulate on computer, and gives smooth transitions between charge, discharge, and recuperation.

Vivian, H. C.↗

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↗

NASA NTRS1994

Elastic-plastic models for multi-site damage

This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.

Actis, Ricardo L.↗

NASA NTRS1995

Transient Studies of G-Induced Capillary Dryout and Rewet

A transient, one-dimensional numerical code is developed to model the liquid motion in an axial groove with square cross section. Axial variation in liquid level, shear stress and heat transfer between the groove wall and the liquid, evaporation and transient body forces are accounted for in the model. Dryout and rewet of the groove are allowed; the front location is determined numerically using conservation of mass and linear extrapolation. Several numerical test results are presented and discussed.

M K Reagan↗

DOE OSTI2025

In‐Plane Anisotropy in van der Waals NiTeSe Ternary Alloy

Abstract The anisotropic properties of materials profoundly influence their electronic, magnetic, optical, and mechanical behaviors and are critical for a wide range of applications. In this study, the anisotropic characteristics of Ni‐based van der Waals materials, specifically NiTe 2 and its alloy NiTeSe, utilizing a combination of comprehensive scanning tunneling microscopy (STM), angle‐resolved photoemission spectroscopy (ARPES), and density functional theory (DFT) calculations, are explored. Unlike 1T‐NiTe 2 , which exhibits trigonal in‐plane symmetry, the substitution of Te with Se in NiTe 2 (resulting in the NiTeSe alloy) induces a pronounced in‐plane anisotropy. This anisotropy is clear in the STM topographs, which reveal a distinct linear order of charge distribution. Corroborating these observations, ARPES measurements and DFT calculations reveal an anisotropic Fermi surface centered at the point, which is notably elongated along the k y direction, leading to directional variations in in‐plane carrier velocities. Consequently, the Fermi velocity is highest along the k x direction where the linear charge distribution aligns in real space and is lowest along the k y direction. These findings offer valuable insights into the tunability of anisotropic properties in ternary transition metal dichalcogenide systems, highlighting their potential applications in the development of anisotropic electronic and optoelectronic devices.

Chemistry↗

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↗

NASA NTRSDate not supplied

Battery Cell-to-Pack Scaling Laws for Electric Aircraft

Battery pack gravimetric energy density is one of the most important, yet often miss-estimated design parameters for sizing all-electric aircraft. Proper accounting for thermal, structural, and operational safety margins are frequently lost when extrapolating performance from the cell level to the aircraft level. This paper summarizes the relevant engineering and certification details needed to better account for the penalties associated when assembling battery packs. The relationship between the cell and pack energy density is not linear, as is often assumed. Furthermore, the relationship varies depending on pack requirements, cell chemistry, and architecture. Parametric, high-fidelity models are used to determine optimal battery pack sizes over a range of conditions to better quantify technology scaling effects.

Battery Electric Aircraft↗

NASA NTRSDate not supplied

Battery Cell-to-Pack Scaling Trends for Electric Aircraft

Battery pack gravimetric energy density is one of the most important, yet often miss- estimated design parameters for sizing all-electric aircraft. Proper accounting for thermal, structural, and operational safety margins are frequently lost when extrapolating performance from the cell level to the aircraft level. This paper summarizes the relevant engineering and certification details needed to better account for the penalties associated when assembling battery packs. The relationship between the cell and pack energy density is not linear, as is often assumed. Furthermore, the relationship varies depending on pack requirements, cell chemistry, and architecture. Parametric, high-fidelity models are used to determine optimal battery pack sizes over a range of conditions to better quantify technology scaling effects.

Battery Electric Aircraft↗

DOE OSTI2026

Topological magnons on the ferromagnetic zigzag lattice

Motivated by the experimental identification of magnetic compounds consisting of zigzag chains, we analyze the band structure topology of magnons in ferromagnets on a zigzag lattice. We account for the general lattice geometry by including spatially anisotropic Heisenberg exchange interactions and by Dzyaloshinskii-Moriya interaction on inversion asymmetric bonds. Within the linear spin-wave theory, we find two magnon branches, whose band structure topology (i.e., Chern numbers) we map out in a comprehensive phase diagram. Notably, besides topologically trivial and gapless phases, we identify topologically nontrivial phases that support chiral edge magnons. We show that these edge states are robust against elastic defect scattering.

Subramanian, Skandan [ORNL] (ORCID:000900087255487↗

DOE OSTI2024

Systematic determination of a material’s magnetic ground state from first principles

Abstract We present a self-consistent method based on first-principles calculations to determine the magnetic ground state of materials, regardless of their dimensionality. Our methodology is founded on satisfying the stability conditions derived from the linear spin wave theory (LSWT) by optimizing the magnetic structure iteratively. We demonstrate the effectiveness of our method by successfully predicting the experimental magnetic structures of NiO, FePS 3 , FeP, MnF 2 , FeCl 2 , and CuO. In each case, we compared our results with available experimental data and existing theoretical calculations reported in the literature. Finally, we discuss the validity of the method and the possible extensions.

Chemistry↗

DOE OSTI2026

Kondo effect in ferromagnetic quantum critical CeRh 6 ⁢Ge 4

The mechanism of a pressure-induced quantum critical point in the heavy fermion ferromagnet CeRh 6 ⁢Ge 4 has attracted interest, as ferromagnetic quantum criticality in a clean itinerant Ce compound is typically avoided. The localized versus itinerant character of the 4⁢𝑓 electrons is a key aspect for understanding this behavior. We investigated the electronic structure of the 4⁢𝑓 shell in CeRh 6 ⁢Ge 4 using core-level photoelectron and x-ray absorption spectroscopy, demonstrating the hybridization of Ce 4⁢𝑓 with the conduction electrons. Linearly polarized x-ray absorption reveals a temperature-dependent linear dichroism consistent with the crystal-electric-field sequence as inferred from the static susceptibility. This dichroism cannot be described by an ionic full-multiplet model alone, but is reproduced by including the Kondo effect within a single-impurity Anderson model in the noncrossing approximation. The Kondo effect mixes higher-lying crystal-field states into a resulting multiorbital ground state with 4⁢𝑓 occupancy, 𝑛 𝑓 ∼ 0.9. Deviations at low temperatures between the measured linear dichroism and calculated dichroism suggest an orbital-dependent Kondo effect. A scenario in which there is a multiorbital ground state and orbital-dependent Kondo hybridization should be a starting point for a model of pressure-induced criticality in CeRh 6 ⁢Ge 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗