Search NASASearch

SEARCH · Search NASA

Results for “linear least squares”

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.

A Case Study of View-Factor Rectification Procedures for Diffuse-Gray Radiation Enclosure Computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Robert P Taylor

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

Remote Guidance of a Mothership-Daughtership Formation for Asteroid Exploration

This work presents an application of Model Predictive Control to maintain long-term stability for a three-satellite formation around the asteroid Apophis via remote guidance. The mothership and two daughterships perform tomographic characterization with a bistatic radar. The mothership estimates the states of all three spacecraft and determines the necessary control to maintain the desired antipodal daughtership states. An Extended Kalman Filter estimates the mothership states. A Batch Least Squares algorithm estimates the daughtership states and performs maneuver reconstruction for the mothership-administered thrust commands under thruster and gravitational uncertainty. The formulation maintains the antipodal angle with growing instability amidst these uncertainties.

Guidance

Satellite Material Contaminant Optical Properties

The Air Force Wright Research and Development Center and the Arnold Engineering Development Center are continuing a program for measuring optical effects of satellite material outgassing products on cryo-optic surfaces. Presented here are infrared (4000 to 700 cm -1 ) transmittance data for contaminant films condensed on a 77 K geranium window. From the transmittance data, the contaminant film refractive and absorptive indices (n, k) were derived using an analytical thin-film interference model with a nonlinear least-squares algorithm. To date 19 materials have been studied with the optical contents determined for 13 of those. The materials include adhesives, paints, composites, films, and lubricants. This program is continuing and properties for other materials will be available in the future.

B E Wood

Determination of Local Experimental Heat-Transfer Coefficients on Combustion Side of an Ammonia-Oxygen Rocket

Local experimental heat-transfer coefficients were measured in the chamber and throat of a 2400-pound-thrust ammonia-oxygen rocket engine with a nominal chamber pressure of 600 pounds per square inch absolute. Three injector configurations were used. The rocket engine was run over a range of oxidant-fuel ratio and chamber pressure. The injector that achieved the best performance also produced the highest rates of heat flux at design conditions. The heat-transfer data from the best-performing injector agreed well with the simplified equation developed by Bartz at the throat region. A large spread of data was observed for the chamber. This spread was attributed generally to the variations of combustion processes. The spread was least evident, however, with the best-performing injector.

Curt H Liebert

Exploration of an Adaptive Routine for Battery Modeling

The purpose of this document is to explore the use of adaptive routines in battery modeling. The adaptive routines consist of real-time state estimators combined with battery parameter model components that are adjusted in real-time as battery data becomes available. Several aspects are explored. It is shown that model parameter identification is possible for simple battery models using available input/output data measurements. The online system identification used is recursive least squares. Model identification may be combined with a state observer such as the extended Kalman filter or the unscented Kalman filter to form an adaptive model combined with state estimation. However, such a combination is found to be problematic due to uncertainty, observability and stability issues. This paper is organized as follows. Section 1 introduces adaptive routines and possible roles they play in battery modeling. In Section 2 real-time parameter identification is described with results based on battery data. Section 3 reviews various state estimators and results using a simple battery model. In Section 4 parameter identification and state estimation are combined to form an adaptive routine. Finally, in Section 5 conclusions are drawn and future work is suggested.

Adaptive

Air-Broadened Line Shapes in the 2n3 R Branch of 12ch4 Between 6014 and 6100 cm-1

Complete and accurate information on line shape parameters of 2?3 methane transitions for air broadening as a function of temperature is critical not only for the correct interpretation of the observed atmospheric spectra but also for the development of a reliable theoretical model. For this reason, we obtained a series of high-resolution, high S/N spectra of high-purity 12CH4 and 12CH4 broadened with dry air at temperatures in the 130 to 295 K range using the Bruker IFS 125HR Fourier transform spectrometer at JPL. Two absorption cells were used in the experiment, a White cell with path length of 13 m for room temperature spectra and a 21 m Herriott cella for cold sample spectra. The 15 spectra used in the analysis consisted of 3 low pressure (0.26 to 2.57 Torr) spectra with pure 12CH4 and 12 air-broadened spectra with total sample pressures of 79-805 Torr and volume mixing ratios of methane between 0.23 and 1%. A multispectrum least-squares fitting techniqueb was employed to fit all 15 spectra simultaneously. Preliminary results for select R(J) manifolds will be presented.

V Malathy Devi

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

Pressurized-Cell Experiment

Pressurized-cell detectors were developed and constructed at the Langley Research Center to obtain a direct measurement of the micrometeoroid puncture hazard to thin structural material. A total of 160 detectors were mounted around the periphery of the rocket motor used on the Explorer XIII satellite. (See fig. IX-1.) Material thicknesses of 0.001, 0.0015, 0.002, 0.0025, and 0.005 inch were used on the detectors. The exposed surface of the test material of each detector was 21.8 square inches and the total exposed surface of the 160 detectors was 24.2 square feet. The average weight of the detectors was 68 grams.

Micrometeoroid

Battery system and method for sensing and balancing the charge state of battery cells

A battery system utilizes a plurality of transformers interconnected with the battery cells. The transformers each have at least one transformer core operable for magnetization in at least a first magnetic state with a magnetic flux in a first direction and a second magnetic state with a magnetic flux in a second direction. The transformer cores retain the first magnetic state and the second magnetic state without current flow through said plurality of transformers. Circuitry is utilized for switching a selected transformer core between the first and second magnetic states to sense voltage and/or balance particular cells or particular banks of cells.

Davies, Francis J.

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

Advanced Space Power PEM Fuel Cell Systems

A model showing mass and heat transfer in proton exchange membrane (PEM) single cells is presented. For space applications, stack operation requiring combined water and thermal management is needed. Advanced hardware designs able to combine these two techniques are available. Test results are shown for membrane materials which can operate with sufficiently fast diffusive water transport to sustain current densities of 300 ma per square centimeter. Higher power density levels are predicted to require active water removal.

PEM

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

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.

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

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.

The Shape Effect: Influence of 1D and 2D Boron Nitride Nanostructures on the Radiation Shielding, Thermal, and Damping Properties of High-Temperature Epoxy Composites

In space exploration, lightweight multifunctional materials capable of shielding neutron radiation, dissipating heat, and providing damping are essential. Polymer composites reinforced with boron nitride (BN) nanomaterials—specifically one-dimensional boron nitride nanotubes (BNNTs) and two-dimensional boron nitride nanoplatelets (BNNPs)—offer promising solutions. This study investigates how BN nanomaterial morphology influences the performance of high-temperature (HT) epoxy composites. We developed ultralightweight, three-dimensional BN foams comprising 1D BNNTs, 2D BNNPs, and hybrid 1D BNNT/2D BNNP structures via freeze-drying, then infiltrated them with HT epoxy to form dense composites. The BNNT foam exhibited the highest neutron radiation shielding, with a mass absorption coefficient of 26.64 cm2 g −1 , outperforming the hybrid foam (18.18 cm 2 g −1 ) and the BNNP foam (11.12 cm 2 g −1 ). A similar trend was observed in the HT epoxy composites; incorporating these foams at least doubled the mass absorption coefficient compared to the neat polymer. In terms of thermal conductivity, the BNNT/BNNP foam-epoxy composite achieved the highest value of 0.34 W m −1 K −1 , a 2.13-fold increase over neat HT epoxy. The BNNT/BNNP foam-epoxy composites also improved by 1.88 and 1.75 times, respectively. Mechanical testing revealed that BNNP foams withstood the highest loads during nanoindentation (3.53 kN), followed by BNNT/BNNP foams (1.93 kN) and BNNT foams (1.56 kN). All BN foam-epoxy composites exhibited enhanced damping properties, with tan δ increasing by at least 30 % compared to neat HT epoxy. These findings elucidate the impact of BN nanomaterial morphology on the multifunctional performance of HT epoxy composites, offering insights for developing high-performance, tailorable materials for demanding environments.

Kazue Orikasa

Cryopump for Large Pulses of Hydrogen

A test involving a flow of hot hydrazine decomposition products at a rate of 13.6 g/s (0.03 lb/s) established the requirement to maintain pressure in an 11.9-m-(39 ft) diameter space chamber below 200 microns. The flow, 2/3 hydrogen and 1/3 nitrogen by volume, continued for several periods ranging from 3 to 15 min. The pressure requirement was necessary to minimize thermal effects of the gas on the test vehicle but was well beyond the capability of the existing facility pumps. Various methods of obtaining additional temporary jump capacity were considered. From these, the slugged-charcoal approach was selected as the quickest and least expensive method to implement.

B C Moore