Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Rendering Three-Dimensional Solar Coronal Structures

An X-ray or EUV image of the corona or chromosphere is a 2D representation of an extended 3D complex for which a general inversion process is impossible. A specific model must be incorporated in order to understand the full 3D structure. We approach this problem by modeling a set of optically-thin 3D plasma flux tubes which we render these as synthetic images. The resulting images allow the interpretation of the X-ray/EUV observations to obtain information on (1) the 3D structure of X-ray images, i.e., the geometric structure of the flux tubes, and on (2) the internal structure using specific plasma characteristics, i.e., the physical structure of the flux tubes. The data-analysis technique uses magnetograms to characterize photospheric magnetic fields and extrapolation techniques to form the field lines. Using a new set of software tools, we have generated 3D flux tube structures around these field lines and integrated the plasma emission along the line of sight to obtain a rendered image. A set of individual flux-tube images is selected by a non-negative least-squares technique to Provide a match with an observed X-ray image. The scheme minimizes the squares of the differences between the synthesized image and the observed image with a non-negative constraint on the coefficients of the brightness of the individual flux-tube loops. The derived images are used to determine the specific photospheric foot points and physical data, i.e., scaling laws for densities and loop lengths. The development has led to Computer efficient integration and display software that is compatible for comparison with observations (e.g., Yohkoh SXT data, NIXT, or EIT). This analysis is important in determining directly the magnetic field configuration, which provides the structure of coronal loops, and indirectly the electric currents or waves, which provide the energy for the heating of the plasma. We have used very simple assumptions (i.e., potential magnetic fields and isothermal corona) to provide an initial test of the techniques before complex models are introduced. We have separated the physical and geometric contributions of the emission for a set of flux tubes and concentrated, in this initial study, on the geometric contributions by making approximations to the physical contributions. The initial results are consistent with the scaling laws derived from Yohkoh SXT data.

Gary, G. Allen↗
NASA NTRS2022

Polarization as a Discriminator of Light-Absorbing Impurities in or Above Snow

This conceptual study presents advanced radiative transfer computations of light polarization originating from a snowpack consisting of nonspherical grains and variable content of light-absorbing impurities, either embedded in the snowpack or (with the same optical properties) lofted above it in the form of atmospheric aerosols. The results highlight the importance of considering shapes other than spherical for the snow grains, which otherwise can lead to non-negligible errors in the retrieval of snow albedo from remote sensing observations. More importantly, it is found that polarimetric measurements provide a means to partition light-absorbing impurities embedded in the snowpack from absorbing aerosols aloft, a task traditionally prohibitive for sensors capable exclusively of measurements of total reflectance. Heritage techniques to obtain snow grain size from shortwave infrared observations of total reflectance are well established, as are those that leverage polarimetric, multiangular observations across the entire optical spectrum to characterize the optical and microphysical properties of atmospheric aerosols. The polarization signatures of near-infrared (e.g., 864 nm) observations carry critical information on snow grain shape. The prospected launch of space-borne polarimeters with proven accuracy, therefore, advocates for the development of data inversion schemes, to boost the accuracy of simultaneous retrievals of atmospheric and surface parameters in the polar and snow-covered regions, critical to climate studies.

light polarization↗
NASA NTRS1995

Accurate orbit determination strategies for the tracking and data relay satellites

The National Aeronautics and Space Administration (NASA) has developed the Tracking and Data Relay Satellite (TDRS) System (TDRSS) for tracking and communications support of low Earth-orbiting satellites. TDRSS has the operational capability of providing 85% coverage for TDRSS-user spacecraft. TDRSS currently consists of five geosynchronous spacecraft and the White Sands Complex (WSC) at White Sands, New Mexico. The Bilateration Ranging Transponder System (BRTS) provides range and Doppler measurements for each TDRS. The ground-based BRTS transponders are tracked as if they were TDRSS-user spacecraft. Since the positions of the BRTS transponders are known, their radiometric tracking measurements can be used to provide a well-determined ephemeris for the TDRS spacecraft. For high-accuracy orbit determination of a TDRSS user, such as the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft, high-accuracy TDRS orbits are required. This paper reports on successive refinements in improved techniques and procedures leading to more accurate TDRS orbit determination strategies using the Goddard Trajectory Determination System (GTDS). These strategies range from the standard operational solution using only the BRTS tracking measurements to a sophisticated iterative process involving several successive simultaneous solutions for multiple TDRSs and a TDRSS-user spacecraft. Results are presented for GTDS-generated TDRS ephemerides produced in simultaneous solutions with the TOPEX/Poseidon spacecraft. Strategies with different user spacecraft, as well as schemes for recovering accurate TDRS orbits following a TDRS maneuver, are also presented. In addition, a comprehensive assessment and evaluation of alternative strategies for TDRS orbit determination, excluding BRTS tracking measurements, are presented.

Oza, D. H.↗
NASA NTRS2019

Observation of Absorbing Aerosols Above Clouds over the South-East Atlantic Ocean from the Geostationary Satellite SEVIRI – Part 1: Method Description and Sensitivity

High temporal resolution observations from satellites have a great potential for studying the impact of biomass burning aerosols and clouds over the South East Atlantic Ocean (SEAO). This paper presents a method developed to retrieve simultaneously aerosol and cloud properties in aerosol above cloud conditions from the geostationary instrument Meteosat Second Generation/Spinning Enhanced Visible and Infrared Imager (MSG/SEVIRI). The above-cloud Aerosol Optical Thickness (AOT), the Cloud Optical Thickness (COT) and the Cloud droplet Effective Radius (CER) are derived from the spectral contrast and the magnitude of the signal measured in three channels in the visible to shortwave infrared region. The impact of the absorption from atmospheric gases on the satellite signal is corrected by applying transmittances calculated using the water vapour profiles from a Met Office forecast model. The sensitivity analysis shows that a 10% error on the humidity profile leads to an 18.5% bias on the above-cloud AOT, which highlights the importance of an accurate atmospheric correction scheme. In situ measurements from the CLARIFY-2017 airborne field campaign are used to constrain the aerosol size distribution and refractive index that is assumed for the aforementioned retrieval algorithm. The sensitivities in the retrieved AOT, COT and CER to the aerosol model assumptions are assessed. Between 09:00-15:00 UTC, an uncertainty of 40% is estimated on the above-cloud AOT, which is dominated by the sensitivity of the retrieval to the single scattering albedo. The absorption AOT is less sensitive to the aerosol assumptions with an uncertainty generally lower than 17% between 09:00-15:00 UTC. Outside of that time range, as the scattering angle decreases, the sensitivity of the AOT and the absorption AOT to the aerosol model increases. The retrieved cloud properties are only weakly sensitive to the aerosol model assumptions throughout the day, with biases lower than 6% on the COT and 3% on the CER. The stability of the retrieval over time is analysed. For observations outside of the backscattering glory region, the time-series of the aerosol and cloud properties are physically consistent, which confirms the ability of the retrieval to monitor the temporal evolution of aerosol above cloud events over the SEAO.

Peers, Fanny↗
DOE OSTI2021

Climatological sensitivities of shallow-cumulus bulk entrainment in continental and oceanic locations

Cumulus entrainment is a complex process that has long challenged conceptual understanding and atmospheric prediction. To investigate this process observationally, two retrievals are used to generate multi-year climatologies of shallow-cumulus bulk entrainment ( ϵ ) at two Atmospheric Radiation Measurement cloud observatories, one in the US southern Great Plains (SGP) and the other in the Azores archipelago in the eastern North Atlantic (ENA). The statistical distributions of ϵ thus obtained, as well as certain environmental and cloud-related sensitivities of ϵ , are consistent with previous findings from large-eddy simulations. The retrieved ϵ robustly increases with cloudlayer relative humidity and decreases in wider clouds and cloud ensembles with larger cloud-base mass fluxes. While ϵ also correlates negatively with measures of cloud-layer vigor (e.g., maximum in-cloud vertical velocity and cloud depth), the extent to which these metrics actually regulate ϵ (or vice-versa) is unclear. Novel sensitivities of ϵ include a robust decrease of ϵ with increasing subcloud wind speed in oceanic flows, as well as a decrease of ϵ with increasing cloud-base mass flux in individual cumuli. A strong land–ocean contrast in ϵ is also found, with median values of 0.5-0.6 km -1 at the continental SGP site and and 1.0-1.1 km -1 at the oceanic ENA site. This trend is associated with drier and deeper cloud layers, along with larger cloud-base mass fluxes, at SGP, all of which favor reduced ϵ . The flow-dependence of retrieved ϵ implies that its various sensitivities should be accounted for in cumulus parameterization schemes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

North Atlantic Oscillation response in GeoMIP experiments G6solar and G6sulfur: why detailed modelling is needed for understanding regional implications of solar radiation management

The realization of the difficulty of limiting global-mean temperatures to within 1.5 or 2.0°C above pre-industrial levels stipulated by the 21st Conference of Parties in Paris has led to increased interest in solar radiation management (SRM) techniques. Proposed SRM schemes aim to increase planetary albedo to reflect more sunlight back to space and induce a cooling that acts to partially offset global warming. Under the auspices of the Geoengineering Model Intercomparison Project, we have performed model experiments whereby global temperature under the high-forcing SSP5-8.5 scenario is reduced to follow that of the medium-forcing SSP2-4.5 scenario. Two different mechanisms to achieve this are employed: the first via a reduction in the solar constant (experiment G6solar) and the second via modelling injections of sulfur dioxide (experiment G6sulfur) which forms sulfate aerosol in the stratosphere. Results from two state-of-the-art coupled Earth system models (UKESM1 and CESM2-WACCM6) both show an impact on the North Atlantic Oscillation (NAO) in G6sulfur but not in G6solar. Both models show a persistent positive anomaly in the NAO during the Northern Hemisphere winter season in G6sulfur, suggesting an increase in zonal flow and an increase in North Atlantic storm track activity impacting the Eurasian continent and leading to high-latitude warming over Europe and Asia. These results are broadly consistent with previous findings which show similar impacts from stratospheric volcanic aerosol on the NAO and emphasize that detailed modelling of geoengineering processes is required if accurate impacts of SRM effects are to be simulated. Differences remain between the two models in predicting regional changes over the continental USA and Africa, suggesting that more models need to perform such simulations before attempting to draw any conclusions regarding potential continental-scale climate change under SRM.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Gravity Wave Variances and Propagation Derived from AIRS Radiances

As the first gravity wave (GW) climatology study using nadir-viewing infrared sounders, 50 Atmospheric Infrared Sounder (AIRS) radiance channels are selected to estimate GW variances at pressure levels between 2-100 hPa. The GW variance for each scan in the cross-track direction is derived from radiance perturbations in the scan, independently of adjacent scans along the orbit. Since the scanning swaths are perpendicular to the satellite orbits, which are inclined meridionally at most latitudes, the zonal component of GW propagation can be inferred by differencing the variances derived between the westmost and the eastmost viewing angles. Consistent with previous GW studies using various satellite instruments, monthly mean AIRS variance shows large enhancements over meridionally oriented mountain ranges as well as some islands at winter hemisphere high latitudes. Enhanced wave activities are also found above tropical deep convective regions. GWs prefer to propagate westward above mountain ranges, and eastward above deep convection. AIRS 90 field-of-views (FOVs), ranging from +48 deg. to −48 deg. off nadir, can detect large-amplitude GWs with a phase velocity propagating preferentially at steep angles (e.g., those from orographic and convective sources). The annual cycle dominates the GW variances and the preferred propagation directions for all latitudes. Indication of a weak two-year variation in the tropics is found, which is presumably related to the Quasi-biennial oscillation (QBO). AIRS geometry makes its out-tracks capable of detecting GWs with vertical wavelengths substantially shorter than the thickness of instrument weighting functions. The novel discovery of AIRS capability of observing shallow inertia GWs will expand the potential of satellite GW remote sensing and provide further constraints on the GW drag parameterization schemes in the general circulation models (GCMs).

AIRS↗
NASA NTRS2019

Observation of absorbing aerosols above clouds over the south-east Atlantic Ocean from the geostationary satellite SEVIRI – Part 1: Method description and sensitivity

High-temporal-resolution observations from satellites have a great potential for studying the impact of biomass burning aerosols and clouds over the south-east Atlantic Ocean (SEAO). This paper presents a method developed to simultaneously retrieve aerosol and cloud properties in aerosol above-cloud conditions from the geostationary instrument Meteosat Second Generation/Spinning Enhanced Visible and Infrared Imager (MSG/SEVIRI). The above-cloud aerosol optical thickness (AOT), the cloud optical thickness (COT) and the cloud droplet effective radius (CER) are derived from the spectral contrast and the magnitude of the signal measured in three channels in the visible to shortwave infrared region. The impact of the absorption from atmospheric gases on the satellite signal is corrected by applying transmittances calculated using the water vapour profiles from a Met Office forecast model. The sensitivity analysis shows that a 10 % error on the humidity profile leads to an 18.5 % bias on the above-cloud AOT, which highlights the importance of an accurate atmospheric correction scheme. In situ measurements from the CLARIFY-2017 airborne field campaign are used to constrain the aerosol size distribution and refractive index that is assumed for the aforementioned retrieval algorithm. The sensitivities in the retrieved AOT, COT and CER to the aerosol model assumptions are assessed. Between 09:00 and 15:00 UTC, an uncertainty of 40 % is estimated on the above-cloud AOT, which is dominated by the sensitivity of the retrieval to the single-scattering albedo. The absorption AOT is less sensitive to the aerosol assumptions with an uncertainty generally lower than 17 % between 09:00 and 15:00 UTC. Outside of that time range, as the scattering angle decreases, the sensitivity of the AOT and the absorption AOT to the aerosol model increases. The retrieved cloud properties are only weakly sensitive to the aerosol model assumptions throughout the day, with biases lower than 6 % on the COT and 3 % on the CER. The stability of the retrieval over time is analysed. For observations outside of the backscattering glory region, the time series of the aerosol and cloud properties are physically consistent, which confirms the ability of the retrieval to monitor the temporal evolution of aerosol above-cloud events over the SEAO.

cloud properties↗
NASA NTRS2022

Unique NISTAR-Based Climate GCM Diagnostics of the Earth’s Planetary Albedo and Spectral Absorption Through Longitudinal Data Slicing

Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar and emitted thermal radiation permit a unique model/data comparison perspective that is not readily available from other satellite data. The key factor is the unique Lissajous orbital viewing geometry from the Lagrangian L1 point, which enables a continuous view of Earth’s sunlit hemisphere. The National Institute of Standards and Technology Advanced Radiometer (NISTAR) is the DSCOVR Mission energy budget instrument, which views the reflected and emitted radiation of the Earth’s sunlit hemisphere by means of single pixel active cavity full-spectrum (Band-A, 0.2–100 μm) and filtered solar wavelength (Band-B, 0.2–4.0 μm; and Band-C, 0.7–4.0 μm) radiometer measurements. An additional solar wavelength photodiode channel (0.3–1.1 μm) provides a calibration reference. The objective of this study is the assessment of climate GCM performance via direct model/data comparisons. Such comparisons are difficult due to quasi-chaotic natural variability present in real-world observational data and in climate GCM simulations. This is where the unique DSCOVR viewing geometry makes possible the longitudinal data slicing methodology for more direct model/data comparison. The key point of the longitudinal slicing approach is that data integration over the entire sunlit hemisphere eliminates the quasi-chaotic meteorological weather-scale noise, while preserving intra-seasonal and planetary-scale variability. The rotation of the Earth that retrieves this climate-style, large-scale longitudinal and seasonal variability. The hemispheric averaging is accomplished automatically in NISTAR measurements with its single-pixel view of the Earth. For climate GCMs, this requires implementing the Sunlit Hemisphere Sampling (SHS) scheme to operate on the GCM run-time output data, utilizing the DSCOVR Satellite Ephemeris data to assure precise viewing geometry between NISTAR measurements and GCM output data, while averaging out the meteorological weather noise. However, GCM generated data are radiative fluxes, while NISTAR (and EPIC) measurements are near-backscattered radiances. Conversing NISTSR measurements into radiative fluxes cannot be accomplished using NISTAR data alone, even with detailed support from conventional satellite data. But the identical viewing geometry of Earth’s sunlit hemisphere, and synergistic analyses of EPIC data make it feasible for this conversion of NISTAR near-backscatter radiances into radiative fluxes.

deep space climate observatory↗
DOE OSTI2024

Ecological Trait-Based Digital Categorization of Microbial Genomes for Denitrification Potential

Microorganisms encode proteins that function in the transformations of useful and harmful nitrogenous compounds in the global nitrogen cycle. The major transformations in the nitrogen cycle are nitrogen fixation, nitrification, denitrification, anaerobic ammonium oxidation, and ammonification. The focus of this report is the complex biogeochemical process of denitrification, which, in the complete form, consists of a series of four enzyme-catalyzed reduction reactions that transforms nitrate to nitrogen gas. Denitrification is a microbial strain-level ecological trait (characteristic), and denitrification potential (functional performance) can be inferred from trait rules that rely on the presence or absence of genes for denitrifying enzymes in microbial genomes. Despite the global significance of denitrification and associated large-scale genomic and scholarly data sources, there is lack of datasets and interactive computational tools for investigating microbial genomes according to denitrification trait rules. Therefore, our goal is to categorize archaeal and bacterial genomes by denitrification potential based on denitrification traits defined by rules of enzyme involvement in the denitrification reduction steps. We report the integration of datasets on genome, taxonomic lineage, ecosystem, and denitrifying enzymes to provide data investigations context for the denitrification potential of microbial strains. We constructed an ecosystem and taxonomic annotated denitrification potential dataset of 62,624 microbial genomes (866 archaea and 61,758 bacteria) that encode at least one of the twelve denitrifying enzymes in the four-step canonical denitrification pathway. Our four-digit binary-coding scheme categorized the microbial genomes to one of sixteen denitrification traits including complete denitrification traits assigned to 3280 genomes from 260 bacteria genera. The bacterial strains with complete denitrification potential pattern included Arcobacteraceae strains isolated or detected in diverse ecosystems including aquatic, human, plant, and Mollusca (shellfish). The dataset on microbial denitrification potential and associated interactive data investigations tools can serve as research resources for understanding the biochemical, molecular, and physiological aspects of microbial denitrification, among others. The microbial denitrification data resources produced in our research can also be useful for identifying microbial strains for synthetic denitrifying communities.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1996

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. We consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise voracity is generated and ejected due to the oscillatory motion of the lid. The present simulations c Implement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University. Jacobson and Reynolds used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds. More recently, Lachowiez and Wlezien are investigating the flow generated by an electro-mechanically driven lid to be used for assertion control in aerodynamic applications. We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin. The viscous effects are taken into account by modifying the strength of the particles, whereas fast multipole schemes employing hundreds of thousands ol'particle's allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies. Estimates of the forces needed to drive the devices will also be presented.

Mansour, N. N.↗
DOE OSTI2024

Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids

This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.

graphics processing units↗
NASA NTRS1993

Post-metamorphic brecciation in type 3 ordinary chondrites

Type 3.1-3.9 ordinary chondrites can be divided into two kinds: those in which the compositions of chondrule silicates are entirely consistent with metamorphism of type 3.0 material, and those in which the computational heterogeneity appears to be too extreme for in situ metamorphism. We present petrologic data for three LL3 chondrites of the second kind--Ngawi, ALH A77278 (both type 3.6), and Hamlet (type 3.9)--and compare these data with results for the first kind of LL3-4 chondrites. Given that chondrules form in the nebula and that metamorphic equilibration occurs in asteroids, our new data imply that Ngawi, A77278, Hamlet, and many other type 3 ordinary chondrites are post-metamorphic breccias containing materials with diverse metamorphic histories; they are not metamorphic rocks or special kinds of 'primitive breccias.' We infer also that metamorphism to type 3.1-3.9 levels produces very friable material that is easily remixed into breccias and lithified by mild shock. Thus, petrologic types and subtypes of chondrites indicate the mean metamorphic history of the ingredients, not the thermal history of the rock. The metamorphic history of individual type 1 or 2 porphyritic chondrules in type 3 breccias is best derived from olivine and pyroxene analyses and the data of McCoy et al. for unbrecciated chondrites. The new chondrule classification schemes of Sears, DeHart et al., appears to provide less information about the original state and metamorphic history of individual porphyritic chondrules and should not replace existing classification schemes.

Scott, E. R. D.↗
NASA NTRS2007

Advanced NASA Earth Science Mission Concept for Vegetation 3D Structure, Biomass and Disturbance

Carbon in forest canopies represents about 85% of the total carbon in the Earth's aboveground biomass (Olson et al., 1983). A major source of uncertainty in global carbon budgets derives from large errors in the current estimates of these carbon stocks (IPCC, 2001). The magnitudes and distributions of terrestrial carbon storage along with changes in sources and sinks for atmospheric C02 due to land use change remain the most significant uncertainties in Earth's carbon budget. These uncertainties severely limit accurate terrestrial carbon accounting; our ability to evaluate terrestrial carbon management schemes; and the veracity of atmospheric C02 projections in response to further fossil fuel combustion and other human activities. Measurements of vegetation three-dimensional (3D) structural characteristics over the Earth's land surface are needed to estimate biomass and carbon stocks and to quantify biomass recovery following disturbance. These measurements include vegetation height, the vertical profile of canopy elements (i.e., leaves, stems, branches), andlor the volume scattering of canopy elements. They are critical for reducing uncertainties in the global carbon budget. Disturbance by natural phenomena, such as fire or wind, as well as by human activities, such as forest harvest, and subsequent recovery, complicate the quantification of carbon storage and release. The resulting spatial and temporal heterogeneity of terrestrial biomass and carbon in vegetation make it very difficult to estimate terrestrial carbon stocks and quantify their dynamics. Vegetation height profiles and disturbance recovery patterns are also required to assess ecosystem health and characterize habitat. The three-dimensional structure of vegetation provides habitats for many species and is a control on biodiversity. Canopy height and structure influence habitat use and specialization, two fundamental processes that modify species richness and abundance across ecosystems. Accurate and consistent 3D measurements of forest structure at the landscape scale are needed for assessing impacts to animal habitats and biodiversity following disturbance.

Ranson, K. Jon↗
NASA NTRS1998

Simulation of Long Lived Tracers Using an Improved Empirically Based Two-Dimensional Model Transport Algorithm

We have developed a new empirically-based transport algorithm for use in our GSFC two-dimensional transport and chemistry model. The new algorithm contains planetary wave statistics, and parameterizations to account for the effects due to gravity waves and equatorial Kelvin waves. As such, this scheme utilizes significantly more information compared to our previous algorithm which was based only on zonal mean temperatures and heating rates. The new model transport captures much of the qualitative structure and seasonal variability observed in long lived tracers, such as: isolation of the tropics and the southern hemisphere winter polar vortex; the well mixed surf-zone region of the winter sub-tropics and mid-latitudes; the latitudinal and seasonal variations of total ozone; and the seasonal variations of mesospheric H2O. The model also indicates a double peaked structure in methane associated with the semiannual oscillation in the tropical upper stratosphere. This feature is similar in phase but is significantly weaker in amplitude compared to the observations. The model simulations of carbon-14 and strontium-90 are in good agreement with observations, both in simulating the peak in mixing ratio at 20-25 km, and the decrease with altitude in mixing ratio above 25 km. We also find mostly good agreement between modeled and observed age of air determined from SF6 outside of the northern hemisphere polar vortex. However, observations inside the vortex reveal significantly older air compared to the model. This is consistent with the model deficiencies in simulating CH4 in the northern hemisphere winter high latitudes and illustrates the limitations of the current climatological zonal mean model formulation. The propagation of seasonal signals in water vapor and CO2 in the lower stratosphere showed general agreement in phase, and the model qualitatively captured the observed amplitude decrease in CO2 from the tropics to midlatitudes. However, the simulated seasonal amplitudes were attenuated too rapidly with altitude in the tropics. Overall, the simulations with the new transport formulation are in substantially better agreement with observations compared with our previous model transport.

Fleming, E. L.↗
NASA NTRS1995

Delayed detonation models for normal and subluminous type Ia sueprnovae: Absolute brightness, light curves, and molecule formation

We compute optical and infrared light curves of the pulsating class of delayed detonation models for Type Ia supernovae (SN Ia's) using an elaborate treatment of the Local Thermodynamic Equilbrium (LTE) radiation transport, equation of state and ionization balance, expansion opacity including the cooling by CO, Co(+), and SiO, and a Monte Carlo gamma-ray deposition scheme. The models have an amount of Ni-56 in the range from approximately or equal to 0.1 solar mass up to 0.7 solar mass depending on the density at which the transition from a deflagration to a detonation occurs. Models with a large nickel production give light curves comparable to those of typical Type Ia supernovae. Subluminous supernovae can be explained by models with a low nickel production. Multiband light curves are presented in comparison with the normally bright event SN 1992bc and the subluminous events Sn 1991bg and SN 1992bo to establish the principle that the delayed detonation paradigm in Chandrasekhar mass models may give a common explosion mechanism accounting for both normal and subluminous SN Ia's. Secondary IR-maxima are formed in the models of normal SN Ia's as a photospheric effect if the photospheric radius continues to increase well after maximum light. Secondary maxima appear later and stronger in models with moderate expansion velocities and with radioactive material closer to the surface. Model light curves for subluminous SN Ia's tend to show only one 'late' IR-maximum. In some delayed detonation models shell-like envelopes form, which consist of unburned carbon and oxygen. The formation of molecules in these envelopes is addressed. If the model retains a C/O-envelope and is subluminous, strong vibration bands of CO may appear, typically several weeks past maximum light. CO should be very weak or absent in normal Sn Ia's.

Hoflich, P.↗
DOE OSTI2025

Determination of Scale Bar for AGHCF Metallography Data

The DOE Nuclear Energy Advanced Reactor Technologies (ART) Fast Reactor Program (FRP) has supported the development of several databases containing information on the safety performance of fast reactors, components, and fuels. This growing collection of legacy experimental data, operating data, and analysis is available online to registered users. Metallography data represents one of the most critical types of post-irradiation examination (PIE) data being collected, organized, and archived in several ART Fast Reactor Databases (https://frdb.ne.anl.gov), including the Metallic Fuels Irradiation & Physics Database (FIPD), Out-of-Pile Transient Database (OPTD), and TREAT (the Transient Reactor Test Facility) Experimental Relational Database (TREXR). These databases contain three principal sets of metallography data. The first set comprises metallography data from Experimental Breeder Reactor-II (EBR-II) and Fast Flux Test Facility (FFTF) irradiated fuel pins examined in the Hot Fuel Examination Facility (HFEF). The second set consists of metallography data from EBR-II irradiated fuel pins examined in the Alpha-Gamma Hot Cell Facility (AGHCF). The third set includes metallography data from transient-tested fuel pins (including both out-of-pile furnace tests and TREAT tests) examined in AGHCF. Since both the second and third sets were generated in AGHCF, they are governed by identical specifications. The metallography data in the databases consist of digital images scanned from either positive or negative photographic films. To analyze the microstructure of a fuel pin, a series of preparatory steps are required, including sectioning, epoxy mounting, mechanical grinding and polishing, and etching. Following sample preparation, specimens are transferred for metallographic examination. The AGHCF and HFEF metallography data were generated using optical microscopes manufactured by Leitz and Bausch and Lomb (B&L). Images were recorded on Polaroid film at preset magnifications. Magnification verification for the Leitz and B&L metallographs was conducted every two months prior to 1989 and at least every six months from 1989 through the conclusion of the IFR program. Magnifications determined from imaging of microslide standards were compared to the instrument settings for magnifications ranging from 50× to 500×. If the magnifications determined from standards deviated from the instrument settings, adjustments were made to the bellows extension until agreement was achieved. The specifications for AGHCF and HFEF legacy metallography data have been established based on available hard-copy and digital records, most of which have been incorporated into the data repositories associated with FIPD, OPTD, and TREXR. Detailed specifications including hard-copy records, digitized records, cutting diagrams and sectioning schemes, high-magnification photographs, photomosaics (composites), information tags, scale bars, and magnification verification procedures can be found in a separate report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2000

GCSS Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison Project, a project of GCSS Working Group on Cirrus Cloud Systems (WG2), involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. The goal of this project is to document and understand the factors resulting in significant inter-model differences. The intent is to foment research leading to model improvement and validation. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 cm/s, respectively. Five models participated. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation formulation, aerosol size, ice crystal growth rate (particularly the deposition coefficient), and water vapor uptake rate are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory data. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 cm/s) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, which is sensitive to the deposition coefficient when ice particles are small, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the nucleated (heterogeneous) ice crystals limits the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗