SEARCH · Search NASA
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.
Lunar depth profiles of conductivity, dielectric constant and magnetic permeability from radar, radiothermic and magnetic measurements, discussing lunar models
Explore the source record for details and available documents.
A predictive temperature estimation technique which can be used to drive a model of the Sunrise/Sunset thermal 'snap' disturbance torque experienced by low Earth orbiting spacecraft is described. The twice per orbit impulsive disturbance torque is attributed to vehicle passage in and out of the Earth's shadow cone (umbra), during which large flexible appendages undergo rapidly changing thermal conditions. Flexible members, in particular solar arrays, experience rapid cooling during umbra entrance (Sunset) and rapid heating during exit (Sunrise). The thermal 'snap' phenomena has been observed during normal on-orbit operations of both the LANDSAT-4 satellite and the Communications Technology Satellite (CTS). Thermal 'snap' has also been predicted to be a dominant source of error for the TOPEX satellite. The fundamental equations used to model the Sunrise/Sunset thermal 'snap' disturbance torque for a typical solar array like structure will be described. For this derivation the array is assumed to be a thin, cantilevered beam. The time varying thermal gradient is shown to be the driving force behind predicting the thermal 'snap' disturbance torque and therefore motivates the need for accurate estimates of temperature. The development of a technique to optimally estimate appendage surface temperature is highlighted. The objective analysis method used is structured on the Gauss-Markov Theorem and provides an optimal temperature estimate at a prescribed location given data from a distributed thermal sensor network. The optimally estimated surface temperatures could then be used to compute the thermal gradient across the body. The estimation technique is demonstrated using a typical satellite solar array.
Technology is being developed to process signals from distributed sensors using distributed computations. These distributed sensors provide a new feedback capability for vibration control that has not been exploited. Additionally, the sensors proposed are of an optical and distributed nature and could be employed with known techniques of distributed optical computation (Fourier optics, etc.) to accomplish the control system functions of filtering and regulation in a distributed computer. This paper extends the traditional digital, optimal estimation and control theory to include distributed sensing and processing for this application. The design model assumes a finite number of modes which make it amenable to empirical determination of the design model via familiar modal-test techniques. The sensors are assumed to be distributed, but a finite number of point actuators are used. The design process is illustrated by application to a Euler beam. A simulation of the beam is used to design an optimal vibration control system that uses a distributed deflection sensor and nine linear force actuators. Simulations are also used to study the influence of design and processing errors on the performance.
Abstract not provided.
The Canadian Fire Weather Index (FWI) system, even though originally developed and calibrated for an upland Jack pine forest, is used globally to estimate fire danger for any fire environment. However, for some environments, such as peatlands, the applicability of the FWI in its current form, is often questioned. In this study, we replaced the original moisture codes of the FWI with hydrological estimates resulting from the assimilation of satellite-based L-band passive microwave observations into a peatland-specific land surface model. In a conservative approach that maintains the integrity of the original FWI structure, the distributions of the hydrological estimates were first matched to those of the corresponding original moisture codes before replacement. The resulting adapted FWI, hereafter called FWIpeat, was evaluated using satellite-based information on fire presence over boreal peatlands from 2010 through 2018. Adapting the FWI with model- and satellite-based hydrological information was found to be beneficial in estimating fire danger, especially when replacing the deeper moisture codes of the FWI. For late-season fires, further adaptations of the fine fuel moisture code show even more improvement due to the fact that late-season fires are more hydrologically driven. The proposed FWIpeat should enable improved monitoring of fire risk in boreal peatlands.
This analysis uses representative utility-scale and distributed wind energy projects to estimate the levelized cost of energy (LCOE) for land-based, offshore, and distributed wind power in the United States. Data and results detailed here are derived from 2021 commissioned plants and representative industry data as well as state-of-the-art modeling capabilities. Modeling is conducted to provide more granular detail on specific cost categories. This study represents the 11th annual installment and is intended to provide insight into current component-level costs as well as a basis for understanding variability in wind energy LCOE across the country.
Estimates of the energetic particle distribution in Jupiter's radiation belts are presented and are compared with previous estimates. Mathematical expressions are developed for the equatorial electron and proton fluxes, shielded electron and proton dose rates, and radiation lifetimes of electronic circuits. It is calculated that a 1 g/sq cm aluminum shield will screen out all protons of energy less than 27.5 MeV, and a 2 g/sq cm shield will screen out protons less than 40.6 MeV. The radiation lifetimes are based on a maximum permissible dose value of 10 million rads, a value 1/2 to 1/3 of the flux at which electronic circuits begin to fail. Estimated increases in lifetimes by using a 3 x 50 orbit instead of a 3 x 3 orbit, and going from 1 to 2 g/sq cm are given.
The magnitudes of the excitation of free precession of the lunar spin axis about the position defined by Cassini's laws, free libration in longitude, and free wobble are determined as a function of meteorite angular momentum relative to the lunar center of mass and the position of impact on the lunar surface. Angular-momentum conservation suffices for the estimates of precession and libration excitation, but a cratering model for the ejecta distribution is necessary for the estimate of the wobble excitation. The simultaneous excitation of free wobble is always associated with the excitation of precession, and the angular amplitude is at least comparable to, but may exceed, that of the induced precession by a factor of 3 or 4. It is possible to excite a free libration in longitude with no first-order excitation of free wobble, but generally, all three free motions are excited simultaneously. The induced libration will nearly always have the largest amplitude. For crater sizes scaled as powers of impact energy, impacts leaving craters as small as a few kilometers in diameter can excite free motions which will ultimately be observable by the lunar laser-ranging experiment.
Nonparametric techniques for probability distribution, probability density, and hazard function estimates for life quality
Traditional watershed modeling often overlooks the role of vegetation dynamics. There is also little quantitative evidence to suggest that increased physical realism of vegetation dynamics in process-based models improves hydrology and water quality predictions simultaneously. In this study, we applied a modified Soil and Water Assessment Tool (SWAT) to quantify the extent of improvements that the assimilation of remotely sensed Leaf Area Index (LAI) would convey to streamflow, soil moisture, and nitrate load simulations across a 16,860 km2 agricultural watershedin the midwestern United States. We modified the SWAT source code to automatically override the model’s built-in semiempirical LAI with spatially distributed and temporally continuous estimates from Moderate Resolution Imaging Spectroradiometer (MODIS). Compared to a “basic” traditional model with limited spatial information, our LAI assimilation model (i) significantly improved daily streamflow simulations during medium-to-low flow conditions, (ii) provided realistic spatial distributions of growing season soil moisture, and (iii) substantially reproduced the long-term observed variability of daily nitrate loads. Further analysis revealed that the overestimation or underestimation of LAI imparted a proportional cascading effect on how the model partitions hydrologic fluxes and nutrient pools. As such, assimilation of MODIS LAI data corrected the model’sLAI overestimation tendency, which led to a proportionally increased rootzone soil moisture and decreased plant nitrogen uptake. With these new findings, our study fills the existing knowledge gap regarding vegetation dynamics in watershed modeling and confirms that assimilation of MODIS LAI data in watershed models can effectively improve both hydrology and water quality predictions.
Natural gas (NG) supply chains contribute substantially to the global energy supply and anthropogenic methane emissions, making them frequent subjects of life cycle assessments (LCAs). To better characterize central tendencies and variability, we systematically reviewed and harmonized published estimates of life cycle greenhouse gas (GHG) emissions from United States NG supply chains. Results informed a streamlined LCA model (SLiNG-GHG: streamlined LCAs of NG-GHGs) that quantifies carbon dioxide and methane from three gates: transmission, distribution, and shipping. Median estimates employing harmonized emission inputs, are 10, 11, and 21 g CO2e/MJ gas (100-year global warming potentials [GWPs]), and 20, 22, and 33 g CO2e/MJ gas (20-year GWPs), delivered to each gate, respectively. Alternatively, inputting available, independent methane measurements, SLiNG-GHG estimates varied from -23% to +316% relative to baseline. Bottom-up inventories used in LCAs tend to underestimate methane compared with measurements. Results underscore the need for open-source, streamlined LCA models that can easily incorporate rapidly evolving measurements for non-experts like investors and regulators.
This report documents a statistical method for the "real-time" characterization of partially observed epidemics. Observations consist of daily counts of symptomatic patients, diagnosed with the disease. Characterization, in this context, refers to estimation of epidemiological parameters that can be used to provide short-term forecasts of the ongoing epidemic, as well as to provide gross information for the time-dependent infection rate. The characterization problem is formulated as a Bayesian inverse problem, and is predicated on a model for the distribution of the incubation period. The model parameters are estimated as distributions using a Markov Chain Monte Carlo (MCMC) method, thus quantifying the uncertainty in the estimates. The method is applied to the COVID-19 pandemic of 2020, using data at the country, provincial (e.g., states) and regional (e.g. county) levels. The epidemiological model includes a stochastic component due to uncertainties in the incubation period. This model-form uncertainty is accommodated by a pseudo-marginal Metropolis-Hastings MCMC sampler, which produces posterior distributions that reflect this uncertainty. We approximate the discrepancy between the data and the epidemiological model using Gaussian and negative binomial error models; the latter was motivated by the over-dispersed count data. For small daily counts we find the performance of the calibrated models to be similar for the two error models. For large daily counts the negative-binomial approximation is numerically unstable unlike the Gaussian error model. Application of the model at the country level (for the United States, Germany, Italy, etc.) generally provided accurate forecasts, as the data consisted of large counts which suppressed the day-to-day variations in the observations. Further, the bulk of the data is sourced over the duration before the relaxation of the curbs on population mixing, and is not confounded by any discernible country-wide second wave of infections. At the state-level, where reporting was poor or which evinced few infections (e.g., New Mexico), the variance in the data posed some, though not insurmountable, difficulties, and forecasts were able to capture the data with large uncertainty bounds. The method was found to be sufficiently sensitive to discern the flattening of the infection and epidemic curve due to shelter-in-place orders after around 90% quantile for the incubation distribution (about 10 days for COVID-19). The proposed model was also used at a regional level to compare the forecasts for the central and north-west regions of New Mexico. Modeling the data for these regions illustrated different disease spread dynamics captured by the model. While in the central region the daily counts peaked in the late April, in the north-west region the ramp-up continued for approximately three more weeks.
MERRA products were used to force an established ocean biogeochemical model to estimate surface carbon inventories and fluxes in the global oceans. The results were compared to public archives of in situ carbon data and estimates. The model exhibited skill for ocean dissolved inorganic carbon (DIC), partial pressure of ocean CO2 (pCO2) and air-sea fluxes (FCO2). The MERRA-forced model produced global mean differences of 0.02% (approximately 0.3 microns) for DIC, -0.3% (about -1.2 (micro) atm; model lower) for pCO2, and -2.3% (-0.003 mol C/sq m/y) for FCO2 compared to in situ estimates. Basin-scale distributions were significantly correlated with observations for all three variables (r=0.97, 0.76, and 0.73, P<0.05, respectively for DIC, pCO2, and FCO2). All major oceanographic basins were represented as sources to the atmosphere or sinks in agreement with in situ estimates. However, there were substantial basin-scale and local departures.
Particle size distributions are used to calculate nucleation rates during new particle formation (NPF). However, the low time resolution and poor sensitivity of instruments deployed to observe particle size distributions can lead to inaccurate estimates of nucleation rates. We calculated nucleation rates in Pittsburgh, PA from particle size distributions measured with the new high time resolution Condensation Particle Counters For Atmospheric Rapid Measurements (CPC FARM) and from gas precursor concentrations measured with a chemical ionization mass spectrometer. The CPC FARM nucleation rates were calculated by examining the appearance rate of particles in the smallest size bin at 1.6 nm. The nucleation rate of 1.6 nm determined every second resulted in positive and negative rates, indicating the influence of inhomogeneous mixing and advection of precursor gases and clusters. These effects are not currently considered when estimating nucleation rates from measurements. The nucleation rates from the CPC FARM were compared to the nucleation rates determined from a transient acid‐base reaction model using measured sulfuric acid concentration, strong base concentrations, and CPC detection efficiencies. The transient form of the model allowed the advection of heterogeneously mixed precursors to influence the nucleation rate. Good agreement between particle and gas‐based nucleation rates is observed during sulfuric acid‐driven NPF events. This study demonstrates that high time resolution measurements of short‐lived and rapidly formed particles are needed to accurately capture nucleation rates, which are likely influenced by advection of nucleation precursors and clusters in heterogeneous environments such as urban areas.
Abstract Large‐scale estimations of flood losses are often based on spatially aggregated inputs. This makes risk assessments vulnerable to aggregation bias, a well‐studied, sometimes substantial outcome in analyses that model fine‐grained spatial phenomena at coarse spatial units. To evaluate this potential in the context of large‐scale flood risk assessments, we use data from a high‐resolution flood hazard model and structure inventory for over 1.3 million properties in Massachusetts and examine how prominent data aggregation approaches affect the magnitude and spatial distribution of flood loss estimates. All considered aggregation approaches rely on aggregate structure inventories but differ in whether flood hazard is also aggregated. We find that aggregating only structure inventories slightly underestimates overall losses (−10% bias), and when flood hazard data is spatially aggregated to even relatively small spatial units (census block), statewide aggregation bias can reach +366%. All aggregation‐based procedures fail to capture the spatial covariation of inputs distributions in the upper tails that disproportionately generate total expected losses. Our findings are robust to several key assumptions, add important context to published risk assessments and highlight opportunities to improve flood loss estimation uncertainty quantification.
When performing the resilience enhancement for distribution networks, there are two obstacles to reliably model the uncertain contingencies: 1) decision-dependent uncertainty (DDU) due to various line hardening decisions, and 2) distributional ambiguity due to limited outage information during extreme weather events (EWEs). Here, to address these two challenges, this paper develops scenario-wise decision-dependent ambiguity sets (SWDD-ASs), where the DDU and distributional ambiguity inherent in EWE-induced contingencies are simultaneously captured for each possible EWE scenario. Then, a two-stage tri-level decision-dependent distributionally robust resilient enhancement (DD-DRRE) model is formulated, whose outputs include the optimal line hardening, distributed generation (DG) allocation, and proactive network reconfiguration strategy under the worst-case distributions in SWDD-ASs. Subsequently, the DD-DRRE model is equivalently recast to a mixed-integer linear programming (MILP)-based master problem and multiple scenario-wise subproblems, facilitating the adoption of a customized column-and-constraint generation (C&CG) algorithm. Finally, case studies demonstrate a remarkable improvement in the out-of-sample performance of our model, compared to its prevailing stochastic and robust counterparts. Moreover, the potential values of incorporating the ambiguity and distributional information are quantitatively estimated, providing a useful reference for planners with different budgets and risk-aversion levels.
Parameterization of lightning NO(x) in global chemical transport models requires a method to specify the temporal and geographic distribution of flash rates, an estimate of NO production per flash, and a method to specify the effective vertical distribution of NO production. Flash rate parameterizations using the upward convective cloud mass flux have been developed and tested for several GCM and assimilated meteorological data sets. This parameterization and others involving the cloud-top height and convective precipitation provide a generally acceptable climatological distribution of flashes after normalization with satellite data. However, considerable improvement in terms of geographic and temporal variability may be achieved in the future using cloud microphysics data available from next-generation GCMs. Cloud-resolving case-study simulations of convective transport and lightning NO production have yielded results which are directly applicable to the NO production per flash and the vertical NO distribution portions of the global model lightning parameterizations. In this work we have used cloud-resolving models (the Goddard Cumulus Ensemble Model (GCE) and MM5) to drive an off-line cloud-scale chemical transport model (CSCTM). The CSCTM, in conjunction with aircraft measurements of NO(x) in thunderstorms and ground-based lightning observations, has been used to constrain the amount of NO produced per flash. Observed lightning flash rates have been incorporated into the CSCTM, and several scenarios of NO production per intracloud (IC) and per cloud-to-ground (CG) flash have been tested in each of several case study storms. The resulting NO(x) mixing ratios are compared with aircraft measurements taken within the storm (typically the anvil region) to determine the most likely NO production scenario. The range of values of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced from observed NO spikes and from anvil flux calculations. Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. Estimates of mean NO production per flash vary by a factor of three from one simulated storm to another. Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. These profiles show a larger percentage of the total lightning NO(x) mass in the free troposphere than those suggested by Pickering et al. (1998).