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 307 records · Page 17

On in situ 'calibration' of shipboard ADCPs

Methods for in situ calibration of acoustic-Doppler current profilers (ADCPs) are considered for measurement of absolute current profiles from a moving ship. Errors are of two types: sensitivity and alignment. Least square error estimates are given for experimental determination of both factors, for use in the 'water track' or 'bottom tracking' mode. Errors in the estimation of either factor may lead to large errors in derived water velocities, although the major contributions of the two factors arise from different sources and are approximately orthogonal to one another.

Joyce, Terrence M.↗

Radiation Pressure Measurements on Micron-Size Individual Dust Grains

Measurements of electromagnetic radiation pressure have been made on individual silica (SiO2) particles levitated in an electrodynamic balance. These measurements were made by inserting single charged particles of known diameter in the 0.2- to 6.82-micron range and irradiating them from above with laser radiation focused to beam widths of approximately 175- 400 microns at ambient pressures particle due to the radiation force is balanced by the electrostatic force indicated by the compensating dc potential applied to the balance electrodes, providing a direct measure of the radiation force on the levitated particle. Theoretical calculations of the radiation pressure with a least-squares fit to the measured data yield the radiation pressure efficiencies of the particles, and comparisons with Mie scattering theory calculations provide the imaginary part of the refractive index of SiO2 and the corresponding extinction and scattering efficiencies.

Abbas, M. M.↗

Adaptive Filtering Using Recurrent Neural Networks

A method for adaptive (or, optionally, nonadaptive) filtering has been developed for estimating the states of complex process systems (e.g., chemical plants, factories, or manufacturing processes at some level of abstraction) from time series of measurements of system inputs and outputs. The method is based partly on the fundamental principles of the Kalman filter and partly on the use of recurrent neural networks. The standard Kalman filter involves an assumption of linearity of the mathematical model used to describe a process system. The extended Kalman filter accommodates a nonlinear process model but still requires linearization about the state estimate. Both the standard and extended Kalman filters involve the often unrealistic assumption that process and measurement noise are zero-mean, Gaussian, and white. In contrast, the present method does not involve any assumptions of linearity of process models or of the nature of process noise; on the contrary, few (if any) assumptions are made about process models, noise models, or the parameters of such models. In this regard, the method can be characterized as one of nonlinear, nonparametric filtering. The method exploits the unique ability of neural networks to approximate nonlinear functions. In a given case, the process model is limited mainly by limitations of the approximation ability of the neural networks chosen for that case. Moreover, despite the lack of assumptions regarding process noise, the method yields minimum- variance filters. In that they do not require statistical models of noise, the neural- network-based state filters of this method are comparable to conventional nonlinear least-squares estimators.

Parlos, Alexander G.↗

Least squares inversion of Stokes profiles in the presence of velocity gradients

The Auer, Heasley and House Stokes inversion procedure in use at High Altitude Observatory is based on the analytic solution of the equation of transfer for polarized light where the representation of the thermodynamic and magnetic structure of the atmosphere is assumed to have a high degree of invariance, namely, a Milne-Eddington (ME) structure with a constant magnetic field. In the presence of invariance breaking gradients the resultant Stokes profiles are represented only approximately, if at all, by analytic forms. The accuracy of the inversion parameters and their significance as measures of actual structure are explored for the ME and the Landman-Finn sunspot models under the effects of velocity gradients. The resulting field parameters are good to a few percent and prove to be insensitive to the errors committed by the use of a ME-representation, but the resulting ME parameters yield a less precise measure of thermal structure.

Skumanich, A.↗

Thermophysical properties of Apollo 14 fines

The vacuum thermal conductivity of lunar fines sample 14163 was measured for the approximate temperature range 100 to 400 deg K. Several sample densities from 1100 to 1800 kg/cu m were used. The temperature dependence of the conductivity was found to be well represented by the relation k = A + BT cubed, which is predicted by elementary theory. The coefficients A and B were obtained by least squares analysis of the data. The thermal diffusivity was calculated for the various densities by using specific heat data from the literature along with the measured conductivities. The results are compared with those obtained for Apollo 11, Apollo 12, and terrestrial basalt samples.

Cremers, C. J.↗

Evaluation of LANDSAT-D Orbit Determination Using a Filter/Smoother (PREFER)

Simulated range and range rate data for five tracking stations were first generated using batch least squares orbit determination (GTDS). Then GTDS was used (in the differential correction mode) to produce a nominal trajectory which was input to PREFER. The GTDS differential correction (DC) run was made using models which differed from those used to produce the simulated data. These model differences were chosen to be fairly realistic approximations to the errors in the models actually used for operational orbit determination. Several different simulation runs were made with different types of model errors in order to determine the sensitivity to these errors. The nominal trajectory and the simulated measurement data were input to PREFER to produce a smoothed ephemeris file. Numerous runs of PREFER were made in which parameters describing the statistics of the model errors were varied. The likelihood function computed by the Kalman filter determined the ""best'' choice of input parameters. There was strong negative correlation between the likelihood function and the errors in the smoothed ephemeris.

Gibbs, B. P.↗

Separation of core and crustal magnetic field sources

Fluid motions in the electrically conducting core and magnetized crustal rocks are the two major sources of the magnetic field observed on or slightly above the Earth's surface. The exact separation of these two contributions is not possible without imposing a priori assumptions about the internal source distribution. Nonetheless models like these were developed for hundreds of years Gauss' method, least squares analysis with a truncated spherical harmonic expansion was the method of choice for more than 100 years although he did not address separation of core and crustal sources, but rather internal versus external ones. Using some arbitrary criterion for appropriate truncation level, we now extrapolate downward core field models through the (approximately) insulating mantle. Unfortunately our view can change dramatically depending on the degree of truncation for describing core sources.

Shure, L.↗

Detection of Non-Symmetrical Damage in Smart Plate-Like Structures

A two-dimensional model for in-plane vibrations of a cantilever plate with a non-symmetrical damage is used in the context of defect identification in materials with piezoelectric ceramic patches bonded to their surface. These patches can act both as actuators and sensors in a self-analyzing fashion, which is a characteristic of smart materials. A Galerkin method is used to approximate the dynamic response of these structures. The natural frequency shifts due to the damage are estimated numerically and compared to experimental data obtained from tests on cantilever aluminum plate-like structures damaged at different locations with defects of different depths. The damage location and extent are determined by an enhanced least square identification method. Efficacy of the frequency shift based algorithms is demonstrated using experimental data.

Blanks, H. T.↗

Variable-permittivity linear inverse problem for the H(sub z)-polarized case

The H(sub z)-polarized inverse problem has rarely been studied before due to the complicated way in which the unknown permittivity appears in the wave equation. This problem is equivalent to the acoustic inverse problem with variable density. We have recently reported the solution to the nonlinear variable-permittivity H(sub z)-polarized inverse problem using the Born iterative method. Here, the linear inverse problem is solved for permittivity (epsilon) and permeability (mu) using a different approach which is an extension of the basic ideas of diffraction tomography (DT). The key to solving this problem is to utilize frequency diversity to obtain the required independent measurements. The receivers are assumed to be in the far field of the object, and plane wave incidence is also assumed. It is assumed that the scatterer is weak, so that the Born approximation can be used to arrive at a relationship between the measured pressure field and two terms related to the spatial Fourier transform of the two unknowns, epsilon and mu. The term involving permeability corresponds to monopole scattering and that for permittivity to dipole scattering. Measurements at several frequencies are used and a least squares problem is solved to reconstruct epsilon and mu. It is observed that the low spatial frequencies in the spectra of epsilon and mu produce inaccuracies in the results. Hence, a regularization method is devised to remove this problem. Several results are shown. Low contrast objects for which the above analysis holds are used to show that good reconstructions are obtained for both permittivity and permeability after regularization is applied.

Moghaddam, M.↗

The spatial distribution of gamma-ray bursts from Pioneer Venus Orbiter

The homogeneity of the burst source population in space with the log N(greater than P sub c)-log P sub c and log (dN/dP sub c)-log P sub c relations is studied on the basis of a sample of 225 observed gamma-ray bursts between September 1978 and July 1988 from the Pioneer Venus Orbiter. The relations N(greater than P sub c) is approximately equal to P sub c exp -1.5 and dN/dP sub c is approximately equal to P sub c exp -2.5 are expected if the gamma bursters are uniformly distributed in space. The Kolmogorov-Smirnov statistic gives a probability of 0.93 that the observed integral distribution N(greater than P sub c) is drawn from the expected uniform distribution P sub c exp -1.5. A least-squares 1 standard deviation fit to the observed differential distribution yields dN/dP sub c of about P sub c -2.42 +/- 0.07. Both are consistent with a uniform distribution in space of the parent population of burst sources.

Chuang, K. W.↗

Maximizing TDRS Command Load Lifetime

The GNC software onboard ISS utilizes TORS command loads, and a simplistic model of TORS orbital motion to generate onboard TORS state vectors. Each TORS command load contains five "invariant" orbital elements which serve as inputs to the onboard propagation algorithm. These elements include semi-major axis, inclination, time of last ascending node crossing, right ascension of ascending node, and mean motion. Running parallel to the onboard software is the TORS Command Builder Tool application, located in the JSC Mission Control Center. The TORS Command Builder Tool is responsible for building the TORS command loads using a ground TORS state vector, mirroring the onboard propagation algorithm, and assessing the fidelity of current TORS command loads onboard ISS. The tool works by extracting a ground state vector at a given time from a current TORS ephemeris, and then calculating the corresponding "onboard" TORS state vector at the same time using the current onboard TORS command load. The tool then performs a comparison between these two vectors and displays the relative differences in the command builder tool GUI. If the RSS position difference between these two vectors exceeds the tolerable lim its, a new command load is built using the ground state vector and uplinked to ISS. A command load's lifetime is therefore defined as the time from when a command load is built to the time the RSS position difference exceeds the tolerable limit. From the outset of TORS command load operations (STS-98), command load lifetime was limited to approximately one week due to the simplicity of both the onboard propagation algorithm, and the algorithm used by the command builder tool to generate the invariant orbital elements. It was soon desired to extend command load lifetime in order to minimize potential risk due to frequent ISS commanding. Initial studies indicated that command load lifetime was most sensitive to changes in mean motion. Finding a suitable value for mean motion was therefore the key to achieving this goal. This goal was eventually realized through development of an Excel spreadsheet tool called EMMIE (Excel Mean Motion Interactive Estimation). EMMIE utilizes ground ephemeris nodal data to perform a least-squares fit to inferred mean anomaly as a function of time, thus generating an initial estimate for mean motion. This mean motion in turn drives a plot of estimated downtrack position difference versus time. The user can then manually iterate the mean motion, and determine an optimal value that will maximize command load lifetime. Once this optimal value is determined, the mean motion initially calculated by the command builder tool is overwritten with the new optimal value, and the command load is built for uplink to ISS. EMMIE also provides the capability for command load lifetime to be tracked through multiple TORS ephemeris updates. Using EMMIE, TORS command load lifetimes of approximately 30 days have been achieved.

Brown, Aaron J.↗

Thermophysical properties of Apollo 14 fines

The vacuum thermal conductivity of lunar fines sample 14163 was measured for the approximate temperature range of 100 to 400 K. Sample densities of 1500 kg/cu m and 1800 kg/cu m were used. The temperature dependence of the conductivity was found to be well represented by the relation k = A + BT-cubed, which is predicted by elementary theory. The coefficients A and B were obtained by least-squares analysis of the data. The thermal diffusivity was calculated for the various densities using specific heat data from the literature along with the measured conductivities. The results are compared with those obtained for Apollo 11, Apollo 12, and terrestrial basalt samples.

Cremers, C. J.↗

Status on Iterative Transform Phase Retrieval Applied to the GBT Data

This slide presentation reviews the use of iterative transform phase retrieval in the analysis of the Green Bank Radio Telescope (GBT) Data. It reviews the NASA projects that have used phase retrieval, and the testbed for the algorithm to be used for the James Webb Space Telescope. It shows the comparison of phase retrieval with an interferometer, and reviews the two approaches used for phase retrieval, iterative transform (ITA) or parametric (non-linear least squares model fitting). The concept of ITA Phase Retrieval is reviewed, and the application to Radio Antennas is reviewed. The presentation also examines the National Radio Astronomy Observatory (NRAO) data from the GBT, and the Fourier model that NRAO uses to analyze the data. The challenge for ITA phase retrieval is reviewed, and the coherent approximation for incoherent data is shown. The validity of the approximation is good for a large tilt. There is a review of the proof of concept of the Phase Review simulation using the input wavefront, and the initial sampling parameters estimate from the focused GBT data.

Dean, Bruce↗

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗

Vibration-rotation intensities of SiO

Dipole moment matrix elements have been computed for a large number of transitions of astrophysical interest for the more abundant isotopes of SiO. The wave functions utilized were obtained from a direct solution of the Schroedinger equation with an accurate RKR potential. The dipole moment function, in the form of a Pade approximant, was chosen to reproduce the experimental measurements near equilibrium, to have the proper united and separated atom limits, and to have the correct long-range asymptotic dependence on internuclear separation. Because of the large number of transitions involved, and to facilitate applications, the squares of the dipole moment matrix elements were fitted by a least-squares procedure to polynomials in v and J. In addition, Einstein A coefficients are given for observed maser transitions and for selected vibrational bands. These latter are compared with previous calculations, and it is concluded that for the higher Delta v transitions, the present results represent a significant improvement.

Tipping, R. H.↗

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗

Testing an Algae-Based Air-Regeneration System

The potential of an air-regeneration system based on the growth of unicellular algae on the surface of porous ceramic tubes was evaluated. The system is fairly robust with respect to environmental conditions and is capable of maintaining algal cultures for up to 365 days. Under standard conditions (50-66 micro mol/sq mm s (PPF), 450 micro mol mol of CO2), mature tubes can remove CO2 at a rate of up to 90 micro mol/sq m min. Under these conditions, approximately 200 square meters of area would be required for each member of the crew. However, the rate of uptake increases with both photon flux and CO2 concentration in accordance with Michaelis-Menton dynamics. An extrapolation to conditions of saturating light and carbon dioxide indicates that the area required can be reduced by a factor of at least 2.5.

Nienow, James↗

Development of MLPG and LBIE Methods for Nonlinear Problems of Fracture

The truly meshless local Petrov-Galerkin (MLPG) method holds a great promise in solving boundary value problems, using a local symmetric weak form as a natural approach. In the present paper, in the context of MLPG and the meshless interpolation of a moving least squares (MLS) type, a method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The sub-domains for the shape functions in the MLS approximation are defined only from the primary nodes, and the secondary nodes use the same sub-domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and/or moved without an additional mesh. The essential boundary conditions can be imposed exactly, and non-convex boundaries can be treated without special techniques. Several numerical examples are presented to illustrate the performance of the present method.

Atluri, Satya N.↗