Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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 217 records · Page 12

Preliminary assessment of the accuracy and precision of TOPEX/POSEIDON altimeter data with respect to the large-scale ocean circulation

TOPEX/POSEIDON sea surface height measurements are examined for quantitative consistency with known elements of the oceanic general circulation and its variability. Project-provided corrections were accepted but are at tested as part of the overall results. The ocean was treated as static over each 10-day repeat cycle and maps constructed of the absolute sea surface topography from simple averages in 2 deg x 2 deg bins. A hybrid geoid model formed from a combination of the recent Joint Gravity Model-2 and the project-provided Ohio State University geoid was used to estimate the absolute topography in each 10-day period. Results are examined in terms of the annual average, seasonal average, seasonal variations, and variations near the repeat period. Conclusion are as follows: the orbit error is now difficult to observe, having been reduced to a level at or below the level of other error sources; the geoid dominates the error budget of the estimates of the absolute topography; the estimated seasonal cycle is consistent with prior estimates; shorter-period variability is dominated on the largest scales by an oscillation near 50 days in spherical harmonics Y(sup m)(sub 1)(theta, lambda) with an amplitude near 10 cm, close to the simplest alias of the M(sub 2) tide. This spectral peak and others visible in the periodograms support the hypothesis that the largest remaining time-dependent errors lie in the tidal models. Though discrepancies attribute to the geoid are within the formal uncertainties of the good estimates, removal of them is urgent for circulation studies. Current gross accuracy of the TOPEX/POSEIDON mission is in the range of 5-10 cm, distributed overbroad band of frequencies and wavenumbers. In finite bands, accuracies approach the 1-cm level, and expected improvements arising from extended mission duration should reduce these numbers by nearly an order of magnitude.

Wunsch, Carl↗

GHRS Observations of Cool, Low-Gravity Stars: The Outer Atmosphere and Wind of the Nearby K Supergiant Lambda Velorum - 5

UV spectra of lambda Velorum taken with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope are used to probe the structure of the outer atmospheric layers and wind and to estimate the mass-loss rate from this K5 lb-II supergiant. VLA radio observations at lambda = 3.6 cm are used to obtain an independent check on the wind velocity and mass-loss rate inferred from the UV observations, Parameters of the chromospheric structure are estimated from measurements of UV line widths, positions, and fluxes and from the UV continuum flux distribution. The ratios of optically thin C II] emission lines indicate a mean chromospheric electron density of log N(sub e) approximately equal 8.9 +/- 0.2 /cc. The profiles of these lines indicate a chromospheric turbulence (v(sub 0) approximately equal 25-36 km/s), which greatly exceeds that seen in either the photosphere or wind. The centroids of optically thin emission lines of Fe II and of the emission wings of self-reversed Fe II lines indicate that they are formed in plasma approximately at rest with respect to the photosphere of the star. This suggests that the acceleration of the wind occurs above the chromospheric regions in which these emission line photons are created. The UV continuum detected by the GHRS clearly traces the mean flux-formation temperature as it increases with height in the chromosphere from a well-defined temperature minimum of 3200 K up to about 4600 K. Emission seen in lines of C III] and Si III] provides evidence of material at higher than chromospheric temperatures in the outer atmosphere of this noncoronal star. The photon-scattering wind produces self-reversals in the strong chromospheric emission lines, which allow us to probe the velocity field of the wind. The velocities to which these self-absorptions extend increase with intrinsic line strength, and thus height in the wind, and therefore directly map the wind acceleration. The width and shape of these self-absorptions reflect a wind turbulence of approximately equal 9-21 km/s. We further characterize the wind by comparing the observations with synthetic profiles generated with the Lamers et al. Sobolev with Exact Integration (SEI) radiative transfer code, assuming simple models of the outer atmospheric structure. These comparisons indicate that the wind in 1994 can be described by a model with a wind acceleration parameter beta approximately 0.9, a terminal velocity of 29-33 km/s, and a mass-loss rate approximately 3 x 10(exp -9) solar M/yr. Modeling of the 3.6 cm radio flux observed in 1997 suggests a more slowly accelerating wind (higher beta) and/or a higher mass-loss rate than inferred from the UV line profiles. These differences may be due to temporal variations in the wind or from limitations in one or both of the models. The discrepancy is currently under investigation.

Carpenter, Kenneth G.↗

Using Bayesian Networks for Candidate Generation in Consistency-based Diagnosis

Consistency-based diagnosis relies heavily on the assumption that discrepancies between model predictions and sensor observations can be detected accurately. When sources of uncertainty like sensor noise and model abstraction exist robust schemes have to be designed to make a binary decision on whether predictions are consistent with observations. This risks the occurrence of false alarms and missed alarms when an erroneous decision is made. Moreover when multiple sensors (with differing sensing properties) are available the degree of match between predictions and observations can be used to guide the search for fault candidates. In this paper we propose a novel approach to handle this problem using Bayesian networks. In the consistency- based diagnosis formulation, automatically generated Bayesian networks are used to encode a probabilistic measure of fit between predictions and observations. A Bayesian network inference algorithm is used to compute most probable fault candidates.

Narasimhan, Sriram↗

X-Ray Detection and Processing Models for Spacecraft Navigation and Timing

The current primary method of deepspace navigation is the NASA Deep Space Network (DSN). High-performance navigation is achieved using Delta Differential One-Way Range techniques that utilize simultaneous observations from multiple DSN sites, and incorporate observations of quasars near the line-of-sight to a spacecraft in order to improve the range and angle measurement accuracies. Over the past four decades, x-ray astronomers have identified a number of xray pulsars with pulsed emissions having stabilities comparable to atomic clocks. The x-ray pulsar-based navigation and time determination (XNAV) system uses phase measurements from these sources to establish autonomously the position of the detector, and thus the spacecraft, relative to a known reference frame, much as the Global Positioning System (GPS) uses phase measurements from radio signals from several satellites to establish the position of the user relative to an Earth-centered fixed frame of reference. While a GPS receiver uses an antenna to detect the radio signals, XNAV uses a detector array to capture the individual xray photons from the x-ray pulsars. The navigation solution relies on detailed xray source models, signal processing, navigation and timing algorithms, and analytical tools that form the basis of an autonomous XNAV system. Through previous XNAV development efforts, some techniques have been established to utilize a pulsar pulse time-of-arrival (TOA) measurement to correct a position estimate. One well-studied approach, based upon Kalman filter methods, optimally adjusts a dynamic orbit propagation solution based upon the offset in measured and predicted pulse TOA. In this delta position estimator scheme, previously estimated values of spacecraft position and velocity are utilized from an onboard orbit propagator. Using these estimated values, the detected arrival times at the spacecraft of pulses from a pulsar are compared to the predicted arrival times defined by the pulsar s pulse timing model. A discrepancy provides an estimate of the spacecraft position offset, since an error in position will relate to the measured time offset of a pulse along the line of sight to the pulsar. XNAV researchers have been developing additional enhanced approaches to process the photon TOAs to arrive at an estimate of spacecraft position, including those using maximum-likelihood estimation, digital phase locked loops, and "single photon processing" schemes that utilize all available time data associated with each photon. Using pulsars from separate, non-coplanar locations provides range and range-rate measurements in each pulsar s direction. Combining these different pulsar measurements solves for offsets in position and velocity in three dimensions, and provides accurate overall navigation for deep space vehicles.

Sheikh, Suneel↗

Uncertainties in the Modelled CO2 Threshold for Antarctic Glaciation

frequently cited atmospheric CO2 threshold for the onset of Antarctic glaciation of approximately780 parts per million by volume is based on the study of DeConto and Pollard (2003) using an ice sheet model and the GENESIS climate model. Proxy records suggest that atmospheric CO2 concentrations passed through this threshold across the Eocene-Oligocene transition approximately 34 million years. However, atmospheric CO2 concentrations may have been close to this threshold earlier than this transition, which is used by some to suggest the possibility of Antarctic ice sheets during the Eocene. Here we investigate the climate model dependency of the threshold for Antarctic glaciation by performing offline ice sheet model simulations using the climate from 7 different climate models with Eocene boundary conditions (HadCM3L, CCSM3, CESM1.0, GENESIS, FAMOUS, ECHAM5 and GISS_ER). These climate simulations are sourced from a number of independent studies, and as such the boundary conditions, which are poorly constrained during the Eocene, are not identical between simulations. The results of this study suggest that the atmospheric CO2 threshold for Antarctic glaciation is highly dependent on the climate model used and the climate model configuration. A large discrepancy between the climate model and ice sheet model grids for some simulations leads to a strong sensitivity to the lapse rate parameter.

simulation↗

Sources, Sinks, and Transatlantic Transport of North African Dust Aerosol: A Multimodel Analysis and Comparison With Remote Sensing Data

This study evaluates model-simulated dust aerosols over North Africa and the North Atlantic from five global models that participated in the Aerosol Comparison between Observations and Models phase II model experiments. The model results are compared with satellite aerosol optical depth (AOD) data from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Sea-viewing Wide Field-of-view Sensor, dust optical depth (DOD) derived from MODIS and MISR, AOD and coarse-mode AOD (as a proxy of DOD) from ground-based Aerosol Robotic Network Sun photometer measurements, and dust vertical distributions/centroid height from Cloud Aerosol Lidar with Orthogonal Polarization and Atmospheric Infrared Sounder satellite AOD retrievals. We examine the following quantities of AOD and DOD: (1) the magnitudes over land and over ocean in our study domain, (2) the longitudinal gradient from the dust source region over North Africa to the western North Atlantic, (3) seasonal variations at different locations, and (4) the dust vertical profile shape and the AOD centroid height (altitude above or below which half of the AOD is located). The different satellite data show consistent features in most of these aspects; however, the models display large diversity in all of them, with significant differences among the models and between models and observations. By examining dust emission, removal, and mass extinction efficiency in the five models, we also find remarkable differences among the models that all contribute to the discrepancies of model-simulated dust amount and distribution. This study highlights the challenges in simulating the dust physical and optical processes, even in the best known dust environment, and stresses the need for observable quantities to constrain the model processes.

Africa↗

Meteorite Ablation at the HYMETS ARC-JET Facility

One of the primary goals of the Asteroid Threat Assessment Program (ATAP) is to understand the ablation phenomena of meteors (e.g., spallation, evaporation, melt flow, etc.). Therefore, a test campaign was executed at the Hypersonic Material Test System (HYMETS) facility on samples of an ordinary chondrite (e.g., Tamdakht) and samples of basalt from the Columbia plateau. High melt viscosities and low ablation rates are predicted for basalt due to higher concentrations of silica and alumina in the matrix phase. However, high-speed video and post-test analysis reveal that basalt experiences greater rates of ablation than Tamdakht. Post-test analysis of emission spectra reveals the presence of volatile species that may explain the discrepancy between model predictions and observed phenomena.

Meteor↗

Magnification, dust, and time-delay constraints from the first resolved strongly lensed Type Ia supernova iPTF16geu

We report lensing magnifications, extinction, and time-delay estimates for the first resolved, multiply imaged Type Ia supernova iPTF16geu, at z = 0.409, using Hubble Space Telescope (HST) observations in combination with supporting ground-based data. Multiband photometry of the resolved images provides unique information about the differential dimming due to dust in the lensing galaxy. Using HST and Keck AO reference images taken after the SN faded, we obtain a total lensing magnification for iPTF16geu of μ = 67.8(+2.6−2.9), accounting for extinction in the host and lensing galaxy. As expected from the symmetry of the system, we measure very short time-delays for the three fainter images with respect to the brightest one: −0.23 ± 0.99, −1.43 ± 0.74, and 1.36 ± 1.07 d. Interestingly, we find large differences between the magnifications of the four supernova images, even after accounting for uncertainties in the extinction corrections: ∆m(1) = −3.88(+0.07−0.06), ∆m(2) = −2.99(+0.09−0.08), ∆m(3) = −2.19(+0.14−0.15),and ∆m(4) = −2.40(+0.14−0.12) mag, discrepant with model predictions suggesting similar image brightnesses. A possible explanation for the large differences is gravitational lensing by substructures, micro- or millilensing, in addition to the large-scale lens causing the image separations. We find that the inferred magnification is insensitive to the assumptions about the dust properties in the host and lens galaxy.

S Dhawan↗

Main and Fringe Field Computations for the Electrostatic Quadrupoles of the Muon g-2 Experiment Storage Ring

We consider semi-infinite electrostatic deflectors with plates of different thickness, including plates with rounded edges, and we calculate their electrostatic potential and field using conformal mappings. To validate the calculations, we compare the fringe fields of these electrostatic deflectors with fringe fields of finite electrostatic capacitors, and we extend the study to fringe fields of adjacent electrostatic deflectors with consideration of electrostatic induction, where field falloffs of semi-infinite electrostatic deflectors are slower than exponential and thus behave differently from most magnetic fringe fields. Building on the success with electrostatic deflectors, we develop a highly accurate and fully Maxwellian conformal mappings method for calculation of main fields of electrostatic particle optical elements. A remarkable advantage of this method is the possibility of rapid recalculations with geometric asymmetries and mispowered plates. We use this conformal mappings method to calculate the multipole terms of the high voltage quadrupole used in the storage ring of the Muon g-2 Experiment (FNALE- 0989). Completing the methodological framework, we present a method for extracting multipole strength falloffs of a particle optical element from a set of Fourier mode falloffs. We calculate the quadrupole strength falloff and its effective field boundary (EFB) for the Muon g-2 quadrupole, which has explained the experimentally measured tunes, while simple estimates based on a linear model exhibited discrepancies up to 2%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dendrite-resolved, full-melt-pool phase-field simulations to reveal non-steady-state effects and to test an approximate model

In this report we study the epitaxial, columnar growth of (multiply oriented) dendrites/cells for a spot melt in a polycrystalline Al–Cu substrate using two-dimensional, phase-field, direct numerical simulations (DNS) at the full-melt-pool scale. Our main objective is to compare the expensive DNS model to a much cheaper but approximate “line” model in which a single-crystal phase-field simulation is confined to a narrow rectangular geometry. To perform this comparison, we develop algorithms that automatically extract quantities of interest (QoIs) from both DNS and line models. These QoIs allow us to quantitatively assess the assumptions in the line model and help us analyze its discrepancy with the DNS model. We consider four sets of heat source parameters, mimicking welding and additive manufacturing conditions, that create a combination shallow and deep melt pools. Our largest DNS simulation used 16K x 14K grid points in space. Our main findings can be summarized as follows. Under AM conditions, the QoIs of line models are in excellent agreement with the full DNS results for both shallow and deep melt pools. Under welding conditions, the primary spacing of the DNS model is smaller than the prediction of line model. We identify a geometric crowding effect that accounts for the discrepancies between the DNS and line models. We propose two potential mechanisms that determine the response of the microstructure to geometric crowding.

36 MATERIALS SCIENCE↗

Resolving the Discrepancies Between Empirical and Rayleigh Charge Limiting Models for Globular Proteins

Starting with the Rayleigh charge limiting model, a slightly different approach is used to account for the well-known discrepancy that exists between the said model and experimental ESI MS data for globular proteins. It is shown using published datasets that for globular proteins, the mass density ρ exhibits a weak second-order dependence on its mass M, according to ρ(M)∝M{sup -α}, α ~ 0.14. A direct equivalence established between ESI MS and x-ray techniques suggests a minimum but critical surface tension of 15.6 ± 5.2 mN/m for the droplet at the liquid-to-gas phase transition point. The packing density factor η for globular proteins is believed to lie between 1 (very tightly packed) and 4.6 (less tight, natively packed). While the Rayleigh charge limiting model has been linked historically to the CRM (J. Chem. Phys. 49:2240–2249, 1968; Anal. Chim. Acta 406:93–104, 2000), this paper does not expressly seek to justify the CRM, but rather uses empirical data and existing knowledge across subfields to help build a consistent picture of ESI MS phenomena that might be difficult to explain otherwise. These results would be useful in molecular dynamics (MD) simulations, understanding liquid-to-gas phase transitions and in opening up new routes for cross-calibration between ESI MS, IM MS, NMR and x-ray crystallography studies. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modeling the High Speed Research Cycle 2B Longitudinal Aerodynamic Database Using Multivariate Orthogonal Functions

The data for longitudinal non-dimensional, aerodynamic coefficients in the High Speed Research Cycle 2B aerodynamic database were modeled using polynomial expressions identified with an orthogonal function modeling technique. The discrepancy between the tabular aerodynamic data and the polynomial models was tested and shown to be less than 15 percent for drag, lift, and pitching moment coefficients over the entire flight envelope. Most of this discrepancy was traced to smoothing local measurement noise and to the omission of mass case 5 data in the modeling process. A simulation check case showed that the polynomial models provided a compact and accurate representation of the nonlinear aerodynamic dependencies contained in the HSR Cycle 2B tabular aerodynamic database.

Morelli, E. A.↗

Magnetic field orientations in Saturn's upper ionosphere inferred from Voyager radio occultations

The radio scintillations observed during occultations of Voyagers 1 and 2 by Saturn are analyzed to determine the morphology of plasma irregularities and hence the magnetic field orientation in Saturn's upper atmosphere. The measurement techniques, the weak scattering theory, and the method used to relate the observed radio scintillations to physical properties of the ionospheric irregularities are briefly described. Results on the spatial characteristics of the irregularities are presented, and the magnetic field orientation in Saturn's ionosphere is inferred. Although the occultation measurements generally confirm the accuracy of the Saturnian magnetic field model of Connerney et al. (1982), it is found that a small adjustment of the coefficients in that model's zonal harmonic expansion would remove the discrepancy between the model predictions and the measurements. A strategy for obtaining improved measurements of Saturn's magnetic field from radio occultation observations of scintillations and Faraday rotation using an orbiting spacecraft is briefly discussed.

Hinson, D. P.↗

Computational Modeling of Two Mars Powered Descent Vehicle Concepts Tested in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not suitable for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). This paper covers computational flowfield predictions completed at wind tunnel conditions and comparisons to the test data. On the blowing HIAD models, the time-averaged pressure inboard of the nozzles was generally well-predicted, especially if the nozzles are canted outward, when the nozzles are located further from the nose. At intermediate CobraMRV thrust coefficients, CFD pressures are more accurately predicted than they are for the HIAD models, largely due to the nozzle locations and pointing directions. Overall, the CFD pressure coefficients were predicted within 0.2 of the steady pressure measurements for all blowing models, with smaller discrepancies at higher HIAD thrust, and larger discrepancies at lower HIAD CobraMRV thrust. All HIAD models were predicted to have a gradually decreasing axial force coefficient as the total thrust increases, in agreement with available pressure sensitive paint data. On models with canted nozzles or with nozzles further from the nose, the force coefficient was slightly higher for a given thrust. On the CobraMRV model, the CFD also shows consistent results between solvers and follows trends revealed in the data; the aerodynamic force coefficient remains near the non-blowing value at a tunnel Mach number of 2.4 regardless of thrust, and increases above that level at a Mach number of 3.5, consistent with the discrete pressure data. CFD analysis at tunnel and flight conditions will continue as flight system designs concepts mature.

Supersonic Retropropulsion↗

Computational Modeling of Two Mars Powered Descent Vehicle Concepts Tested in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not suitable for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). This paper covers computational flowfield predictions completed at wind tunnel conditions and comparisons to the test data. On the blowing HIAD models, the time-averaged pressure inboard of the nozzles was generally well-predicted, especially if the nozzles are canted outward, when the nozzles are located further from the nose. At intermediate CobraMRV thrust coefficients, CFD pressures are more accurately predicted than they are for the HIAD models, largely due to the nozzle locations and pointing directions. Overall, the CFD pressure coefficients were predicted within 0.2 of the steady pressure measurements for all blowing models, with smaller discrepancies at higher HIAD thrust, and larger discrepancies at lower HIAD CobraMRV thrust. All HIAD models were predicted to have a gradually decreasing axial force coefficient as the total thrust increases, in agreement with available pressure sensitive paint data. On models with canted nozzles or with nozzles further from the nose, the force coefficient was slightly higher for a given thrust. On the CobraMRV model, the CFD also shows consistent results between solvers and follows trends revealed in the data; the aerodynamic force coefficient remains near the non-blowing value at a tunnel Mach number of 2.4 regardless of thrust, and increases above that level at a Mach number of 3.5, consistent with the discrete pressure data. CFD analysis at tunnel and flight conditions will continue as flight system designs concepts mature.

Supersonic Retropropulsion↗

An evaluation of the Tsyganenko magnetic field model

A data set of more than 22,000 vector averages of the magnetosphere magnetic field over 0.5 R(E) regions is used to evaluate Tsyganenko's 1982 and 1987 magnetospheric magnetic field models. The magnetic field predicted by the model in various regions is compared to observations to find systematic discrepancies which future models might address. While agreement is generally good, discrepancies are noted which include: (1) a lack of adequate field line stretching in the tail and ring current regions; (2) an inability to predict weak enough fields in the polar cusps; and (3) a deficiency of Kp as a predictor of the field configuration.

Fairfield, D. H.↗

Radiative Heating in MSL Entry: Comparison of Flight Heating Discrepancy to Ground Test and Predictive Models

During the recent entry of the Mars Science Laboratory (MSL), the heat shield was equipped with thermocouple stacks to measure in-depth heating of the thermal protection system (TPS). When only convective heating was considered, the derived heat flux from gauges in the stagnation region was found to be underpredicted by as much as 17 W/sq cm, which is significant compared to the peak heating of 32 W/sq cm. In order to quantify the contribution of radiative heating phenomena to the discrepancy, ground tests and predictive simulations that replicated the MSL entry trajectory were performed. An analysis is carried through to assess the quality of the radiation model and the impact to stagnation line heating. The impact is shown to be significant, but does not fully explain the heating discrepancy.

MSL↗

Moving closer to experimental level materials property prediction using AI

Abstract While experiments and DFT-computations have been the primary means for understanding the chemical and physical properties of crystalline materials, experiments are expensive and DFT-computations are time-consuming and have significant discrepancies against experiments. Currently, predictive modeling based on DFT-computations have provided a rapid screening method for materials candidates for further DFT-computations and experiments; however, such models inherit the large discrepancies from the DFT-based training data. Here, we demonstrate how AI can be leveraged together with DFT to compute materials properties more accurately than DFT itself by focusing on the critical materials science task of predicting “formation energy of a material given its structure and composition”. On an experimental hold-out test set containing 137 entries, AI can predict formation energy from materials structure and composition with a mean absolute error (MAE) of 0.064 eV/atom; comparing this against DFT-computations, we find that AI can significantly outperform DFT computations for the same task (discrepancies of $$>0.076$$ > 0.076 eV/atom) for the first time.

36 MATERIALS SCIENCE↗