Search NASA⌕ Search

SEARCH · Search NASA

Results for “Estimation method”

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,621 records · Page 90

Methodologies to determine forces on bones and muscles of body segments during exercise, employing compact sensors suitable for use in crowded space vehicles

Work under this grant was carried out by the author and by a graduate research assistant. An instrumented bicycle ergometer was implemented focusing on the stated objective: to estimate the forces exerted by each muscle of the feet, calf, and thigh of an individual while bicycling. The sensors used were light and compact. These were probes to measure muscle EMG activity, miniature accelerometers, miniature load sensors, and small encoders to measure angular positions of the pedal. A methodology was developed and implemented to completely describe the kinematics of the limbs using data from the sensors. This work has been published as a Master's Thesis by the Graduate student supported by the grant. The instrumented ergometer along with the sensors and instrumentation were tested during a KC-135 Zero-Gravity flight in July, 1994. A complete description of the system and the tests performed have been published as a report submitted to NASA Johnson Space Center. The data collected during the KC-135 flight is currently being processed so that a kinematic description of the bicycling experiment will be soon determined. A methodology to estimate the muscle forces has been formulated based on previous work. The methodology involves the use of optimization concepts so that the individual muscle forces that represent variables in dynamic equations of motion may be estimated. Optimization of a criteria (goal) function such as minimization of energy will be used along with constraint equations defined by rigid body equations of motion. Use of optimization principles is necessary, because the equations of motion alone constitute an indeterminate system of equations with respect to the large amount of muscle forces which constitute the variables in these equations. The number of variables is reduced somewhat by using forces measured by the load cells installed on the pedal. These load cells measure pressure and shear forces on the foot. The author and his collaborators at NASA and at the University of Alabama, Tuscaloosa, are continuing the work of reducing the experimental data from the KC-135 flight, and the implementation of the optimization methods to estimate muscle forces. As soon as results from these efforts are available, they will be published in reputable journals. Results of this work will impact studies addressing bone density loss and development of countermeasures to minimize bone loss in zero gravity conditions. By analyzing muscle forces on Earth and in Space during exercise, scientists could eventually formulate new exercises and machines to help maintain bone density. On Earth, this work will impact studies concerning arthritis, and will provide the means to study possible exercise countermeasures to minimize arthritis problems.

Figueroa, Fernando↗

Numerical Study and Validation of Melting and Solidification in PCM Embedded Heat Exchangers with Straight Tube

Latent heat thermal energy storage (LHTES) systems have shown great potential to enable reliable use of renewable energy and load shifting. LHTES offer high storage density and release energy at near constant temperature because of its use of phase change materials (PCMs). The cylindrical PCM heat exchangers (PCMHX) are one of the most used technologies due to their simplicity. Numerical models for such PCMHX enable engineers to estimate their performance for different design parameters and operating conditions without having to test them all. However, modeling the phase change phenomena can be challenging. To better understand the difficulties involving accurate modeling of PCMHX, a cylindrical latent storage unit filled with PCM and water as in-tube heat transfer fluid (HTF) is numerically investigated. This paper presents a study based on a 2D-axisymmetric model of a straight tube embedded in PCM in a cylindrical container. CFD is used to study the charging (melting) and discharging (solidification) phenomena. The models are validated against experimental and numerical data from the literature. The predicted local PCM temperature profile over time agrees within 2K compared to the experimental values. The paper also presents a simple method to estimate the melting and solidification phase change temperature range from limited data provided by PCM manufacturers.

PCM, melting, solidification, CFD, heat exchanger↗

Exploring Telematics Big Data for Truck Platooning Opportunities

NREL completed a temporal and geospatial analysis of telematics data to estimate the fraction of platoonable miles traveled by class 8 tractor trailers currently in operation. This paper discusses the value and limitations of very large but low time-resolution data sets, and the fuel consumption reduction opportunities from large scale adoption of platooning technology for class 8 highway vehicles in the US based on telematics data. The telematics data set consist of about 57,000 unique vehicles traveling over 210 million miles combined during a two-week period. 75% of the total fuel consumption result from vehicles operating in top gear, suggesting heavy highway utilization. The data is at a one-hour resolution, resulting in a significant fraction of data be uncategorizable, yet significant value can still be extracted from the remaining data. Multiple analysis methods to estimate platoonable miles are discussed. Results indicate that 63% of total miles driven at known hourly-average speeds happens at speeds amenable to platooning. When also considering availability of nearby partner vehicles, results indicate 55.7% of all classifiable miles driven were platoonable. Analysis also address the availability of numerous partners enabling platoons greater than 2 trucks and the percentage of trucks that would be required to be equipped with platooning equipment to realize more than 50% of the possible savings.

33 ADVANCED PROPULSION SYSTEMS↗

Monte Carlo control loops for cosmic shear cosmology with DES Year 1 data

Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. Here, we present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum $C_ℓ$ and the redshift distribution $n(z)$ of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of $C_ℓ$ and $n(z)$ curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of $Ω_m$, $σ_8$, and intrinsic alignment amplitude $A_{IA}$, defined as $S_8D_{IA}=σ_8(Ω_m/0.3)^{0.5}D_{IA}$, where $D_{IA}=1-0.11(A_{IA}-1)$. We find $S_8D_{IA}=0.8954_{-0.039}^{+0.054}$, where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Air Pathway Dose Modeling for the E-Area Low-Level Waste Facility

The US Department of Energy (DOE) Order 435.1 performance assessment (PA) process (USDOE 1999) prescribes a performance objective (10 mrem y -1 ) for evaluating atmospheric releases of radionuclides from DOE low-level waste (LLW) disposal facilities. The potential dose to an individual from exposure to radionuclides released into the atmosphere from LLW disposals can be estimated by application of radionuclide -specific dose-release factors (DRFs) (mrem Ci -1 ) to estimated flux rates (Ci y -1 ) at a particular time and location. In Revision 0 of this report, Lee (2006) calculated DRFs for potential atmospheric releases of 15 volatile radionuclides from seven ELLWF disposal units in support of the 2008 ELLWF PA (WSRC 2008). The mainframe version of the US Environmental Protection Agency’s (EPA) dose model CAP88 was used in the Revision 0 assessment. In Revision 1 of this report (Dixon and Minter 2017), the site-specific Savannah River National Laboratory (SRNL) atmospheric dose models MAXDOSE-SR version 2013 (Stone and Jannik 2013a) and MAXINE version 2017 (Bell 2017) were used to calculate DRFs for ten radionuclides. This work (Revision 2) supersedes Revision 0 (Lee 2006) and Revision 1 (Dixon and Minter 2017) of this report. The Revision 2 method for estimating new DRFs for the ELLWF disposal units is described in this report.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimation of Rotary Stability Derivatives at Subsonic and Transonic Speeds

The first part of this paper pertains to the estimation of subsonic rotary stability derivatives of wings. The unsteady potential flow problem is solved by a superposition of steady flow solutions. Numerical results for the damping coefficients of triangular wings are presented as functions of aspect ratio and Mach number, and are compared with experimental results over the Mach number range 0 to 1. In the second part, experimental results are used. to point out a close correlation between the nonlinear variations with angle of attack of the static pitching-moment curve slope and the damping-in-pitch coefficient. The underlying basis for the correlation is found as a result of an analysis in which the indicial function concept and. the principle of super-position are adapted to apply to the nonlinear problem. The form of the result suggests a method of estimating nonlinear damping coefficients from results of static wind-tunnel measurements.

Tobak, Murray↗

Applying Gaussian Process Machine Learning and Modern Probabilistic Programming to Satellite Data to Infer CO 2 Emissions

Satellite data provides essential insights into the spatiotemporal distribution of CO 2 concentrations. However, many atmospheric inverse models fail to adequately incorporate the spatial and temporal correlations inherent in satellite observations and often lack rigorous methods for estimating parameters like spatial length scales. We introduce an inference model that processes the spatiotemporal covariance in satellite data and estimates hyperparameters such as covariance length scales. Our approach uses the Gaussian process (GP) machine learning (ML) and modern probabilistic programming languages (PPLs) to perform atmospheric inversions of emissions from satellite data. We develop a GP ML inversion system based on modern PPLs and the GEOS-Chem chemical transport model, simulating atmospheric CO 2 concentrations corresponding to the Orbiting Carbon Observatory-2/3 (OCO-2/3) data for July 2020. In our supervised learning framework, we treat the GEOS-Chem simulated data set as the target, with predictors derived by scaling the target with sector-specific factors hidden from the GP machine. Our results show that the GP model, combined with GPU-enabled PPLs, effectively retrieves true emission scaling factors and infers noise levels concealed within the data. This suggests that our method could be applied over larger areas with more complex covariance structures, enabling comprehensive analysis of the spatiotemporal patterns observed in OCO-2/3 and similar satellite data sets.

54 ENVIRONMENTAL SCIENCES↗

Water use in a young Pinus taeda bioenergy plantation: Effect of intensive management on stand evapotranspiration

Abstract The increasing demand for plant‐derived bioenergy is projected to expand tree plantations with intensive silviculture and improved tree genetics. These silvicultural practices result in faster stand development and canopy closure, which may also influence the systems' water dynamics. Here, we studied the evapotranspiration (ET) of a young (5 years old) intensively managed loblolly pine ( Pinus taeda ) stand and investigated the components of ET to determine its contribution to overall water use. We also compared ET with plantations that received less intensive management to determine whether our stand used more water. We used the eddy covariance method to estimate ecosystem‐level total ET (ET EC ), while plot‐level estimates of ET (ET P ) were obtained via soil lysimeters, sap flow sensors, and throughfall collectors, enabling measurement of the components of ET. Soil evaporation (Es) was the largest component of ET P (36%) over the course of the study, while transpiration and canopy interception accounted for 27% and 22%, respectively. Es decreased with stand development, while transpiration and canopy interception increased. Leaf area index (LAI) and precipitation were the most significant factors controlling ET and its components. Compared to previous studies in different sites that have similar age but lower LAI, our stand had higher water use. This high water use in the early stages of stand development was primarily due to high Es before the canopy was fully developed. While there are potential sources of uncertainty when comparing ET EC and the component fluxes in ET P , results from the two methods were not significantly different. This study had the advantage of using multiple methods to understand and verify the component processes that contribute to ET. Therefore, we recommend that multiple measurement techniques be used in the long‐term observation of ET, and in particular for the evaluation of the impact that intensively managed forests have on water resources in the southeastern United States.

54 ENVIRONMENTAL SCIENCES↗

How Noninnocent Spectator Species Improve the Oxygen Reduction Activity of Single-Atom Catalysts: Microkinetic Models from First-Principles Calculations

Graphene-based single-atom catalysts are promising alternatives to platinum-based catalysts for fuel cell applications. Different transition metals have been screened using electronic structure methods by estimating onset potentials from the most endergonic elementary reaction step. Here, we calculate onset potentials for the oxygen reduction reaction on metal atoms embedded in N-substituted graphene di-vacancies by virtue of first-principles-informed microkinetic analysis. We find that for more oxophilic metals (Cr, Fe, Mn, and Ru), purely thermodynamic models systematically underestimate onset potentials. Furthermore, the oxophilic metals (Cr, Fe, Mn, and Ru) are oxidized under reaction conditions, leading to an increase in activity compared to their reduced state. Importantly, coadsorbed O m H n species actively participate in the reaction, which requires a dynamic treatment of spectator species. These findings highlight the limitations of thermodynamic analyses for electrocatalytic processes, which commonly assume the same oxidation state for each metal, and show that deviations between computational and experimental onset potentials cannot be solely attributed to the shortcomings of the electronic structure methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of Model Based Parameter Estimation for RCS Frequency Response Calculations Using Method of Moments

An implementation of the Model Based Parameter Estimation (MBPE) technique is presented for obtaining the frequency response of the Radar Cross Section (RCS) of arbitrarily shaped, three-dimensional perfect electric conductor (PEC) bodies. An Electric Field Integral Equation (EFTE) is solved using the Method of Moments (MoM) to compute the RCS. The electric current is expanded in a rational function and the coefficients of the rational function are obtained using the frequency derivatives of the EFIE. Using the rational function, the electric current on the PEC body is obtained over a frequency band. Using the electric current at different frequencies, RCS of the PEC body is obtained over a wide frequency band. Numerical results for a square plate, a cube, and a sphere are presented over a bandwidth. Good agreement between MBPE and the exact solution over the bandwidth is observed.

Reddy, C. J.↗

Information depth in backscattered electron microscopy of nanoparticles within a solid matrix

Highlights: • Backscattered electron microscopy was used for reliable determination of nanoparticle density and sizes within Al alloys • Backscattered electron information depth was determined as a function of nanoparticle size from Monte Carlo simulations • Unlike transmission electron microscopy, the method is non-destructive and allows analyzing large sample areas • Good agreement with transmission electron microscopy data shows that our method is suitable for nanocomposite materials - Abstract: Measuring the dimensions and number density of nanoparticles dispersed in a solid matrix is usually accomplished via transmission electron microscopy (TEM) which suffers from high cost, low throughput, and small analytical volume. In comparison, scanning backscattered electron microscopy is inexpensive, requires little sample preparation, and allows for the analysis of large sample areas. However, the information depth is usually not known precisely and depends on several factors such as the composition of the nanoparticles and the matrix as well as the size of the nanoparticles, hindering the reconstruction of the actual size distribution and three-dimensional number density. Here we present a method to estimate the information depth for spherical nanoparticles of different sizes in order to accurately determine size distribution and number density. The approach is based on Monte Carlo simulation of electron trajectories in the material and analysis of the obtained backscattered electron signal-to-noise-ratio. Our experimental results are compared to those obtained via TEM and good agreement is demonstrated; this shows that TEM can be replaced by scanning electron microscopy for studying nanocomposites in many cases.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Scintillation and cherenkov photon counting detectors with analog silicon photomultipliers for TOF-PET

Objective. Standard signal processing approaches for scintillation detectors in positron emission tomography (PET) derive accurate estimates for 511 keV photon time of interaction and energy imparted to the detection media from aggregate characteristics of electronic pulse shapes. The ultimate realization of a scintillation detector for PET is one that provides a unique timestamp and position for each detected scintillation photon. Detectors with these capabilities enable advanced concepts for three-dimensional (3D) position and time of interaction estimation with methods that exploit the spatiotemporal arrival time kinetics of individual scintillation photons. Approach. In this work, we show that taking into consideration the temporal photon emission density of a scintillator, the channel density of an analog silicon photomultiplier (SiPM) array, and employing fast electronic readout with digital signal processing, a detector that counts and timestamps scintillation photons can be realized. To demonstrate this approach, a prototype detector was constructed, comprising multichannel electronic readout for a bismuth germanate (BGO) scintillator coupled to an SiPM array. Main Results. In proof-of-concept measurements with this detector, we were able to count and provide unique timestamps for 66% of all optical photons, where the remaining 34% (two-or-more-photon pulses) are also independently counted, but each photon bunch shares a common timestamp. We show this detector concept can implement 3D positioning of 511 keV photon interactions and thereby enable corrections for time of interaction estimators. The detector achieved 17.6% energy resolution at 511 keV and 237 ± 10 ps full-width-at-half-maximum coincidence time resolution (CTR) (fast spectral component) versus a reference detector. We outline the methodology, readout, and approach for achieving this detector capability in first-ever, proof-of-concept measurements for scintillation photon counting detector with analog silicon photomultipliers. Significance. The presented detector concept is a promising design for large area, high sensitivity TOF-PET detector modules that can implement advanced event positioning and time of interaction estimators, which could push state-of-the-art performance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Exploiting Measurement Uncertainty Estimation in Evaluation of GOES-R ABI Image Navigation Accuracy Using Image Registration Techniques

In evaluating GOES-R Advanced Baseline Imager (ABI) image navigation quality, upsampled sub-images of ABI images are translated against downsampled Landsat 8 images of localized, high contrast earth scenes to determine the translations in the East-West and North-South directions that provide maximum correlation. The native Landsat resolution is much finer than that of ABI, and Landsat navigation accuracy is much better than ABI required navigation accuracy and expected performance. Therefore, Landsat images are considered to provide ground truth for comparison with ABI images, and the translations of ABI sub-images that produce maximum correlation with Landsat localized images are interpreted as ABI navigation errors. The measured local navigation errors from registration of numerous sub-images with the Landsat images are averaged to provide a statistically reliable measurement of the overall navigation error of the ABI image. The dispersion of the local navigation errors is also of great interest, since ABI navigation requirements are specified as bounds on the 99.73rd percentile of the magnitudes of per pixel navigation errors. However, the measurement uncertainty inherent in the use of image registration techniques tends to broaden the dispersion in measured local navigation errors, masking the true navigation performance of the ABI system. We have devised a novel and simple method for estimating the magnitude of the measurement uncertainty in registration error for any pair of images of the same earth scene. We use these measurement uncertainty estimates to filter out the higher quality measurements of local navigation error for inclusion in statistics. In so doing, we substantially reduce the dispersion in measured local navigation errors, thereby better approximating the true navigation performance of the ABI system.

navigation performance↗

Air Pathway Dose Modeling for the E-Area Low-Level Waste Facility

The US Department of Energy (DOE) Order 435.1 performance assessment (PA) process (USDOE 1999) prescribes a performance objective (10 mrem y -1 ) for evaluating atmospheric releases of radionuclides from DOE low-level waste (LLW) disposal facilities. The potential dose to an individual from exposure to radionuclides released into the atmosphere from LLW disposals can be estimated by application of radionuclide-specific dose-release factors (DRFs) (mrem Ci -1 ) to estimated flux rates (Ci y -1 ) at a particular time and location. In Revision 0 of this report, Lee (2006) calculated DRFs for potential atmospheric releases of 15 volatile radionuclides from seven ELLWF disposal units in support of the 2008 ELLWF PA (WSRC 2008). The mainframe version of the US Environmental Protection Agency’s (EPA) dose model CAP88 was used in the Revision 0 assessment. In Revision 1 of this report (Dixon and Minter 2017), the site-specific Savannah River National Laboratory (SRNL) atmospheric dose models MAXDOSE-SR version 2013 (Stone and Jannik 2013a) and MAXINE version 2017 (Bell 2017) were used to calculate DRFs for ten radionuclides. Revision 2 provided updated DRFs using the current version of the US Environmental Protection Agency’s (EPA) dose model CAP88-PC Version 4.1. After Revision 2 was finalized, DRFs were requested for Ar-37, Ar-39, Kr-83m, and Hg-206. Revision 3 provides the additional DRFs calculated using the same methods as used in Revision 2 of this document. The method for estimating new DRFs for the ELLWF disposal units is described in this report.

61 RADIATION PROTECTION AND DOSIMETRY↗

Cladding Profilometry Analysis of Experimental Breeder Reactor-II Metallic Fuel Pins with HT9, D9, and SS316 Cladding

BISON finite element method fuel performance simulations were conducted using an existing automated process that couples the Fuels Irradiation & Physics Database (FIPD) and the Integral Fast Reactor Materials Information System database by writing input files and comparing the BISON output to post-irradiation fuel pin profilometry measurements contained within the databases. The importance of this work is to demonstrate the ability to benchmark fuel performance metallic fuel models within BISON using Experimental Breeder Reactor-II fuel pin data for a number of similar pins, while building off previous modeling efforts. Changes to the generic BISON input file include implementing pin specific axial power and flux profiles, pin specific fluences, frictional contact, and irradiation-induced volumetric swelling models for cladding. A statistical analysis of irradiation-induced volumetric swelling models for HT9, D9, and SS316 was performed for experiments X421/X421A, X441/X441A, and X486. Between these three experiments, there were 174 post-irradiation examination (PIE) profilometries used for validating the swelling models presented using a standard error of the estimate (SEE) method. Implementation of the volumetric swelling models for D9 and SS316 claddings was found to have a significant impact on the BISON profilometry simulated, where HT9 clad pins had an insignificant change due to low fluence values. BISON profilometry simulated for HT9, D9, and SS316 fuel pins agreed with PIE profilometry measurements, with assembly SEE values being 4.4 × 10−3 for X421A, 2.0 × 10−3 for X441A, and 2.8 × 10−3 for X486. D9 clad pins in X421/X421A had the highest SEE values, which is due to the BISON simulated profilometry being shifted axially. While this work accomplished its purpose to demonstrate the modeling of multiple fuel pins from the databases to help validate models, the results suggest that the continued development of metallic fuel models is necessary for qualifying new metallic fuel systems to better capture some physical performance phenomena, such as the hot pressing of U-Pu-Zr and the fuel cladding chemical interaction.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗

Lighting the dark molecular gas: H{sub 2} as a direct tracer

Robust knowledge of molecular gas mass is critical for understanding star formation in galaxies. The H{sub 2} molecule does not emit efficiently in the cold interstellar medium, hence the molecular gas content of galaxies is typically inferred using indirect tracers. At low metallicity and in other extreme environments, these tracers can be subject to substantial biases. We present a new method of estimating total molecular gas mass in galaxies directly from pure mid-infrared rotational H{sub 2} emission. By assuming a power-law distribution of H{sub 2} rotational temperatures, we can accurately model H{sub 2} excitation and reliably obtain warm (T ≳ 100 K) H{sub 2} gas masses by varying only the power law’s slope. With sensitivities typical of Spitzer/IRS, we are able to directly probe the H{sub 2} content via rotational emission down to ∼80 K, accounting for ∼15% of the total molecular gas mass in a galaxy. By extrapolating the fitted power-law temperature distributions to a calibrated single lower cutoff temperature, the model also recovers the total molecular content within a factor of ∼2.2 in a diverse sample of galaxies, and a subset of broken power-law models performs similarly well. In ULIRGs, the fraction of warm H{sub 2} gas rises with dust temperature, with some dependency on α {sub CO}. In a sample of five low-metallicity galaxies ranging down to 12+log[O/H]=7.8, the model yields molecular masses up to ∼100× larger than implied by CO, in good agreement with other methods based on dust mass and star formation depletion timescale. This technique offers real promise for assessing molecular content in the early universe where CO and dust-based methods may fail.

79 ASTRONOMY AND ASTROPHYSICS↗