Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sum of squares”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Analysis of Sting Balance Calibration Data Using Optimized Regression Models

Calibration data of a wind tunnel sting balance was processed using a search algorithm that identifies an optimized regression model for the data analysis. The selected sting balance had two moment gages that were mounted forward and aft of the balance moment center. The difference and the sum of the two gage outputs were fitted in the least squares sense using the normal force and the pitching moment at the balance moment center as independent variables. The regression model search algorithm predicted that the difference of the gage outputs should be modeled using the intercept and the normal force. The sum of the two gage outputs, on the other hand, should be modeled using the intercept, the pitching moment, and the square of the pitching moment. Equations of the deflection of a cantilever beam are used to show that the search algorithm s two recommended math models can also be obtained after performing a rigorous theoretical analysis of the deflection of the sting balance under load. The analysis of the sting balance calibration data set is a rare example of a situation when regression models of balance calibration data can directly be derived from first principles of physics and engineering. In addition, it is interesting to see that the search algorithm recommended the same regression models for the data analysis using only a set of statistical quality metrics.

Ulbrich, Norbert↗
DOE OSTI2025

Using the optimal combined index weight ratio to improve the probability of anomaly detection in big area additive manufacturing

Big Area Additive Manufacturing (BAAM) of composites requires significant time, energy, and material, so it is critical to reduce production inefficiencies to make functional parts without multiple iterations. Statistical process control coupled with Principal Component Analysis (PCA) is a powerful technique that provides a quick, computationally inexpensive, and intuitive way for operators to detect defects that form in a manufacturing process without massive datasets. Recently, a combined index that is a weighted sum of the Hotelling's T 2 and squared residual error statistics has been proposed that can be monitored in one chart, improving interpretation accuracy and simplicity. However, the literature does not offer a formal method to optimise the weights. Here, we introduce two new approaches to the traditional weight selection approach using simulated and BAAM image data. Approach 1 uses a theoretically motivated optimum inspired by probabilistic principal component analysis. Approach 2 systematically varies the ratio of the weights to find the optimum. We show that approach 1 delivers optimal anomaly detection performance in select cases while approach 2 fares better in practice. Surprisingly, we also show that choosing a more complex PCA model has a minimal negative impact on anomaly detection performance compared to a more simplistic model.

3-dimensional printing↗
NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗
NASA NTRS1996

First-Order System Least Squares for Velocity-Vorticity-Pressure Form of the Stokes Equations, with Application to Linear Elasticity

In this paper, we study the least-squares method for the generalized Stokes equations (including linear elasticity) based on the velocity-vorticity-pressure formulation in d = 2 or 3 dimensions. The least squares functional is defined in terms of the sum of the L(exp 2)- and H(exp -1)-norms of the residual equations, which is weighted appropriately by by the Reynolds number. Our approach for establishing ellipticity of the functional does not use ADN theory, but is founded more on basic principles. We also analyze the case where the H(exp -1)-norm in the functional is replaced by a discrete functional to make the computation feasible. We show that the resulting algebraic equations can be uniformly preconditioned by well-known techniques.

Cai, Zhiqiang↗
DOE OSTI2020

Hydrodynamic electron transport near charge neutrality

We develop the theory of hydrodynamic electron transport in a long-range disorder potential for conductors in which the underlying electron liquid lacks Galilean invariance. For weak disorder, we express the transport coefficients of the system in terms of the intrinsic kinetic coefficients of the electron liquid and the correlation function of the disorder potential. We apply these results to analyze the doping and temperature dependence of transport coefficients of graphene devices. In this work, we show that at charge neutrality, long-range disorder increases the conductivity of the system above the intrinsic value. The enhancement arises from the predominantly vortical hydrodynamic flow caused by local deviations from charge neutrality. Its magnitude is inversely proportional to the shear viscosity of the electron liquid and scales as the square of the disorder correlation radius. This is qualitatively different from the situation away from charge neutrality. In that case, the flow is predominantly potential and produces negative viscous contributions to the conductivity, which are proportional to the sum of shear and bulk viscosities and inversely proportional to the square of disorder correlation radius.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1973

Aerodynamic influence coefficient method using singularity splines.

A new numerical formulation with computed results, is presented. This formulation combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of the loading function methods. The formulation employs a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfies all of the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise (termed 'spline'). Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral.

Mercer, J. E.↗
NASA NTRS1974

Aerodynamic influence coefficient method using singularity splines

A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.

Mercer, J. E.↗
NASA NTRS1979

Close-coupling calculations with an INDOX/1s static potential, semiclassical exchange, and a semi-empirical polarisation potential for electron-CO2 elastic scattering and rotational excitation

A semiempirical molecular-orbital method for modelling the effective potential for electron-molecule scattering is applied to elastic scattering and rotational excitation of CO2 at 20 eV impact energy. Agreement with experiment is reasonably good. The calculated rotationally summed integral cross section is 67.8 a sub 0 squared.

Onda, K.↗
NASA NTRS2011

Perimetric Complexity of Binary Digital Images: Notes on Calculation and Relation to Visual Complexity

Perimetric complexity is a measure of the complexity of binary pictures. It is defined as the sum of inside and outside perimeters of the foreground, squared, divided by the foreground area, divided by 4p . Difficulties arise when this definition is applied to digital images composed of binary pixels. In this paper we identify these problems and propose solutions. Perimetric complexity is often used as a measure of visual complexity, in which case it should take into account the limited resolution of the visual system. We propose a measure of visual perimetric complexity that meets this requirement.

Watson, Andrew B.↗
NASA NTRS1974

A redetermination of the galactic H I half-thickness and a discussion of some dynamical consequences

We have redetermined the half-thickness of the neutral hydrogen layer of our Galaxy over a significant portion of the disk (R, distance from the galactic center between 1.4 and 16.8 kpc). Four high-resolution 21-cm surveys were employed, and corrections were made for the effects of finite optical depth in the 21-cm line. We find that the half-thickness is approximately constant for R from 4.5 to 10 kpc and averages 260 pc for 0 to 90 deg long and 230 pc for 270 to 360 deg long. Using this thickness in conjunction with a recent galactic mass model, we are able to determine the quantity Q over a large portion of the galactic disk, where Q-squared is the ratio at the galactic plane of the sum of the gas, magnetic, and cosmic-ray pressures to the gas density.

Jackson, P. D.↗
NASA NTRS1995

Ulysses observations of latitude gradients in the heliospheric magnetic field

Several parameters measured by Ulysses as it traveled southward to heliographic latitudes of -50 deg are presented and analyzed. The radial component of the magnetic field, averaged over 5 deg latitude increments and extrapolated back to 1 AU, is found to agree with baseline measurements provided by IMP-8. There is little, if any, evidence of a latitude gradient, a result consistent with the dominance of the magnetic field associated with the heliospheric current sheet and with recent models which include the effect of the current sheet as well as of source surface fields. Thus far, the spiral angle agrees with the Parker spiral assuming a rate of rotation of the field lines at the Sun equal to the equatorial value. No evidence is seen of either a change in rotation rate with latitude or an unwinding of the spiral as suggested by a recent analysis. Hourly variances in the field magnitude and in the sum of the variances in the components, normalized to the square of the observed field strenght, show the former to be independent of latitude while the latter shows a strong increase with latitude. These two observations are shown to be associated with Alfven waves that are continuously present at high latitudes. The waves have large amplitudes, extend to long periods, and have important implications for galactic cosmic rays and the solar wind.

Smith, E. J.↗
DOE OSTI2024

Stochastic average model methods

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

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Effect of contrast on the active control of a moving line

In many passive visual tasks, human perceptual judgments are contrast dependent. To explore whether these contrast dependencies of visual perception also affect closed-loop manual control tasks, we examined visuomotor performance as humans actively controlled a moving luminance-defined line over a range of contrasts. Four subjects were asked to use a joystick to keep a horizontal line centered on a display as its vertical position was perturbed by a sum of sinusoids under two control regimes. The total root mean square (RMS) position error decreased quasi-linearly with increasing log contrast across the tested range (mean slope across subjects: -8.0 and -7.7% per log2 contrast, for the two control regimes, respectively). Frequency-response (Bode) plots showed a systematic increase in open-loop gain (mean slope: 1.44 and 1.30 dB per log2 contrast, respectively), and decrease in phase lag with increasing contrast, which can be accounted for by a decrease in response time delay (mean slope: 32 and 40 ms per log2 contrast, respectively). The performance data are well fit by a Crossover Model proposed by McRuer and Krendel, which allowed us to identify both visual position and motion cues driving performance. This analysis revealed that the position and motion cues used to support manual control under both control regimes appear equally sensitive to changes in stimulus contrast. In conclusion, our data show that active control of a moving visual stimulus is as dependent on contrast as passive perception and suggest that this effect is attributed to a shared contrast sensitivity early in the visual pathway, before any specialization for motion processing.

Contrast Sensitivity/physiology↗
NASA NTRS2016

Photosynthetic Efficiency of Northern Forest Ecosystems Using a MODIS-Derived Photochemical Reflectance Index (PRI)

This study evaluates a direct remote sensing approach from space for the determination of ecosystem photosynthetic light use efficiency (LUE), through measurement of vegetation reflectance changes expressed with the Photochemical Reflectance Index (PRI). The PRI is a normalized difference index based on spectral changes at a physiologically active wavelength (approximately 531 nanometers) as compared to a reference waveband, and is only available from a very few satellites. These include the two Moderate-Resolution Imaging Spectroradiometers (MODIS) on the Aqua and Terra satellites each of which have a narrow (10-nanometer) ocean band centered at 531 nanometers. We examined several PRI variations computed with candidate reference bands, since MODIS lacks the traditional 570-nanometer reference band. The PRI computed using MODIS land band 1 (620-670 nanometers) gave the best performance for daily LUE estimation. Through rigorous statistical analyses over a large image collection (n equals 420), the success of relating in situ daily tower-derived LUE to MODIS observations for northern forests was strongly influenced by satellite viewing geometry. LUE was calculated from CO2 fluxes (moles per moles of carbon absorbed quanta) measured at instrumented Canadian Carbon Program flux towers in four Canadian forests: a mature fir site in British Columbia, mature aspen and black spruce sites in Saskatchewan, and a mixed deciduous/coniferous forest site in Ontario. All aspects of the viewing geometry had significant effects on the MODIS-PRI, including the view zenith angle (VZA), the view azimuth angle, and the displacement of the view azimuth relative to the solar principal plane, in addition to illumination related variables.Nevertheless, we show that forward scatter sector views (VZA, 16 degrees-45 degrees) provided the strongest relationships to daily LUE, especially those collected in the early afternoon by Aqua (r squared = 0.83, RMSE (root mean square error) equals 0.003 moles per moles of carbon absorbed quanta). Nadir (VZA, 0 degrees plus or minus 15 degrees) and backscatter views (VZA, −16 degrees to −45 degrees) had lower performance in estimating LUE (nadir: r squared approximately equal to 0.62-0.67; backscatter: r squared approximately equal to 0.54-0.59) and similar estimation error (RMSE equals 0.004-0.005).When directional effects were not considered, only a moderately successful MODIS-PRI vs. LUE relationship (r squared equals 0.34, RMSE equals 0.007) was obtained in the full dataset (all views & sites, both satellites), but site-specific relationships were able to discriminate between coniferous and deciduous forests. Overall, MODIS-PRI values from Terra (late morning) were higher than those from Aqua (early afternoon), before/after the onset of diurnal stress responses expressed spectrally. Therefore, we identified ninety-two Terra-Aqua "same day" pairs, for which the sum of Terra morning and Aqua afternoon MODIS-PRI values (PRI (sub sum) using all available directional observations was linearly correlated with daily tower LUE (r squared equals 0.622, RMSE equals 0.013) and independent of site differences or meteorological information. Our study highlights the value of off-nadir directional reflectance observations, and the value of pairing morning and afternoon satellite observations to monitor stress responses that inhibit carbon uptake in Canadian forest ecosystems. In addition, we show that MODIS-PRI values, when derived from either: (i) forward views only, or (ii) Terra/Aqua same day (any view) combined observations, provided more accurate estimates of tower-measured daily LUE than those derived from either nadir or backscatter views or those calculated by the widely used semi-operational MODIS GPP model (MOD17) which is based on a theoretical maximum LUE and environmental data. Consequently, we demonstrate the importance of diurnal as well as off-nadir satellite observations for detecting vegetation physiological processes.

Middleton, E. M.↗
DOE OSTI2023

Calculation of Dynamical Response Functions Using a Bound-State Method

Abstract We investigate a method to extract response functions (dynamical polarisabilities) directly from a bound-state approach applied to calculations of perturbation-induced reactions. The use of a square-integrable basis leads to a response in the form of a sum of $$\delta $$ δ functions. We integrate this over energy and fit a smooth function to the resulting stepwise-continuous one. Its derivative gives the final approximation to the physical response function. We show that the method reproduces analytical results where known, and analyse the details for a variety of models. We apply it to some simple models, using the stochastic variational method as the numerical method. Albeit we find that this approach, and other numerical techniques, have some difficulties with the threshold behavior in coupled-channel problems with multiple thresholds, its stochastic nature allows us to extract robust results even for such cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1976

A quasi-Newton approach to optimization problems with probability density constraints

A quasi-Newton method is presented for minimizing a nonlinear function while constraining the variables to be nonnegative and sum to one. The nonnegativity constraints were eliminated by working with the squares of the variables and the resulting problem was solved using Tapia's general theory of quasi-Newton methods for constrained optimization. A user's guide for a computer program implementing this algorithm is provided.

Tapia, R. A.↗
NASA NTRS1986

Comparison of the STI NIPIP tracking dynamics identification with the on-line Fourier analyzer DFA results including a time varying case

The Non-Intrusive Pilot Identification Procedure (NIPIP) recently developed at STI and described at the 1981 Annual Manual was used to identify operators who were compensatory tracking a sub-critical-instability task. NIPIP uses a time domain least squares procedure converting to frequency domain coefficients. The forcing function was a sum of sinusoids supplied by the STI Mark II Describing Function Analyzer, which computes on-line Fourier coefficients of the operator's error/input describing function. The resulting open-loop and operator dynamics computed by each procedure are compared, and they are shown to be reasonably close when there is reasonable power in the error signal at the measurement frequencies. A special run was made in which the operator abruptly reduced gain within 1 sec, and the ability of the NIPIP to identify this step time variation in the operator is illustrated.

Jex, H. R.↗
DOE OSTI2024

Online and Offline Identification of False Data Injection Attacks in Battery Sensors Using a Single Particle Model

The cells in battery energy storage systems are monitored, protected, and controlled by battery management systems whose sensors are susceptible to cyberattacks. False data injection attacks (FDIAs) targeting batteries’ voltage sensors affect cell protection functions and the estimation of critical battery states like the state of charge (SoC). Inaccurate SoC estimation could result in battery overcharging and over discharging, which can have disastrous consequences on grid operations. This paper proposes a three-pronged online and offline method to detect, identify, and classify FDIAs corrupting the voltage sensors of a battery stack. To accurately model the dynamics of the series-connected cells a single particle model is used and to estimate the SoC, the unscented Kalman filter is employed. FDIA detection, identification, and classification was accomplished using a tuned cumulative sum (CUSUM) algorithm, which was compared with a baseline method, the chi-squared error detector. Online simulations and offline batch simulations were performed to determine the effectiveness of the proposed approach. Throughout the batch simulations, the CUSUM algorithm detected attacks, with no false positives, in 99.83% of cases, identified the corrupted sensor in 97% of cases, and determined if the attack was positively or negatively biased in 97% of cases.

25 ENERGY STORAGE↗