Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “least square approximations”

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 289 records · Page 16

Regularization via f -Divergence: An Application to Multi-Oxide Spectroscopic Analysis

In this paper, we explore the application of convolutional neural networks (CNNs) for predicting the chemical composition of complex geologic samples in a simulated Martian atmospheric environment. Specifically, we aim to characterize oxide weight percentages (wt.%) of rock samples analyzed by remote Laser-Induced Breakdown Spectroscopy (LIBS), framing the problem as a multi-target regression task . Neural networks trained on LIBS spectra are prone to overfitting due to high spectral complexity, limited labeled data, and measurement noise. While regularization is critical for improving generalization, common methods (e.g., ℓ 2 regularization) impose constraints not directly tied to data distribution properties. We propose a novel regularization method based on a specific ƒ-divergence induced by a graph-based estimator, designed to constrain the distributional discrepancy between predictions and targets. This regularizer serves a dual purpose: (a) mitigating overfitting by enforcing a constraint on the distributional difference between predictions and noisy targets, and (b) acting as an auxiliary loss that penalizes large divergences. To enable backpropagation, we develop a differentiable approximation of this particular ƒ-divergence, making the method feasible for neural networks. Experiments on ChemCam and SuperCam LIBS calibration spectra show that mathematical equation-divergence regularization outperforms or matches standard regularization methods (ℓ 1 , ℓ 2 , dropout) and the classical baseline, partial least squares (PLS). Combining ƒ-divergence regularization with standard regularization yields further performance gains, indicating that distributional regularization is useful in this context giving a promising direction for robust model training in planetary science applications. Source code is publicly available at Klein and Li (2025), https://doi.org/10.11578/dc.20250530.7.

58 GEOSCIENCES↗

Algorithms

The implementation of the algorithms used in the flight program to approximate elementary functions and mathematical procedures was checked. This was done by verifying that at least one, and in most cases, more than one function computed through the use of the algorithms was calculated properly. The following algorithms were checked: sine-cosine, arctangent, natural logarithm, square root, inverse square root, as well as the vector dot and cross products.

Source record↗

Kyoto Meteor Radar Observation in November, 1981 and May, 1982

The Kyoto meteor radar a pulsed Doppler radar and operates at 31.57 MHz ith peak power of 10 kW. An interferometer is adopted for arrival angle measurement. Height resolution is approximately 1 to 3 km depending on both azimuth and elevation angles. Data is analyzed according to an algorithm approved by GROVES. All data in height range 84 to 104 km are separated into five layers, whose altitude width is 4 km. In each layer, least square fitting is performed for ten parameters, which are amplitudes of northward/eastward mean wind and sin/cos component of diurnal/semidiurnal tide. Reliable range for each parameter is also calculated.

Tsuda, T.↗

Fuzzy Modeling and Parallel Distributed Compensation for Aircraft Flight Control from Simulated Flight Data

A method is described that combines fuzzy system identification techniques with Parallel Distributed Compensation (PDC) to develop nonlinear control methods for aircraft using minimal a priori knowledge, as part of NASA’s Learn-to-Fly initiative. A fuzzy model was generated with simulated flight data, and consisted of a weighted average of multiple linear time invariant state-space cells having parameters estimated using the equation-error approach and a least-squares estimator. A compensator was designed for each subsystem using Linear Matrix Inequalities (LMI) to guarantee closed-loop stability and performance requirements. This approach is demonstrated using simulated flight data to automatically develop a fuzzy model and design control laws for a simplified longitudinal approximation of the F-16 nonlinear flight dynamics simulation. Results include a comparison of flight data with the estimated fuzzy models and simulations that illustrate the feasibility and utility of the combined fuzzy modeling and control approach.

Weinstein, Rose↗

Bayesian Monte Carlo Evaluation Framework for Imperfect Nuclear Data

Bayesian evaluation of resolved resonance region (RRR) nuclear data has historically been carried out using the generalized least squares (GLS) formalism, as implemented in, e.g., SAMMY. We have recently developed a prototype of Bayesian Monte Carlo (BMC) evaluation framework, implemented using a Markov Chain Monte Carlo (MCMC) method with a Metropolis-Hastings (MH) acceptance criterion. This was done in order to remove the approximations underlying the conventional GLS evaluations, namely, the linear approximation, and the approximation that all probability density functions (PDFs) are of the normal kind. Recent works by others have used similar stochastic approaches to quantify cross section uncertainties from ENDF evaluated co-variances, and/or, from integral benchmark data, but those have not been conceived as an evaluation framework like the one presented here.

97 MATHEMATICS AND COMPUTING↗

Neural networks for function approximation in nonlinear control

Two neural network architectures are compared with a classical spline interpolation technique for the approximation of functions useful in a nonlinear control system. A standard back-propagation feedforward neural network and a cerebellar model articulation controller (CMAC) neural network are presented, and their results are compared with a B-spline interpolation procedure that is updated using recursive least-squares parameter identification. Each method is able to accurately represent a one-dimensional test function. Tradeoffs between size requirements, speed of operation, and speed of learning indicate that neural networks may be practical for identification and adaptation in a nonlinear control environment.

Linse, Dennis J.↗

Bayesian Monte Carlo Evaluation Framework for Imperfect Data and Models [Abstract]

Nuclear data evaluation methods conventionally make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate for minimization of a cost function (even for non-linear models); and that both the model (of, e.g., neutron cross section) and experimental data (including their covariance data) are perfect. These assumptions are inherent to the well-known generalized linear least squares (GLLS) minimization method commonly used for evaluations of resolved resonance region (RRR) neutron cross sections. However, these assumptions are often not justified due to the presence of non-normal PDFs, non-linear models (e.g. R -matrix formalism), and inherent imperfections in data and models (e.g. discrepant data sets, discrepancies between the previous evaluation and newly measured data, or imperfect covariance data). We remove the said assumptions in a mathematical framework of Bayes’ theorem, and implement it using the Metropolis-Hastings Monte Carlo method. Parameters of a new kind are introduced to parameterize inherent imperfections, e.g. , any discrepancies between the theoretical model and measured data. These new parameters enable evaluators to quantify their expert judgement about any discrepancies or imperfections in a reproducible manner. We demonstrate the framework with an ongoing evaluation of 233 U in the eV region using the ENDF-B/VIII library and transmission data measured by Guber, et al. , and compare the posterior parameters to those obtained by conventional evaluation methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Self‐Potential Tomography Preconditioned by Particle Swarm Optimization—Application to Monitoring Hyporheic Exchange in a Bedrock River

Abstract A self‐potential (SP) data‐inversion algorithm was developed and tested on an analytical model of electrical‐potential profile data attributed to single and multiple polarized electrical sources. The developed algorithm was then validated by an application to SP‐monitoring field data measured on the floodplain of East Fork Poplar Creek, Oak Ridge, Tennessee, to image electrical sources in areas conducive to preferential flow into the flood plain from the bedrock‐lined riverbed. The algorithm combined stochastic source‐localization by particle‐swarm‐optimization (PSO) of electrical sources characterized by simplified geometries with source tomography by regularized weighted least‐squares minimization of a quadratic objective function. Prior information was incorporated by preconditioning the tomography algorithm by PSO results. Variable percentages of random noise were added to analytical‐model data to evaluate the algorithm performance. Results indicated that true parameters of single‐source models were inverted and approximated with small residual error, whereas inversion of analytical‐model data representing multiple electrical sources accurately approximated the locations of the sources but miscalculated some parameters because of the non‐uniqueness of the inverse‐model solution. Source tomography applied to analytical model data during testing produced a spatially continuous parameter field that identified the locations of point‐scale synthetic dipole sources of electrical current flow with varying degrees of accuracy depending on the prior information incorporated into the tomography. When applied to SP‐monitoring field data, the algorithm imaged electrical sources within a known fault that intersects the bedrock riverbed and flood plain of East Fork Poplar Creek and depicted dynamic electrical conditions attributed to hyporheic exchange.

54 ENVIRONMENTAL SCIENCES↗

Time delay occultation data of the Helios spacecrafts and preliminary analysis for probing the solar corona

S-band time delay measurements were collected from the spacecrafts Helios A and B during three solar occultations in 1975/76 within heliocentric distances of about 3 and 215 solar radii in terms of range, Doppler frequency shift, and electron content. A description is given concerning some characteristic features of the methods of measurement and data processing. Typical data sets are discussed to probe the electron density distribution near the sun (west and east limb as well) including the outer and extended corona. Steady-state and dynamical aspects of the solar corona are presented and compared with earth-bound K-coronagraph measurements. Using a weighted least squares estimation 3 parameters of an average coronal electron density profile are derived in a preliminary analysis to yield electron densities of about 130 billion; 100 million; 7 million/cu m at r?3; 65; 215 solar radii. Transient phenomena are discussed and a velocity of propagation v approximately 900 km/s for plasma ejecta from a solar flare is determined from an extraordinary set of Helios B electron content measurements on April 30/May 1, 1976.

Edenhofer, P.↗

Total ionospheric electron content calibration using SERIES GPS satellite data

The current status of the Deep Space Network advanced systems research into ionospheric calibration techniques, based on Global Positioning System (GPS) data is described. A GPS-based calibration system is planned to replace the currently used Faraday rotation method by 1989. The SERIES receiver system determines the differential group delay of signals transmitted at two different carrier frequencies. This differential delay includes an ionospheric component and a GPS transmitter offset. The transmitter offsets are different for each GPS satellite. Tests were conducted to assess the effect of the offsets on the ionospheric calibration accuracy. From the obtained data, the total electron content and GPS transmitter offsets are calculated by a least squares estimation method employing a local model of total ionospheric electron content. The end product is an estimation of the total ionospheric content for an arbitrary line-of-sight direction. For the presented polynomial fitting technique, the systematic error due to mismodeling is estimated to be approximately 6 x 10 to the 16th power el/sq m, while the formal error is approximately 2 x 10 to the 16th power el/sq m. The final goal is an error of 3 x 10 to the 16th power el/sq m (approximately 0.7 ns at 2.3 GHz).

Lanyi, G.↗

Pressure and higher-order spectra for homogeneous isotropic turbulence

The spectra of the pressure, and other higher-order quantities including the dissipation, the enstrophy, and the square of the longitudinal velocity derivative are computed using data obtained from direct numerical simulation of homogeneous isotropic turbulence at Taylor-Reynolds numbers R(sub lambda) in the range 38 - 170. For the pressure spectra we find reasonable collapse in the dissipation range (of the velocity spectrum) when scaled in Kolmogorov variables and some evidence, which is not conclusive, for the existence of a k(exp -7/3) inertial range where k = absolute value of K, is the modulus of the wavenumber. The power spectra of the dissipation, the enstrophy, and the square of the longitudinal velocity derivative separate in the dissipation range, but appear to converge together in the short inertial range of the simulations. A least-squares curve-fit in the dissipation range for one value of R(sub lambda) = 96 gives a form for the spectrum of the dissipation as k(exp 0)exp(-Ck eta), for k(eta) greater than 0.2, where eta is the Kolmogorov length and C is approximately equal to 2.5.

Pullin, D. I.↗

Investigation of multiple scattering effects in aerosols

The results are presented of investigations on the various aspects of multiple scattering effects on visible and infrared laser beams transversing dense fog oil aerosols contained in a chamber (4' x 4' x 9'). The report briefly describes: (1) the experimental details and measurements; (2) analytical representation of the aerosol size distribution data by two analytical models (the regularized power law distribution and the inverse modified gamma distribution); (3) retrieval of aerosol size distributions from multispectral optical depth measurements by two methods (the two and three parameter fast table search methods and the nonlinear least squares method); (4) modeling of the effects of aerosol microphysical (coagulation and evaporation) and dynamical processes (gravitational settling) on the temporal behavior of aerosol size distribution, and hence on the extinction of four laser beams with wavelengths 0.44, 0.6328, 1.15, and 3.39 micrometers; and (5) the exact and approximate formulations for four methods for computing the effects of multiple scattering on the transmittance of laser beams in dense aerosols, all of which are based on the solution of the radiative transfer equation under the small angle approximation.

Deepak, A.↗

A numerical scheme for the identification of hybrid systems describing the vibration of flexible beams with tip bodies

A cubic spline based Galerkin-like method is developed for the identification of a class of hybrid systems which describe the transverse vibration to flexible beams with attached tip bodies. The identification problem is formulated as a least squares fit to data subject to the system dynamics given by a coupled system of ordnary and partial differential equations recast as an abstract evolution equation (AEE) in an appropriate infinite dimensional Hilbert space. Projecting the AEE into spline-based subspaces leads naturally to a sequence of approximating finite dimensional identification problems. The solutions to these problems are shown to exist, are relatively easily computed, and are shown to, in some sense, converge to solutions to the original identification problem. Numerical results for a variety of examples are discussed.

Rosen, I. G.↗

An efficient approach to ARMA modeling of biological systems with multiple inputs and delays

This paper presents a new approach to AutoRegressive Moving Average (ARMA or ARX) modeling which automatically seeks the best model order to represent investigated linear, time invariant systems using their input/output data. The algorithm seeks the ARMA parameterization which accounts for variability in the output of the system due to input activity and contains the fewest number of parameters required to do so. The unique characteristics of the proposed system identification algorithm are its simplicity and efficiency in handling systems with delays and multiple inputs. We present results of applying the algorithm to simulated data and experimental biological data In addition, a technique for assessing the error associated with the impulse responses calculated from estimated ARMA parameterizations is presented. The mapping from ARMA coefficients to impulse response estimates is nonlinear, which complicates any effort to construct confidence bounds for the obtained impulse responses. Here a method for obtaining a linearization of this mapping is derived, which leads to a simple procedure to approximate the confidence bounds.

Non-NASA Center↗

Stress corrosion cracking susceptibility of 18 Ni maraging steel

The stress corrosion cracking (SCC) resistance of 18Ni maraging steel (grades 200, 250, 300, and 350) was determined in 3.5 percent salt (NaCl) solution, synthetic sea water, high humidity, and outside MSFC atmosphere. All grades of the maraging steel were found to be susceptible to SCC in varying degrees according to their strengths, with the lowest strength steel (grade 200) being the least susceptible and the highest strength steel (grade 350), the most susceptible to SCC. The SCC resistance of 250 grade maraging steel was also evaluated in salt and salt-chromate solutions using fracture mechanics techniques. The threshold value, K sub SCC, was found to be approximately 44 MN/sq m square root m, (40 ksi square root in.) or 40 percent of the K sub Q value.

Humphries, T. S.↗

Lunar Rotation and the Lunar Interior

Variations in rotation and orientation of the Moon are sensitive to solid-body tidal dissipation, dissipation due to relative motion at the fluid-core/ solid-mantle boundary, and tidal Love number k2. There is weaker sensitivity to flattening of the core-mantle boundary (CMB) and fluid core moment of inertia. Accurate Lunar Laser Ranging (LLR) measurements of the distance from observatories on the Earth to four retroreflector arrays on the Moon are sensitive to lunar rotation and orientation variations and tidal displacements. Past solutions using the LLR data have given results for dissipation due to solid-body tides and fluid core plus Love number. Past detection of CMB flattening has been marginal but is improving, while direct detection of the core moment has not yet been achieved. Three decades of Lunar Laser Ranging (LLR) data are analyzed using a weighted least-squares approach. The lunar solution parameters include dissipation at the fluid-core/solid-mantle boundary, tidal dissipation, dissipation-related coefficients for rotation and orientation terms, potential Love number k2, a correction to the constant term in the tilt of the equator to the ecliptic which is meant to approximate the influence of core-mantle boundary flattening, and displacement Love numbers h2 and l2. Several solutions, with different combinations of solution parameters and constraints, are considered.

Williams, J. G.↗

Morphology of Mars North Polar Ice Cap

The northern ice cap of Mars consists of a parabolic dome centered within 13 km of the pole, plus an arm-like ridge extending from the dome between about 135 and 225 east. Chasma Boreale lies between the dome and the extended ridge. The base of the dome is approximately elliptical with a major axis of 1100 km along the 90 east to 270 east direction and minor axis of 700 km along zero east to 180 deg. The heights of the dome and the extended ridge are respectively 2900 inches and 1700 inches above the surrounding basin. Least-squares fitting of a parabola through height profiles of the dome along longitudes 90 deg to 270 deg and zero deg to 180 deg gives an elliptic-paraboloid equation for the dome: Z(m) = 2800 - [(X-x)(exp 2)/113.6] - [(Y-y)(exp 2)/50.3], where X is the 90 deg to 270 deg axis, x = 9.90 km, y = 13.32 km, and the slightly-different fitted heights for the two axes are averaged. The center of the dome is shifted 13.32 km from the pole along zero deg longitude and 9.90 km along 90 deg longitude. Typical mean surface slopes on the ice cap are the order of 1/100 (0.6 deg), A small central portion of the cap, about 100 km by 200 km, extends in elevation about 200 inches above the parabolic shape of the cap. Additional information is contained in the original extended abstract.

Zwally, H. J.↗

A two-spacecraft test of a single spacecraft method of estimating shock normals

By assuming the validity of a subset of the Rankine-Hugoniot conservation relations for interplanetary (IP) shocks in an isotropic medium it has been demonstrated, in principle, that improved shock normals can be calculated by using a least-squares technique on combined magnetic field and plasma data from a single spacecraft. The scheme uses those six conservation relations not involving pressure and temperature. This paper deals with a test of the scheme by examining in detail a shock across which the magnetic field changed direction by a small amount (approximately 10 deg). On January 26, 1968 at about 1430 UT this shock was observed by the plasma and magnetic field instruments in Explorers 33 and 35. The spacecraft were 76.6 and 56.9 R sub E sunward of the earth, respectively (and 43.5 R sub E from each other), and therefore well outside the earth's bow shock region, a necessary condition for a valid test.

Lepping, R. P.↗