Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “least-squares analysis”

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 73 records · Page 4

A machine learning framework for accurate and robust analysis of radiation detector pulses

The microscopic properties of atomic nuclei are used to study various scientific questions. They are essential for understanding the fundamental forces of nature and the chemical evolution of the universe. Detecting decay radiation from radioactive nuclei makes it possible to probe these fundamental nuclear properties. Detector waveform traces may contain additional information about the radiation. Generally, advanced signal processing techniques are needed to extract this additional information, often involving fitting the waveform with model response functions using non-linear least-squares optimization with second-order gradient methods. While this is a powerful technique, it is also computationally expensive, leading to slow processing time, which scales with the volume of data. To address this problem, we have developed a machine learning (ML) approach that infers the characteristics of traces from a model detector response function. In particular, we are interested in classifying whether a single recorded trace consists of one or two pulse constituents and estimating the pulse parameters. Furthermore, our proposed ML method can precisely extract the pulses’ parameters, such as energy and timing information, and accurately classify the pulse multiplicity of a trace. Unlike non-learning-based approaches, our ML approach uses neural networks that are significantly faster at inference, as they do not require any optimization during this stage.

Curve fitting↗

Measurement-informed Dynamic Aggregation of Distribution Systems

This paper proposes a measurement-informed dynamic aggregation methodology in order to create equivalent representations of distribution systems that are compatible with large-scale transmission analysis. By optimizing an equivalent feeder parameters using time-series measurements of active power, reactive power, and voltage at the Point of Interconnection (POI), the approach yields simplified yet dynamically accurate equivalents. Implemented in PSCAD with models of photovoltaic–battery systems, three-phase motors, and static loads, the method employs hybrid differential evolution and bounded least-squares optimization laying the foundation for for real-time state estimation and optimized sensor placement in distribution networks.

Ahmed, Kazi Ishrak [University of Tennessee, Knoxv↗

Analysis of 12-700 micron reflectance spectra of three optical black samples

Normal incidence, specular reflectance spectra of three samples of optical black coatings have been measured in 10 filter passbands at effective wavelengths between 12 and 700 microns. Strong absorption, due to amorphous silicate material in the coatings, is found in the 15-50 micron region of each spectrum. Reflecting-layer model calculations are fitted by a nonlinear least-squares routine to the spectra at wavelengths beyond the absorption band to determine far-infrared optical parameters. It is found that, with a similar wavelength dependence, the absorption coefficient of a new black coating (ECP 2200) made by the Minnesota Mining and Manufacturing Company (3M) is less than 1/10 that of the discontinued coating, 3M Black Velvet Nextel. It is also shown that the extinction effect of a very rough substrate is negligible compared with the effects of increased upper surface roughness and coating thickness attendant upon fully covering a very rough substrate.

Smith, S. M.↗

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗

Analysis of estimation algorithms for autonomous navigation with TDRSS data

Estimation techniques for onboard orbit determination using Tracking and Data Relay Satellite System (TDRSS) data are investigated. The two user satellite orbits studied are similar to Landsat-D (near-circular, 700 kilometers altitude, near-polar inclination). The following estimation algorithms are identified as candidates for use in autonomous navigation: (1) the extended Kalman filter with process noise, (2) the EKF with consider parameters, (3) the sequential Kalman filter with consider parameters, and (4) the batch least-squares differential correction technique. The candidate estimators are evaluated with respect to their performance with both baseline and worst TDRSS measurement errors and tracking configurations.

Dunham, J.↗

Photo-orthotropic elasticity - A new technique for stress analysis of composites

This paper deals with the application of photoelastic methods to glass fiber reinforced plastics. The basic principles governing the photoelastic response in birefringent orthotropic composites are briefly reviewed. Then, three important aspects of photo-orthotropic elasticity, namely, photoelastic calibration, separation of principal stresses, and application to problems of engineering importance, are considered in detail. A new calibration specimen is proposed for determining the three independent photoelastic constants: an orthotropic half-plane subjected to a concentrated edge load. From the isochromatic response, the three stress-fringe values can be obtained by a least-squares procedure based on the theoretical stress distribution. Next, four methods of separating the principal stresses in birefringent composite models are proposed: the method of drilling small holes at the points of interest, the method of oblique incidence, the method of combining the transmitted photoelastic response with the reflected photoelastic response, and the method of birefringence dispersion.

Prabhakaran, R.↗

The effect of data analysis techniques on the interpretation of wide-angle longwave radiation measurements

Three different data analysis techniques - shape factor, parameter estimation, and deconvolution - have been applied to the same set of satellite radiation measurements to determine their effect on the estimated radiation field. The measurements are from a wide-angle, horizon-to-horizon, nadir-pointing sensor. The shape factor technique reduces each measurement to a radiant exitance at the top of the atmosphere by simple division by a constant. The parameter estimation technique processes all measurements together as a batch and defines the radiant exitance as a least-squares fit to the data. The deconvolution technique takes advantage of the fact that spherical harmonics are the eigenfunctions of the measurement operator. All three techniques are derived, and their assumptions, advantages and disadvantages are discussed. Their results are compared globally, zonally, regionally and on a spatial spectrum basis. All three techniques give comparable results for global parameters; however, results on a regional scale were quite different. The standard deviations of the regional differences in radiant exitance varied from 7.4 to 13.5 W/sq m. Of the three techniques, the parameter estimation technique produced the best regional results and is the choice of the author.

Green, R. N.↗

Simultaneous elastic and photoelastic calibration of birefringent orthotropic model materials

Transmission photoelastic analysis of composite models has been shown in recent years to be a viable stress analysis method for composite structures. However, difficulties in the fabrication of the model materials require the development of efficient mechanical and optical calibration procedures. In this paper, the methods for measuring elastic and photoelastic constants of composites are briefly reviewed. Then a new method, utilizing a relatively small calibration specimen, is described. In this method, electrical resistance strain gauges are mounted at selected points of a half-plane model and photoelastic measurements and strain readings are simultaneously taken. The data are processed by a least-squares procedure to yield the elastic and photoelastic constants. Comparison of results obtained by the proposed method with the results from conventional tests shows good agreement.

Prabhakaran, R.↗

Seasat orbit refinement for altimetry application

This paper describes the use of stochastic differential correction models in refining the Seasat orbit based on post-flight analysis of tracking data. The objective is to obtain orbital-height precision that is commensurate with the inherent Seasat altimetry data precision level of 10 cms. Local corrections to a mean ballistic arc, perturbed principally by atmospheric drag variations and local gravitational anomalies, are obtained by the introduction of stochastic dynamical models in conjunction with optimal estimation/smoothing techniques. Assessment of the resulting orbit with 'ground truth' provided by Seasat altimetry data shows that the orbital height precision is improved by 32% when compared to a conventional least-squares solution using the same data set. The orbital height precision realized by employing stochastic differential correction models is in the range of 73 cms to 208 cms rms.

Mohan, S. N.↗

Correction factors for on-line microprobe analysis of multielement alloy systems

An on-line correction technique was developed for the conversion of electron probe X-ray intensities into concentrations of emitting elements. This technique consisted of off-line calculation and representation of binary interaction data which were read into an on-line minicomputer to calculate variable correction coefficients. These coefficients were used to correct the X-ray data without significantly increasing computer core requirements. The binary interaction data were obtained by running Colby's MAGIC 4 program in the reverse mode. The data for each binary interaction were represented by polynomial coefficients obtained by least-squares fitting a third-order polynomial. Polynomial coefficients were generated for most of the common binary interactions at different accelerating potentials and are included. Results are presented for the analyses of several alloy standards to demonstrate the applicability of this correction procedure.

Unnam, J.↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

An analysis of the shortening and lengthening of the phase of satellite-to-satellite (SST) data that passed within 40-700 km above the earth surface during its ATS-6 to GEOS-3 to ATS-6 path resulted in refractivity vs height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbations from the moon, sun, and a 15th-order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occulation period and a model ionosphere was estimated by a least-squares adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the two methods were caused by orbital errors, which propagated into the solution.

Liu, A. S.↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

When the radio link between two satellites (GEOS-3/ATS-6) is intercepted by the earth's ionosphere and neutral atmosphere, a change in the Doppler frequency results. Travel through the atmosphere causes the Doppler phase to be advanced in the ionosphere's portion and retarded in the neutral portion of the atmosphere. Analysis of the shortening and lengthening of the phase of the Satellite-to-Satellite Tracking (SST) data that passed within 40-700 km above the earth's surface during its ATS-6 to GEOS-3 to ATS-6 path, caused by the atmosphere, results in refractivity versus height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbation from the Moon, Sun and a 15th order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occultation period and a model ionosphere was estimated by a least-square adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the 2 methods were caused by orbital errors, which propagated into the solution. The SST data yielded refractive index profiles in a novel economical manner because no additional or special on-board equipment were required.

Liu, A. S.↗

Center of mass estimation in closed vortices - A verification in principle and practice

The problem of tracking closed mesoscale vortices using center of mass estimation techniques is studied. Three estimators are evaluated using data from a warm core Gulf Stream ring. The comparisons show that a method based on the intersection of perpendicular bisectors and one using a least-squares fit of a conic section perform comparably. The perpendicular bisector algorithm is used in conjunction with a Gaussian ring model and a star-shaped survey pattern to produce an expected error curve as a function of vortex translation, survey speed and vortex size. For typical ring parameters, center estimation is usually possible to within + or - 5 km. The feasibility of using differing data sets to construct a history of ring motion based on a coordinate system moving with the ring is also investigated. In this way, the validity of using satellite-derived data and drifter trajectories to estimate the center of mass of a mesoscale feature is assessed. The results of the analysis demonstrate that the location of the deeper structure of the ring and the surface expression are sufficiently well correlated to permit dynamically relevant calculations based on surface measurements. It is shown that satellite-derived data can be used to approximate the center of mass trajectory to within the error in the individual center estimates for the period analyzed. The Lagrangian-drifter-derived centers are offset from the center of mass trajectory in a manner consistent with kinematic arguments.

Hooker, S. B.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗

Rapid monitoring of fermentations: a feasibility study on biological 2,3-butanediol production

2,3-butanediol (2,3-BDO) is an economically important platform chemical that can be produced by the fermentation of sugars using an engineered strain of Zymomonas mobilis . These fermentations require continuous monitoring and modification of fermentation conditions to maximize 2,3-BDO yields and minimize the production of the undesired coproducts glycerol and acetoin. Because of the time required for sampling and off-line chromatographic measurement of fermentation samples, the ability of fermentation scientists to modify fermentation conditions in a timely manner is limited. The goal of this study was to test if near-infrared spectroscopy (NIRS) along with multivariate statistics could reduce the time needed for this analysis and enable real-time monitoring and control of the fermentation. In this work we developed partial least squares (PLS) calibration models to predict the concentrations of glucose, xylose, 2,3-BDO, acetoin, and glycerol in fermentations via NIRS using two different spectrometers and two different spectroscopy modalities. We first evaluated the feasibility of rapid NIRS monitoring through experiments where we measured the signals from each analyte of interest and built NIRS-based PLS models using spectra from synthetic samples containing uncorrelated concentrations of these analytes. All analytes showed unique spectral signatures, and this initial modeling showed that all analytes could be detected simultaneously. We then began work with samples from laboratory fermentation experiments and tested the feasibility of regression model development across two spectral collection modalities (at-line and on-line) and two instruments: a laboratory-grade instrument and a low-cost instrument with a more limited spectral range. All modalities showed promise in the ability to monitor Z. mobilis fermentations of glucose and xylose to 2,3-BDO. The low-cost instrument displayed a lower signal-to-noise ratio than the laboratory-grade instrument, which led to comparatively lower performance overall, but still provided sufficient accuracy to monitor fermentation trends. While the ease of use of on-line monitoring systems was favored as compared to at-line systems due to the lack of sampling required and potential for automated process control, we observed some decrease in performance due to the additional complexity of the sample matrix. We have demonstrated that NIRS combined with multivariate analysis can be used for at-line and on-line monitoring of the concentrations of glucose, xylose, 2,3-BDO, acetoin, and glycerol during Z. mobilis fermentations. The decrease in signal-to-noise ratio when using a low-cost spectrometer led to greater prediction error than the laboratory-grade spectrometer for at-line monitoring. The on-line monitoring modality showed great promise for real time process control via NIRS.

09 BIOMASS FUELS↗

Tracking seasonal variability in plant traits from spaceborne PRISMA and NEON AOP across forest types and ecoregions

Plant traits serve as critical indicators of how plants adapt to environmental changes and influence ecosystem functions. While airborne hyperspectral remote sensing effectively maps plant traits through detailed reflectance properties, it is limited by cost and scale, making large-scale and temporal studies challenging. The recently launched spaceborne hyperspectral imager, PRecursore IperSpettrale della Missione Applicativa (PRISMA), offers frequent, large scale and high-fidelity observations on a spatial resolution of 30 m and a revisit time of around 29 days, making it suitable for large-scale seasonal trait mapping. However, their potential remains largely unexplored. This study developed a multi-stage framework by leveraging the PRISMA spaceborne hyperspectral data and National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) hyperspectral data to investigate the seasonal dynamics of four key plant traits — chlorophyll content, carotenoid content, equivalent water thickness, and nitrogen content — across eleven NEON sites representing diverse forest types and ecoregions in the contiguous U.S. Our results demonstrated that PRISMA hyperspectral data can reliably track seasonal variability in plant traits, achieving overall R 2 values ranging from 0.78 to 0.88 and normalized root mean square error (NRMSE) values ranging from 5.4% to 8.4% for the four traits. Seasonal patterns revealed bell-shaped trajectories for chlorophyll and carotenoids, while equivalent water thickness decreased steadily across most sites, driven by structural changes during leaf maturation and senescence. Nitrogen content exhibited less pronounced seasonal variation but followed expected nutrient resorption patterns. Analysis of environmental drivers showed that seasonal variability is primarily controlled by solar radiation and day length in northern sites, vapor pressure in semi-arid regions, and temperature in mid-southeastern sites. Spatial variability, meanwhile, was primarily driven by soil properties, particularly during the peak growing season. However, the influence of soil variables slightly declines toward the end of the season at several sites, as climatic factors become more prominent. This study highlights the capability of PRISMA, and potentially other similar spaceborne hyperspectral data for large-scale, time-series plant trait mapping and provides valuable insights into the interactions between plant traits and environmental factors. In conclusion, these findings contribute to advancing our understanding of plant functional ecology and improving predictions of ecosystem responses to environmental changes.

Environmental drivers↗