Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter 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 703 records · Page 39

Are light curve classification metrics good proxies for SN Ia cosmological constraining power?

Context. When selecting a light curve classifier for use as part of a photometric supernova Ia (SN Ia) cosmological analysis, it is common to make decisions based on metrics of classification performance, such as the contamination within the photometrically classified SN Ia sample, rather than a measure of cosmological constraining power. If the former is an appropriate proxy for the latter, this practice would eliminate the computational expense of a full cosmology forecast in the analysis pipeline design process. Aims. This study tests the assumption that light curve classification metrics are an appropriate proxy for cosmology metrics. Methods. We emulated photometric SN Ia cosmology light curve samples with controlled contamination rates of individual contaminant classes and evaluated each of them under a set of classification metrics. We then derived cosmological parameter constraints from all samples under two common analysis approaches and quantified the impact of contamination by each contaminant class on the resulting cosmological parameter estimates. Results. We observe that cosmology metrics are sensitive to both the contamination rate and the class of the contaminating population, whereas the classification metrics are shown to be insensitive to the latter. Conclusions. Based on these findings, we discourage any exclusive reliance on light curve classification-based metrics for analysis design decisions, which (counterintuitively) include but are not limited to the classifier choice. Instead, we recommend optimising science analysis pipeline design choices using a metric of the information gained about the physical parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

Deconvolution estimation theory applied to Nimbus 6 ERB data

It is pointed out that the ERB (Earth Radiation Budget) Experiment aboard the Nimbus 6 spacecraft has provided nearly 3 years of data thus far from its wide field of view (WFOV) radiometers. Each data point is an integral of the irradiance from all points within the field of view of the WFOV sensor, which is an approximately 60 deg diameter circular region on the earth. House (1972) proposed that the data, being a convolution of the flux field at the top of the atmosphere, could be convoluted so as to enhance the resolution. The problem was solved by Smith and Green (1975-76) for the case of earth emitted radiation. A parameter estimation approach to the deconvolution problem was formulated. A description is presented of the deconvolution estimation concept and the results obtained by its application to the Nimbus 6 ERB WFOV data for earth emitted radiation for August 1975.

Green, R. N.↗

Estimating tropical-forest density profiles from multibaseline interferometric SAR

Vertical profiles of forest density are potentially robust indicators of forest biomass, fire susceptibility and ecosystem function. Tropical forests, which are among the most dense and complicated targets for remote sensing, contain about 45% of the world's biomass. Remote sensing of tropical forest structure is therefore an important component to global biomass and carbon monitoring. This paper shows preliminary results of a multibasline interfereomtric SAR (InSAR) experiment over primary, secondary, and selectively logged forests at La Selva Biological Station in Costa Rica. The profile shown results from inverse Fourier transforming 8 of the 18 baselines acquired. A profile is shown compared to lidar and field measurements. Results are highly preliminary and for qualitative assessment only. Parameter estimation will eventually replace Fourier inversion as the means to producing profiles.

profiles↗

Input design for identification of aircraft stability and control derivatives

An approach for designing inputs to identify stability and control derivatives from flight test data is presented. This approach is based on finding inputs which provide the maximum possible accuracy of derivative estimates. Two techniques of input specification are implemented for this objective - a time domain technique and a frequency domain technique. The time domain technique gives the control input time history and can be used for any allowable duration of test maneuver, including those where data lengths can only be of short duration. The frequency domain technique specifies the input frequency spectrum, and is best applied for tests where extended data lengths, much longer than the time constants of the modes of interest, are possible. These technqiues are used to design inputs to identify parameters in longitudinal and lateral linear models of conventional aircraft. The constraints of aircraft response limits, such as on structural loads, are realized indirectly through a total energy constraint on the input. Tests with simulated data and theoretical predictions show that the new approaches give input signals which can provide more accurate parameter estimates than can conventional inputs of the same total energy. Results obtained indicate that the approach has been brought to the point where it should be used on flight tests for further evaluation.

Gupta, N. K.↗

On the VLBI-satellite laser ranging Iron Triangle intercomparison experiment

Simulations were performed to determine the suitability of the proposed station locations. The criterion was a comparison among the possible station configuration of the standard deviations of baseline and Earth rotation parameters estimated from a least squares covariance analysis. Only the relative magnitudes of the standard deviations were addressed in the analysis. Only random errors were assumed and no provision was made for systematic effects. The Iron Triangle, consisting of the stations at Westford (Haystack), Massachusetts, Owens Valley, California, and Ft. Davis, Texas, was regarded as the structure of the proposed network with options to incorporate either the Goldstone, California, or Green Bank, West Virginia, station or both. It was decided to include the Richmond, Florida, station in the analysis since it offered more North-South separation and therefore could strengthen the geometry of the network especially in the recovery of Earth rotation parameters.

Bock, Y.↗

An iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions, Addendum

New results and insights concerning a previously published iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions were discussed. It was shown that the procedure converges locally to the consistent maximum likelihood estimate as long as a specified parameter is bounded between two limits. Bound values were given to yield optimal local convergence.

Peters, B. C., Jr.↗

Inference of dust opacities for the 1977 Martian great dust storms from Viking Lander 1 pressure data

The tidal heating components for the dusty Martian atmosphere are computed based on dust optical parameters estimated from Viking Lander imaging data, and used to compute the variation of the tidal surface pressure components at the Viking Lander sites as a function of season and the total vertical extinction optical depth of the atmosphere. An atmospheric tidal model is used which is based on the inviscid, hydrostatic primitive equations linearized about a motionless basic state the temperature of which varies only with height, and the profiles of the tidal forcing components are computed using a delta-Eddington approximation to the radiative transfer equations. Comparison of the model results with the observed variations of surface pressure and overhead dust opacity at the Viking Lander 1 site reveal that the dust opacities and optical parameters derived from imaging data are roughly representative of the global dust haze necessary to reproduce the observed surface pressure amplitudes, with the exception of the model-inferred asymmetry parameter, which is smaller during the onset of a great storm. The observed preferential enhancement of the semidiurnal tide with respect to the diurnal tide during dust storm onset is shown to be due primarily to the elevation of the tidal heating source in a very dusty atmosphere.

Zurek, R. W.↗

Performance Analysis on Fault Tolerant Control System

In a fault tolerant control (FTC) system, a parameter varying FTC law is reconfigured based on fault parameters estimated by fault detection and isolation (FDI) modules. FDI modules require some time to detect fault occurrences in aero-vehicle dynamics. In this paper, an FTC analysis framework is provided to calculate the upper bound of an induced-L(sub 2) norm of an FTC system with existence of false identification and detection time delay. The upper bound is written as a function of a fault detection time and exponential decay rates and has been used to determine which FTC law produces less performance degradation (tracking error) due to false identification. The analysis framework is applied for an FTC system of a HiMAT (Highly Maneuverable Aircraft Technology) vehicle. Index Terms fault tolerant control system, linear parameter varying system, HiMAT vehicle.

Shin, Jong-Yeob↗

GVIPS Models and Software

Two reports discuss, respectively, (1) the generalized viscoplasticity with potential structure (GVIPS) class of mathematical models and (2) the Constitutive Material Parameter Estimator (COMPARE) computer program. GVIPS models are constructed within a thermodynamics- and potential-based theoretical framework, wherein one uses internal state variables and derives constitutive equations for both the reversible (elastic) and the irreversible (viscoplastic) behaviors of materials. Because of the underlying potential structure, GVIPS models not only capture a variety of material behaviors but also are very computationally efficient. COMPARE comprises (1) an analysis core and (2) a C++-language subprogram that implements a Windows-based graphical user interface (GUI) for controlling the core. The GUI relieves the user of the sometimes tedious task of preparing data for the analysis core, freeing the user to concentrate on the task of fitting experimental data and ultimately obtaining a set of material parameters. The analysis core consists of three modules: one for GVIPS material models, an analysis module containing a specialized finite-element solution algorithm, and an optimization module. COMPARE solves the problem of finding GVIPS material parameters in the manner of a design-optimization problem in which the parameters are the design variables.

Arnold, Steven M.↗

Controls for space structures

Assembly and operation of large space structures (LSS) in orbit will require robot-assisted docking and berthing of partially-assembled structures. These operations require new solutions to the problems of controls. This is true because of large transient and persistent disturbances, controller-structure interaction with unmodeled modes, poorly known structure parameters, slow actuator/sensor dynamical behavior, and excitation of nonlinear structure vibrations during control and assembly. For on-orbit assembly, controllers must start with finite element models of LSS and adapt on line to the best operating points, without compromising stability. This is not easy to do, since there are often unmodeled dynamic interactions between the controller and the structure. The indirect adaptive controllers are based on parameter estimation. Due to the large number of modes in LSS, this approach leads to very high-order control schemes with consequent poor stability and performance. In contrast, direct model reference adaptive controllers operate to force the LSS to track the desirable behavior of a chosen model. These schemes produce simple control algorithms which are easy to implement on line. One problem with their use for LSS has been that the model must be the same dimension as the LSS - i.e., quite large. A control theory based on the command generator tracker (CGT) ideas of Sobel, Mabins, Kaufman and Wen, Balas to obtain very low-order models based on adaptive algorithms was developed. Closed-loop stability for both finite element models and distributed parameter models of LSS was proved. In addition, successful numerical simulations on several LSS databases were obtained. An adaptive controller based on our theory was also implemented on a flexible robotic manipulator at Martin Marietta Astronautics. Computation schemes for controller-structure interaction with unmodeled modes, the residual mode filters or RMF, were developed. The RMF theory was modified to compensate slow actuator/sensor dynamics. These new ideas are being applied to LSS simulations to demonstrate the ease with which one can incorporate slow actuator/sensor effects into our design. It was also shown that residual mode filter compensation can be modified for small nonlinearities to produce exponentially stable closed-loop control. A theory for disturbance accommodating controllers based on reduced order models of structures was developed, and stability results for these controllers in closed-loop with large-scale finite element models of structures were obtained.

Balas, Mark↗

Estimating short-period dynamics using an extended Kalman filter

An extended Kalman filter (EKF) is used to estimate the parameters of a low-order model from aircraft transient response data. The low-order model is a state space model derived from the short-period approximation of the longitudinal aircraft dynamics. The model corresponds to the pitch rate to stick force transfer function currently used in flying qualities analysis. Because of the model chosen, handling qualities information is also obtained. The parameters are estimated from flight data as well as from a six-degree-of-freedom, nonlinear simulation of the aircraft. These two estimates are then compared and the discrepancies noted. The low-order model is able to satisfactorily match both flight data and simulation data from a high-order computer simulation. The parameters obtained from the EKF analysis of flight data are compared to those obtained using frequency response analysis of the flight data. Time delays and damping ratios are compared and are in agreement. This technique demonstrates the potential to determine, in near real time, the extent of differences between computer models and the actual aircraft. Precise knowledge of these differences can help to determine the flying qualities of a test aircraft and lead to more efficient envelope expansion.

Bauer, Jeffrey E.↗

Computational control of flexible aerospace systems

The main objective of this project is to establish a distributed parameter modeling technique for structural analysis, parameter estimation, vibration suppression and control synthesis of large flexible aerospace structures. This report concentrates on the research outputs produced in the last two years. The main accomplishments can be summarized as follows. A new version of the PDEMOD Code had been completed based on several incomplete versions. The verification of the code had been conducted by comparing the results with those examples for which the exact theoretical solutions can be obtained. The theoretical background of the package and the verification examples has been reported in a technical paper submitted to the Joint Applied Mechanics & Material Conference, ASME. A brief USER'S MANUAL had been compiled, which includes three parts: (1) Input data preparation; (2) Explanation of the Subroutines; and (3) Specification of control variables. Meanwhile, a theoretical investigation of the NASA MSFC two-dimensional ground-based manipulator facility by using distributed parameter modeling technique has been conducted. A new mathematical treatment for dynamic analysis and control of large flexible manipulator systems has been conceived, which may provide an embryonic form of a more sophisticated mathematical model for future modified versions of the PDEMOD Codes.

Sharpe, Lonnie, Jr.↗

Feasibility Study of SDAS Instrumentation's Ability to Identify Mobile Launcher (ML)/Crawler-Transporter (CT) Modes During Rollout Operations

The Space Launch System (SLS) and its Mobile Launcher (ML) will be transported to the launch pad via the Crawler-Transporter (CT) system. Rollout (i.e., transportation) loads produce structural loads on the integrated SLS/Orion Multi-Purpose Crew Vehicle (MPCV) launch vehicle which are of a concern with respect to fatigue. As part of the risk reduction process and in addition to the modal building block test approach that has been adopted by the SLS Program, acceleration data will be obtained during rollout for use in modal parameter estimation. There are several occurrences where the ML/CT will be transported either into the Vertical Assembly Building (VAB) or to the launch pad and back without the SLS stack as part of the Kennedy Space Center (KSC) Exploration Ground Systems (EGS) Integrated Test and Checkout (ITCO). NASA KSC EGS has instrumentation installed on both the ML and CT to record data during rollout, at the launch pad, and during liftoff. The EGS instrumentation on the ML, which includes accelerometers, is referred to as the Sensor Data Acquisition System (SDAS). The EGS instrumentation on the CT, which also includes accelerometers, is referred to as the CT Data Acquisition System (CTDAS). The forces and accelerations applied to the ML and CT during a rollout event will be higher than any of the planned building block modal tests. This can be very beneficial in helping identify nonlinear behavior in the structure. Developing modal parameters from the same test hardware in multiple boundary conditions and under multiple levels of excitation is a key step in developing a well correlated FEM. The purpose of this study was three fold. First, determine the target modes of the ML/CT in its rollout configuration. Second, determine if the test degrees of freedom (DOF) corresponding to the layout of the SDAS/CTDAS accelerometers (i.e. position and orientation) is sufficient to identify the target modes. Third, determine if the Generic Rollout Forcing Functions (GRFF's) is sufficient for identifying the ML/CT target modes accounting for variations in CT speed, modal damping, and sensor/ambient background noise levels. The finding from the first part of this study identified 28 target modes of the ML/CT rollout configuration based upon Modal Effective Mass Fractions (MEFF) and engineering judgement. The finding from the second part of this study showed that the SDAS/CTDAS accelerometers (i.e. position and orientation) are able to identify a sufficient number of the target modes to support model correlation of the ML/CT FEM. The finding from the third part of this study confirms the GRFFs sufficiently excite the ML/CT such that varying quantities of the defined target modes should be able to be extracted when utilizing an Experimental Modal Analysis (EMA) Multi-Input Multi-Output (MIMO) analysis approach. An EMA analysis approach was used because Operational Modal Analysis (OMA) tools were not available and the GRFFs were sufficiently uncorrelated. Two key findings from this third part of the study are that the CT speed does not show a significant impact on the ability to extract the modal parameters and that keeping the ambient background noise observed at each accelerometer location at or below 30 µgrms is essential to the success of this approach.

Winkel, James P.↗

Catalog-level blinding on the bispectrum for DESI-like galaxy surveys

We evaluate the performance of the catalog-level blind analysis technique (blinding) presented in Brieden et al. (2020) in the context of a fixed template power spectrum and bispectrum analysis. This blinding scheme, which is tailored for galaxy redshift surveys similar to the Dark Energy Spectroscopic Instrument (DESI), has two components: the so-called “AP blinding” (concerning the dilation parameters α$_{∥}$, α$_{⊥}$) and “RSD blinding” (redshift space distortions, affecting the growth rate parameter f). Through extensive testing, including checks for the RSD part in cubic boxes, the impact of AP blinding on mocks with realistic survey sky coverage, and the implementation of a full AP+RSD blinding pipeline, our analysis demonstrates the effectiveness of the technique in preserving the integrity of cosmological parameter estimation when the analysis includes the bispectrum statistic. We emphasize the critical role of sophisticated — and difficult to accidentally unblind — blinding methods in precision cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

ECLEIRS: Exact conservation law embedded identification of reduced states for parameterized nonlinear conservation laws from sparse and noisy data

Multi-query applications such as parameter estimation, uncertainty quantification and design optimization for parameterized partial differential equation (PDE) systems are expensive. While reduced/latent state dynamics approaches for parameterized PDEs offer a viable alternative, these approaches rely on high-quality data and struggle with highly sparse spatiotemporal noisy measurements typically obtained from experiments. Furthermore, there is no guarantee that these models satisfy governing physical conservation laws. In this article, we propose a reduced state dynamics approach, referred to as ECLEIRS, that embeds exact conservation in the solution and flux representation by utilizing a space-time divergence-free neural network formulation. We compare ECLEIRS with other reduced state dynamics approaches, those that do not enforce any physical constraints and those with physics-informed loss functions, for three shock-propagation problems: 1-D advection, 1-D Burgers and 2-D Euler equations. In conclusion, the numerical experiments conducted in this study demonstrate that ECLEIRS provides the most accurate prediction of dynamics for unseen parameters even in the presence of highly sparse and noisy data.

97 MATHEMATICS AND COMPUTING↗

Stroboscopic Goggles as a Countermeasure for Dynamic Visual Acuity and Landing Sickness in Crewmembers Returning from Long-Duration Spaceflight

Long-term exposure to microgravity causes sensorimotor adaptations that result in functional deficits upon returning to a gravitational environment. At landing, the vestibular system and the central nervous system, responsible for coordinating head and eye movements via the vestibulo-occular reflex (VOR), are adapted to microgravity and must re-adapt to the Earth's gravitational environment. This re-adaptation causes decrements in gaze control and dynamic visual acuity, with astronauts reporting oscillopsia and blurred vision. These effects are caused by retinal slip, or the inability to keep an image focused on their retina, which is thought to drive motion sickness symptoms experienced upon landing. Retinal slip can be estimated by dynamic visual acuity (DVA); visual acuity while in motion. Peters et al. (2011) find that DVA is worsened in astronauts by an average of 0.75 eye-chart lines one day after landing. Previously, the use of stroboscopic goggles has shown to be effective in minimizing motion sickness symptoms due to retinal slip (Reschke et al. 2007). In this study, we simulated the decrement in DVA caused by sensorimotor re-adaptation by using minifying lenses and then testing the efficacy of stroboscopic goggles in preventing retinal slip and improving DVA. Dynamic visual acuity is assessed using an oscillating chair developed in the Neuroscience Laboratory at JSC. This chair is motor-driven and oscillates vertically at 2 Hz with a vertical displacement of +/- 2 cm to simulate the vertical translations that occur while walking. As the subject is being oscillated, they are asked to discern the direction of Landolt-C optotypes of varying sizes and record their direction using a gamepad. The visual acuity thresholds are determined using an algorithm that alters the size of the optotype based on the previous responses of the subject using a forced-choice best parameter estimation that is able to rapidly converge on the threshold value. Visual acuity thresholds were determined both for static (seated) and dynamic (oscillating) conditions. Dynamic visual acuity is defined as the difference between the dynamic and static conditions. We found that healthy subjects (n=20) have a significantly impaired DVA while wearing the minifying lenses, demonstrating that the VOR is in an adaptive state and retinal slip is occurring. When subjects' acuity was tested wearing the stroboscopic goggles with the minifying lenses, there was no significant difference in their DVA compared to their baseline DVA. This suggests that stroboscopic goggles are preventing retinal slip and would function as an efficient countermeasure for VOR adaptations and thus help mitigate landing sickness symptoms experienced by long-duration crewmembers. These goggles might also be used to counter blurred vision (caused by retinal slip) experienced by crewmembers during launch where the vehicle vibrations are greatest. The stroboscopic effect could be built into a section of their head mounted displays on the visor of their helmets to be used in these high vibration situation if a mission critical task is necessary.

Rosenberg, M. J. F.↗

Assimilation of Precipitation Information Retrieved from Passive Microwave Sensors

Estimates of surface precipitation and the associated vertical latent heating structure can be retrieved from space-borne passive microwave sensors such as the TRMM Microwave Imager (TMI) and Special Sensor Microwave/Imager (SSM/I). Studies have shown that assimilation of microwave rain rates can improve the quality of global analysis and forecast. Numerical weather prediction (NWP) centers are beginning to make operational use of these data. In the next few years, there will be a gradual increase in microwave rain products available from operational and research satellites, culminating to a target constellation of 9 satellites to provide global rain measurements every 3 hours with the proposed Global Precipitation Measurement (GPM) mission to be launched around 2007. Realizing the full potential of these observations in data assimilation will require continued advances in retrieval algorithms, assimilation techniques, and model physics. I present results from two research efforts at NASA GSFC. The first is on the assimilation of tropical precipitation using temperature/moisture tendency corrections within a continuous variational assimilation framework to improve analyses and forecasts produced by the Goddard Earth Observing System (GEOS) global data assimilation system. Results show that assimilation of TMI and SSM/I rain rates improves not only precipitation and moisture, but also the related clouds, radiation energy fluxes, and large-scale circulations in GEOS analyses. The improved analyses also yield better short-range forecasts, quantitative precipitation forecast (QPF) threat scores, and storm track predictions. The second part consists of results from exploratory experiments on variational assimilation of convective and stratiform latent heating profiles within the general framework of model parameter estimation as a way to identify systematic model errors and optimize physical paremeterization schemes in global forecast and climate models.

Hou, Arthur Y.↗

Fast, differentiable, and extensible big bang nucleosynthesis package

Here, we introduce light isotope nucleosynthesis with JAX (LINX), a new differentiable public big bang nucleosynthesis code designed for fast parameter estimation. By leveraging JAX, LINX achieves both speed and differentiability, enabling the use of Bayesian inference, including gradient-based methods. We discuss the formalism used in LINX for rapid primordial elemental abundance predictions and give examples of how LINX can be used. When combined with differentiable cosmic microwave background power spectrum emulators, LINX can be used for joint cosmic microwave background and big bang nucleosynthesis analyses without requiring extensive computational resources, including on personal hardware.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗