Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function Estimation”

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 829 records · Page 46

Assimilation of SMAP and ASCAT Soil Moisture Retrievals into the JULES Land Surface Model Using the Local Ensemble Transform Kalman Filter

A land data assimilation system is developed to merge satellite soil moisture retrievals into the Joint U.K. Land Environment Simulator (JULES) land surface model (LSM) using the Local Ensemble Transform Kalman Filter (LETKF). The system assimilates microwave soil moisture retrievals from the Soil Moisture Active Passive (SMAP) radiometer and the Advanced Scatterometer (ASCAT) after bias correction based on cumulative distribution function fitting. The soil moisture assimilation estimates are evaluated with ground-based soil moisture measurements over the continental U.S. for five consecutive warm seasons (May–September of 2015–2019). The result shows that both SMAP and ASCAT retrievals improve the accuracy of soil moisture estimates. Especially, the SMAP single-sensor assimilation experiment shows the best performance with the increase of temporal anomaly correlation by ΔR ~ 0.05 for surface soil moisture and ΔR ~ 0.03 for root-zone soil moisture compared with the LSM simulation without satellite data assimilation. SMAP assimilation is more skillful than ASCAT assimilation primarily because of the greater skill of the assimilated SMAP retrievals compared to the ASCAT retrievals. The skill improvement also depends significantly on the region; the higher skill improvement in the western U.S. compared to the eastern U.S. is explained by the Kalman gain in the two experiments. Additionally, the regional skill differences in the single-sensor assimilation experiments are attributed to the number of assimilated observations. Finally, the soil moisture assimilation estimates provide more realistic land surface information than model-only simulations for the 2015 and the 2016 western U.S. droughts, suggesting the advantage of using satellite soil moisture retrievals in the current drought monitoring system.

SMAP↗

First Attempt of Global-scale Assimilation of Subdaily Scale Soil Moisture Estimates from CYGNSS and SMAP into a Land Surface Model

Soil moisture performs a key function in the hydrologic process and understanding the global-scale water cycle. However, estimations of soil moisture taken from current sun-synchronous orbit (SSO) satellites are limited in that they are neither spatially nor temporally continuous. This limitation creates discontinuous soil moisture observation from space and hampers our understanding of the fundamental processes that control the surface hydrologic cycle across both time and space domains. Here, we propose to use frequent soil moisture observations from NASA’s constellation of eight micro-satellites called the Cyclone Global Navigation Satellite System(CYGNSS) together with the Soil Moisture Active Passive (SMAP) to assimilate subdaily-scale soil moisture intoa land surface model(LSM). Our results, which are based on triple collocation analysis(TCA), show how current scientific advances in satellite systems can fill previous gaps in soil moisture observations in subdaily scale bypast observations, and eventually adds value to improvements in global scale soil moisture estimates in LSMs. Overall, TCA-based fractional mean square errors (fMSE) of LSM soil moisture are improved by 61% with the synergetic assimilation of CYGNSS data with SMAP soil moisture observations. However, assimilating satellite-based soil moisture over dense vegetation areas can degrade the performance of LSMs as these areas propagate erroneous soil moisture information to LSMs. To our knowledge, this study isthe first global assimilation of GNSS-based soil moisture observations in land surface models.

Hyunglok Kim↗

Monte Carlo Tree Search Methods for the Earth-Observing Satellite Scheduling Problem

This work explores on-board planning for the single spacecraft, multiple ground station Earth-observing satellite scheduling problem through artificial neural network function approximation of state–action value estimates generated by Monte Carlo tree search (MCTS). An extensive hyperparameter search is conducted for MCTS on the basis of performance, safety, and downlink opportunity utilization to determine the best hyperparameter combination for data generation. A hyperparameter search is also conducted on neural network architectures. The learned behavior of each network is explored, and each network architecture’s robustness to orbits and epochs outside of the training distributions is investigated. Furthermore, each algorithm is compared with a genetic algorithm, which serves to provide a baseline for optimality. MCTS is shown to compute near-optimal solutions in comparison to the genetic algorithm. The state–action value networks are shown to match or exceed the performance of MCTS in six orders of magnitude less execution time, showing promise for execution on board spacecraft.

Adam P. Herrmann↗

Modeling, estimation and identification methods for static shape determination of flexible structures

This paper outlines methods for modeling, identification and estimation for static determination of flexible structures. The shape estimation schemes are based on structural models specified by (possibly interconnected) elliptic partial differential equations. The identification techniques provide approximate knowledge of parameters in elliptic systems. The techniques are based on the method of maximum-likelihood that finds parameter values such that the likelihood functional associated with the system model is maximized. The estimation methods are obtained by means of a function-space approach that seeks to obtain the conditional mean of the state given the data and a white noise characterization of model errors. The solutions are obtained in a batch-processing mode in which all the data is processed simultaneously. After methods for computing the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the related estimation error is conducted. In addition to outlining the above theoretical results, the paper presents typical flexible structure simulations illustrating performance of the shape determination methods.

Rodriguez, G.↗

Parameter Estimation of Actuators for Benchmark Active Control Technology (BACT) Wind Tunnel Model with Analysis of Wear and Aerodynamic Loading Effects

This report describes the development of transfer function models for the trailing-edge and upper and lower spoiler actuators of the Benchmark Active Control Technology (BACT) wind tunnel model for application to control system analysis and design. A simple nonlinear least-squares parameter estimation approach is applied to determine transfer function parameters from frequency response data. Unconstrained quasi-Newton minimization of weighted frequency response error was employed to estimate the transfer function parameters. An analysis of the behavior of the actuators over time to assess the effects of wear and aerodynamic load by using the transfer function models is also presented. The frequency responses indicate consistent actuator behavior throughout the wind tunnel test and only slight degradation in effectiveness due to aerodynamic hinge loading. The resulting actuator models have been used in design, analysis, and simulation of controllers for the BACT to successfully suppress flutter over a wide range of conditions.

Waszak, Martin R.↗

Estimating phase integrals - A generalization of Russell's law

An attempt is made to demonstrate, by means of Hapke's (1981) photometric function, that a simple and reliable method exists for the estimation of phase integrals from limited higher phase angle measurements. This method is a generalization of the approximation first proposed by Russell (1916) and more recently treated by Veverka (1971). It is shown that this generalization of Russell's law can employ observations anywhere in the range of phase angles from 40 to 90 deg; optimum estimates of q will be obtained if the data near 70 deg are obtainable.

Verbiscer, Anne J.↗

An optical study of the faint end of the stellar luminosity function

We implement a new method by which to study the faint end of the field star luminosity function. The method relies on deep, multicolor photometry of fields projected against highly obscured, nearby molecular clouds. The clouds act as nearby opaque screens and delimit a well-defined survey volume which is in principle free of the problem of distinguishing nearby, intrinsically faint dwarf stars from more distant red giants. This study is based upon deep photographic and CCD photometry at optical (V, R, I) bandpasses toward the most highly obscured portions of the Taurus and Ophiuchus molecular clouds. The total volume delimited by the clouds is approximately 200 cu pc. Within this region our survey is complete for all stars brighter than M(sub V) = 16-17 mag; at R and I, the survey is complete down to the lowest mass stars capable of sustaining core hydrogen burning. We estimate the faint end of the field star luminosity function for the composite Taurus and Ophiuchus foreground sample and find that it resembles the local luminosity function down to M(sub V) approx. 16. At still fainter magnitudes we find more stars than do photometric parallex studies of the polar regions. This difference widens dramatically if even the simplest correction for incompleteness is applied to our data. We therefore tentatively conclude that the luminosity function rises beyond M(sub V) approx. 16; even if we discard our attempts to correct for incompleteness in the faintest magnitude bins, the luminosity function at least remains flat for the lowest mass stars. Our provisional finding that the luminosity function rises beyond its well-known peak at M(sub V) approx. 12-13, implies that the initial mass function (IMF) probably rises beyond the turnover point associated with this peak. Even if our most conservative estimate for the faint end of the luminosity function is used-in which no corrections are made for incompleteness-the IMF must at least remain flat down to the edge of the hydrogen-burning main sequence.

Jarrett, T. H.↗

Adjoint-Based Mesh Adaptation for the Sonic Boom Signature Loudness

The mesh adaptation functionality of FUN3D is utilized to obtain a mesh optimized to calculate sonic boom ground signature loudness. During this process, the coupling between the discrete-adjoints of the computational fluid dynamics tool FUN3D and the atmospheric propagation tool sBOOM is exploited to form the error estimate. This new mesh adaptation methodology will allow generation of suitable meshes adapted to reduce the estimated errors in the ground loudness, which is an optimization metric employed in supersonic aircraft design. This new output-based adaptation could allow new insights into meshing for sonic boom analysis and design, and complements existing output-based adaptation techniques such as adaptation to reduce estimated errors in off-body pressure functional. This effort could also have implications for other coupled multidisciplinary adjoint capabilities (e.g., aeroelasticity) as well as inclusion of propagation specific parameters such as prevailing winds or non-standard atmospheric conditions. Results are discussed in the context of existing methods and appropriate conclusions are drawn as to the efficacy and efficiency of the developed capability.

Rallabhandi, Sriram K.↗

Computer program reduces calculation time of normal response functions

FORTRAN 2 computer program rapidly calculates parameters of maximum likelihood estimates from sensitivity experiment data populations. The program uses the Newton-Raphson iterative procedure to calculate the mean and standard deviation of portions of the cumulative normal response function.

Alexander, M. J.↗

Characterization, parameter estimation, and aircraft response statistics of atmospheric turbulence

A nonGaussian three component model of atmospheric turbulence is postulated that accounts for readily observable features of turbulence velocity records, their autocorrelation functions, and their spectra. Methods for computing probability density functions and mean exceedance rates of a generic aircraft response variable are developed using nonGaussian turbulence characterizations readily extracted from velocity recordings. A maximum likelihood method is developed for optimal estimation of the integral scale and intensity of records possessing von Karman transverse of longitudinal spectra. Formulas for the variances of such parameter estimates are developed. The maximum likelihood and least-square approaches are combined to yield a method for estimating the autocorrelation function parameters of a two component model for turbulence.

Mark, W. D.↗

Temporal analysis of the October 1989 proton flare using computerized anatomical models

The GOES-7 time history data of hourly averaged integral proton fluxes at various particle kinetic energies are analyzed for the solar proton event that occurred between October 19 and 29, 1989. By analyzing the time history data, the dose rates which may vary over many orders of magnitude in the early phases of the flare can be estimated as well as the cumulative dose as a function of time. Basic transport calculations are coupled with detailed body organ thickness distributions from computerized anatomical models to estimate dose rates and cumulative doses to 20 critical body organs. For a 5-cm-thick water shield, cumulative skin, eye, and blood-forming-organ dose equivalents of 1.27, 1.23, and 0.41 Sv, respectively, are estimated. These results are approximately 40-50 percent less than the widely used 0- and 5-cm slab dose estimates. The risk of cancer incidence and mortality are also estimated for astronauts protected by various water shield thicknesses.

Simonsen, L. C.↗

An approach to the analysis of performance of quasi-optimum digital phase-locked loops.

An approach to the analysis of performance of quasi-optimum digital phase-locked loops (DPLL's) is presented. An expression for the characteristic function of the prior error in the state estimate is derived, and from this expression an infinite dimensional equation for the prior error variance is obtained. The prior error-variance equation is a function of the communication system model and the DPLL gain and is independent of the method used to derive the DPLL gain. Two approximations are discussed for reducing the prior error-variance equation to finite dimension. The effectiveness of one approximation in analyzing DPLL performance is studied.

Polk, D. R.↗

Blocking reduction of Landsat Thematic Mapper JPEG browse images using optimal PSNR estimated spectra adaptive postfiltering

Two representative sample images of Band 4 of the Landsat Thematic Mapper are compressed with the JPEG algorithm at 8:1, 16:1 and 24:1 Compression Ratios for experimental browsing purposes. We then apply the Optimal PSNR Estimated Spectra Adaptive Postfiltering (ESAP) algorithm to reduce the DCT blocking distortion. ESAP reduces the blocking distortion while preserving most of the image's edge information by adaptively postfiltering the decoded image using the block's spectral information already obtainable from each block's DCT coefficients. The algorithm iteratively applied a one dimensional log-sigmoid weighting function to the separable interpolated local block estimated spectra of the decoded image until it converges to the optimal PSNR with respect to the original using a 2-D steepest ascent search. Convergence is obtained in a few iterations for integer parameters. The optimal logsig parameters are transmitted to the decoder as a negligible byte of overhead data. A unique maxima is guaranteed due to the 2-D asymptotic exponential overshoot shape of the surface generated by the algorithm. ESAP is based on a DFT analysis of the DCT basis functions. It is implemented with pixel-by-pixel spatially adaptive separable FIR postfilters. PSNR objective improvements between 0.4 to 0.8 dB are shown together with their corresponding optimal PSNR adaptive postfiltered images.

Linares, Irving↗

An Ab Initio Study of CuCO

Modified coupled-pair functional (MCPF) calculations and coupled cluster singles and doubles calculations, which include a perturbational estimate of the connected triples [CCSD(T)], yield a bent structure for CuCO, thus, supporting the prediction of a nonlinear structure based on density functional (DF) calculations. Our best estimate for the binding energy is 4.9 +/- 1.4 kcal/mol; this is in better agreement with experiment (6.0 +/- 1.2 kcal/mol) than the DF approach which yields a value (19.6 kcal/mol) significantly larger than experiment.

Bauschlicher, Charles W., Jr.↗

Diagnostics of sources of disturbances and distribution of vibrations over the width of a tape in tape-feed mechanisms

Disturbances created by certain assemblies and components of tape-feed mechanisms (TFM) and acting on a moving magnetic tape are studied. The method, based on elements of digital logic, is established by stress-strain diagrams of the longitudinal deformations and vibrations across the width of a magnetic tape. Experimental studies were carried out for determination of the functional relationships of longitudinal deformations in a section of magnetic tape to the magnitude of roller play, tension vibrations rate of movement, and elasticity of magnetic tapes. A block diagram of the measurements is shown. Appropriate digital computer algorithms and programs were proposed for statistical analysis of the data obtained. Estimates of the mathematical expectation, dispersion, intercorrelation function, energy spectral density, and distribution pattern of the random process values were calculated.

Kenstavichyus, A. B. B.↗

Minimum silicon wafer thickness for ID wafering

An analytical model, based on fracture mechanics analysis, is proposed for estimating the minimum wafer thickness as a function of the diameter requirement for solar cells. The conditions under which the model can be applied are discussed with reference to the critical flaw size, the applied force, and the width of the side support. It is shown that the equivalent cantilever force applied during ID slicing can be estimated from the wafering mechanical yield data. The width of the wafer side support was found to be a significant factor in controlling the minimum allowable wafer thickness during slicing. Wafer side support width requirements were found to increase with decreasing wafer thickness.

Chen, C. P.↗

Attitude and Trajectory Estimation Using Earth Magnetic Field Data

The magnetometer has long been a reliable, inexpensive sensor used in spacecraft momentum management and attitude estimation. Recent studies show an increased accuracy potential for magnetometer-only attitude estimation systems. Since the Earth's magnetic field is a function of time and position, and since time is known quite precisely, the differences between the computer and measured magnetic field components, as measured by the magnetometers throughout the entire spacecraft orbit, are a function of both the spacecraft trajectory and attitude errors. Therefore, these errors can be used to estimate both trajectory and attitude. Traditionally, satellite attitude and trajectory have been estimated with completely separate system, using different measurement data. Recently, trajectory estimation for low earth orbit satellites was successfully demonstrated in ground software using only magnetometer data. This work proposes a single augmented extended Kalman Filter to simultaneously and autonomously estimate both spacecraft trajectory and attitude with data from a magnetometer and either dynamically determined rates or gyro-measured body rates.

Deutschmann, Julie↗

Atmospheric attenuation for correlated satellite communication ground sites

Link margin estimates are calculated for 30/20 GHz satellite communication systems employing closely-spaced (4 to 100 km) 'dual diversity' switched ground sites. The link margin estimates are based on a new analysis in which the bivariate rain attenuation density function for two correlated ground sites is modelled by an exponential density function. The results of the exponential density function analysis enable derivation of a direct relation between rain margin estimates and probability of exceedance (link availability). Margins typically in the range 2 to 12 dB are calculated for various ground site separations and summarized for seven city locations and five satellite orbit positions over the U.S. The results can be easily extended to other EHF satellite frequencies and other locations.

Christopher, P.↗