Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interval analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Single-mode projection filters for modal parameter identification for flexible structures

Single-mode projection filters are developed for eigensystem parameter identification from both analytical results and test data. Explicit formulations of these projection filters are derived using the orthogonal matrices of the controllability and observability matrices in the general sense. A global minimum optimization algorithm is applied to update the filter parameters by using the interval analysis method. The updated modal parameters represent the characteristics of the test data. For illustration of this new approach, a numerical simulation for the MAST beam structure is shown by using a one-dimensional global optimization algorithm to identify modal frequencies and damping. The projection filters are practical for parallel processing implementation.

Huang, Jen-Kuang↗
NASA NTRS1992

Interval prediction in structural dynamic analysis

Methods for assessing the predictive accuracy of structural dynamic models are examined with attention given to the effects of modal mass, stiffness, and damping uncertainties. The methods are based on a nondeterministic analysis called 'interval prediction' in which interval variables are used to describe parameters and responses that are unknown. Statistical databases for generic modeling uncertainties are derived from experimental data and incorporated analytically to evaluate responses. Covariance matrices of modal mass, stiffness, and damping parameters are propagated numerically in models of large space structures by means of three methods. The test data tend to fall within the predicted intervals of uncertainty determined by the statistical databases. The present findings demonstrate the suitability of using data from previously analyzed and tested space structures for assessing the predictive accuracy of an analytical model.

Hasselman, Timothy K.↗
NASA NTRS1988

Modal identification using single-mode projection filters and comparison with ERA and MLE results

The Single-Mode Projection Filter (SPF) is a newly developed algorithm for eigensystem parameter identification from both analytical results and test data. The SPF is formulated with a single mode only and practical for parallel processing implementation. Explicit formulations of SPF are derived for the multi-input multi-output (MIMO) system by using the orthogonal matrices of the controllability and observability matrices in the general sense. The modal parameters of SPF are initially obtained from an analytical model in modal space. The experimental data are then processed through SPF to update its modal parameters and to minimize a cost function defined by the norm of an error matrix. The updated modal parameters represent the characteristics of the test data. A two-dimensional global minimum optimization algorithm is developed and applied for the filter update by using the interval analysis method. The SPF is developed based on a single-mode subsystem and identifies only one modal frequency and one modal damping within a specified region. For an n-modes structure, n SPF can be implemented for parallel processing to reduce the computational burden. The SPF is applied to analyze the simulated data for the MAST beam structure. The estimated modal parameters are comparable to those from the Eigensystem Realization Algorithm (ERA) and repeated modal frequencies are identified. The modal analysis of the Spacecraft Control Laboratory Experiment (SCOLE) data is also performed by using the ERA and the Maximum Likelihood Estimate (MLE). The result shows that the first five modal frequencies are very close from ERA and MLE. However, there are slight disparities in the damping rates and the computational burdens are quite different among these two algorithms.

Huang, Jen-Kuang↗
NASA NTRS1994

High charge state carbon and oxygen ions in Earth's equatorial quasi-trapping region

Observations of energetic (1.5 - 300 keV/e) medium-to-high charge state (+3 less than or equal to Q less than or equal to +7) solar wind origin C and O ions made in the quasi-trapping region (QTR) of Earth's magnetosphere are compared to ion trajectories calculated in model equatorial magnetospheric magnetic and electric fields. These comparisons indicate that solar wind ions entering the QTR on the nightside as an energetic component of the plasma sheet exit the region on the dayside, experiencing little or no charge exchange on the way. Measurements made by the CHarge Energy Mass (CHEM) ion spectrometer on board the Active Magnetospheric Particle Tracer Explorer/Charge Composition Explorer (AMPTE/CCE) spacecraft at 7 less than L less than 9 from September 1984 to January 1989 are the source of the new results contained herein: quantitative long-term determination of number densities, average energies, energy spectra, local time distributions, and their variation with geomagnetic disturbance level as indexed by Kp. Solar wind primaries (ions with charge states unchanged) and their secondaries (ions with generally lower charge states produced from primaries in the magnetosphere via charge exchange)are observed throughout the QTR and have distinctly different local time variations that persist over the entire 4-year analysis interval. During Kp larger than or equal to 3 deg intervals, primary ion (e.g., O(+6)) densities exhibit a pronounced predawn maximum with average energy minimum and a broad near-local-noon density minimum with average energy maximum. Secondary ion (e.g., O(+5)) densities do not have an identifiable predawn peak, rather they have a broad dayside maximum peaked in local morning and a nightside minimum. During Kp less than or equal to 2(-) intervals, primary ion density peaks are less intense, broader in local time extent, and centered near midnight, while secondary ion density local time variations diminish. The long-time-interval baseline helps to refine and extend previous observations; for example, we show that ionospheric contribution to O(+3)) is negligible. Through comparison with model ion trajectories, we interpret the lack of pronounced secondary ion density peaks colocated with the primary density peaks to indicate that: (1) negligible charge exchange occurs at L greater than 7, that is, solar wind secondaries are produced at L less than 7, and (2) solar wind secondaries do not form a significant portion of the plasma sheet population injected into the QTR. We conclude that little of the energetic solar wind secondary ion population is recirculated through the magnetosphere.

Christon, S. P.↗
NASA NTRS1995

Analysis of the heliospheric current sheet at Earth's orbit and model comparisons

IMP 8 magnetic field data for the first half of the year 1994, i.e., for about 6 solar rotations, are analyzed around regions of sector boundary crossings with the purpose of obtaining both gross- and fine-scale characteristics of the related heliospheric current sheets separating the observed sectors. For purposes of estimating the attitudes of the normals to the sector boundaries. analysis intervals (sometimes 30 min or more in length) allowing the field to fully complete an excursion of about 180 deg were used in the study, which consisted of variance analyses of the field within those intervals. The resulting boundary normals were analyzed and compared to known (generic) models of projected heliospheric current sheets and to a coronal field model for the same time period. One of the most outstanding features of the resulting ensemble of estimated boundary normals for this period is that they strongly prefer low inclinations, indicating that the observations do not support a 1 AU model that predicts a current sheet whose surface is approximately parallel with the sun's equator, such as the 'sombrero' model. They instead support a model that predicts a relatively high inclination current sheet at 1 AU. Also the normals assume a surprisingly large range of longitudes, somewhat favoring those consistent with a Parker model (45 deg and 225 deg) and/or radial alignment (0 deg and 180 deg). These boundary structures, as defined, are shown typically to be as broad as several hundred proton gyroradii, but having embedded within them very thin structures associated with stronger currents. Such thin structures have normals usually differing markedly from the gross boundary. For some crossings there are indications of a wave-like structure in the current sheet as it passed the spacecraft.

Lepping, R. P.↗
NASA NTRS2002

The effect of exercise and rest duration on the generation of venous gas bubbles at altitude

BACKGROUND: Decompression, as occurs with aviators and astronauts undergoing high altitude operations or with deep-sea divers returning to surface, can cause gas bubbles to form within the organism. Pressure changes to evoke bubble formation in vivo during depressurization are several orders of magnitude less than those required for gas phase formation in vitro in quiescent liquids. Preformed micronuclei acting as "seeds" have been proposed, dating back to the 1940's. These tissue gas micronuclei have been attributed to a minute gas phase located in hydrophobic cavities, surfactant-stabilized microbubbles, or arising from musculoskeletal activity. The lifetimes of these micronuclei have been presumed to be from a few minutes to several weeks. HYPOTHESIS: The greatest incidence of venous gas emboli (VGE) will be detected by precordial Doppler ultrasound with depressurization immediately following lower extremity exercise, with progressively reduced levels of VGE observed as the interval from exercise to depressurization lengthens. METHODS: In a blinded cross-over design, 20 individuals (15 men, 5 women) at sea level exercised by performing knee-bend squats (150 knee flexes over 10 min, 235-kcal x h(-1)) either at the beginning, middle, or end of a 2-h chair-rest period without an oxygen prebreathe. Seated subjects were then depressurized to 6.2 psia (6,706 m or 22,000 ft altitude equivalent) for 120 min with no exercise performed at altitude. RESULTS: Of the 20 subjects with VGE in the pulmonary artery, 10 demonstrated a greater incidence of bubbles with exercise performed just prior to depressurization, compared with decreasing bubble grades and incidence as the interval of rest increased prior to depressurization. No decompression illness was reported. CONCLUSIONS: There is a significant increase in decompression-induced bubble formation at 6.2 psia when lower extremity exercise is performed just prior to depressurization as compared with longer rest intervals. Analysis indicated that micronuclei half-life is on the order of an hour under these hypobaric conditions.

Controlled Clinical Trial↗
NASA NTRS1992

ICE/ISEE plasma wave data analysis

The interval reported on, from Jan. 1990 to Dec. 1991, has been one of continued processing and archiving of ICE plasma wave (pw) data and transition from analysis of ISEE 3 and ICE cometary data to ICE data taken along its cruise trajectory, where coronal mass ejections are the focus of attention. We have continued to examine with great interest the last year of ISEE 3's precomet phase, when it spent considerable time far downwind from Earth, recording conditions upstream, downstream, and across the very weak, distant flank bow shock. Among other motivations was the apparent similarity of some shock and post shock structures to the signatures of the bow wave surrounding comet Giacobini-Zinner, whose ICE-phase data was revisited. While pursuing detailed, second-order scientific inquiries still pending from the late ISEE 3 recordings, we have also sought to position ourselves for study of CME's by instituting a data processing format new to the ISEE 3/ICE pw detector. Processed detector output has always been summarized and archived in 24-hour segments, with all pw channels individually plotted and stacked one above the next down in frequency, with each channel calibrated separately to keep all data patterns equally visible in the plots, regardless of gross differences in energy content at the various frequencies. Since CME's, with their preceding and following solar wind plasmas, can take more than one day to pass by the spacecraft, a more condensed synoptic view of the pw data is required to identify, let alone assess, CME characteristics than has been afforded by the traditional routines. This requirement is addressed in a major new processing initiative in the past two years. Besides our own ongoing and fresh investigations, we have cooperated, within our resources, with studies conducted extramurally by distant colleagues irrespective of the phase of the ISEE 3/ICE mission under scrutiny. The remainder of this report summarizes our processing activities, our investigations, both internal and cooperative, our scientific results, and our publication activity.

Greenstadt, E. W.↗
NASA NTRS2016

Advanced Interval Management: A Benefit Analysis

This document is the final report for the NASA Langley Research Center (LaRC)- sponsored task order 'Possible Benefits for Advanced Interval Management Operations.' Under this research project, Architecture Technology Corporation performed an analysis to determine the maximum potential benefit to be gained if specific Advanced Interval Management (AIM) operations were implemented in the National Airspace System (NAS). The motivation for this research is to guide NASA decision-making on which Interval Management (IM) applications offer the most potential benefit and warrant further research.

Timer, Sebastian↗
NASA NTRS1986

Atropine unmasks bed-rest effect - A spectral analysis of cardiac interbeat intervals

Heart rate spectral data obtained for 10 male subjects between 35-49 years following orthostatic tolerance testing with lower body negative pressure prebed rest and after 7-10 days of bed rest, while on placebo and after intravenous atropine are analyzed. Comparison of the spectral atropine rms for subjects prebed rest and after bed rest reveal a decrease from 63 + or - 24 ms to 40 + or - 23 ms. It is observed that heart rate interval variability for subjects after bed rest and with atropine is reduced; the heart rate at bed rest with atropine is increased from 70.4 + or - 12.4 beats/min prebed rest to 83.7 + or - 18.9 beats/min; and the exercise tolerance time for subjects in the atropine prebed-rest phase (658 + or - 352 s) is higher than the bed-rest phase (505 + or - 252 s). It is noted that bed rest impairs the cardiovascular capacity to adaptively modulate physiological responses, atropine exposes bed-rest deconditioning effects, and spectral analysis is useful for studying the effects of bed-rest deconditioning on cardiac dynamics.

Goldberger, Ary L.↗
NASA NTRS1978

Average properties of X-ray burst sources

Total integrated fluxes, maximum fluxes, and spectra of X-ray bursts from 10 sources are analyzed, and a uniformity in observed properties is detected. If it is assumed that the peak luminosity in a burst is a standard candle, the data suggest that (1) the effective black body radius observed in bursts is the same for all burst sources; (2) the tails of X-ray bursts are emitted by cooling surfaces whose size is the same for all burst sources (to within + or - 20%); and (3) the total energy emitted from X-ray bursts is the same for all burst sources (to within + or - 30%). This analysis applies to Type I bursts which are emitted by all burst sources and repeat at comparatively long time intervals. The analysis does not apply to Type II bursts, which repeat on time scales of seconds to minutes.

Van Paradijs, J.↗
NASA NTRS1982

Stationarity of magnetohydrodynamic fluctuations in the solar wind

Solar wind research and studies of charged particle propagation often assume that the interplanetary magnetic field represents a stationary random process. The extent to which ensemble averages of the solar wind magnetic fields follow the asymptotic behavior predicted by the ergodic theorem was investigated. Several time periods, including a span of nearly two years, are analyzed. Data intervals which span many solar rotations satisfy the conditions of weak stationarity if the effects of solar rotation are included in the asymptotic analysis. Shorter intervals which include a small integral number of interplanetary sectors also satisfy weak stationarity. The results are illustrated using magnetometer data from the ISEE-3, Voyager and IMP spacecraft.

Matthaeus, W. H.↗
NASA NTRS2011

Assessment of Radiative Heating Uncertainty for Hyperbolic Earth Entry

This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.

Johnston, Christopher O.↗
DOE OSTI2025

Middle and Late-Holocene coastal environments of far Southeastern Russia as inferred from palynological and diatom data

Existing discontinuous palynological records from coastal and river valley exposures, with varying quality of radiocarbon control, suggested that the regional Holocene climate in southern areas of the Russian Far East was characterized by as many as 10 fluctuations in temperature and/or precipitation. In this study, palynological data from Zerkalnoye Lake, located on the western coast of the Sea of Japan, indicate a gradual decline in temperature through the Middle and Late-Holocene with a Holocene thermal maximum between 8800 and 5500 cal yr BP. Three wetter than present intervals, marked by an increase in Pinus koraiensis, occurred c. 3600–3500 cal yr BP, 2340–2050 cal yr BP, and 1830–1800 cal yr BP. The Zerkalnoye record shows the dominance of Quercus-broadleaf forests during the Middle and Early Holocene, although Quercus shows a gradual decrease as climate cooled during this interval. Diatom analysis of the Zerkalnoye sediments documents that the site was a shallow bay or coastal lagoon until c. 3540 cal yr BP. After that time, a freshwater lake was established, which had variable marine influences probably caused by sea level changes following deglaciation. The diatom data indicate a cool water interval between 2900 and 2580 cal yr BP, also noted in other sites in the region. Changes in the basin’s depositional environment do not affect the palynological record, indicating such sites can provide reliable paleovegetational and paleoclimatic records. In conclusion, the discrepancies of the various regional paleoclimatic scenarios indicate the need for the further collection of continuous records from coastal to alpine zones in this region of northeastern Asia.

Environmental sciences↗
NASA NTRS1986

Periodicity of the earth's magnetic reversals

Results are presented from an attempt to perform a relatively comprehensive analysis of the evidence for a periodicity, with harmonics, of the observed regular reversals of the earth's magnetic field. The database considered covers 296 reversals over the past 165 Myr. Histograms with bins 1 Myr apart reveal only 30 Myr reversal patterns. The reversal dates are fitted to a linear periodic function and a spectrum is computed for the residuals at the adopted dates. The possible presence of multiple periodicities is evaluated and over various time intervals. The analysis shows that a recently observed 15 Myr periodicity is probably a harmonic of the 29.5-30.5 Myr period. The calculations do not confirm an inherent magnetic reversal property of the earth. The reversals may arise from tectonic events or from impacts from extraterrestrial objects.

Stothers, R. B.↗
NASA NTRS2011

An Empirical State Error Covariance Matrix for Batch State Estimation

State estimation techniques serve effectively to provide mean state estimates. However, the state error covariance matrices provided as part of these techniques suffer from some degree of lack of confidence in their ability to adequately describe the uncertainty in the estimated states. A specific problem with the traditional form of state error covariance matrices is that they represent only a mapping of the assumed observation error characteristics into the state space. Any errors that arise from other sources (environment modeling, precision, etc.) are not directly represented in a traditional, theoretical state error covariance matrix. Consider that an actual observation contains only measurement error and that an estimated observation contains all other errors, known and unknown. It then follows that a measurement residual (the difference between expected and observed measurements) contains all errors for that measurement. Therefore, a direct and appropriate inclusion of the actual measurement residuals in the state error covariance matrix will result in an empirical state error covariance matrix. This empirical state error covariance matrix will fully account for the error in the state estimate. By way of a literal reinterpretation of the equations involved in the weighted least squares estimation algorithm, it is possible to arrive at an appropriate, and formally correct, empirical state error covariance matrix. The first specific step of the method is to use the average form of the weighted measurement residual variance performance index rather than its usual total weighted residual form. Next it is helpful to interpret the solution to the normal equations as the average of a collection of sample vectors drawn from a hypothetical parent population. From here, using a standard statistical analysis approach, it directly follows as to how to determine the standard empirical state error covariance matrix. This matrix will contain the total uncertainty in the state estimate, regardless as to the source of the uncertainty. Also, in its most straight forward form, the technique only requires supplemental calculations to be added to existing batch algorithms. The generation of this direct, empirical form of the state error covariance matrix is independent of the dimensionality of the observations. Mixed degrees of freedom for an observation set are allowed. As is the case with any simple, empirical sample variance problems, the presented approach offers an opportunity (at least in the case of weighted least squares) to investigate confidence interval estimates for the error covariance matrix elements. The diagonal or variance terms of the error covariance matrix have a particularly simple form to associate with either a multiple degree of freedom chi-square distribution (more approximate) or with a gamma distribution (less approximate). The off diagonal or covariance terms of the matrix are less clear in their statistical behavior. However, the off diagonal covariance matrix elements still lend themselves to standard confidence interval error analysis. The distributional forms associated with the off diagonal terms are more varied and, perhaps, more approximate than those associated with the diagonal terms. Using a simple weighted least squares sample problem, results obtained through use of the proposed technique are presented. The example consists of a simple, two observer, triangulation problem with range only measurements. Variations of this problem reflect an ideal case (perfect knowledge of the range errors) and a mismodeled case (incorrect knowledge of the range errors).

Frisbee, Joseph H., Jr.↗
DOE OSTI2022

Robust trajectory-constrained frequency control for microgrids considering model linearization error

Grid supportive modes integrated within inverter-based resources can improve the frequency response of renewable-rich microgrids. The synthesis of grid supportive modes to guarantee frequency trajectory constraints under a predefined disturbance set is challenging but essential. To tackle this challenge, a numerical optimal control (NOC)-based control synthesis methodology is proposed. Without loss of generality, a wind-diesel fed microgrid is studied, where we aim to design grid supportive functions in the wind turbine. In the control design, linearized models are used, and the linearization-induced errors are quantitatively analyzed by reachability and interval arithmetics and represented in the form of interval uncertainties. Then, the NOC problem can be formulated into a robust mixed-integer linear program. The control structure is strategically configured into two levels to realize online deployment. The proposed control is verified on the modified 33-node microgrid with a full-order three-phase nonlinear model in Simulink. In conclusion, the simulation results show the effectiveness of the proposed control paradigm and the necessity of considering linearization-induced uncertainty.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Monte Carlo Dropout Uncertainty Quantification of Long Short-Term Memory Autoencoder Anomaly Detection in a Liquid Sodium Cold Trap

Advanced high-temperature fluid reactors, such as sodium-cooled fast reactors (SFRs) and molten salt–cooled reactors (MSCRs), require coolant purification systems to prevent fluid contamination and local freezing that can lead to plugging. Liquid sodium purification can be achieved with a cold trap, where the sodium temperature is reduced to a near-freezing point to precipitate out impurities. Automation of monitoring of the cold trap performance with machine learning algorithms can aid in early detection of incipient anomalies. An efficient approach to loss-of-coolant–type anomaly detection in a cold trap monitored with more than two dozen thermal-hydraulic sensors consists of a long short-term memory (LSTM) autoencoder. This work develops the uncertainty quantification of the LSTM autoencoder performance for cold trap anomaly detection using the Monte Carlo (MC) dropout method. The MC dropout methodology creates a distribution of sister distributions that all slightly differ from each other because of random neurons being turned off for testing. The variances of the sister network distributions are used to make an uncertainty interval. Our analysis shows that the uncertainty in the autoencoder performance is largest near the peak of the anomaly signal. Using the MC dropout method, we investigate the uncertainty in the anomaly detection with missing sensor inputs. This capability allows the reactor operator to evaluate resilience of the anomaly detection system and to make informed decisions about continuity of operation in the event of sensor failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗