Search NASA⌕ Search

SEARCH · Search NASA

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

At least 145 records · Page 8

Signal Prediction With Input Identification

A novel coding technique is presented for signal prediction with applications including speech coding, system identification, and estimation of input excitation. The approach is based on the blind equalization method for speech signal processing in conjunction with the geometric subspace projection theory to formulate the basic prediction equation. The speech-coding problem is often divided into two parts, a linear prediction model and excitation input. The parameter coefficients of the linear predictor and the input excitation are solved simultaneously and recursively by a conventional recursive least-squares algorithm. The excitation input is computed by coding all possible outcomes into a binary codebook. The coefficients of the linear predictor and excitation, and the index of the codebook can then be used to represent the signal. In addition, a variable-frame concept is proposed to block the same excitation signal in sequence in order to reduce the storage size and increase the transmission rate. The results of this work can be easily extended to the problem of disturbance identification. The basic principles are outlined in this report and differences from other existing methods are discussed. Simulations are included to demonstrate the proposed method.

Juang, Jer-Nan↗

STS-8 bet results

The final Best Estimate Trajectory (BET) products, i.e., the reconstructed trajectory, the Extended BET, AEROBET and MMLE input files, generated for the eighth NASA Space Shuttle flight are documented. The reconstructed trajectory (inertial BET) for this Challenger flight, the first night landing is discussed. State (position, velocity, and attitude) plus three accelerometer scale factors were determined from fitting the Guam S-band data, seven C-band passes, and pseudo Doppler and altimeter during rollout on Runway 22. The anchor epoch utilized for the batch weighted-least-squares determination was Sept. 5, 1983 7h1m50s.0 (25310 GMT seconds). The spacecraft altitude at epoch is approx. 617 kft. IMU2 data were selected for the reconstruction.

Findlay, J. T.↗

Magnetic field structure of interplanetary magnetic clouds at 1 AU

Interplanetary magnetic clouds emerge as a feature of the solar wind at 1 AU, exhibiting enhanced field strength and lower plasma temperature and density than the surrounding plasma. A least-squares program has been developed which fits magnetic field data within a cloud, while estimating such cloud properties as its size, maximum field strength, and axis inclination. The results obtained from a study of 12 clouds observed at 1 AU point to a probable cloud axis direction within 15 deg of the ecliptic plane and about 100 deg from the sun's direction, when projected into the ecliptic plane. A wide variety of orientations is observed; some extend to 80 deg from the ecliptic.

Lepping, R. P.↗

East-west inclination of large-scale photospheric magnetic fields

Sixteen years of Wilcox Solar Observatory (WSO) magnetogram data have been studied to determine the solar cycle variation and latitude dependence of the east-west inclination of photospheric magnetic field lines. East-west inclination is here defined as the angle between a field line and its local radial vector, as projected onto the plane of the latitude and line of sight. Inclination is determined by a least-squares fit of observed magnetic fields to a simple projection model, and is found to depend on polarity and to change with the solar cycle. Leading and following polarities are tipped towards each by about 9 deg and have an overall net tilt in the direction of rotation (to the west) of 0.6 deg. New cycles are seen to begin at high latitudes and to grow through the lower latitudes over approximately 5 years, providing evidence for an extended cycle length of 16-18 years.

Shrauner, J. A↗

Analysis of altimetry over inland seas

Satellite-based altimetric data taken by GEOS-3 and SEASAT over the Black Sea and Caspian Sea are analyzed and a least squares collocation technique is used to predict the geoid undulation on a .25-degree by .25-degree grid and to transform these geoid undulations to free air gravity anomalies. This project entailed processing satellite altimeter data over inland seas for recovery of area mean gravity information. Gravity information in this area of the world is not readily available, so the possibility of obtaining it from the processing of altimeter observations is attractive. The principal objective was to complete and extend analyses done in a previous study, verify those results, and document the results and techniques. A secondary objective was to improve the algorithms and results, if possible. The approach used involved editing geoid height data to remove overland data; evaluating geoid height differences at crossover points; removing orbit errors from geoid heights using crossover differences; gridding geoid height data at .25-degree by .25-degree intervals; and estimating the gravity anomalies from gridded geoid heights using the collocation technique.

Au, A. Y.↗

Uncertainty Analysis of Inertial Model Attitude Sensor Calibration and Application with a Recommended New Calibration Method

Statistical tools, previously developed for nonlinear least-squares estimation of multivariate sensor calibration parameters and the associated calibration uncertainty analysis, have been applied to single- and multiple-axis inertial model attitude sensors used in wind tunnel testing to measure angle of attack and roll angle. The analysis provides confidence and prediction intervals of calibrated sensor measurement uncertainty as functions of applied input pitch and roll angles. A comparative performance study of various experimental designs for inertial sensor calibration is presented along with corroborating experimental data. The importance of replicated calibrations over extended time periods has been emphasized; replication provides independent estimates of calibration precision and bias uncertainties, statistical tests for calibration or modeling bias uncertainty, and statistical tests for sensor parameter drift over time. A set of recommendations for a new standardized model attitude sensor calibration method and usage procedures is included. The statistical information provided by these procedures is necessary for the uncertainty analysis of aerospace test results now required by users of industrial wind tunnel test facilities.

Tripp, John S.↗

Standard remote manipulator system docking target augmentation for automated docking

A docking target is provided for use in automated docking of a first vehicle on which the target is located. The target comprises a pair of laterally extending arm portions lying in substantially the same plane and a central post extending outwardly from the plane of the arm portions. At least three reflectors are located on the target. Two of the reflectors are located at the outboard ends of the arms portions and another reflector is located at the end of the central post. In an important embodiment, the reflectors comprise individual pieces of retroreflective tape. The reflectors, when viewed from the front of the target, are aligned along the longitudinal center line of the target, and can take a number of different shapes including circular or square.

Dabney, Richard W.↗

Chemical kinetic studies of atmospheric reactions using tunable diode laser spectroscopy

IR absorption using tunable diode laser spectroscopy provides a sensitive and quantitative detection method for laboratory kinetic studies of atmospheric trace gases. Improvements in multipass cell design, real time signal processing, and computer controlled data acquisition and analysis have extended the applicability of the technique. We have developed several optical systems using off-axis resonator mirror designs which maximize path length while minimizing both the sample volume and the interference fringes inherent in conventional 'White' cells. Computerized signal processing using rapid scan (300 kHz), sweep integration with 100 percent duty cycle allows substantial noise reduction while retaining the advantages of using direct absorption for absolute absorbance measurements and simultaneous detection of multiple species. Peak heights and areas are determined by curve fitting using nonlinear least square methods. We have applied these techniques to measurements of: (1) heterogeneous uptake chemistry of atmospheric trace gases (HCl, H2O2, and N2O5) on aqueous and sulfuric acid droplets; (2) vapor pressure measurements of nitric acid and water over prototypical stratospheric aerosol (nitric acid trihydrate) surfaces; and (3) discharge flow tube kinetic studies of the HO2 radical using isotopic labeling for product channel and mechanistic analysis. Results from each of these areas demonstrate the versatility of TDL absorption spectroscopy for atmospheric chemistry applications.

Worsnop, Douglas R.↗

Towards a Third-Order Accurate, Second-Derivative-Free, Shock-Capturing Finite-Volume Method for Hypersonic Flows on Tetrahedral Grids

In this paper, we report progress in the development of a third-order accurate, second-derivative-free, shock-capturing finite-volume solver for three-dimensional unstructured grids. The method is economical in the sense that the computation and storage of second derivatives are not required for third-order accuracy. It is based on point-valued numerical solutions stored at cells, gradients computed and stored at nodes, and an efficient projected-derivative formula that eliminates the need for second derivatives in a quadratic solution interpolation. The projected-derivative formula is also used to eliminate second derivatives from a high-order flux quadrature formula, so that it can be implemented conveniently in the form of a numerical flux at a face center plus a correction term. Similarly, a high-order source quadrature formula can also be implemented in the form of a cell-center point evaluation plus a similar correction term. These features make it relatively straightforward to extend an existing second-order finite-volume code to third-order. This paper reports progress of implementing the method in the NASA VULCAN-CFD code and discusses the implementation of a high-order accurate limiter for shock capturing.

Weighted Least-Squares↗

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗

Radar altimeter waveform modeled parameter recovery

Satellite-borne radar altimeters include waveform sampling gates providing point samples of the transmitted radar pulse after its scattering from the ocean's surface. Averages of the waveform sampler data can be fitted by varying parameters in a model mean return waveform. The theoretical waveform model used is described as well as a general iterative nonlinear least squares procedures used to obtain estimates of parameters characterizing the modeled waveform for SEASAT-1 data. The six waveform parameters recovered by the fitting procedure are: (1) amplitude; (2) time origin, or track point; (3) ocean surface rms roughness; (4) noise baseline; (5) ocean surface skewness; and (6) altitude or off-nadir angle. Additional practical processing considerations are addressed and FORTRAN source listing for subroutines used in the waveform fitting are included. While the description is for the Seasat-1 altimeter waveform data analysis, the work can easily be generalized and extended to other radar altimeter systems.

Source record↗

Wind-tunnel investigation of six shielded total-pressure tubes at high angles of attack subsonic speeds

The effect of inclination of the air stream on the measured pressures of six shielded total-pressure tubes (all designed for end-mounting on a horizontal boom) has been determined for an angle-of-attack range of approximately -30° to 65°. The tests were conducted at a Mach number of 0.26 in the Langley stability tunnel and at Mach numbers of 0.50, 0.70, 0.90, and 0.95 in the Langley high-speed 7- by 10-foot tunnel. The results of the tests showed that curved venturi entries are less sensitive to inclination of the air stream than the conical entry of the standard Kiel design. Of the three curved entries tested, the tube having the bluntest entry proved the least sensitive to inclination. The critical angle (that is, the sngle.at which the total-pressure error equals 1 percent of the impact pressure) for this tube was ±63.2° at a Mach number of 0.26. The tests also showed that the critical angle of a shielded tube could be extended to higher positive angles of attack by means of a slant profile. The critical angle of a tube having a conical entry with a 10° slant profile, for example, was found to average 9° greater than that of a similar tube having a square profile. Although the critical angle of this tube at angles of yaw was less than that at angles of attack, the critical angle in yaw was equal to that of the tube with the square profile. Tests of the effect of varying the position of the total-pressure probe in the shield showed that, in general, no advantage was to be gained by changing from the position used in the standard Kiel design. Limited tests of the effect of varying the internal diameter of the shield showed that an increase in the diameter of the throat corresponding an increase in the internal area of 50 percent resulted in a relatively small increase (2.5°) in the critical angle. The effect of Mach number on the sensitivity of the tubes to inclination was shown to be appreciable. An increase in Mach number from 0.26 to 0.90, for example, resulted in a decrease in the critical angle of most of the tubes of about 4°.

Walter R Russell↗

Tracking and Data Relay Satellite (TDRS) Orbit Estimation Using an Extended Kalman Filter

Alternatives to the Tracking and Data Relay Satellite (TDRS) orbit estimation procedure were studied to develop a technique that both produces more reliable results and is more amenable to automation than the prior procedure. The Earth Observing System (EOS) Terra mission has TDRS ephemeris prediction 3(sigma) requirements of 75 meters in position and 5.5 millimeters per second in velocity over a 1.5-day prediction span. Meeting these requirements sometimes required reruns of the prior orbit determination (OD) process, with manual editing of tracking data to get an acceptable solution. After a study of the available alternatives, the Flight Dynamics Facility (FDF) began using the Real-Time Orbit Determination (RTOD(Registered TradeMark)) Kalman filter program for operational support of TDRSs in February 2007. This extended Kalman filter (EKF) is used for daily support, including within hours after most thrusting, to estimate the spacecraft position, velocity, and solar radiation coefficient of reflectivity (C(sub R)). The tracking data used are from the Bilateration Ranging Transponder System (BRTS), selected TDRS System (TDRSS) User satellite tracking data, and Telemetry, Tracking, and Command (TT&C) data. Degraded filter results right after maneuvers and some momentum unloads provided incentive for a hybrid OD technique. The results of combining EKF strengths with the Goddard Trajectory Determination System (GTDS) Differential Correction (DC) program batch-least-squares solutions, as recommended in a 2005 paper on the chain-bias technique, are also presented.

Ward, Douglas T.↗

Separation of lithospheric, external, and core components of the south polar geomagnetic field at satellite altitudes

We present a new approach to producing scalar Magsat magnetic anomaly maps based on correlation coefficient filtering and the use of almost all of the available orbits. Our method differs from earlier techniques with respect to the following: (1) Passes are selected based on their variance properties rather than planetary indices such as Kp. (2) The core field model is least squares fit to individual passes and subsequently removed instead of substracting the model directly. This technique replaces band pass filtering or polynomial trend removal methods. (3) Each selected pass is sorted geographically and by local time, placed into one of four different altitude bands, and correlation coefficient filtered with the two adjacent passes. The filtering is the second step toward isolating the static lithospheric signal from the more dynamic external field signals. (4) Least squares collocation is used to grid the correlated passes; subsequently, the dawn and dusk maps are also correlation filtered providing another step toward removal of external fields. (5) The four resultant total field maps are continued to a common altitude and again correlation filtered for the static lithospheric anomalies. (6) The filtered results are then averaged together to provide a new total field map of the lithosphere south of 40 deg S latitude. Our total field map differs from previous efforts over the crustal blocks of West Antarctica. We obtained a positive anomaly over Edward VII Peninsula, extending into the Byrd subglacial basin and obtained a negative anomaly over the Ellsworth Mountains and parts of the Byrd subglacial basin. Also, a positive anomaly extending from the Ross Sea to offshore Wilkes Land is present in our map; however, this feature is absent in other maps. Positive anomalies marking the Weddell Sea in previous efforts are not present in our map. Prominent external field anomalies in the quadrant offshore of Wilkes Land are present in all previous efforts; however, these signatures are reduced in our total field map. This map significantly advances efforts for the identification and interpretation of magnetic anomalies of the south polar lithosphere.

Alsdorf, Douglas E.↗

Possible Quantum Paramagnetism in Compressed Sr 2 IrO 4

The result of compression on the magnetic ground state of Sr 2 IrO 4 is studied with x-ray resonant techniques in the diamond anvil cell. The weak interlayer exchange coupling between square-planar 2D IrO 2 layers is easily modified upon compression, with a crossover between magnetic structures around 7 GPa mimicking the effect of an applied magnetic field at ambient pressure. Higher pressures drive an order-disorder magnetic phase transition with no magnetic order detected above 17–20 GPa. The persistence of strong exchange interactions between J eff = 1 / 2 magnetic moments within the insulating IrO 2 layers up to at least 35 GPa points to a highly frustrated magnetic state in compressed Sr 2 IrO 4 , opening the door for realization of novel quantum paramagnetic phases driven by extended 5 d orbitals with entangled spin and orbital degrees of freedom.

36 MATERIALS SCIENCE↗

Stochastic average model methods

We consider the solution of finite-sum minimization problems, such as those appearing in nonlinear least-squares or general empirical risk minimization problems. We are motivated by problems in which the summand functions are computationally expensive and evaluating all summands on every iteration of an optimization method may be undesirable. Here we present the idea of stochastic average model (SAM) methods, inspired by stochastic average gradient methods. SAM methods sample component functions on each iteration of a trust-region method according to a discrete probability distribution on component functions; the distribution is designed to minimize an upper bound on the variance of the resulting stochastic model. We present promising numerical results concerning an implemented variant extending the derivative-free model-based trust-region solver POUNDERS, which we name SAM-POUNDERS.

97 MATHEMATICS AND COMPUTING↗

Quantification of Rare Earth Elements in the Parts Per Million Range: A Novel Approach in the Application of Laser-Induced Breakdown Spectroscopy

This work extends a previous percentage level concentration study of the optical emission spectra for six rare earth elements, europium (Eu), gadolinium (Gd), lanthanum (La), praseodymium (Pr), neodymium (Nd), and samarium (Sm), along with the transition metal, yttrium (Y) using laser-induced breakdown spectroscopy (LIBS). The concentration of these six rare earth elements and yttrium has been attempted for the first time systematically down to parts per million (ppm) concentration levels ranging from 30 to 300 ppm. In this study, the authors have developed multivariate models for each element capable of predicting concentration with acceptable to excellent levels of accuracy. Additionally, partial least squares regression coefficients were used to identify key spectral features able to be used in this lower concentration regime. This study has demonstrated that it is conceivable to quantify the six rare earth elements along with yttrium at low concentrations in the parts per million levels.

59 BASIC BIOLOGICAL SCIENCES↗

LightForce Photon-Pressure Collision Avoidance: Updated Efficiency Analysis Utilizing a Highly Parallel Simulation Approach

This paper provides an updated efficiency analysis of the LightForce space debris collision avoidance scheme. LightForce aims to prevent collisions on warning by utilizing photon pressure from ground based, commercial off the shelf lasers. Past research has shown that a few ground-based systems consisting of 10 kilowatt class lasers directed by 1.5 meter telescopes with adaptive optics could lower the expected number of collisions in Low Earth Orbit (LEO) by an order of magnitude. Our simulation approach utilizes the entire Two Line Element (TLE) catalogue in LEO for a given day as initial input. Least-squares fitting of a TLE time series is used for an improved orbit estimate. We then calculate the probability of collision for all LEO objects in the catalogue for a time step of the simulation. The conjunctions that exceed a threshold probability of collision are then engaged by a simulated network of laser ground stations. After those engagements, the perturbed orbits are used to re-assess the probability of collision and evaluate the efficiency of the system. This paper describes new simulations with three updated aspects: 1) By utilizing a highly parallel simulation approach employing hundreds of processors, we have extended our analysis to a much broader dataset. The simulation time is extended to one year. 2) We analyze not only the efficiency of LightForce on conjunctions that naturally occur, but also take into account conjunctions caused by orbit perturbations due to LightForce engagements. 3) We use a new simulation approach that is regularly updating the LightForce engagement strategy, as it would be during actual operations. In this paper we present our simulation approach to parallelize the efficiency analysis, its computational performance and the resulting expected efficiency of the LightForce collision avoidance system. Results indicate that utilizing a network of four LightForce stations with 20 kilowatt lasers, 85% of all conjunctions with a probability of collision Pc > 10 (sup -6) can be mitigated.

conjunction analysis↗