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 955 records · Page 53

Resistive hose modes in tokamak runaway electron beams

Beams of energetic runaway electrons are generated during disruptions in tokamaks, and fluid models are used to study their effects on macroscale dynamics. Linear computations of a massless, runaway electron beam coupled to MHD plasma show that resistive hose instabilities grow faster than tearing modes at large resistivity. Eigenvalue results with reduced models of the resistive hose instability are compared with results from the full MHD and beam system, showing that the resistive hose decouples from any plasma response. An estimate of plasma temperature at which growth of the resistive hose dominates tearing for post-disruption DIII-D plasma parameters is in a physically relevant regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control of linear uncertain systems utilizing mismatched state observers

The control of linear continuous dynamical systems is investigated as a problem of limited state feedback control. The equations which describe the structure of an observer are developed constrained to time-invarient systems. The optimal control problem is formulated, accounting for the uncertainty in the design parameters. Expressions for bounds on closed loop stability are also developed. The results indicate that very little uncertainty may be tolerated before divergence occurs in the recursive computation algorithms, and the derived stability bound yields extremely conservative estimates of regions of allowable parameter variations.

Goldstein, B.↗

The mass function of Seyfert 1 nuclei

The first mass function of Seyfert 1 nuclei is derived from optical spectra of the complete CfA sample of Seyfert galaxies by estimating the mass for each object from a dynamical relation. An independent estimate is also derived using a complete infrared-selected sample. The two mass functions are indistinguishable. The mean mass of Seyfert 1 nuclei is about 2 x 10 to the 7th solar masses, and the integrated mass density is about 6 x 10 to the 11th solar masses/cu Gpc. This is approximately two orders of magnitude less than the value inferred from the energetics associated with quasar counts. A careful analysis of the various parameters and assumptions involved suggests that this large difference is not due to systematic errors in the determinations. Therefore, the bulk of mass related to the accretion processes connected with past quasar activity does not reside in Seyfert 1 nuclei. Instead, the remnants of past activity must be present in a much larger number of galaxies, and a one-to-one relation between distant and local active galactic nuclei seems then to be excluded.

Padovani, P.↗

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗

Multiple rings around Wolf-Rayet evolution

We present optical narrow-band imaging of multiple rings existing around galactic Wolf-Rayet (WR) stars. The existence of multiple rings of material around Wolf-Rayet stars clearly illustrates the various phases of evolution that massive stars go through. The objects presented here show evidence of a three stage evolution. O stars produce an outer ring with the cavity being partially filled by ejecta from a red supergiant of luminous blue variable phase. A wind from the Wolf-Rayet star then passes into the ejecta materials. A simple model is presented for this three stage evolution. Using observations of the size and dynamics of the rings allows estimates of time scales for each stage of the massive star evolution. These are consistent with recent theoretical evolutionary models. Mass estimates for the ejecta, from the model presented, are consistent with previous ring nebula mass estimates from IRAS data, showing a number of ring nebulae to have large masses, most of which must in be in the form of neutral material. Finally, we illustrate how further observations will allow the determination of many of the parameters of the evolution of massive stars such as total mass loss, average mass loss rates, stellar abundances, and total time spent in each evolutionary phase.

Marston, A. P.↗

Results of a Parametric Study on 10,000 lb. and 50,000 lb. Permanent Magnetic Actuators

The purpose of this study was to generate design data and complete dynamic performance estimates for a high performance permanent magnet actuator. The basic configuration selected for analysis is an axisymmetric Nd-B-Fe permanent magnet actuator capable of providing force in one direction along its major axis. The actuator consisted of two main axisymmetric components separated by an air gap. The design was optimized for each value of force, gap and magnetic field to yield minimum weight and maximum lift to weight ratio. The basic conclusion is that, within parameters considered, the 10,000 lb. and 50,000 lb. actuators are lightweight and compact. As expected for most permanent magnet devices, the smaller ones have higher lift to eight ratios.

Stekly, Z. J. J.↗

In-Flight System Identification

A method is proposed and studied whereby the system identification cycle consisting of experiment design and data analysis can be repeatedly implemented aboard a test aircraft in real time. This adaptive in-flight system identification scheme has many advantages, including increased flight test efficiency, adaptability to dynamic characteristics that are imperfectly known a priori, in-flight improvement of data quality through iterative input design, and immediate feedback of the quality of flight test results. The technique uses equation error in the frequency domain with a recursive Fourier transform for the real time data analysis, and simple design methods employing square wave input forms to design the test inputs in flight. Simulation examples are used to demonstrate that the technique produces increasingly accurate model parameter estimates resulting from sequentially designed and implemented flight test maneuvers. The method has reasonable computational requirements, and could be implemented aboard an aircraft in real time.

Morelli, Eugene A.↗

Efficient and Collision-Free Anticipative Cruise Control in Randomly Mixed Strings

Connected intelligent vehicle following can improve safety and efficiency compared to today's road transport, but real traffic in the near future will not provide an ideal setting for its deployment. Unconnected human-driven vehicles will follow variable behavior patterns that combine with differences in dynamic capability to create heterogeneous scenarios. In this paper, connected automated vehicles employ model predictive control for following in traffic that may include both heavy and passenger vehicles at quasi-random positions. Adding further realism, some quasi-randomly mixed vehicles are completely unconnected and equipped with the reactive intelligent driver model using pseudorandom parameters. A mixed-integer quadratic programming formulation adapts the predictive algorithm to diverse powertrain operating point constraints. The preceding vehicle's control input is estimated from velocity and brake light observations and used to probabilistically generate a preview for the ego vehicle. Here, a terminal constraint designed using particle kinematics prevents collisions due to shortsightedness. Results are simulated at various heavy and predictive vehicle concentrations in the presence of packet loss. Linear fuel economy improvements between 1.4% and 1.9% per 10% increase in predictive vehicle penetration rate are shown.

33 ADVANCED PROPULSION SYSTEMS↗

Predicting the Earth encounters of (99942) Apophis

Arecibo delay-Doppler measurements of (99942) Apophis in 2005 and 2006 resulted in a five standard-deviation trajectory correction to the optically predicted close approach distance to Earth in 2029. The radar measurements reduced the volume of the statistical uncertainty region entering the encounter to 7.3% of the pre-radar solution, but increased the trajectory uncertainty growth rate across the encounter by 800% due to the closer predicted approach to the Earth. A small estimated Earth impact probability remained for 2036. With standard-deviation plane-of-sky position uncertainties for 2007-2010 already less than 0.2 arcsec, the best near-term ground-based optical astrometry can only weakly affect the trajectory estimate. While the potential for impact in 2036 will likely be excluded in 2013 (if not 2011) using ground-based optical measurements, approximations within the Standard Dynamical Model (SDM) used to estimate and predict the trajectory from the current era are sufficient to obscure the difference between a predicted impact and a miss in 2036 by altering the dynamics leading into the 2029 encounter. Normal impact probability assessments based on the SDM become problematic without knowledge of the object's physical properties; impact could be excluded while the actual dynamics still permit it. Calibrated position uncertainty intervals are developed to compensate for this by characterizing the minimum and maximum effect of physical parameters on the trajectory. Uncertainty in accelerations related to solar radiation can cause between 82 and 4720 Earth-radii of trajectory change relative to the SDM by 2036. If an actionable hazard exists, alteration by 2-10% of Apophis' total absorption of solar radiation in 2018 could be sufficient to produce a six standard-deviation trajectory change by 2036 given physical characterization; even a 0.5% change could produce a trajectory shift of one Earth-radius by 2036 for all possible spin-poles and likely masses. Planetary ephemeris uncertainties are the next greatest source of systematic error, causing up to 23 Earth-radii of uncertainty. The SDM Earth point-mass assumption introduces an additional 2.9 Earth-radii of prediction error by 2036. Unmodeled asteroid perturbations produce as much as 2.3 Earth-radii of error. We find no future small-body encounters likely to yield an Apophis mass determination prior to 2029. However, asteroid (144898) 2004 VD17, itself having a statistical Earth impact in 2102, will probably encounter Apophis at 6.7 lunar distances in 2034, their uncertainty regions coming as close as 1.6 lunar distances near the center of both SDM probability distributions.

Orbit determination↗

Investigation of Aerodynamic Matching Sensitivities for Conceptual Design

A methodology for aerodynamic matching was implemented in the Higher fidelity Conceptual Design and structural optimization (HCDstruct) tool to support transonic flutter simulations. The Enhanced Factor Correction Technique was implemented in HCDstruct to correct the Nastran doublet lattice method using data from corresponding computational fluid dynamics simulations. In this paper, sensitivities associated with this aerodynamic matching technique were investigated using HCDstruct and a Boeing 737-200-like aircraft configuration. Specifically, sensitivities of the aerodynamic matching method were investigated with respect to changes in flight conditions, such as changes in Mach number and angle of attack, and also in key design variables for the main wing, such as changes in aspect ratio and sweep angle. The variations in flight conditions and design parameters were mapped to estimated prediction error in sectional lift for the main wing, which helps to demonstrate the extent to which a user may expect a given aerodynamic correction matrix to provide suitable aerodynamic predictions for a realistic aircraft wing.

Jesse R. Quinlan↗

Investigation of Aerodynamic Matching Sensitivities for Conceptual Design

A methodology for aerodynamic matching was implemented in the Higher fidelity Conceptual Design and structural optimization (HCDstruct) tool to support transonic flutter simulations. The Enhanced Factor Correction Technique was implemented in HCDstruct to correct the Nastran doublet lattice method using data from corresponding computational fluid dynamics simulations. In this paper, sensitivities associated with this aerodynamic matching technique were investigated using HCDstruct and a Boeing 737-200-like aircraft configuration. Specifically, sensitivities of the aerodynamic matching method were investigated with respect to changes in flight conditions, such as changes in Mach number and angle of attack, and also in key design variables for the main wing, such as changes in aspect ratio and sweep angle. The variations in flight conditions and design parameters were mapped to estimated prediction error in sectional lift for the main wing, which helps to demonstrate the extent to which a user may expect a given aerodynamic correction matrix to provide suitable aerodynamic predictions for a realistic aircraft wing.

Jesse Ray Quinlan↗

Frequency Domain Quasi Maximum Likelihood Identification of Low Order Aeroservoelastic Models from Flight-Test Data

In this paper a quasi maximum likelihood method for estimating a low order model of a flexible vehicle has been developed and demonstrated. The quasi maximum likelihood method uses a large number of sensors to estimate the parameters and covariance in a method consistent with the maximum likelihood filter-error method. The cost function has been defined in a frequency domain to improve computational efficiency and to enable estimation of an unstable plant. The method has been demonstrated using flight data from the National Aeronautics and Space Administration X-56A Multi-Utility Technology Testbed flex wing tests, which include unstable flutter modes. The method was able to effectively estimate the frequency and damping of the dynamics of the aircraft and generate a transfer function that is useful for control system evaluations.

Jeffrey Ouellette↗

Cutting out the Middleman: Calibrating and Validating a Dynamic Vegetation Model (ED2-PROSPECT5) Using Remotely Sensed Surface Reflectance

Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant processes as well as understanding land–atmosphere feedbacks. In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty. Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified. Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state. We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types. Using a single common set of parameters across all sites, the calibrated model was able to accurately reproduce surface reflectance for sites with highly varied forest composition and structure. However, the calibrated model's predictions of leaf area index (LAI) were less robust, capturing only 46 % of the variability in the observations. Comparing the ED2 radiative transfer model with another two-stream soil–leaf–canopy radiative transfer model commonly used in remote sensing studies (PRO4SAIL) illustrated structural errors in the ED2 representation of direct radiation backscatter that resulted in systematic underestimation of reflectance. In addition, we also highlight that, to directly compare with a two-stream radiative transfer model like EDR, we had to perform an additional processing step to convert the directional reflectance estimates of AVIRIS to hemispherical reflectance (also known as “albedo”). In future work, we recommend that vegetation models add the capability to predict directional reflectance, to allow them to more directly assimilate a wide range of airborne and satellite reflectance products. We ultimately conclude that despite these challenges, using dynamic vegetation models to predict surface reflectance is a promising avenue for model calibration and validation using remote sensing data.

Alexey N Shiklomanov↗

Bottomonium suppression and elliptic flow using Heavy Quarkonium Quantum Dynamics

We introduce a framework called Heavy Quarkonium Quantum Dynamics (HQQD) which can be used to compute the dynamical suppression of heavy quarkonia propagating in the quark-gluon plasma using real-time in-medium quantum evolution. Using HQQD we compute large sets of real-time solutions to the Schrödinger equation using a realistic in-medium complex-valued potential. We sample 2 million quarkonia wave packet trajectories and evolve them through the QGP using HQQD to obtain their survival probabilities. The computation is performed using three different HQQD model parameter sets in order to estimate our systematic uncertainty. After taking into account final state feed down we compare our results to existing experimental data for the suppression and elliptic flow of bottomonium states and find that HQQD predictions are good agreement with available data for R AA as a function of N part and p T collected at \( \sqrt{s_{\mathrm{NN}}} \) = 5 . 02 TeV. In the case of v 2 for the various states, we find that the path-length dependence of ϒ(1 s ) suppression results in quite small v 2 for ϒ(1 s ). Our prediction for the integrated elliptic flow for ϒ(1 s ) in the 10 – 90% centrality class, which now includes an estimate of the systematic error, is v 2 [ϒ(1 s )] = 0 . 003 ± 0 . 0007 ± \( {}_{0.0013}^{0.0006} \) . We also find that, due to their increased suppression, excited bottomonium states have a larger elliptic flow. Based on this observation we make predictions for v 2 [ϒ(2 s )] and v 2 [ϒ(3 s )] as a function of centrality and transverse momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fluctuation-driven excess noise near superconducting phase transition

We discuss intrinsic mechanisms of nonequilibrium excess noise in superconducting devices and transition edge sensors. In particular, we present an overview of fluctuation-driven contributions to the current noise in the vicinity of the superconducting transition. We argue that sufficiently close to the critical temperature fluctuations of conductivity may become correlated provided that the rate of quasiparticle relaxation is slow as compared to dynamics of superconducting fluctuations. In this regime, fluctuations of conductivity adiabatically follow the fluctuations of the electron distribution. This leads to a substantial enhancement of current noise. Here, the corresponding spectral power density of noise has a Lorentzian shape in the frequency domain while its magnitude scales proportionally to the inelastic relaxation time. It also sensitively depends on the dephasing and Ginzburg–Landau timescales. Further estimates suggest that this mechanism dominates over the conventional temperature fluctuations in the same range of parameters. To describe these effects microscopically, we use the nonequilibrium Keldysh technique in the semiclassical approximation of superconductivity with Boltzmann–Langevin random forces to account for correlations of fluctuations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Morphotype-resolved characterization of microalgal communities in a nutrient recovery process with ARTiMiS flow imaging microscopy

Microalgae-driven nutrient recovery represents a promising technology for phosphorus removal from wastewater while simultaneously generating biomass that can be valorized to offset treatment costs. As full-scale processes come online, system parameters including biomass composition must be carefully monitored to optimize performance and prevent culture crashes. In this study, flow imaging microscopy (FIM) was leveraged to characterize microalgal community composition in near real-time at a full-scale municipal wastewater treatment plant (WWTP) in Wisconsin, USA, and population and morphotype dynamics were examined to identify relationships between water chemistry, biomass composition, and system performance. Two FIM technologies, FlowCam and ARTiMiS, were evaluated as monitoring tools. ARTiMiS provided a more accurate estimate of total system biomass, and estimates derived from particle area as a proxy for biovolume yielded better approximations than particle counts. Deep learning classification models trained on annotated image libraries demonstrated equivalent performance between FlowCam and ARTiMiS, and convolutional neural network (CNN) classifiers proved significantly more accurate when compared to feature table-based dense neural network (DNN) models. Across a two-year study period, Scenedesmus spp. appeared most important for phosphorus removal, and were negatively impacted by elevated temperatures and increase in nitrite/nitrate concentrations. Chlorella and Monoraphidium also played an important role in phosphorus removal. For both Scenedesmus and Chlorella, smaller morphological types were more often associated with better system performance, whereas larger morphotypes likely associated with stress response(s) correlated with poor phosphorus recovery rates. Furthermore, these results demonstrate the potential of FIM as a critical technology for high-resolution characterization of industrial microalgal processes.

59 BASIC BIOLOGICAL SCIENCES↗

Small Body GN and C Research Report: G-SAMPLE - An In-Flight Dynamical Method for Identifying Sample Mass [External Release Version]

G-SAMPLE is an in-flight dynamical method for use by sample collection missions to identify the presence and quantity of collected sample material. The G-SAMPLE method implements a maximum-likelihood estimator to identify the collected sample mass, based on onboard force sensor measurements, thruster firings, and a dynamics model of the spacecraft. With G-SAMPLE, sample mass identification becomes a computation rather than an extra hardware requirement; the added cost of cameras or other sensors for sample mass detection is avoided. Realistic simulation examples are provided for a spacecraft configuration with a sample collection device mounted on the end of an extended boom. In one representative example, a 1000 gram sample mass is estimated to within 110 grams (95% confidence) under realistic assumptions of thruster profile error, spacecraft parameter uncertainty, and sensor noise. For convenience to future mission design, an overall sample-mass estimation error budget is developed to approximate the effect of model uncertainty, sensor noise, data rate, and thrust profile error on the expected estimate of collected sample mass.

in-flight sample-mass identification↗