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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Unsteady Aerodynamic Force Sensing from Strain Data

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm.

aerodynamic load sensing↗
NASA NTRS1960

Determination of Satellite Orbits from Radar Data

An optimum method for determining satellite orbits from radar data is presented in this report. Offering a good combination of speed and accuracy, the method makes use of orbit inclination and orbit elements in the plane, and proceeds with a differential correction of the orbit elements. Rapid, accurate methods of computing orbit elements are required to predict satellite positions for acquisition by other radars at points later along the orbit. In some cases the data are limited to a single pass over the observing station. The dynamical method is described in detail, and its accuracy is compared with those of two other methods: the purely geometrical, and the least-squares geometrical. By this optimum method the computing time, including the differential correction time, is 1 minute. Without differential correction, the rough determination takes from 10 to 20 seconds with approximately 5 miles positional uncertainty.

Cahill, W. F.↗
NASA NTRS2005

New Algorithms for Estimating Spacecraft Position Using Scanning Techniques for Deep Space Network Antennas

As more and more nonlinear estimation techniques become available, our interest is in finding out what performance improvement, if any, they can provide for practical nonlinear problems that have been traditionally solved using linear methods. In this paper we examine the problem of estimating spacecraft position using conical scan (conscan) for NASA's Deep Space Network antennas. We show that for additive disturbances on antenna power measurement, the problem can be transformed into a linear one, and we present a general solution to this problem, with the least square solution reported in literature as a special case. We also show that for additive disturbances on antenna position, the problem is a truly nonlinear one, and we present two approximate solutions based on linearization and Unscented Transformation respectively, and one 'exact' solution based on Markov Chain Monte Carlo (MCMC) method. Simulations show that, with the amount of data collected in practice, linear methods perform almost the same as MCMC methods. It is only when we artificially reduce the amount of collected data and increase the level of noise that nonlinear methods show significantly better accuracy than that achieved by linear methods, at the expense of more computation.

Metropolis-Hastings (MH)↗
NASA NTRS2010

Monte Carlo Modeling of Sodium in Mercury's Exosphere During the First Two MESSENGER Flybys

We present a Monte Carlo model of the distribution of neutral sodium in Mercury's exosphere and tail using data from the Mercury Atmospheric and Surface Composition Spectrometer (MASCS) on the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft during the first two flybys of the planet in January and September 2008. We show that the dominant source mechanism for ejecting sodium from the surface is photon-stimulated desorption (PSD) and that the desorption rate is limited by the diffusion rate of sodium from the interior of grains in the regolith to the topmost few monolayers where PSD is effective. In the absence of ion precipitation, we find that the sodium source rate is limited to approximately 10(exp 6) - 10(exp 7) per square centimeter per second, depending on the sticking efficiency of exospheric sodium that returns to the surface. The diffusion rate must be at least a factor of 5 higher in regions of ion precipitation to explain the MASCS observations during the second MESSENGER f1yby. We estimate that impact vaporization of micrometeoroids may provide up to 15% of the total sodium source rate in the regions observed. Although sputtering by precipitating ions was found not to be a significant source of sodium during the MESSENGER flybys, ion precipitation is responsible for increasing the source rate at high latitudes through ion-enhanced diffusion.

Burger, Matthew H.↗
NASA NTRS2012

Chandra Observations of Galaxy Zoo Mergers: Frequency of Binary Active Nuclei in Massive Mergers

We present the results from a Chandra pilot study of 12 massive mergers selected from Galaxy Zoo. The sample includes major mergers down to a host galaxy mass of 10(sup 11) solar mass that already have optical AGN signatures in at least one of the progenitors. We find that the coincidences of optically selected active nuclei with mildly obscured (N(sub H) less than or approximately 1.1 x 10(exp 22) per square centimeter) X-ray nuclei are relatively common (8/12), but the detections are too faint (less than 40 counts per nucleus; f(sub 2-10 keV) less than or approximately 1.2 x 10(exp -13) ergs per second per square centimeter) to separate starburst and nuclear activity as the origin of the X-ray emission. Only one merger is found to have confirmed binary X-ray nuclei, though the X-ray emission from its southern nucleus could be due solely to star formation. Thus, the occurrences of binary AGN in these mergers are rare (0-8%), unless most merger-induced active nuclei are very heavily obscured or Compton thick.

Teng, Stacy H.↗
NASA NTRS1985

Venus Gravity: Global Field Results

The Pioneer Venus Orbiter has provided a vast amount of gravity data, since March 1979. High resolution results have revealed the high correlation between topography and gravity. These data were acquired at relatively low spacecraft altitudes (150 km) where atmospheric effects and high frequency variations were significant and modelling with global spherical harmonics was difficult. During 1982 data were acquired over a complete Venus rotation where the lowest altitudes were 1000 km or more, thus removing atmospherics and high frequency gravity effects. These data that have been reduced to produce a tenth degree and order spherical harmonic model of the global gravity field of Venus. The reduction technique used a least squares where 78 independent arcs of data were combined in a simultaneous inversion. The tenth degree and order gravity field solution has 117 parameters describing the spherical harmonic coefficients. This provides approximately 1800 km feature resolution. These coefficients are presented as a geoid map.

Mottinger, N. A.↗
NASA NTRS1988

Analysis of positron lifetime spectra in polymers

A new procedure for analyzing multicomponent positron lifetime spectra in polymers was developed. It requires initial estimates of the lifetimes and the intensities of various components, which are readily obtainable by a standard spectrum stripping process. These initial estimates, after convolution with the timing system resolution function, are then used as the inputs for a nonlinear least squares analysis to compute the estimates that conform to a global error minimization criterion. The convolution integral uses the full experimental resolution function, in contrast to the previous studies where analytical approximations of it were utilized. These concepts were incorporated into a generalized Computer Program for Analyzing Positron Lifetime Spectra (PAPLS) in polymers. Its validity was tested using several artificially generated data sets. These data sets were also analyzed using the widely used POSITRONFIT program. In almost all cases, the PAPLS program gives closer fit to the input values. The new procedure was applied to the analysis of several lifetime spectra measured in metal ion containing Epon-828 samples. The results are described.

Singh, Jag J.↗
NASA NTRS2005

Data Assimilation and Data Fusion for Planetary Atmospheres

The overarching goal of this Cooperative Agreement was to develop a model and procedures for the data assimilation of planetary spacecraft atmospheric observations. Data assimilation - in its application to weather analysis and prediction - is the process of finding an initial state of the meteorological variables (winds, temperatures, pressures, etc.) of an atmosphere, which, when propagated forward in time using a deterministic general circulation model, reproduces all of the available observations over that time to within the measurement and computational errors. With this definition, data assimilation is seen to be a natural extension of well-known least-squares minimization techniques. The primary complication arises from the scale of the problem: For the Martian atmosphere with the available nadir-viewing Thermal Emission Spectrometer data from Mars Global Surveyor, approximately 1,000,000 individual measurements of channel radiances (in the 15-micrometer region, where these radiances relate directly to the surface and atmospheric temperature) were made per day. A suitable general circulation model for dealing with this data set has on the order of 20,000 independent variables. After some spatial and temporal averaging of the data - which provides a necessary statistical estimate of the representativeness of the measurements, a crucial issue in data assimilation - the problem reduces in scale to the solution of approximately 50,000 equations for the 20,000 variables.

Houben, Howard↗
NASA NTRS1994

Isostatic compensation of equatorial highlands on Venus

Spherical harmonic models for Venus' global topography and gravity incorporating Magellan data are used to test isostatic compensation models in five 30 deg x 30 deg regions representative of the main classes of equatorial highlands. The power spectral density for the harmonic models obeys a power-law scaling with spectral slope Beta approximately 2 (Brown noise) for the topography and Beta approximately 3 (Kaula's law) for the geoid, similar to what is observed for Earth. The Venus topography spectrum has lower amplitudes than Earth's which reflects the dominant lowland topography on Venus. Observed degree geoid to topography ratios (GTRs) on Venus are significantly smaller than degree GTRs for uncompensated topography, indicative of substantial compensation. Assuming a global Airy compensation, most of the topography is compensated at depths greater than 100 km, suggesting a thick lithosphere on Venus. For each region considered we obtain a regional degree of compensation C from a linear regression of Bouguer anomaly versus Bouguer gravity data. Geoid anomaly (N) versus topography variation (h) data for each sample were compared, in the least-squares sense, to theoretical correlations for Pratt, Airy, and thermal thinning isostasy models yielding regional GTR, zero-elevation crustal thickness (H), and zero elevation thermal lithosphere thickness (y(sub L(sub 0)), respectively. We find the regional compensation to be substantial (C approximately 52-80%), and the h, N data correlations in the chosen areas can be explained by isostasy models applicable on the Earth and involving variations in crustal thickness (Airy) and/or lithospheric (thermal thinning) thickness. However, a thick crust and lithosphere (y(sub L(sub 0)) approximately 300 km) must be assumed for Venus.

Kucinskas, Algis B.↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗
NASA NTRS1993

An improved calibration technique for wind tunnel model attitude sensors

Aerodynamic wind tunnel tests at NASA Langley Research Center (LaRC) require accurate measurement of model attitude. Inertial accelerometer packages have been the primary sensor used to measure model attitude to an accuracy of +/- 0.01 deg as required for aerodynamic research. The calibration parameters of the accelerometer package are currently obtained from a seven-point tumble test using a simplified empirical approximation. The inaccuracy due to the approximation exceeds the accuracy requirement as the misalignment angle between the package axis and the model body axis increases beyond 1.4 deg. This paper presents the exact solution derived from the coordinate transformation to eliminate inaccuracy caused by the approximation. In addition, a new calibration procedure is developed in which the data taken from the seven-point tumble test is fit to the exact solution by means of a least-squares estimation procedure. Validation tests indicate that the new calibration procedure provides +/- 0.005-deg accuracy over large package misalignments, which is not possible with the current procedure.

Tripp, John S.↗
NASA NTRS1996

Exposure History of Lunar Meteorites Queen Alexandra Range 93069 and 94269

Cosmic-ray produced C-14 (t(sub 1/2) = 5730 years), 36Cl (3.01 x 10(exp 5 years), Al-26 (7.05 x 10(exp 5 years), and Be-10 (1.5 x 10(exp 6 years) in the recently discovered lunar meteorites Queen Alexandra Range 93069 (QUE 93069) and 94269 (QUE 94269) were measured by accelerator mass spectrometry. The abundance pattern of these four cosmogenic radionuclides and of noble gases indicates QUE 93069 and QUE 94269 were a paired fall and were exposed to cosmic rays near the surface of the Moon for at least several hundred million years before ejection. After the meteorite was launched from the Moon, where it had resided at a depth of 65-80 g/cm square, it experienced a short transition time, approximately 20-50 ka, before colliding with the Earth. The terrestrial age of the meteorite is 5-10 ka. Comparison ofthe cosmogenic nuclide concentrations in QUE 93069/94269 and MAC 88104/88105 clearly shows that these meteorites were not ejected by a common event from the Moon.

Nishiizumi, K.↗
NASA NTRS2004

The Sombrero: An Edge-On Spiral: Halo Dynamics, A Massive Dark Object, and XRBS

We present an analysis of the XMM-Newton and Chandra ACIS-S observations of the LINER nucleus of the Sombrero galaxy, and we discuss possible explanations for its very Sub-Eddington luminosity by complementing the X-ray results with high angular resolution observations in other bands. The X-ray investigation shows a hard (gamma=1.89) and moderately absorbed (NH=1.8x1021 per square centimeter) nuclear source of 1.5x1040 ergs per second in the 2-10 keV band, surrounded by hot gas at a temperature of approximately 0.6 keV. The bolometric nuclear luminosity is at least approximately 200 times lower than expected if mass accreted on the supermassive black hole, which Hubble Space Telescope shows to reside at the center of this galaxy, at the rate predicted by the spherical and adiabatic Bondi accretion theory and with the high radiative efficiency of a standard accretion disk. The low luminosity, coupled to the observed absence of Fe K emission in the nuclear spectrum, indicates that such a disk is not present. This nucleus also differs from bright unobscured active galactic nuclei in the lack of high flux variability and prominent broad H emission. However, it is also too faint for the predictions of simple radiatively inefficient accretion taking place at the Bondi rate; it could instead be too radio bright for radiatively inefficient accretion that includes strong mass outflows or convection. This discrepancy could be solved by the possible presence of nuclear radio jets. An alternative explanation of the low luminosity, in place of radiative inefficiency, could be unsteady accretion.

Fabbiano, G.↗
DOE OSTI2022

Nonlinear dimension reduction for surrogate modeling using gradient information

We introduce a method for the nonlinear dimension reduction of a high-dimensional function $u:{\mathbb{R}}^d\rightarrow{\mathbb{R}}$, $d\gg 1$. Our objective is to identify a nonlinear feature map $g:{\mathbb{R}}^d\rightarrow{\mathbb{R}}^m$, with a prescribed intermediate dimension $m\ll d$, so that $u$ can be well approximated by $f\circ g$ for some profile function $f:{\mathbb{R}}^m\rightarrow{\mathbb{R}}$. We propose to build the feature map by aligning the Jacobian $\nabla g$ with the gradient $\nabla u$, and we theoretically analyze the properties of the resulting $g$. Once $g$ is built, we construct $f$ by solving a gradient-enhanced least squares problem. Our practical algorithm uses a sample $\{{\textbf{x}}^{(i)},u({\textbf{x}}^{(i)}),\nabla u({\textbf{x}}^{(i)})\}_{i=1}^N$ and builds both $g$ and $f$ on adaptive downward-closed polynomial spaces, using cross validation to avoid overfitting. We numerically evaluate the performance of our algorithm across different benchmarks, and explore the impact of the intermediate dimension $m$. We show that building a nonlinear feature map $g$ can permit more accurate approximation of $u$ than a linear $g$, for the same input data set.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Projection-based multifidelity linear regression for data-scarce applications

Surrogate modeling for systems with high-dimensional quantities of interest remains challenging, particularly when training data are costly to acquire. This work develops multifidelity methods for multiple-input multiple-output linear regression targeting data-limited applications with high-dimensional outputs. Multifidelity methods integrate many inexpensive low-fidelity model evaluations with limited, costly high-fidelity evaluations. We introduce two projection-based multifidelity linear regression approaches with linear and nonlinear features that leverage principal component basis vectors for dimensionality reduction and combine multifidelity data through: (i) a direct data augmentation using low-fidelity data, and (ii) a data augmentation incorporating explicit linear corrections between low-fidelity and high-fidelity data. The data augmentation approaches combine high-fidelity and low-fidelity data into a unified training set and train the linear regression model through weighted least squares with fidelity-specific weights. We introduce a proximity-based weighting scheme with automatic weight selection strategy through cross-validation. Here, the proposed multifidelity linear regression methods are demonstrated on approximating the surface pressure field of a hypersonic vehicle in flight and the temperature field on an aircraft disc braking system. In an ultra low-data regime of no more than twelve high-fidelity samples, multifidelity linear regression achieves approximately 2% – 12% improvement in median accuracy and a higher R 2 score relative to single-fidelity methods at comparable computational cost.

data augmentation↗
NASA NTRS1983

A statistical study of large electric field events in the earth's magnetotail

A statistical study of large electric field events, measured on ISEE 1 at geocentric distances of 3 to 23 earth radii, has been performed in order to determine the occurrence frequency and spatial distribution of such events. Large electric fields were seen at all altitudes, with the lowest altitude events occurring on auroral field lines and the higher altitude events occurring on the high latitude boundary of the plasma sheet. The largest events had amplitudes that were approximately proportional to the square root of the local magnetic field strength, in agreement with expectations of a simple mapping model. The apparent concentration of big field events near local midnight is caused, at least in part, by a nonuniform frequency of satellite crossings through the plasma sheet boundary as a function of local time. Thus, the spatial distribution and amplitudes of the large field events are not inconsistent with the mapping of low-altitude electrostatic shocks to the plasma sheet boundary along equipotential magnetic field lines.

Levin, S.↗
NASA NTRS1986

Modified superposition: A simple time series approach to closed-loop manual controller identification

Single-channel pilot manual control output in closed-tracking tasks is modeled in terms of linear discrete transfer functions which are parsimonious and guaranteed stable. The transfer functions are found by applying a modified super-position time series generation technique. A Levinson-Durbin algorithm is used to determine the filter which prewhitens the input and a projective (least squares) fit of pulse response estimates is used to guarantee identified model stability. Results from two case studies are compared to previous findings, where the source of data are relatively short data records, approximately 25 seconds long. Time delay effects and pilot seasonalities are discussed and analyzed. It is concluded that single-channel time series controller modeling is feasible on short records, and that it is important for the analyst to determine a criterion for best time domain fit which allows association of model parameter values, such as pure time delay, with actual physical and physiological constraints. The purpose of the modeling is thus paramount.

Biezad, D. J.↗