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.
The CPU allocation was primarily used for nucleosyn thesis calculations using LANL’s PRISM code, radiative transfer simulations using LANL’s SuperNu code, and parameter inference of the radiative transfer simulations. Below, I document all relevant publications where Inti tutional Computing (IC) resources were used in relation to this allocation and their scientific impacts.
Abstract Space‐based lightning imagers have shown that complex cloud scenes consisting of multiple tall convective features, anvil clouds, and warm boundary cloud layers are illuminated by lightning in many different ways depending on where the lightning occurs and how energetic it is. Modifications to the optical lightning signals from radiative transfer in the cloud medium can lead to reductions in detection efficiency and location accuracy for these instruments and can also cause some of the optical signals that are detected to have unexpected spatial energy distributions. In this study, we perform Monte Carlo radiative transfer simulations of optical lightning emissions in clouds with complex 3‐D geometries to shed some light on the origins of certain irregular radiance patterns that have been recorded from orbit. We show that diffuse reflections off nearby cloud faces can explain lightning signals in nonelectrified clouds, tall clouds can result in poor optical transmission and suppressed radiances that could lead to missed events and that particularly favorable viewing conditions can cause otherwise normal lightning to produce a superbolt that is orders of magnitude brighter than the same flash seen from a different angle.
Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.
In this study we incorporate a new grid of kilonova simulations produced by the Monte Carlo radiative transfer code SuperNu in an inference pipeline for astrophysical transients and evaluate their performance. These simulations contain four different two-component ejecta morphology classes. We analyze follow-up observational strategies by Vera Rubin Observatory in optical and James Webb Space Telescope (JWST) in mid-infrared (MIR). Our analysis suggests that, within these strategies, it is possible to discriminate between the four different morphologies only when late-time JWST observations in MIR are available. We conclude that follow-ups by the new Vera Rubin Observatory alone are not sufficient to determine ejecta morphology. Additionally, we make comparisons between surrogate models based on radiative transfer simulation grids by SuperNu and POSSIS, by analyzing the historic kilonova AT2017gfo that accompanied the gravitational wave event GW170817. We show that both SuperNu and POSSIS models provide similar fits to photometric observations but their qualitative interpretations differ.
Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations.
Radiative transfer models (RTMs) describe how light is absorbed, scattered, and transmitted within leaves by simulating mechanistic light propagation processes. The PROSPECT model is based on measurable parameters (the leaf biochemical content) and a non-measurable parameter (the leaf anatomical structure represented by the leaf structure parameter (N)). The effect of N on the optical properties of leaves has been investigated through a number of local and global sensitivity analyses. Other studies have directly evaluated the effect of the leaf anatomical structure on spectral reflectance, particularly in the near infrared region. However, the relationship between N and the anatomical structure is unclear. Here, in this study, we leveraged eLeaf, a ray tracing-based 3D rice leaf simulator, to establish relationships between leaf anatomical features and spectral properties, enabling us to replace N in the PROSPECT-4 model with measurable leaf anatomical parameters and develop the RSPECT model. The leaf thickness at minor vein, leaf thickness at bulliform cells, mesophyll thickness at minor vein, and distance between two minor veins could be used to predict N effectively. The RSPECT model achieved spectral simulation accuracy comparable to PROSPECT-4 and was more suitable for parameter inversion of the physical and chemical properties of rice leaves, with relative root mean square errors of 7.4% for chlorophyll content, 5.6% for equivalent water thickness, and 7.5% for dry matter content. In conclusion, RSPECT improves radiative transfer modeling by integrating measurable anatomical features and provides a framework for extending this approach to other vegetation types.
Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Radiation model was effectively calibrated by modifying radiation model parameters for methanol. Computing integrated buoyancy flux, entrainment rate, and turbulent kinetic energy allowed for evaluation of less typical quantities in this validation study. Quantities were compared with experimental data or correlations and generally showed agreement. TKE needed fine mesh to be computed accurately. Predicted flame height less sensitive to variations in mixture fraction than temperature. Mixture fraction is a preferable threshold variable for this application.
This final report details all work performed on the project by both project partners. This project provided support to properly couple RTE+RRTMGP, a high-performance broadband radiation code, within DOE’s Energy Exascale Earth System Model (E3SM). RTE+RRTMGP is a successor to the RRTMG radiation code, which has been widely accepted for its speed and accuracy by the global modeling community, and has been in use in the NCAR CESM for many years and was implemented in the initial version of E3SM. However, the computational cost of RRTMG remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, issues that were rectified by the development of RTE+RRTMGP. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. One focus of the project was to enhance the code’s optimization on the limited number of emerging computing systems on which the model is expected be used, including Many Integrated Core (MIC) architectures and Graphics Processing Unit (GPU) hardware. We also developed additional capabilities for RTE+RRTMGP that E3SM scientists identified as important for the planned applications of the model. The result of our project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.
The "Reduced Speed of Light" (RSL) approximation is commonly used to speed up radiative transfer calculations in cosmological simulations. However, it has been shown previously that the RSL approximation leads to photon non-conservation when the radiation field is rapidly evolving in time. I show that these missing photons can be counted exactly for some numerical schemes. Adding them back into a simulation, however, is a much harder task. I show one example of such a scheme, which achieves sub-percent accuracy on simple tests. Unfortunately, the scheme performs much worse on semi-realistic simulations of cosmic reionization, leading to a faster overlap and significant errors in the point-wise comparison of the RSL radiation field with the reference simulation that maintains the full speed of light for the radiative transfer.
As part of the Advanced Simulation and Computing Verification and Validation (ASCVV) program, a 0.3-m diameter hydrocarbon pool fire with multiple fuels was modeled and simulated. In the study described in this report, systematic examination was performed on the radiation model used in a series of coupled Fuego/Nalu simulations. A calibration study was done with a medium-scale methanol pool fire and the effect of calibration traced throughout the radiation model. This analysis provided a more detailed understanding of the effect of radiation model parameters on each other and on other quantities in the simulations. Heptane simulation results were also examined using this approach and possible areas for further improvement of the models were identified. The effect of soot on radiative losses was examined by comparing heptane and methanol results.
Abstract Several mechanisms have been proposed to explain why the isotope ratios of precipitation vary in space and time and why they correlate with other climate variables like temperature and precipitation. Here, we argue that this behavior is best understood through the lens of radiative transfer, which treats the depletion of atmospheric vapor transport by precipitation as analogous to the attenuation of light by absorption or scattering. Building on earlier work by Siler et al., we introduce a simple model that uses the equations of radiative transfer to approximate the two-dimensional pattern of the oxygen isotope composition of precipitation ( δ p ) from monthly mean hydrologic variables. The model accurately simulates the spatial and seasonal variability in δ p within a state-of-the-art climate model and permits a simple decomposition of δ p variability into contributions from gradients in evaporation and the length scale of vapor transport. Outside the tropics, δ p is mostly controlled by gradients in evaporation, whose dependence on temperature explains the positive correlation between δ p and temperature (i.e., the temperature effect). At low latitudes, δ p is mostly controlled by gradients in the transport length scale, whose inverse relationship with precipitation explains the negative correlation between δ p and precipitation (i.e., the amount effect). This suggests that the temperature and amount effects are both mostly explained by the variability in upstream rainout, but they reflect distinct mechanisms governing rainout at different latitudes. Significance Statement The isotopic composition of precipitation has long been used to make inferences about past climates based on its observed relationship with precipitation in the tropics and with temperature at higher latitudes. These relationships—known as the “amount effect” and “temperature effect,” respectively—have been attributed to many different mechanisms, most of which are thought to operate at either high or low latitudes but not both. Here, we present a unified framework for interpreting the isotope variability that can explain the latitude dependence of the temperature and amount effects despite making no distinction between high and low latitudes. Although our results are generally consistent with certain interpretations of the amount effect, they suggest that the temperature effect is widely misunderstood.
This study introduces a new technique to compute direct normal irradiance (DNI) for improving the National Solar Radiation Data Base (NSRDB). A finite-surface integration algorithm is developed to compute solar radiation in differential solid angles and efficiently infer its contribution to a surface perpendicular to the solar direction. A lookup table of cloud bi-directional transmittance distribution function (BTDF) is developed by use of the discrete ordinates radiative transfer (DISORT) model for possible solar and observing directions and various cloud optical and microphysical properties. In each solar incident direction, DNI is given by the cloud BTDFs from approximately 200 differential solid angles. The simulated DNI is calibrated and evaluated using surface observations at the National Renewable Energy Laboratory's (NREL's) Solar Energy Research Laboratory (SRRL) and the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) facility.
Kilonovae are the electromagnetic transients created by the radioactive decay of freshly synthesized elements in the environment surrounding a neutron star merger. To study the fundamental physics in these complex environments, kilonova modeling requires, in part, the use of radiative transfer simulations. The microphysics involved in these simulations results in high computational cost, prompting the use of emulators for parameter inference applications. Utilizing a training set of 22 248 high-fidelity simulations (composed of 412 unique ejecta parameter combinations evaluated at 54 viewing angles), we use a neural network to efficiently train on existing radiative transfer simulations and predict light curves for new parameters in a fast and computationally efficient manner. Our neural network can generate millions of new light curves in under a minute. We discuss our emulator's degree of off-sample reliability and parameter inference of the AT2017gfo observational data. Finally, we discuss tension introduced by multiband inference in the parameter inference results, particularly with regard to the neural network's recovery of viewing angle. Published by the American Physical Society 2024
Microwave cloud polarized observations have shown the potential to improve precipitation retrievals since they are linked to the orientation and shape of ice habits. Stratiform clouds show larger brightness temperature (TB) polarization differences (PDs), defined as the vertically polarized TB (TBV) minus the horizontally polarized TB (TBH), with ~10 K PD values at 89 GHz due to the presence of horizontally aligned snowflakes, while convective regions show smaller PD signals, as graupel and/or hail in the updraft tend to become randomly oriented. The launch of the global precipitation measurement (GPM) microwave imager (GMI) has extended the availability of microwave polarized observations to higher frequencies (166 GHz) in the tropics and midlatitudes, previously only available up to 89 GHz. This study analyzes one year of GMI observations to explore further the previously reported stable relationship between the PD and the observed TBs at 89 and 166 GHz, respectively. The latitudinal and seasonal variability is analyzed to propose a cloud scattering polarization parameterization of the PD-TB relationship, capable of reconstructing the PD signal from simulated TBs. Given that operational radiative transfer (RT) models do not currently simulate the cloud polarized signals, this is an alternative and simple solution to exploit the large number of cloud polarized observations available. Finally, the atmospheric radiative transfer simulator (ARTS) is coupled with the weather research and forecasting (WRF) model, in order to apply the proposed parameterization to the RT simulated TBs and hence infer the corresponding PD values, which show to reproduce the observed GMI PDs well.
Attenuation of solar radiation between the receiver and the heliostat field in concentrated solar power (CSP) tower plants can reduce the overall system performance significantly. The attenuation varies strongly with time and the average attenuation at different sites might also vary strongly from each other. If no site specific attenuation data is available, the optimal plant design cannot be determined and rough estimations of the attenuation effect are required leading to high uncertainties of yield analysis calculations. The attenuation is caused mainly by water vapor content and aerosol particles in the lower atmospheric layer above ground. Although several on-site measurement systems have been developed during recent years, attenuation data sets are usually not available to be included during the plant project development. An Atmospheric Attenuation (AATTENUATION) model to derive the atmospheric transmittance between a heliostat and receiver on the basis of common direct normal irradiance (DNI), temperature, relative humidity, and barometric pressure measurements was developed and validated by the authors earlier. The model allows the accurate estimation of attenuation for sites with low attenuation and gives an estimation of the attenuation for less clear sites. However, the site-dependent coefficients of the AATTENUATION model had to be developed individually for each site of interest, which required time-consuming radiative transfer simulations, considering the exact location and altitude, as well as the pre-dominant aerosol type at the location. This strongly limited the application of the model despite its typically available input data. In this manuscript, a look-up table (LUT) is presented which enables the application of the AATTENUATION model at the site of interest without the necessity to perform the according complex radiative transfer calculations for each site individually. This enables the application of the AATTENUATION model for virtually all resource assessments for tower plants and in an operational mode in real time within plant monitoring systems around the world. The LUT also facilitates the generation of solar attenuation maps on the basis of long-term meteorological data sets which can be considered during resource assessment for CSP tower plant projects. The LUTs are provided together with this manuscript as supplementary files. The LUT for the AATTENUATION model was developed for a solar zenith angle (SZA) grid of 1°, an altitude grid of 100 m, 7 different standard aerosol types and the standard AFGL atmospheres for mid-latitudes and the tropics. The LUT was tested against the original version of the AATTENUATION model at 4 sites in Morocco and Spain, and it was found that the additional uncertainty introduced by the application of the LUT is negligible. With the information of latitude, longitude, altitude above mean sea level, DNI, relative humidity (RH), ambient temperature (Tair), and barometric pressure (bp), the attenuation can be now derived easily for each site of interest.
Very young (t ≲ 10 Myr) stars possess strong magnetic fields that channel ionized gas from the interiors of their circumstellar disks to the surface of the star. Upon impacting the stellar surface, the shocked gas recombines and emits hydrogen spectral lines. To characterize the density and temperature of the gas within these accretion streams, we measure equivalent widths of Brackett (Br) 11–20 emission lines detected in 1101 APOGEE spectra of 326 likely pre-main-sequence accretors. For sources with multiple observations, we measure median epoch-to-epoch line strength variations of 10% in Br11 and 20% in Br20. We also fit the measured line ratios to predictions of radiative transfer models by Kwan & Fischer. We find characteristic best-fit electron densities of n e = 10 11 –10 12 cm -3 , and excitation temperatures that are inversely correlated with electron density (from T ~ 5000 K for n e ~ 10 12 cm -3 to T ~ 12,500 K at n e ~ 10 11 cm -3 ). These physical parameters are in good agreement with predictions from modeling of accretion streams that account for the hydrodynamics and radiative transfer within the accretion stream. We also present a supplementary catalog of line measurements from 9733 spectra of 4255 Brackett emission-line sources in the APOGEE Data Release 17 data set.
Cold-air outbreak (CAO) clouds in the Arctic are commonly mixed-phase (MP); however, the partitioning of the amount of ice and water in CAO clouds and precipitation is not always well observed. Understanding how cloud phases partition as a function of cloud life cycle is important for predicting snowfall rates, convective life cycle, and intensity at weather timescales. The partitioning into liquid versus ice also has radiative impacts that are consequential for climate. These concerns motivated the incorporation of an airborne G-band vapor radiometer (GVR) into an National Science Foundation-supported aircraft campaign named the Cold Air outbreak Experiment in the Sub-Arctic Region (CAESAR). The GVR is an upward-pointing passive microwave radiometer using four frequencies centered around the 183.31 GHz water vapor absorption line, displaced by +- 1, 3, 7, and 14 GHz. For context, The U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility operates a surface-based GVR at its North Slope of Alaska (NSA) site. The same GVR has been used previously for a field campaign in the southeast Pacific, where an offset was noticed between brightness temperatures (Tbs) measured under clear skies compared to those calculated from a radiative transfer model. To account for any calibration offsets to the Tbs, a request was made to DOE to allow the GVR to operate at ARM’s Southern Great Plains (SGP) observatory and enable comparisons between its measurements and those available at SGP. This request was granted, titled ‘UMGVR_CAL’, short for ‘UMGVR_Calibration’. From October 30 to November 10, 2023, the GVR was deployed to the ARM SGP site to take advantage of their regular, nearby radiosonde launches under clear-sky conditions. The latter were determined using the SGP total sky imagery data. The GVR brightness temperatures in these clear-sky conditions were compared to those calculated by a radiative transfer code (PAMTRA) based on the SGP radiosondes. During the campaign, four suitable clear-sky episodes could be used for the GVR calibration. While few in number, these proved to be enough to satisfy our goal.