Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model comparison”

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 1,009 records · Page 56

Data Analysis and Theoretical Studies of the Upper Mesosphere and Lower Thermosphere

Three separate tasks were proposed under this award. The first involved extending our continuing study of electrodynamical feedback between the thermosphere/ionosphere and the magnetosphere. The second was a model-experiment comparison study of global dynamics and the third was a 'spectral energetics' analysis of tidal dissipation and energy exchange mechanisms. The Earth's mesosphere and lower-thermosphere/ionosphere (MLTI), between approximately 60 and 180 km altitude, is the most poorly understood region of the Earth's atmosphere, primarily because of its relative inaccessibility. This lack of knowledge has been widely recognized and has provided important scientific rationale for the upcoming NASA TIMED mission. While the data gathered during the TIMED era will revolutionize our understanding of the MLTI region, much work can be done prior to the mission, both to develop data-analysis and modeling techniques and to study the more limited relevant experimental data from previous missions. The grant reported on here continues and extends an existing successful program of scientific research into the energetics, dynamics and electrodynamics of the MLTI, using available theoretical and data analysis tools.

Burns, Alan↗

Compositional vs. thermal buoyancy and the evolution of subducted lithosphere

We formulate 2-D Cartesian finite element models that explore the fate of compositionally defined lithosphere as it encounters a viscosity increase at the boundary between the upper and lower mantle. Subducted lithosphere is represented as a cold, stiff, layered composite of denser eclogite underlain by more buoyant harzburgite. Slabs impinging on a lower mantle 30 and 100 times more viscous than the upper mantle thicken and fold strongly as they penetrate the lower mantle. Approximately a factor of two thickening occurs via pure shear just above the discontinuity, with additional enhancement due to folding by over a factor of two. No separation of the individual slab components occurs at the discontinuity, and direct comparison with models in which compositional buoyancy is explicitly ignored indicates that slab evolution is largely controlled by the thermal buoyancy. These results are at odds with hypotheses about slab evolution in which the compositional buoyancy contributions lead to component separation and the formation of slab megaliths or a compositionally layered upper mantle.

Gaherty, James B.↗

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗

Sensitivity Studies of Dust Ice Nuclei Effect on Cirrus Clouds with the Community Atmosphere Model CAM5

In this study the effect of dust aerosol on upper tropospheric cirrus clouds through heterogeneous ice nucleation is investigated in the Community Atmospheric Model version 5 (CAM5) with two ice nucleation parameterizations. Both parameterizations consider homogeneous and heterogeneous nucleation and the competition between the two mechanisms in cirrus clouds, but differ significantly in the number concentration of heterogeneous ice nuclei (IN) from dust. Heterogeneous nucleation on dust aerosol reduces the occurrence frequency of homogeneous nucleation and thus the ice crystal number concentration in the Northern Hemisphere (NH) cirrus clouds compared to simulations with pure homogeneous nucleation. Global and annual mean shortwave and longwave cloud forcing are reduced by up to 2.0+/-0.1Wm (sup−2) (1 uncertainty) and 2.4+/-0.1Wm (sup−2), respectively due to the presence of dust IN, with the net cloud forcing change of −0.40+/-0.20W m(sup−2). Comparison of model simulations with in situ aircraft data obtained in NH mid-latitudes suggests that homogeneous ice nucleation may play an important role in the ice nucleation at these regions with temperatures of 205-230 K. However, simulations overestimate observed ice crystal number concentrations in the tropical tropopause regions with temperatures of 190- 205 K, and overestimate the frequency of occurrence of high ice crystal number concentration (greater than 200 L(sup-1) and underestimate the frequency of low ice crystal number concentration (less than 30 L(sup-1) at NH mid-latitudes. These results highlight the importance of quantifying the number concentrations and properties of heterogeneous IN (including dust aerosol) in the upper troposphere from the global perspective.

Liu, Xiaohong↗

Interpretation of satellite airglow observations during the March 22, 1979, magnetic storm, using the coupled ionosphere-thermosphere model developed at University College, London

The University of California, Berkeley, extreme ultraviolet spectrometer aboard the U.S. Air Force STP 78-1 satellite measured emission features in the Earth's dayglow due to neutral and ionized species in the atmosphere, in the 35 to 140-nm range. The spectrometer was operating between March 1979 and March 1980, including the period of the magnetic storm on March 22, 1979. Some of these measurements are interpreted using the predictions of the three-dimensional time-dependent coupled ionosphere-thermosphere model developed at University College, London. The observations show a reduction in the atomic oxygen 130.4-nm airglow emission at high northern latitudes following the storm. Model simulations show that this reduction in 130.4-nm emission is associated with an increase in the O2/O ratio. Analysis of model results using electron transport and radiative transport codes show that the brightness of 130.4-nm emission at high latitudes due to resonantly scattered sunlight is approximately twice that due to photoelectron impact excitation. However, the observed decrease in the brightness at high northern latitudes is mainly due to a change in the photoelectron impact source, which contributes approximately 75% of the total, as well as its multiple scattering component; for the photoelectron impact source at 70 deg latitude and 200 km altitude, the reduction in multiple scattering is 1.5 times greater than the reduction in the initial excitation. The reduction in the airglow emission is visible only in the norther n hemisphere because the south pole was not sunlit over the storm period. The comparison of model results with observations suggests that 130.4-nm emission may be useful as a tracer for global changes in the concentration of atomic energy.

Parish, H. F.↗

A Method for Obtaining High Frequency, Global, IR-Based Convective Cloud Tops for Studies of the TTL

Models of varying complexity that simulate water vapor and clouds in the Tropical Tropopause Layer (TTL) show that including convection directly is essential to properly simulating the water vapor and cloud distribution. In boreal winter, for example, simulations without convection yield a water vapor distribution that is too uniform with longitude, as well as minimal cloud distributions. Two things are important for convective simulations. First, it is important to get the convective cloud top potential temperature correctly, since unrealistically high values (reaching above the cold point tropopause too frequently) will cause excessive hydration of the stratosphere. Second, one must capture the time variation as well, since hydration by convection depends on the local relative humidity (temperature), which has substantial variation on synoptic time scales in the TTL. This paper describes a method for obtaining high frequency (3-hourly) global convective cloud top distributions which can be used in trajectory models. The method uses rainfall thresholds, standard IR brightness temperatures, meteorological temperature analyses, and physically realistic and documented corrections IR brightness temperature corrections to derive cloud top altitudes and potential temperatures. The cloud top altitudes compare well with combined CLOUDSAT and CALIPSO data, both in time-averaged overall vertical and horizontal distributions and in individual cases (correlations of .65-.7). An important finding is that there is significant uncertainty (nearly .5 km) in evaluating the statistical distribution of convective cloud tops even using lidar. Deep convection whose tops are in regions of high relative humidity (such as much of the TTL), will cause clouds to form above the actual convection. It is often difficult to distinguish these clouds from the actual convective cloud due to the uncertainties of evaluating ice water content from lidar measurements. Comparison with models show that calculated cloud top altitudes are generally higher than those calculated by global analyses (e.g., MERRA). Interannual variability in the distribution of convective cloud top altitudes is also investigated.

hydration↗

Parameterization of planetary wave breaking in the middle atmosphere

A parameterization of planetary wave breaking in the middle atmosphere has been developed and tested in a numerical model which includes governing equations for a single wave and the zonal-mean state. The parameterization is based on the assumption that wave breaking represents a steady-state equilibrium between the flux of wave activity and its dissipation by nonlinear processes, and that the latter can be represented as linear damping of the primary wave. With this and the additional assumption that the effect of breaking is to prevent further amplitude growth, the required dissipation rate is readily obtained from the steady-state equation for wave activity; diffusivity coefficients then follow from the dissipation rate. The assumptions made in the derivation are equivalent to those commonly used in parameterizations for gravity wave breaking, but the formulation in terms of wave activity helps highlight the central role of the wave group velocity in determining the dissipation rate. Comparison of model results with nonlinear calculations of wave breaking and with diagnostic determinations of stratospheric diffusion coefficients reveals remarkably good agreement, and suggests that the parameterization could be useful for simulating inexpensively, but realistically, the effects of planetary wave transport.

Garcia, Rolando R.↗

Terrestrial Carbon Sinks in the Brazilian Amazon and Cerrado Region Predicted from MODIS Satellite Data and Ecosystem Modeling

A simulation model based on satellite observations of monthly vegetation cover from the Moderate Resolution Imaging Spectroradiometer (MODIS) was used to estimate monthly carbon fluxes in terrestrial ecosystems of Brazilian Amazon and Cerrado regions over the period 2000-2004. Net ecosystem production (NEP) flux for atmospheric CO2 in the region for these years was estimated. Consistently high carbon sink fluxes in terrestrial ecosystems on a yearly basis were found in the western portions of the states of Acre and Rondonia and the northern portions of the state of Par a. These areas were not significantly impacted by the 2002-2003 El Nino event in terms of net annual carbon gains. Areas of the region that show periodically high carbon source fluxes from terrestrial ecosystems to the atmosphere on yearly basis were found throughout the state of Maranhao and the southern portions of the state of Amazonas. As demonstrated though tower site comparisons, NEP modeled with monthly MODIS Enhanced Vegetation Index (EVI) inputs closely resembles the measured seasonal carbon fluxes at the LBA Tapajos tower site. Modeling results suggest that the capacity for use of MODIS Enhanced Vegetation Index (EVI) data to predict seasonal uptake rates of CO2 in Amazon forests and Cerrado woodlands is strong.

Potter, C.↗

Influence of Aerosol Heating on the Stratospheric Transport of the Mt. Pinatubo Eruption

On June 15th, 1991 the eruption of Mt. Pinatubo (15.1 deg. N, 120.3 Deg. E) in the Philippines injected about 20 Tg of sulfur dioxide in the stratosphere, which was transformed into sulfuric acid aerosol. The large perturbation of the background aerosol caused an increase in temperature in the lower stratosphere of 2-3 K. Even though stratospheric winds climatological]y tend to hinder the air mixing between the two hemispheres, observations have shown that a large part of the SO2 emitted by Mt. Pinatubo have been transported from the Northern to the Southern Hemisphere. We simulate the eruption of Mt. Pinatubo with the Goddard Earth Observing System (GEOS) version 5 global climate model, coupled to the aerosol module GOCART and the stratospheric chemistry module StratChem, to investigate the influence of the eruption of Mt. Pinatubo on the stratospheric transport pattern. We perform two ensembles of simulations: the first ensemble consists of runs without coupling between aerosol and radiation. In these simulations the plume of aerosols is treated as a passive tracer and the atmosphere is unperturbed. In the second ensemble of simulations aerosols and radiation are coupled. We show that the set of runs with interactive aerosol produces a larger cross-equatorial transport of the Pinatubo cloud. In our simulations the local heating perturbation caused by the sudden injection of volcanic aerosol changes the pattern of the stratospheric winds causing more intrusion of air from the Northern into the Southern Hemisphere. Furthermore, we perform simulations changing the injection height of the cloud, and study the transport of the plume resulting from the different scenarios. Comparisons of model results with SAGE II and AVHRR satellite observations will be shown.

Aquila, Valentina↗

Backscatter Modeling at 2.1 Micron Wavelength for Space-Based and Airborne Lidars Using Aerosol Physico-Chemical and Lidar Datasets

Space-based and airborne coherent Doppler lidars designed for measuring global tropospheric wind profiles in cloud-free air rely on backscatter, beta from aerosols acting as passive wind tracers. Aerosol beta distribution in the vertical can vary over as much as 5-6 orders of magnitude. Thus, the design of a wave length-specific, space-borne or airborne lidar must account for the magnitude of 8 in the region or features of interest. The SPAce Readiness Coherent Lidar Experiment under development by the National Aeronautics and Space Administration (NASA) and scheduled for launch on the Space Shuttle in 2001, will demonstrate wind measurements from space using a solid-state 2 micrometer coherent Doppler lidar. Consequently, there is a critical need to understand variability of aerosol beta at 2.1 micrometers, to evaluate signal detection under varying aerosol loading conditions. Although few direct measurements of beta at 2.1 micrometers exist, extensive datasets, including climatologies in widely-separated locations, do exist for other wavelengths based on CO2 and Nd:YAG lidars. Datasets also exist for the associated microphysical and chemical properties. An example of a multi-parametric dataset is that of the NASA GLObal Backscatter Experiment (GLOBE) in 1990 in which aerosol chemistry and size distributions were measured concurrently with multi-wavelength lidar backscatter observations. More recently, continuous-wave (CW) lidar backscatter measurements at mid-infrared wavelengths have been made during the Multicenter Airborne Coherent Atmospheric Wind Sensor (MACAWS) experiment in 1995. Using Lorenz-Mie theory, these datasets have been used to develop a method to convert lidar backscatter to the 2.1 micrometer wavelength. This paper presents comparison of modeled backscatter at wavelengths for which backscatter measurements exist including converted beta (sub 2.1).

Srivastava, V.↗

Measurements of soot formation and hydroxyl concentration in near critical equivalence ratio premixed ethylene flame

The testing and development of existing global and detailed chemical kinetic models for soot formation requires measurements of soot and radical concentrations in flames. A clearer understanding of soot particle inception relies upon the evaluation and refinement of these models in comparison with such measurements. We present measurements of soot formation and hydroxyl (OH) concentration in sequences of flat premixed atmospheric-pressure C2H4/O2/N2 flames and 80-torr C2H4/O2 flames for a unique range of equivalence ratios bracketting the critical equivalence ratio (phi(sub c)) and extending to more heavily sooting conditions. Soot volume fraction and number density profiles are measured using a laser scattering-extinction apparatus capable of resolving a 0.1 percent absorption. Hydroxyl number density profiles are measured using laser-induced fluorescence (LIF) with broadband detection. Temperature profiles are obtained from Rayleigh scattering measurements. The relative volume fraction and number density profiles of the richer sooting flames exhibit the expected trends in soot formation. In near-phi(sub c) visibility sooting flames, particle scattering and extinction are not detected, but an LIF signal due to polycyclic aromatic hydrocarbons (PAH's) can be detected upon excitation with an argon-ion laser. A linear correlation between the argon-ion LIF and the soot volume fraction implies a common mechanistic source for the growth of PAH's and soot particles. The peak OH number density in both the atmospheric and 80-torr flames declines with increasing equivalence ratio, but the profile shape remains unchanged in the transition to sooting, implying that the primary reaction pathways for OH remain unchanged over this transition. Chemical kinetic modeling is demonstrated by comparing predictions using two current reaction mechanisms with the atmospheric flame data. The measured and predicted OH number density profiles show good agreement. The predicted benzene number density profiles correlate with the measured trends in soot formation, although anomalies in the benzene profiles for the richer and cooler sooting flames suggest a need for the inclusion of benzene oxidation reactions.

Inbody, Michael Andrew↗

Assimilation of GOES Land Surface Data into a Mesoscale Models

A technique has been developed for assimilating Geostationary Operational Environmental Satellite (GOES)-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The assimilation technique has been applied to the Oklahoma-Kansas region during the spring-summer 2000 time period when dynamic changes in vegetation cover occur. In April, central Oklahoma is characterized by large NDVI associated with winter wheat while surrounding areas are primarily rangeland with lower NDVI. In July the vegetation pattern reverses as the central wheat area changes to low NDVI due to harvesting and the surrounding rangeland is greener than it was in April. The goal of this study is to determine if assimilating satellite land surface data can improve simulation of the complex spatial distribution of surface energy and water fluxes across this region. The PSU/NCAR NM5 V3 system is used in this study. The grid configuration consists of a 36-km CONUS domain and a 12-km nest over the area of interest. Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint indicates that assimilation of the satellite data results reduces both the bias and RMSE for both state variables. In addition, comparison of model data with ARM/CART EBBR flux observations reveals that the assimilation technique adjusts the bowen ratio in a realistic fashion.

Lapenta, William M.↗

On the Quality of the Nimbus 7 LIMS Version 6 Ozone for Studies of the Middle Atmosphere

The Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) radiance profile dataset of 1978/79 was reconditioned and reprocessed to Version 6 (V6) profiles of temperature and species that are improved significantly over those from Version 5 (V5). The LIMS V6 dataset was archived for public use in 2002. Improvements for its ozone include: (1) a more accurate accounting for instrument and spacecraft motion effects in the radiances, (2) the use of better spectroscopic line parameters for its ozone forward model, (3) retrievals of all its scans, (4) more accurate and compatible temperature versus pressure profiles (or T(p)), which are needed for the registration of the ozone radiances and for the removal of temperature effects from them, and (5) a better accounting for interfering species in the lower stratosphere. The retrieved V6 ozone profiles extend from near cloud top altitudes to about 80 km and from 64S to 84N latitude with better sampling along the orbit than for the V5 dataset. Calculated estimates of the single-profile precision and accuracy are provided for the V6 ozone from this study. Precision estimates based on the data themselves are of order 3% or better from 1 to 30 hPa. Estimates of total systematic error for a single profile are hard to generalize because the separate sources of error may not all be of the same sign and they depend somewhat on the atmospheric state. It is estimated that the V6 zonal mean ozone distributions are accurate to within 9% to 7% from 50 hPa to 3 hPa, respectively. Effects of a temperature bias can be significant and may be present at 1 to 2 hPa though. There may be ozone biases of order 10% at those levels due to possible biases of up to +2 K, but there is no indication of a similar problem elsewhere in the stratosphere. Simulation studies show that the LIMS retrievals are also underestimating slightly the small amplitudes of the atmospheric temperature tides, which affect its retrieved day/night ozone differences. There are small biases in the middle to lower stratosphere for the ascending versus descending node LIMS ozone, due principally to not accounting for the asymmetric weighting of its radiances across the tangent layer. The estimates of total accuracy were assessed by comparing the daily zonal mean LIMS ozone distributions against those from the Nimbus 7 SBUV Version 8 (V8) dataset for the same period. Generally, the LIMS V6 ozone agrees well with SBUV, except perhaps in the tropical lower stratosphere where the LIMS ozone is less. Still, the accuracy for LIMS V6 ozone in the lower stratosphere is improved over that found for LIMS V5, as indicated by several LIMS comparisons with ECC ozonesonde profiles. The LIMS V6 ozone is considered especially suitable for detailed studies of large-scale stratospheric processes above the 100-hPa level. Comparison of diurnal, photochemical model calculations with the monthly-averaged, upper stratospheric ozone observed from LIMS V6 indicates only a slight ozone deficit for the model at about 2 hPa. However, that deficit exhibits little to no seasonal variation and is in good agreement with similar model comparisons for a seasonal time series of ozone obtained with ground-based microwave instruments. Because the LIMS V6 ozone in the lower stratosphere has improved accuracy and sampling versus that of V5, it should now be possible to conduct quantitative studies of ozone transport and chemistry for the northern hemisphere, polar stratospheric winter of 1978/79 a time period when the catalytic loss of ozone due to reactive chlorine should not have been a major factor for the Arctic region.

Remsberg, Ellis↗

A seasonal model of the Saturnian upper troposphere Comparison with Voyager infrared measurements

In the present extension of Cess and Caldwell's (1979) seasonal climate model to the upper troposphere of Saturn, the ring-modulated latitudinal dependence of the insolation, the ring thermal emission, the oblateness and orbital eccentricity of the planet, and the latitudinal variation of the internal heat flux, are taken into account. While these calculations agree with the temperature-latitude profiles retrieved from Voyager IRIS measurements above 0.2-bar pressure level atmospheric strata, the model fails to predict the retrieved temperature-latitude profiles below the 0.3-bar level. This discrepancy may be due to the existence of clouds at these levels.

Bezard, B.↗

A comparison of the radio data and model calculations of Jupiter's synchrotron radiation. I - The high energy electron distribution in Jupiter's inner magnetosphere. II - East-west asymmetry in the radiation belts as a function of Jovian longitude

A comparison has been made between detailed model calculations of Jupiter's synchrotron radiation and the radio data at wavelengths of 6, 21, and 50 cm. The calculations were performed for a Jovian longitude of 200 deg and were based on the multipole field configurations as derived from the Pioneer data. The electron distribution in the inner magnetosphere was derived as a function of energy, pitch angle, and spatial coordinates. In addition, the hot region or east-west asymmetry in the radiation belts is investigated. It is suggested that this asymmetry is due to the combined effect of an overabundance of electrons at jovicentric longitudes of 240-360 deg and the existence of a dusk-to-dawn directed electric field over the inner magnetosphere generated by the wind system in the upper atmosphere.

De Pater, I.↗