Search NASA⌕ Search

SEARCH · Search NASA

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.

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

23 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Optical properties of morphologically complex black carbon aerosols: Effects of coatings

Optical properties are computed for fractal-like aggregate black carbon (BC) aerosols coated with different substances. Two models are used for these aerosols: (i) the coated aggregate model (model I), where coating is added to the voids and surroundings of BC fractal-like aggregates; and (ii) the closed-cell aggregate model (model II), where coating is added concentrically to each monomer of BC fractal-like aggregates. Our results favor choosing the coated aggregate model (model I) to simulate scattering and absorption by coated BC aerosols because this model is morphologically more realistic, and because this model yields mass absorption cross section (MAC) and backscattering linear depolarization ratio (LDR) values that are consistent with field measurements. Moreover the corresponding computed degree of linear polarization (DoLP) and LDR values are very sensitive to changes in the coating volume fraction (fvol) and the coating refractive index (m). With the same absorbing BC core, the MAC value increases steadily with increasing fvol and m values. For example, using a density of 1.8 g/cm 3 and BC refractive index 1.95 + i0.79, the calculated MAC values for uncoated BC aggregates range from 6.2 to 6.8 m 2 /g at 0.55 µm. When coating material is applied to the BC aggregates, the calculated MAC values for model I particles increase to between 9.8 and 13.2 m 2 /g (depending on m) when fvol = 87.5%. The backscattering LDR values also tend to increase with the increasing m values for the shapes and sizes considered in our study. For model I particles, the backscattering LDR values span a wide range of 4.2–27.8% at a wavelength of 0.35 µm at fvol = 87.5% when m increases from 1.33 to 1.55. Our results are relevant to analyses of polarimetric and lidar observations of smoke particles, especially when these particles undergo hygroscopic growth.

Black carbon aerosols↗
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↗
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↗
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↗
NASA NTRS1998

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Isaac Elishakoff↗
NASA NTRSDate not supplied

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

Accurate prediction of aerothermodynamic loads in thermal non-equilibrium flows requires precise modeling of internal energy exchange. While previous direct simulation Monte Carlo frameworks have successfully implemented discrete rotational energy models for diatomic species, the treatment of polyatomic molecules has traditionally relied on continuous energy assumptions that break down at low temperatures and neglect critical high-temperature corrections. This study extends the discrete rotational energy models of Boyd and Gimelshein to fully encompass polyatomic molecules. The proposed framework implements quantized rotational energy level sampling for linear, spherical, and symmetric/asymmetric top rotors. Crucially, the model incorporates centrifugal distortion to address the limitations of the rigid-rotor assumption at hypersonic temperatures, and accounts for nuclear spin parity, which dictates the permissible rotational states and macroscopic specific heats at low temperatures. The model is verified through equilibrium sampling procedures, demonstrating agreement with theoretical quantum Boltzmann distributions and accurately reproducing thermophysical properties across a wide range of temperatures.

DSMC↗
NASA NTRSDate not supplied

Softening the Gap between Wöhler and Paris – New Approaches for Fatigue Analysis –

Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.

cohesive elements↗
NASA NTRS2026

Parabolic Deformable Mirror Project Development and Device Characterization

Existing state-of-the-art wavefront sensing and control setups for the imaging of exoplanets utilize two high-actuator count flat deformable mirrors (DMs) to control both phase and amplitude aberrations. One of the objectives of the Exoplanet Spectroscopy (ExoSpec) project at Goddard Space Flight Center (GSFC) is to improve upon this regime through the development of parabolic deformable mirrors (PDMs), which can address several of the key limitations in existing setups. Specifically, the need for flat DMs to be located in conjugate pupil planes drives unique packaging challenges that can be overcome by allowing the off-axis imaging elements to be controllable as well. In addition, simulations show that the use of deformable powered elements will increase the overall controllable bandwidth, up to 35% over 5-12 lambda/D. To pursue this effort, the ExoSpec group has begun testing a first-generation parabolic DM, which is already informing a second generation design. This paper will provide a comprehensive update on its characterization, including influence function measurements, linearity analysis, closed-loop wavefront control performance, and thermal stability assessment.

spectroscopy↗
NASA NTRSDate not supplied

Stationary Testing of a Pulse Tube Cryocooler Designed to Rotate at 6800 RPM with the Superconducting Rotor of an Electric Machine

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Most of these electric machines employ superconductors on their rotor, thereby requiring cryogenic cooling of the rotating system. An attractive solution is to conductively cool the rotor using a cryocooler that rotates with the rotor. NASA has been developing a 50 K pulse tube cryocooler for its High Efficiency Megawatt Motor (HEMM) that can operate while rotating at 6,800 rpm. The design of this cryocooler and testing of its linear motor have been discussed in prior work. This presentation presents the first measurements of the complete HEMM cryocooler. The design and results of a stationary (non-rotating) test of the cryocooler are described.

electrified aircraft propulsion↗
NASA NTRS1985

Summary of recent NASA propeller research

Advanced high speed propellers offer large performance improvements for aircraft that cruise in the Mach 0.7 to 0.8 speed regime. At these speeds, studies indicate that there is a 15 to near 40 percent block fuel savings and associated operating cost benefits for advanced turboprops compared to equivalent technology turbofan powered aircraft. Recent wind tunnel results for five eight to ten blade advanced models are compared with analytical predictions. Test results show that blade sweep was important in achieving net efficiencies near 80 percent at Mach 0.8 and reducing nearfield cruise noise about 6 dB. Lifting line and lifting surface aerodynamic analysis codes are under development and some results are compared with propeller force and probe data. Also, analytical predictions are compared with some initial laser velocimeter measurements of the flow field velocities of an eight bladed 45 swept propeller. Experimental aeroelastic results indicate that cascade effects and blade sweep strongly affect propeller aeroelastic characteristics. Comparisons of propeller nearfield noise data with linear acoustic theory indicate that the theory adequately predicts nearfield noise for subsonic tip speeds, but overpredicts the noise for supersonic tip speeds.

Daniel C Mikkelson↗
NASA NTRS2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗
NASA NTRSDate not supplied

Predicting Team Functioning in Long Term Space Missions Using Acoustic and Linguistic Measures

Maintaining optimal team functioning is critical for long-duration space exploration missions, yet traditional monitoring methods, such as self-reports and wearable sensors, often impose operational burdens or suffer from bias. This paper investigates a non-intrusive speech-based artificial intelligence (AI) framework to predict degradations in team functioning using data from the Human Exploration Research Analog (HERA) of the U.S. National Aeronautics and Space Administration (NASA). Using acoustic features, linguistic descriptors, and semantic embeddings, we evaluate static non-linear and temporal machine learning models to predict both objective (task accuracy) and subjective (self-reported efficacy and cohesion) team functioning outcomes. Results indicate that temporal models outperform static approaches, with prediction of objective task accuracy in Team Interaction Battery (TIB) improving from near chance to 71%. Self-reported outcomes, including team efficacy and cohesion, are predicted more reliably than task performance, achieving balanced accuracies of up to 85.56% and 78.12%, respectively, and are found to be most strongly associated with acoustic features. In a second interdependent task, the MMSEV–EVA, accuracies of up to 78% are achieved using temporal models with acoustic features. Furthermore, incorporating just 1–2 days of team-specific historical data systematically improved performance, and acoustic markers from informal pre-task interactions provided modest predictive gains. Finally, while automated preprocessing yielded viable accuracy, humancorrected data provided moderate performance gains, though transcription error rates did not significantly correlate with model performance. These findings highlight the potential of speech as a passive, high-fidelity monitoring tool for autonomous habitats.

Temporal modeling↗
NASA NTRS2026

In-Flight Calibration of ESA Hera’s Hyperscout-H Imager

ESA’s Hera space mission is on its way to the mission target, the binary asteroid (65803) Didymos. HyperScout-H, one of the instruments onboard Hera, is a hyperspectral imager operating in the visible and near-infrared regions between 0.65 and 0.95 μm. HyperScout-H will enable a detailed assessment of the composition of both objects, Didymos and its satellite Dimorphos, the characterization of space weathering effects, and the possible presence of exogenous material on their surfaces. To monitor instrument functionality, calibration exposures are acquired regularly. This article describes the in-flight calibrations carried out for HyperScout-H during the commissioning and cruise phases. Bias and dark exposures, as well as stellar field observations, were acquired several times after launch. We update the calibration data and monitor instrument performance in the space environment. In addition, images of Earth and Moon were acquired from distances of 1.5 × 10 6 to 2.0 × 10 6 km, and Mars and its satellite were imaged during the flyby. In five images, the surface of Mars fills the entire field of view, enabling cross-validation of HyperScout-H results with those reported by other Mars missions. We characterize the detector under in-flight operational conditions. The calibration data indicate that the bias pattern is stable, the dark current remains negligible for short exposures, and the detector response is highly linear. We quantify the field-of-view alignment and geometric distortion, and evaluate the point spread function based on the stellar field observations. Stellar observations and Mars swing-by data provide updated radiometric calibration constants, suggesting that in-flight conditions have slightly modified the detector’s spectral response. In-flight calibrations are essential to ensure data quality and reliability.

Image Processing↗