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.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Stress intensities for cracks emanating from pin-loaded holes

A series of stress freezing photoelastic experiments were conducted on large plates containing central holes with cracks emanating from the edge formed by the intersection of the hole with the plate surface. Loads were applied through rigid pins with neat fits in the holes. Stress-intensity factors (SIF) were estimated by a computer assisted least squares analysis of the photoelastic data taken from slices near the points of intersection of the flaw border with the hole boundary and the plate surface. Results indicate that the local mode of loading changes from Mode 1 near the hole boundary to mixed mode near the plate surface. The analysis is extended to include mixed mode loading, and results are compared with an existing approximate theory.

Smith, C. W.↗

NASA NTRS1992

Computer program for calculating and fitting thermodynamic functions

A computer program is described which (1) calculates thermodynamic functions (heat capacity, enthalpy, entropy, and free energy) for several optional forms of the partition function, (2) fits these functions to empirical equations by means of a least-squares fit, and (3) calculates, as a function of temperture, heats of formation and equilibrium constants. The program provides several methods for calculating ideal gas properties. For monatomic gases, three methods are given which differ in the technique used for truncating the partition function. For diatomic and polyatomic molecules, five methods are given which differ in the corrections to the rigid-rotator harmonic-oscillator approximation. A method for estimating thermodynamic functions for some species is also given.

Mcbride, Bonnie J.↗

NASA NTRS1994

Parallel Preconditioning for CFD Problems on the CM-5

Up to today, preconditioning methods on massively parallel systems have faced a major difficulty. The most successful preconditioning methods in terms of accelerating the convergence of the iterative solver such as incomplete LU factorizations are notoriously difficult to implement on parallel machines for two reasons: (1) the actual computation of the preconditioner is not very floating-point intensive, but requires a large amount of unstructured communication, and (2) the application of the preconditioning matrix in the iteration phase (i.e. triangular solves) are difficult to parallelize because of the recursive nature of the computation. Here we present a new approach to preconditioning for very large, sparse, unsymmetric, linear systems, which avoids both difficulties. We explicitly compute an approximate inverse to our original matrix. This new preconditioning matrix can be applied most efficiently for iterative methods on massively parallel machines, since the preconditioning phase involves only a matrix-vector multiplication, with possibly a dense matrix. Furthermore the actual computation of the preconditioning matrix has natural parallelism. For a problem of size n, the preconditioning matrix can be computed by solving n independent small least squares problems. The algorithm and its implementation on the Connection Machine CM-5 are discussed in detail and supported by extensive timings obtained from real problem data.

Simon, Horst D.↗

NASA NTRS2000

Preliminary development and evaluation of an algae-based air regeneration system

The potential of air regeneration system based on the growth of microalgae on the surface of porous ceramic tubes is evaluated. The algae have been maintained in the system for extended periods, up to 360 days. Preliminary measurements of the photosynthetic capacity have been made for Chlorella vulgaris (UTEX 259), Neospongiococcum punctatum (UTEX 786), Stichococcus sp., and Gloeocapsa sp. Under standard test conditions (photosynthetic photon flux approximately 66 micromoles m-2 s-1, initial CO2 concentration approximately 450 micromoles mol-1), mature tubes remove up to 0.2 micromoles of CO2 per tube per minute. The rate of removal increases with photon flux up to at least 225 micromoles m-2 s-1 (PPF); peak rates of 0.35 micromoles of CO2 per tube per minute have been achieved with Chlorella vulgaris. These rates correspond to between 120 and 210 micromoles of CO2 removed per square meter of projected area per minute.

Non-NASA Center↗

DOE OSTI2022

High-throughput characterization, correlation, and mapping of leaf photosynthetic and functional traits in the soybean ( Glycine max ) nested association mapping population

We report photosynthesis is a key target to improve crop production in many species including soybean [Glycine max (L.) Merr.]. A challenge is that phenotyping photosynthetic traits by traditional approaches is slow and destructive. There is proof-of-concept for leaf hyperspectral reflectance as a rapid method to model photosynthetic traits. However, the crucial step of demonstrating that hyperspectral approaches can be used to advance understanding of the genetic architecture of photosynthetic traits is untested. To address this challenge, we used full-range (500–2,400 nm) leaf reflectance spectroscopy to build partial least squares regression models to estimate leaf traits, including the rate-limiting processes of photosynthesis, maximum Rubisco carboxylation rate, and maximum electron transport. In total, 11 models were produced from a diverse population of soybean sampled over multiple field seasons to estimate photosynthetic parameters, chlorophyll content, leaf carbon and leaf nitrogen percentage, and specific leaf area (with R 2 from 0.56 to 0.96 and root mean square error approximately <10% of the range of calibration data). We explore the utility of these models by applying them to the soybean nested association mapping population, which showed variability in photosynthetic and leaf traits. Genetic mapping provided insights into the underlying genetic architecture of photosynthetic traits and potential improvement in soybean. Notably, the maximum Rubisco carboxylation rate mapped to a region of chromosome 19 containing genes encoding multiple small subunits of Rubisco. We also mapped the maximum electron transport rate to a region of chromosome 10 containing a fructose 1,6-bisphosphatase gene, encoding an important enzyme in the regeneration of ribulose 1,5-bisphosphate and the sucrose biosynthetic pathway. The estimated rate-limiting steps of photosynthesis were low or negatively correlated with yield suggesting that these traits are not influenced by the same genetic mechanisms and are not limiting yield in the soybean NAM population. Leaf carbon percentage, leaf nitrogen percentage, and specific leaf area showed strong correlations with yield and may be of interest in breeding programs as a proxy for yield. This work is among the first to use hyperspectral reflectance to model and map the genetic architecture of the rate-limiting steps of photosynthesis.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1986

Rain cell size statistics derived from radar observations at Wallops Island, Virginia

An investigation of two-dimensional rain cell size statistics has analyzed regression relations relating radar-determined rain rates to disdrometer data. This has yielded least-squares fits of radar reflectivity factors, and the application of a contouring program has generated 22,000 contours in which each isopleth belongs to predefined rain-rate intervals. An abundance of total and cell contours were observed belonging to all rain-rate categories. Both the computed number distributions and the conditional cumulative distributions as a function of contour diameter were found to be represented with good approximation by given exponential functions.

Goldhirsh, J.↗

NASA NTRS1997

Registration assisted mosaic generation

This paper presents a general strategy for assembling mosaics from numerous individual images where uncertainty exists in the position and orientation of those images. Both of the presented applications relate to remotely operated camera platforms, the first being the Galileo solid state imaging (SSI) camera presently in orbit around Jupiter, and the second being the Imager for Mars Pathfinder (IMP) stereo camera on Mars. A basic strategy in both applications is to determine the correct relative camera pointing followed by direct map projection of the images. It is assumed that approximate camera pointing exists sufficient to locate adjacent images and to place initial tiepoints within reach of the correlator. Spatial correlation is used to fix tiepoints whose initial locations are predicted by the camera pointing. We use either an fast fourier transform (fft) algorithm or a variant of Gruen's scheme permitting limited image rotation and skew. The Gruen correlator has three hierarchical modes: 1) A classical spatial least squares correlation on integral pixel boundaries used when rotation is small. 2) An annealing non-deterministic search used when rotations are unknown. A simplex deterministic search used for the end game. The correlation operation can be performed either interactively or autonomously. The final camera pointing solution relies upon a simplex downhill search in 2n or 3n dimensions where n is the number of images comprising the mosaic and the objective function to be minimized is the disagreement between tiepoint locations predicted from the camera pointing with those observed by the correlator. For Galileo the 3n unknowns are euler angles defining camera pointing in planet coordinates, and for Mars Pathfinder they are 2n unknowns representing commanded azimuth and elevation in the Lander coordinate system.

Lorre, Jean J.↗

DOE OSTI2022

Ion energy distribution in an electron beam ion trap inferred from simulations of the trapped ion cloud

For this work, we have inferred the energy distribution of trapped ions in an electron beam ion trap (EBIT) from simulations of the spatial distribution of Fe 13+ ions and a comparison with measured visible light images of the ion cloud. We simulated the cloud of Fe 13+ ions by computing ion trajectories in the EBIT for different ion energy distributions used to initialize the trajectories. We then performed a least-squares fit to infer the ion energy distribution that best reproduced the measured ion cloud. These best-fit distributions were typically non-Maxwellian. For electron beam energies of 395–475 eV and electron beam currents of 1–9 mA, we find that the average ion energy is in the range of 10–300 eV. We also find that the average ion energy increases with increasing beam current approximately as $\langle$E$\rangle$ ≈ 25 I e eV, where I e is the electron beam current in mA. We have also compared our results to Maxwell-Boltzmann-distribution ion clouds. We find that our best-fit non-thermal distributions have an $\langle$E$\rangle$ that is less than half that of the T from the best-fit Maxwell-Boltzmann distributions ($\langle$E$\rangle$/q)/T=0.41 ± 0.05.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1994

Mesoscale roughness of Venus

The global distribution of multi-kilometer (approx. 9 km) length scale 'roughness' (hereafter mesoscale roughness or MR) on Venus can be estimated from the Magellan global altimetry dataset (GxDR) and then compared with MR data derived for Earth from 5' ETOP5 data and for Mars (from USGS Mars DTM dataset). The mesoscale roughness parameter (MR) represents the RMS variance in meters of the actual planetary surface topography relative to the best fitting tangent plane defined on the basis of a 3x3 pixel sliding window. The best-fit plane was computed using a least-squares solution which minimizes delta H, the sum of the squares of the differences between the 9 local elevation values (H(sub i)), and the elevation of best-fit plane at the same grid location. Using the best-fit plane and delta H, we have computed the RMS 'roughness' var(delta R), where this parameter is always minimized on the basis of its calculation using least squares. We have called this 'ruggedness' parameter the Mesoscale Roughness (MR) because it is directly related to the high-frequency variance of topography after mesoscale slopes and tilts (i.e., for Venus, the baseline over which MR is computed (dx) is approx. 8.8 km and dx for Earth is approx. 9.3 km) are removed. As such, MR represents the degree to which a planetary surface is more rugged than approximately 10 km scale facets or tilts. It should not be confused with the radar 'RMS Roughness' parameter computed at 0.1 to 10 m length scales on the basis of the Magellan radar altimeter echo. We will use our MR parameter to investigate the global ruggedness properties of Venus as they relate to geological provinces and in comparison with the spatial pattern of MR for Earth and Mars.

Garvin, J. B.↗

NASA NTRS1994

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↗

NASA NTRS1977

Undulation and anomaly estimation using Geos-3 altimeter data without precise satellite orbits

Geos-3 altimeter data was used to obtain an approximate geoid undulation which was contaminated by long wavelength errors caused primarily by altimeter bias and orbit error. This error was reduced by fitting, with a low degree polynomial, the raw undulation data to the undulations implied by the GEM 7 potential coefficients in an adjustment process that included conditions on tracks that cross. These adjusted undulations were used to construct a geoid map in the Geos-3 calibration area using a least squares filter to remove the remaining noise in the undulations.

Rummel, R.↗

NASA NTRS2000

Baseline Calibration of the GeoSAR Interferometric Mapping Instrument

GeoSAR is a dual frequency, X-Band and P-Band, interferometric mapping radar flown on a Gulfstream II with a nominal flying altitude of 10000 and a swath width of 10 km on both left and right sides of the aircraft. The radar is designed to generate high resolution and high accuracy digital elevation maps both above and below the vegetation canopy. In order to generate the high accuracy digital elevations models it is necessary to have extremely precise knowledge of the interferometric baseline, approximately .1 mm for the X-band system and 3 mm for the P-band system. Baseline knowledge to this level of accuracy is achieved through a combination of active baseline measurement using a Laser Baseline Metrology System (LBMS) and baseline calibration using an accurately surveyed array of corner reflectors. In this paper, we describe the LBMS system used to measure the dynamic component of the interferometric baseline and the least squares procedure used to estimate all three components of the interferometric baseline based on differences between calculated and surveyed corner reflector locations.

Hensley, Scott↗

NASA NTRS2009

The Satellite Passive-Microwave Record of Sea Ice in the Ross Sea Since Late 1978

Satellites have provided us with a remarkable ability to monitor many aspects of the globe day-in and day-out and sea ice is one of numerous variables that by now have quite substantial satellite records. Passive-microwave data have been particularly valuable in sea ice monitoring, with a record that extends back to August 1987 on daily basis (for most of the period), to November 1970 on a less complete basis (again for most of the period), and to December 1972 on a less complete basis. For the period since November 1970, Ross Sea sea ice imagery is available at spatial resolution of approximately 25 km. This allows good depictions of the seasonal advance and retreat of the ice cover each year, along with its marked interannual variability. The Ross Sea ice extent typically reaches a minimum of approximately 0.7 x 10(exp 6) square kilometers in February, rising to a maximum of approximately 4.0 x 10(exp 6) square kilometers in September, with much variability among years for both those numbers. The Ross Sea images show clearly the day-by-day activity greatly from year to year. Animations of the data help to highlight the dynamic nature of the Ross Sea ice cover. The satellite data also allow calculation of trends in the ice cover over the period of the satellite record. Using linear least-squares fits, the Ross Sea ice extent increased at an average rate of 12,600 plus or minus 1,800 square kilometers per year between November 1978 and December 2007, with every month exhibiting increased ice extent and the rates of increase ranging from a low of 7,500 plus or minus 5,000 square kilometers per year for the February ice extents to a high of 20,300 plus or minus 6,100 kilometers per year for the October ice extents. On a yearly average basis, for 1979-2007 the Ross Sea ice extent increased at a rate of 4.8 plus or minus 1.6 % per decade. Placing the Ross Sea in the context of the Southern Ocean as a whole, over the November 1978-December 2007 period the Ross Sea had the highest rate of increase in sea ice coverage of any of five standard divisions of the Southern Ocean, although the Weddell Sea, Indian Ocean, and Western Pacific Ocean all also had sea ice increases, while only the Bellingshausen/Smundsen Seas experienced overall sea ice decreases. Overall, the Southern Ocean sea ice cover increased at an average rate of 10,800 plus or minus 2,500 square kilometers per year between November 1978 and December 2007, with every month showing positive values although with some of these values not being statistically significant. The sea ice increase since November 1978 was preceded by a sharp decrease in Southern Ocean ice coverage in the 1970's and is in marked contrast to the decrease in Arctic sea ice coverage that has occurred both in the period since November 1978 and since earlier in the 1970's. On a yearly average bases, for 1979-2007 the Southern Ocean sea ice extent increased at a rate of 1.0 plus or minus 0.4% per decade, whereas the Arctic ice extent decreased at the much greater rate of 4.0 plus or minus 0.4 percent per decade (closer to the % per decade rate of increase in the Ross Sea). Considerable research is ongoing to explain the differences.

Parkinson, Claire L.↗

NASA NTRS1972

Determination of stability derivatives from flight data using a Newton-Raphson minimization technique

A modified Newton-Raphson or quasilinearization minimization technique for determining stability derivatives from flight data was developed and compared with simple-equations, analog-matching, least-squares, and Shinbrot methods of analysis. For the data analyzed, the solutions computed by using the estimates obtained from the Newton-Raphson technique fit the data and determined coefficients adequately. A further modification to include a priori information was found to be useful. A model statistically similar to the flight data was analyzed using the same methods (excluding analog matching), and the Newton-Raphson technique was found to yield superior estimates. An approximate Cramer-Rao bound was compared with the error covariance matrix of the model and was found to provide information about the reliability of the individual estimates obtained. The technique was successfully applied to data obtained from a light airplane, a large supersonic airplane, and a lifting body vehicle. It was shown that the reliability of the estimates of a given coefficient obtained from these vehicles depends upon the data analyzed.

Iliff, K. W.↗

NASA NTRS1977

The analysis of GEOS-3 altimeter data in the Tasman and Coral seas

A technique was developed for preprocessing GEOS-3 altimetry data to establish a model of the regional sea surface. The algorithms developed models for a 35,000,000 sq km area with an internal precision of + or - 1 m. There were discrepancies between the sea surface model so obtained and GEM6 based geoid profiles with wavelengths of approximately 2500 km and amplitudes of up to 5 m in this region. The amplitudes were smaller when compared with GEM10-based geoid determinations. However, the comparison of 14 pairs of overlapping passes in the region indicated altimeter resolution of the + or - 25 cm level if the wavelength corresponding to the Nyquist frequency were 30 km. The spectral analysis of such comparisons indicated the existence of significant signal strength in the discrepancies after least squares fitting, with wavelengths in excess of 200 km.

Mather, R. S.↗

NASA NTRS2021

Performance Assessment of a Western Range 915-MHz Doppler Radar Wind Profiler

Please find enclosed the performance assessment for a Western Range 915-MHz Doppler Radar Wind Profiler (DRWP). This assessment was conducted following the installation of the new DRWP systems at the Western and Eastern Ranges. DRWP profiles were compared to Automated Meteorological Profiling System balloons launched from Vandenberg Space Force Base. Balloon data was collected from both normal synoptic and launch operations as well as some releases at the DRWP site. These results will include errors due to the spatial separation of winds. Root Mean Square wind speed deltas between the balloons and DRWP system mainly ranged between 1.0 and 2.3 m/s. Mean wind direction deltas mainly ranged between 10 and 28 degrees for wind speeds greater than or equal to 4 m/s. Data availability at discrete altitudes ranged from 20% to 94% with the least amount of available profiles being in the lower altitude ranges between 100 m and 500 m. The effective vertical resolution was determined to be approximately 250 m. The results from this analysis were requested by the United States Space Force for use of the acceptance of the DRWP system and for internal EV44 reference.

B J Barbre↗

DOE OSTI2023

CP decomposition for tensors via alternating least squares with QR decomposition

The CP tensor decomposition is used in applications such as machine learning and signal processing to discover latent low-rank structure in multidimensional data. Computing a CP decomposition via an alternating least squares (ALS) method reduces the problem to several linear least squares problems. The standard way to solve these linear least squares subproblems is to use the normal equations, which inherit special tensor structure that can be exploited for computational efficiency. However, the normal equations are sensitive to numerical ill-conditioning, which can compromise the results of the decomposition. In this paper, we develop versions of the CP-ALS algorithm using the QR decomposition and the singular value decomposition, which are more numerically stable than the normal equations, to solve the linear least squares problems. Our algorithms utilize the tensor structure of the CP-ALS subproblems efficiently, have the same complexity as the standard CP-ALS algorithm when the input is dense and the rank is small, and are shown via examples to produce more stable results when ill-conditioning is present. Our MATLAB implementation achieves the same running time as the standard algorithm for small ranks, and we show that the new methods can obtain lower approximation error.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2018

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗