Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling 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 253 records · Page 14

AVHRR Surface Temperature and Narrow-Band Albedo Comparison with Ground Measurements for the Greenland Ice Sheet

An ice-surface temperature retrieval algorithm for the Greenland ice sheet was developed using NOAA 11 thermal radiances from channels 4 and 5. Temperature, pressure and humidity profiles, cloud observations and skin temperatures from the Swiss Federal Institute of Technology (ETH) camp, located at the equilibrium line altitude at 49 deg17 min W, 69 deg 34 min N, were used in the LOWTRAN 7 model. Through a statistical analysis of daily clear sky profiles, the coefficients that correct for the atmospheric effects were determined for the ETH-Camp field season (May to August). Surface temperatures retrieved by this method were then compared against the in situ observations with a maximum difference of 0.6 K. The NOAA 11 narrow-band planetary albedo values for channels 1 and 2 were calculated using pre-launch calibration coefficients. Scattering and absorption by the atmosphere were modelled with LOWTRAN 7. Then, narrow-band albedo values for the AVHRR visible and near infrared channels were compared with in situ high resolution spectral reflectance measurements. In the visible band (580-680 nm), AVHRR-derived narrow-band albedo and the in situ measurements corrected with radiative transfer model LOWTRAN 7 showed a difference of less than 2%. For the near infrared channel (725-1100 nm) the difference between the measured and modelled narrow-band albedo was 14%. These discrepancies could be either the result of inaccurate aerosol scattering modelling (lack of the in situ observation), or the result of sensor drift due to degradation.

Haefliger, M.↗

Inadequacy of the statistical many-body model for molecular infrared absorption in the far wings

It is known that the infrared absorption coefficient of CO2 in the far wings is far less than what is derived from the assumption of a Lorentzian absorption profile. Sorgen and Ben-Aryeh (1972) tried to show that by incorporating the statistical many-body model into impact theory, the discrepancy between the Lorentzian absorption profile and the experimental infrared absorption coefficient can be removed. This paper attempts to show that the high frequency or short-time property of the dipole-dipole interaction cannot be adequately explained by the many-body model; in other words, the statistical model as conventionally understood is not applicable to infrared spectral lines in the far wings.

Fowler, B. W.↗

Frequency based localization of structural discrepancies

The intent of modal analysis is to develop a reliable model of a structure by working with the analytical and experimental modal properties of frequency, damping and mode shape. In addition to identifying these modal properties, it would be desirable to determine spatially which parts of the structure are modelled poorly or well. It is shown how the pattern of discrepancies in the analytical and experimental test values for the pole and the driving point zero frequencies of a structure can be linked to discrepancies in the mass or stiffness of the structural elements. The success of the procedure depends on the numerical conditioning of a modal reference matrix. Strategies to insure adequate numerical conditioning require a formulation which avoids geometric and energy storage symmetries of the structure, and ignores structural elements which contribute negligibly small potential or kinetic energy to the excited modes. Physical insight into the numerical conditioning problem is provided by a numerical example and by localization of a mass discrepancy in a real structure based on lab tests.

Shepard, G. D.↗

Grain Destruction in a Supernova Remnant Shock Wave

Dust grains are sputtered away in the hot gas behind shock fronts in supernova remnants, gradually enriching the gas phase with refractory elements. We have measured emission in C IV (lambda)1550 from C atoms sputtered from dust in the gas behind a non-radiative shock wave in the northern Cygnus Loop. Overall, the intensity observed behind the shock agrees approximately with predictions from model calculations that match the Spitzer 24 micron and the X-ray intensity profiles. Thus these observations confirm the overall picture of dust destruction in SNR shocks and the sputtering rates used in models. However, there is a discrepancy in that the CIV intensity 10'' behind the shock is too high compared to the intensities at the shock and 25'' behind it. Variations in the density, hydrogen neutral fraction and the dust properties over parsec scales in the pre- shock medium limit our ability to test dust destruction models in detail.

dust↗

A critical analysis of CLO and O3 in the mid-latitude stratosphere

The Caltech one-dimensional photochemical model is used to analyze important questions concerning the upper stratospheric O3 and the key HO(x), NO(x), and ClO(x) free radicals. The model is described and first order effects of the inclusion of diffuse radiation in a spherical rather than plane parallel atmosphere at solar zenith angles close to 90 deg are assessed. A comparison is made between photochemical theory and observations for upper stratospheric O3, where local photochemistry rather than dynamics should control the O3 abundance. The photochemical equilibrium relation for odd oxygen is studied, including the important radicals involves. It is concluded that a significant model ozone deficit exists in the upper stratosphere. Possible causes for such a discrepancy are investigated in light of current model and experimental uncertainties. In particular, the observational ratio of atomic oxygen to zone appears to disagree with the model results.

Froidevaux, L.↗

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗

The use of transition heights for the representation of ion composition

Several characteristic transition heights can be found in ionospheric ion-composition profiles. These are the boundaries between the regions dominated by light ions, atomic oxygen ions, molecular ions, and cluster ions, respectively. Ion-composition modeling can benefit from the use of these transition heights. Special emphasis is given to the 'upper' transition height H sub T (light ions to atomic oxygen) and to the 'lower' transition height h sub t (atomic oxygen to molecular ions). Transition-height models deduced from rocket and satellite measurements are compared with the heights predicted by the International Reference Ionosphere. Considerable discrepancies are found between the different models as well as between the models and independent measurements.

Bilitza, Dieter↗

A critical evaluation of semianalytic methods in the study of centrally heated, unresolved, infrared sources

We critically evaluate current methods of analysis in infrared (IR) astronomy and investigate the conditions under which these semianalytic methods are reliable. Specifically we examine the usual assumptions of homogeneities in dust density and temperature, and neglect of opacity effects when applied to internally heated, unresolved IR sources. To accomplish this, a series of radiation transport models for these sources have been constructed. The model results are treated as observed quantities and analyzed to derive the source parameters, using simple semianalytic methods. The discrepancies between the derived and actual model parameters can then be attributed to the limitations of the analysis methods and provide a measure of their reliability. Applying this approach to centrally heated, unresolved IR sources, we have studied in detail the following diagnostic problems: (1) determination of dust mass from monochromatic and integrated luminosities; (2) estimation of dust temperature distribution from color temperatures derived from the flux spectrum; and (3) determination of the empirical grain emissivity law (opacity function) for both continuum and spectral features.

Doty, Steven D.↗

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.↗

Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line

Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.

Chilldown↗

Validation Experiments for Turbulent Separation over an Axisymmetric Body of Revolution

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges to validating numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires data from advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed measurements of the flow field undergoing a smooth, adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states progressing from fully attached, to incipient separation, and finally to small-scale separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configurations of the axisymmetric model host a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation, as well as a slightly larger adverse pressure gradient, inducing a small-scale region of turbulent separation. Experiments include steady pressure measurements and 2-D PIV to provide the preliminary dataset for simulation studies, which examine the effect of variable grid domains and RANS turbulence models. This is done in an effort to understand and evaluate the critical variability between solutions for the present model configurations. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and slight variability in reference parameters. Ongoing work will focus on higher fidelity experimental campaigns to obtain surface flow visualizations and Stereo-Particle Image Velocimetry (SPIV) to deliver higher spatial resolution of the three-dimensional flow field to aid turbulence modelers.

Validation↗

Modeling of the 22 July 2009 Storm

The magnetic storm on 22 July 2009 is the greatest storm observed since the summer of 2008 when TWINS began its stereo imaging of the magnetosphere, On 22 July 2009, the Dst dropped to nearly -80 nT at 07:00 and 10:00 UT. During the main phase and at the peak of the storm, TWINS 1 and 2 were near apogee and moving from pre-dawn to post-dawn local time. The energetic neutral atom (ENA) Imagers on the 2 spacecraft captured the storm intensification and the formation of the partial ring current. The development of this storm has been simulated using the Comprehensive Ring Current Model (CRCM) to understand and interpret the observed signatures. The CRCM reproduced the double-dip in the Dst index and the simulated ENA flux intensities agree very well with the TWINS images. However, the peak of ion flux predicted by the model is always eastward of the observed maximum by TWINS. This discrepancy posts a challenge to reexamine the physical models employed in the CRCM.

Fok, Mei-Ching↗

Quantitative calculations of helium ion escape fluxes from the polar ionospheres

An attempt is made to reduce the discrepancy between the earlier theoretical model and the recent experimental observations of H(+) outward fluxes for winter and summer hemispheres by improving the theoretical model. The helium photoionization cross sections used are accurate to 10%, the latest solar EUV fluxes measured by the Atmosphere Explorer satellites are incorporated, and the most recent MSIS model of the neutral atmosphere is contained in the model. A range of conditions covering solar cycle, seasonal, and geomagnetic conditions are studied. The results show a maximum H(+) escape flux of 1.4 times 10 to the 7th per sq cm per sec for solar maximum, winter, low magnetic activity conditions, which is within the scatter of the measured fluxes. The computed summer H(+) escape flux is a factor of 20 lower than the winter value. Possible reasons for the slight discrepancy between theory and experiment in summer are discussed.

Raitt, W. J.↗

Vibrational kinetics in CO electric discharge lasers - Modeling and experiments

A model of CO laser vibrational kinetics is developed, and predicted vibrational distributions are compared with measurements. The experimental distributions were obtained at various flow locations in a transverse CW discharge in supersonic (M = 3) flow. Good qualitative agreement is obtained in the comparisons, including the prediction of a total inversion at low discharge current densities. The major area of discrepancy is an observed loss in vibrational energy downstream of the discharge which is not predicted by the model. This discrepancy may be due to three-dimensional effects in the experiment which are not included in the model. Possible kinetic effects which may contribute to vibrational energy loss are also examined.

Stanton, A. C.↗

Coulomb Stress Accumulation along the San Andreas Fault System

Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 x 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1 - 10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.

Smith, Bridget↗

Bennu Shape Model Validation Methods and Results

The OSIRIS-REx Independent Shape Modeling Team, staffed by the Jet Propulsion Laboratory (JPL), has used stereophotoclinometry (SPC) to produce a three-dimensional shape model of Bennu. The SPC process is informed by (but is separate from) orbit determination, so it is important to ensure that the resulting shape model is consistent with all available spacecraft tracking data. Specifically, checking the shape model’s consistency with LIDAR measurements can illuminate any discrepancies, because LIDAR measurements are highly correlated with the shape model. This study focuses on the JPL experience with these LIDAR measurements and the greater context of the shape model validation process.

Lauretta, Dante S.↗

Measurement of Low Amounts of Precipitable Water Vapor Using Ground-Based Millimeterwave Radiometry

Extremely dry conditions characterized by amounts of precipitable water vapor (PWV) as as 1-2 mm commonly occur in high-latitude regions during the winter months. While such atmospheres carry only a few percent of the latent heat energy compared to tropical atmospheres, the effects of low vapor amounts on the polar radiation budget - both directly through modulation of longwave radiation and indirectly through the formation of clouds - are considerable. Accurate measurements of precipitable water vapor (PWV) during such dry conditions are needed to improve polar radiation models for use in understanding and predicting change in the climatically sensitive polar regions. To this end, the strong water vapor absorption at 183.310 GHz provides a unique means of measuring low amounts of PWV. Weighting function analysis, forward model calculations based upon a 7-year radiosonde dataset, and retrieval simulations consistently predict that radiometric measurements made using several millimeter-wavelength (MMW) channels near the 183 GHz line, together with established microwave (MW) measurements at the 22.235 GHz water vapor line and -3 1 GHz atmospheric absorption window can be used to determine within 5% uncertainty the full range of PWV expected in the Arctic. This unique collective capability stands in spite of accuracy limitations stemming from uncertainties due to the sensitivity of the vertical distribution of temperature and water vapor at MMW channels. In this study the potential of MMW radiometry using the 183 GHz line for measuring low amounts of PWV is demonstrated both theoretically and experimentally. The study uses data obtained during March 1999 as part of an experiment conducted at the Department of Energy s Cloud and Radiation Testbed (CART) near Barrow, Alaska. Several radiometers from both NOAA and NASA were deployed during the experiment to provide the first combined MMW and MW ground-based data set during dry arctic conditions. Single-channel retrievals of PWV were performed using the MW and MMW data. Discrepancies in the retrieved values were found to be consistent with differences observed between measured brightness temperatures (TBs) and forward-modeled TBs based on concurrent radiosonde profiles. These discrepancies are greater than can be explained by measurement error alone and are attributed to absorption model uncertainty. We discuss here the measurements, retrieval technique, and line model discrepancies along with difficulties and potential of MMW/MW PWV measurement.

Racette, Paul E.↗

Estimation of the accuracy of dynamic flight-determined coefficients

This paper discusses means of assessing the accuracy of maximum likelihood parameter estimates obtained from dynamic flight data. The commonly used analytical predictors of accuracy are compared from both statistical and simplified geometric standpoints. Emphasizing practical considerations, such as modeling error, the accuracy predictions are evaluated with real and simulated data. Improved computations of the Cramer-Rao bound to correct large discrepancies caused by colored noise and modeling error are presented. This corrected Cramer-Rao bound is the best available analytical predictor of accuracy. Engineering judgement, aided by such analytical tools, is the final arbiter of accuracy estimation.

Maine, R. E.↗