Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 883 records · Page 49

Evaluation of a load cell model for dynamic calibration of the rotor systems research aircraft

The Rotor Systems Research Aircraft uses load cells to isolate the rotor/transmission system from the fuselage. An analytical model of the relationship between applied rotor loads and the resulting load cell measurements is derived by applying a force-and-moment balance to the isolated rotor/transmission system. The model is then used to estimate the applied loads from measured load cell data, as obtained from a ground-based shake test. Using nominal design values for the parameters, the estimation errors, for the case of lateral forcing, were shown to be on the order of the sensor measurement noise in all but the roll axis. An unmodeled external load appears to be the source of the error in this axis.

Duval, R. W.↗

A nudged hybrid analysis and modeling approach for realtime wake-vortex transport and decay prediction

Here, we put forth a long short-term memory (LSTM) nudging framework for the enhancement of reduced order models (ROMs) of fluid flows utilizing noisy measurements for air traffic improvements. Toward emerging applications of digital twins in aviation, the proposed approach allows for constructing a realtime predictive tool for wake-vortex transport and decay systems. We build on the fact that in realistic application, there are uncertainties in initial and boundary conditions, model parameters, as well as measurements. Moreover, conventional nonlinear ROMs based on Galerkin projection (GROMs) suffer from imperfection and solution instabilities, especially for advection-dominated flows with slow decay in the Kolmogorov n-width. In the presented LSTM nudging (LSTM-N) approach, we fuse forecasts from a combination of imperfect GROM and uncertain state estimates, with sparse Eulerian sensor measurements to provide more reliable predictions in a dynamical data assimilation framework. We illustrate our concept by solving the two-dimensional vorticity transport equation. We investigate the effects of measurements noise and state estimate uncertainty on the performance of the LSTM-N behavior. We also demonstrate that it can sufficiently handle different levels of temporal and spatial measurement sparsity, and offer a huge potential in developing next-generation digital twin technologies for aerospace applications.

42 ENGINEERING↗

Utilization of Hydroturbines in Wastewater Treatment Plants

Wastewater treatment plants (WWTPs) are a significant energy consumer, yet there are several opportunities for implementing on-site power generation systems. Within the treatment process, the high flow rate of effluent is produced and discharged to a nearby water body by gravity. Thus, hydroturbines can be utilized to generate power in such an application due to a difference in elevation and high flow rate. This paper presents a case study of introducing a hydroturbine in a WWTP in Wisconsin and evaluating the power output in addition to determining the energy savings. The WWTP considered in this study has an effluent flow rate of 190 MGD (million gallons per day) and elevation difference of 3 m (10 ft) between the final stage of treatment and the discharge point. Based on the parameters above; hubless rim-drive Kaplan type hydroturbine (RDT) is the optimal choice to be used in such an application. The RDT is designed and optimized by using in-house code. A computational fluid dynamics (CFD) software is applied to evaluate the performance of the proposed model, and the system is simulated through homer software to validate the results generated by the CFD. The expected savings is estimated to be 1564 MWh/yr.

Energy & Fuels↗

Computational Assessment of Hydrodynamic Loads on Rockeries for River Bank Protection

Colorado experienced an extreme flooding event in September of 2013 that ravaged canyon roadways southeast of Estes Park, CO. Central Federal Lands was asked to rebuild SH 34 and CO 43 after the floods. To minimize river encroachments in narrow canyons, roadway fill embankments must be over steepened. Using a rockery to stabilize the roadway embankment rather than traditional vertical, concrete structures was the solution. However, no analytical techniques were available for evaluating the hydraulic stability of such structures. Argonne National Laboratory was employed to estimate the hydrodynamic forces by using advanced, computational fluid dynamics (CFD). An investigation was conducted on a section of a model rockery comprised of large boulders exposed to a variety of hydraulic flow conditions. A repeatable geometric section of a typical rockery structure was modeled and used to create boundaries for the CFD model. The test case matrix covers computer simulations of various rockery and channel geometries. The geometrical parameters taken into account are: rockery height, outside face batter angle of the rockery, rockery base width, river bed width, and the hydraulic parameters are: flow depth, flow velocity, and hydraulic angle of attack. Two failure mechanisms of rockery structures, sliding and overturning, are considered. The existing design equations for both modes are modified by addition of hydrodynamic force terms. The new terms were calculated for a parametric set of geometry and condition variation using the CFD model. Moreover, formulas for dimensionless hydraulic force coefficients are proposed that can be used in the design of new rockeries or in evaluation of existing dry-stack rockeries in river environments.

42 ENGINEERING↗

Boundary layer roll circulations during FIRE

The probable mechanism underlying the development of boundary layer roll circulations are studied using wind and temperature profiles measured by the National Center for Atmospheric Research (NCAR) Electra during the stratocumulus phase of the First ISCCP Regional Experiment (FIRE). The expected, or preferred, roll orientations, horizontal wavelengths, and propagation periods are determined by finding the minimum values of the dynamic and thermodynamic forcing parameters, which here are the eddy Reynolds number (Re) and moist Rayleigh number (Ra sub m). These minimum values depend on the height z sub T of the capping temperature inversion and on the values of the Fourier coefficients of the background height-dependent vector wind profile. As input to our nonlinear spectral model, descent and ascent runs by the Electra provide for initial estimates of the inversion height and the wind profiles. In the first phase of the investigation presented here, a mechanism is said to be a probable contributor to the development of roll circulations within the stratocumulus-topped boundary layer if the modeled roll orientation and wavelengths agree with their observed values. Preliminary results using the 14-coefficient model of Haack-Hirschberg (1988) are discussed for the 7 July 1987 Electra Mission 188-A (Flight 5). This mission was flown across a sharp cloud boundary that was within a LANDSAT/SPOT scene. The stratocumulus deck was relatively solid in the eastern part of the scene, while there was a rapid decrease in cloud cover to scattered cumulus clouds aligned in streets to the west. These cloud streets were oriented nearly parallel to the mean wind direction in the layer, which was approximately 340 degrees. The hypothesis that roll circulations occurred in both the relatively clear and the cloudy regions is investigated using as model input a descent profile obtained in the relatively clear air and an ascent profile obtained in the cloudy air. Initial results for the clear air case are that the pure inflection point mode is not possible and the pure thermal mode was oriented 35 degrees to the right of the mean wind direction. The origin of this unacceptably large discrepancy between the observed and modeled results will be investigated further and the conclusions reported at the next FIRE workshop.

Shirer, Hampton N.↗

Evaluation of the Event Driven Phenology Model Coupled with the VegET Evapotranspiration Model Through Comparisons with Reference Datasets in a Spatially Explicit Manner

A new model coupling scheme with remote sensing data assimilation was developed for estimation of daily actual evapotranspiration (ET). The scheme represents a mix of the VegET, a physically based model to estimate ET from a water balance, and an event driven phenology model (EDPM), where the EDPM is an empirically derived crop specific model capable of producing seasonal trajectories of canopy attributes. In this experiment, the scheme was deployed in a spatially explicit manner within the croplands of the Northern Great Plains. The evaluation was carried out using 2007-2009 land surface forcing data from the North American Land Data Assimilation System (NLDAS) and crop maps derived from remotely sensed data of NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). We compared the canopy parameters produced by the phenology model with normalized difference vegetation index (NDVI) data derived from the MODIS nadir bi-directional reflectance distribution function (BRDF) adjusted reflectance (NBAR) product. The expectations of the EDPM performance in prognostic mode were met, producing determination coefficient (r2) of 0.8 +/-.0.15. Model estimates of NDVI yielded root mean square error (RMSE) of 0.1 +/-.0.035 for the entire study area. Retrospective correction of canopy dynamics with MODIS NDVI brought the errors down to just below 10% of observed data range. The ET estimates produced by the coupled scheme were compared with ones from the MODIS land product suite. The expected r2=0.7 +/-.15 and RMSE = 11.2 +/-.4 mm per 8 days were met and even exceeded by the coupling scheme0 functioning in both prognostic and retrospective modes. Minor setbacks of the EDPM and VegET performance (r2 about 0.5 and additional 30 % of RMSR) were found on the peripheries of the study area and attributed to the insufficient EDPM training and to spatially varying accuracy of crop maps. Overall the experiment provided sufficient evidence of soundness and robustness of the EDPM and VegET coupling scheme, assuring its potential for spatially explicit applications.

Kovalskyy, V.↗

Particle Filtering for Model-Based Anomaly Detection in Sensor Networks

A novel technique has been developed for anomaly detection of rocket engine test stand (RETS) data. The objective was to develop a system that postprocesses a csv file containing the sensor readings and activities (time-series) from a rocket engine test, and detects any anomalies that might have occurred during the test. The output consists of the names of the sensors that show anomalous behavior, and the start and end time of each anomaly. In order to reduce the involvement of domain experts significantly, several data-driven approaches have been proposed where models are automatically acquired from the data, thus bypassing the cost and effort of building system models. Many supervised learning methods can efficiently learn operational and fault models, given large amounts of both nominal and fault data. However, for domains such as RETS data, the amount of anomalous data that is actually available is relatively small, making most supervised learning methods rather ineffective, and in general met with limited success in anomaly detection. The fundamental problem with existing approaches is that they assume that the data are iid, i.e., independent and identically distributed, which is violated in typical RETS data. None of these techniques naturally exploit the temporal information inherent in time series data from the sensor networks. There are correlations among the sensor readings, not only at the same time, but also across time. However, these approaches have not explicitly identified and exploited such correlations. Given these limitations of model-free methods, there has been renewed interest in model-based methods, specifically graphical methods that explicitly reason temporally. The Gaussian Mixture Model (GMM) in a Linear Dynamic System approach assumes that the multi-dimensional test data is a mixture of multi-variate Gaussians, and fits a given number of Gaussian clusters with the help of the wellknown Expectation Maximization (EM) algorithm. The parameters thus learned are used for calculating the joint distribution of the observations. However, this GMM assumption is essentially an approximation and signals the potential viability of non-parametric density estimators. This is the key idea underlying the new approach.

Solano, Wanda↗

Adiabatic Gravitational Waveform Model for Compact Objects Undergoing Quasi-Circular Inspirals Into Rotating Massive Black Holes

We present bhpwave: a new python-based, open-source tool for generating the gravitational waveforms of stellar-mass compact objects undergoing quasicircular inspirals into rotating massive black holes. These binaries, known as extreme-mass-ratio inspirals (EMRIs), are exciting mHz gravitational wave sources for future space-based detectors such as the Laser Interferometer Space Antenna (LISA). Relativistic models of EMRI gravitational wave signals are necessary to unlock the full scientific potential of mHz detectors, yet few open-source EMRI waveform models exist. Thus we built bhpwave, which uses the adiabatic approximation from black hole perturbation theory to rapidly construct gravitational waveforms based on the leading-order inspiral dynamics of the binary. In this work, we present the theoretical and numerical foundations underpinning bhpwave. We also demonstrate how bhpwave can be used to assess the impact of EMRI modeling errors on LISA gravitational wave data analysis. In particular, we find that for retrograde orbits and slowly spinning black holes we can mismodel the gravitational wave phasing by as much as ∼10 radians without significantly biasing EMRI parameter estimation.

Zachary Nasipak↗

The GEODYN Precision Orbit Determination Software

GEODYN-II is NASA’s state-of-the-art geodetic parameter estimation and precision orbit determination (POD) system, developed at Goddard Space Flight Center and maintained within the Geodesy and Geophysics Laboratory. GEODYN supports POD and geodetic parameter estimation for both Earth orbiting and planetary/planetary body missions. It implements a comprehensive set of high-fidelity force models for integrating satellite orbits. GEODYN supports approximately 100 measurement types, including support for Global Navigation Satellite System (GNSS), Satellite Laser Ranging (SLR), Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS), Very Long Base Interferometry (VLBI), landmark (camera data), landmark crossovers, radar and laser direct altimeter ranging, and dynamic altimetry crossovers.

Joseph B Nicholas↗

Correlation Filtering of Modal Dynamics using the Laplace Wavelet

Wavelet analysis allows processing of transient response data commonly encountered in vibration health monitoring tasks such as aircraft flutter testing. The Laplace wavelet is formulated as an impulse response of a single mode system to be similar to data features commonly encountered in these health monitoring tasks. A correlation filtering approach is introduced using the Laplace wavelet to decompose a signal into impulse responses of single mode subsystems. Applications using responses from flutter testing of aeroelastic systems demonstrate modal parameters and stability estimates can be estimated by correlation filtering free decay data with a set of Laplace wavelets.

Freudinger, Lawrence C.↗

Estimating lithospheric properties at Atla Regio, Venus

Magellan spehrical harmonic gravity and topography models are used to estimate lithospheric properties at Alta Regio, Venus, a proposed hotspot with dynamic support from mantle plume(s). Global spherical harmonic and local representations of the gravity field share common properties in the Atla region interms of their spectral behavior over a wavelength band from approximately 2100 to approximately 700 km. The estimated free-air admittance spectrum displays a rather featureless long-wavelength portion followed by a sharp rise at wavelengths shorter than about 1000 km. This sharp rise requires significant flexural support of short-wavelength structures. The Bouguer coherence also displays a sharp drop in this wavelength band, indicating a finite flexural rigidity of the lithosphere. A simple model for lithospheric loading from above and below is introduced (D. W. Forsyth, 1985) with four parameters: f, the ratio of bottom loading to top loading; z(sub m), crustal thickness; z(sub l) depth to bottom loading source; and T(sub e) elastic lithosphere thickness. A dual-mode compensation model is introduced in which the shorter wavelengths (lambda approximately less than 1000 km) might be explained best by a predominance of top loading by the large shield volcanoes Maat Mons, Ozza Mons, and Sapas Mons, and the longer wavelengths (lambda approximately greater than 1500 km) might be explained best by a deep depth of compensation, possibly representing bottom loading by a dynamic source. A Monte Carlo inversion technique is introduced to thoroughly search out the four-space of the model parameters and to examine parameter correlation in the solutions. Venus either is a considerabe deficient in heat sources relative to Earth, or the thermal lithosphere is overthickened in response to an earlier episode of significant heat loss from the planet.

Phillips, Roger J.↗

Quantum criticality using a superconducting quantum processor

Quantum criticality emerges from the collective behavior of many interacting quantum particles, often at the transition between different phases of matter. It is one of the cornerstones of condensed matter physics, which we access on noisy intermediate-scale (NISQ) quantum devices by leveraging a dynamically driven phenomenon. Here we probe the critical properties of the one-dimensional quantum Ising model on a programmable superconducting quantum chip via a Kibble-Zurek process, obtain scaling laws, and estimate critical exponents despite inherent sources of errors on the hardware. In addition, we investigate how the improvement of NISQ computers (more qubits, less noise) will consolidate the computation of those universal physical properties. A one-parameter noise model captures the effect of imperfections and reproduces the experimental data. Its systematic study reveals that the noise, analogously to temperature, induces a new length scale in the system. We introduce and successfully verify modified scaling laws, directly accounting for the noise without any prior knowledge. It makes data analyses for extracting physical properties transparent to noise. By understanding how imperfect quantum hardware modifies the genuine properties of quantum states of matter, we enhance the power of NISQ processors considerably for addressing quantum criticality and potentially other phenomena and algorithms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fed-DeepONet: Stochastic Gradient-Based Federated Training of Deep Operator Networks

The Deep Operator Network (DeepONet) framework is a different class of neural network architecture that one trains to learn nonlinear operators, i.e., mappings between infinite-dimensional spaces. Traditionally, DeepONets are trained using a centralized strategy that requires transferring the training data to a centralized location. Such a strategy, however, limits our ability to secure data privacy or use high-performance distributed/parallel computing platforms. To alleviate such limitations, in this paper, we study the federated training of DeepONets for the first time. That is, we develop a framework, which we refer to as Fed-DeepONet, that allows multiple clients to train DeepONets collaboratively under the coordination of a centralized server. To achieve Fed-DeepONets, we propose an efficient stochastic gradient-based algorithm that enables the distributed optimization of the DeepONet parameters by averaging first-order estimates of the DeepONet loss gradient. Then, to accelerate the training convergence of Fed-DeepONets, we propose a moment-enhanced (i.e., adaptive) stochastic gradient-based strategy. Finally, we verify the performance of Fed-DeepONet by learning, for different configurations of the number of clients and fractions of available clients, (i) the solution operator of a gravity pendulum and (ii) the dynamic response of a parametric library of pendulums.

Moya, Christian↗

Digital active control law synthesis for aeroservoelastic systems

This paper presents a formulation for synthesis of digital active control laws for aeroservoelastic systems, which are typically modeled by large order equations in order to accurately represent the rigid and flexible body modes, unsteady aerodynamic forces, actuator dynamics, and gust spectra. The control law is expected to satisfy multiple design requirements on the dynamic loads, responses, actuator deflection and rate limitations, as well as maintain certain stability margins, yet should be simple enough to be implemented by an onboard digital microprocessor. The synthesis procedure minimizes a linear quadratic Gaussian type cost function, by updating selected free parameters of the control law, while satisfying a set of inequality constraints on the design loads, responses and stability margins. A stable classical control law or an estimator based full or reduced order control law can be used as an initial design starting point. The gradients of the cost function and the constraints, with respect to the digital control law design variables are derived analytically, to facilitate rapid convergence. Selected design responses can be treated as constraints instead of lumping them into the cost function, in order to satisfy individual root-mean-square load and response limitations. Constraints are also imposed on the minimum singular value requirements for stability robustness improvement.

Mukhopadhyay, Vivekananda↗

Geopotential Model Improvement Using POCM_4B Dynamic Ocean Topography Information: PGM2000A

The two-year mean (1993-1994) Dynamic Ocean Topography (DOT) field implied by the POCM_4B circulation model was used to develop normal equations for DOT, in a surface spherical harmonic representation. These normal equations were combined with normal equations from satellite tracking data, surface gravity data, and altimeter data from TOPEX/Poseidon and ERS-1. Several least-squares combination solutions were developed in this fashion, by varying parameters such as the maximum degree of the estimated DOT and the relative weights of the different data. The solutions were evaluated in terms of orbit fit residuals, GPS/Leveling-derived undulations, and independent DOT information from in situ WOCE hydrographic data. An optimal solution was developed in this fashion which was originally presented at the 1998 EGS meeting in Nice, France. This model, designated here PGM2000A, maintains the orbit and land geoid modeling performance of EGM96, while improving its marine geoid modeling capability. In addition, PGM2000A's error spectrum is considerably more realistic than those of other contemporary gravitational models and agrees well with the error spectrum of EGM96. We will present the development and evaluation of PGM2000A, with particular emphasis on the weighting of the DOT information implied by POCM_4B. We will also present an inter-comparison of PGM2000A with the GRIM5-C1 and TEG-4 models. Directions for future work and problematic areas will be identified.

Pavlis, N. K.↗

Incremental States for Precise On-Orbit Relative Knowledge in Formation Flight

High precision and close proximity formation flight is an enabling technology for future space missions and requires an on-board relative navigation capability that is accurate to the mm-level and robust to formation parameters. Common estimation techniques linearize the entire formation about one spacecraft’s position, resulting in degraded accuracy due to linearization errors when separation distances become large. Additionally, methods which decouple absolute and relative state estimates usually require ad-hoc methods to incorporate the estimates together. This work discusses an alternative ”incremental” formulation of the relative navigation problem which is in variant to formation size, robust to coupling between absolute and relative dynamics, and can undergo similarity transformations to smoothly incorporate either absolute or relative information without numerical issues. A specific example of this architecture is presented in the context of formation navigation using Carrier-Differential Global Positioning System measurements, and is compared to a traditional leader-linearized filter. In the presence of accurate measurements, the incremental architecture is shown to reduce linearization errors and mean estimate errors by up to three orders of magnitude without the incorporation of new sensor information.

Seubert, Carl↗

Galaxy Cluster Mass Reconstruction Project – IV. Understanding the effects of imperfect membership on cluster mass estimation

Here, the primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters with precise spectroscopic redshifts. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples, respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated–true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

79 ASTRONOMY AND ASTROPHYSICS↗

Paramagnon dispersion and damping in doped Na x Ca 2 – x CuO 2 Cl 2

Here, using resonant inelastic x-ray scattering, we measure the paramagnon dispersion and damping of undoped, antiferromagnetic Ca 2 CuO 2 Cl 2 as well as doped, superconducting Na x Ca 2–x CuO 2 Cl 2 . Our estimation of the spin-exchange parameter and width of the paramagnon peak at the zone boundary X=(0.5, 0) confirms that no simple relation can be drawn between these parameters and the critical temperature T c . Consistently with other cuprate compounds, we show that upon doping there is a slight softening at (0.25,0) but not at the zone boundary X. In combination with these measurements we perform calculations of the dynamical spin structure factor of the one-band Hubbard model using cluster dynamical mean-field theory. The calculations are in excellent agreement with the experiment in the undoped case, both in terms of energy position and width. While the increase in width is also captured upon doping, the dynamical spin structure factor shows a sizable softening at X, which provides insightful information on the length-scale of the spin fluctuations in doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗