Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative transfer simulations”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

C HIMERA : A Massively Parallel Code for Core-collapse Supernova Simulations

Herein we provide a detailed description of the C HIMERA code, a code developed to model core collapse supernovae (CCSNe) in multiple spatial dimensions. The CCSN explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the fundamental physics of the neutrino–stellar core weak interactions and the equations of state of all stellar core constituents—particularly, the nuclear equation of state associated with core nucleons, both free and bound in nuclei. C HIMERA , by incorporating detailed neutrino transport, realistic neutrino–matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with CCSNe. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. C HIMERA is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Impact of Radiative Transfer at Reduced Spectral Resolution in Large‐Eddy Simulations of Convective Clouds

Abstract Many radiative transfer schemes approximate the spectral integration over ∼10 5 to ∼10 6 wavelengths with correlated k ‐distributions methods that typically require only 10 1 –10 2 spectral integration points ( g ‐points). The exact number of g ‐points is then chosen as an optimal balance between computational costs and accuracy, normally assessed in terms of a number of radiative quantities. How this radiative accuracy propagates to simulation accuracy, however, is not straightforward. In this study, we therefore explore the sensitivity of cloud properties in large‐eddy simulations (LES) to the accuracy of radiative fluxes and heating rates. We first generate smaller sets of g ‐points from existing k ‐distributions by repeatedly combining adjacent g ‐points while maintaining the highest possible accuracy on a chosen set of radiative metrics. Next, we perform three sets of LES with varying cloud—radiation coupling pathways, and therefore different requirements for the accuracy of the radiative transfer computations, to investigate how these smaller and thus less accurate k ‐distributions affect simulation characteristics. The decrease in radiative accuracy with 3–4 times smaller k ‐distributions results in biases in cloud properties that are relative small compared to their temporal fluctuations. These results show potential for speeding up radiative transfer computations in cloud‐resolving models by reducing the resolved spectral detail. However, more statistically converged simulations and a wider set of case studies is required to fully assess the robustness of our results.

Meteorology & Atmospheric Sciences↗
NASA NTRS1965

Radiant heating simulation.

Radiation energy transfer rates characteristic of environments for advanced entry missions, examining laboratory simulation techniques, particularly plasma arc generators

PLASMA GENERATOR↗
NASA NTRS1979

Transfer of solar radiation in the atmosphere

Current needs in the field of atmospheric radiative transfer include the development and validation of fast but accurate calculation techniques for use in climate modeling. Progress toward meeting these needs can be made by application of existing computational techniques and use of observational data for radiation transfer across the earth-atmosphere system boundaries. In this study radiative transfer model simulations are compared to satellite observations of reflected solar radiation. Calculations are presented from both detailed and approximate transfer models for the solar radiation reflected from the earth-atmosphere system. These calculations are compared with the measurements made by the Earth Radiation Budget (ERB) experiment on Nimbus 6 in August 1975. Results show good agreement between simulations and observations.

Suttles, J. T.↗
NASA NTRS2014

Optimization of a Radiative Transfer Forward Operator for Simulating SMOS Brightness Temperatures over the Upper Mississippi Basin, USA

The Soil Moisture and Ocean Salinity (SMOS) satellite mission is routinely providing global multi-angular observations of brightness temperature (TB) at both horizontal and vertical polarization with a 3-day repeat period. The assimilation of such data into a land surface model (LSM) may improve the skill of operational flood forecasts through an improved estimation of soil moisture (SM). To accommodate for the direct assimilation of the SMOS TB data, the LSM needs to be coupled with a radiative transfer model (RTM), serving as a forward operator for the simulation of multi-angular and multi-polarization top of atmosphere TBs. This study investigates the use of the Variable Infiltration Capacity (VIC) LSM coupled with the Community Microwave Emission Modelling platform (CMEM) for simulating SMOS TB observations over the Upper Mississippi basin, USA. For a period of 2 years (2010-2011), a comparison between SMOS TBs and simulations with literature-based RTM parameters reveals a basin averaged bias of 30K. Therefore, time series of SMOS TB observations are used to investigate ways for mitigating these large biases. Specifically, the study demonstrates the impact of the LSM soil moisture climatology in the magnitude of TB biases. After CDF matching the SM climatology of the LSM to SMOS retrievals, the average bias decreases from 30K to less than 5K. Further improvements can be made through calibration of RTM parameters related to the modeling of surface roughness and vegetation. Consequently, it can be concluded that SM rescaling and RTM optimization are efficient means for mitigating biases and form a necessary preparatory step for data assimilation.

SMOS↗
DOE OSTI2025

Angular-spatial hp -adaptivity for radiative transfer with discontinuous Galerkin spectral element methods

Radiative transfer is important for many science and engineering applications, and numerical simulations of radiative transfer can be challenging. For instance, the radiation field is seven-dimensional – three spatial, two angular, one wavelength, and one temporal – and often features steep gradients. Therefore, memory usage is a key issue. To reduce memory, some past work has investigated the use of adaptive mesh refinement (AMR), typically for either the spatial or angular coordinate, and typically for only h -adaptivity. Here, we propose the use of AMR for the spatial and angular coordinates together, and the use of h - and p -adaptivity together as hp -AMR for the potential for further memory savings. We implemented the proposed method for several test cases in two spatial and one angular dimension, with the discontinuous Galerkin spectral element method. These test cases featured highly anisotropic angular radiation, with or without steep spatial gradients. Our primary findings from these test cases were: (1) Angular hp -adaptivity can deliver the radiation solution with the same accuracy as, and with much less computational memory than, uniform angular h - or p -refinements, or angular h -adaptivity alone. This is most obvious when the incoming radiation is highly anisotropic, in which case the savings can be orders of magnitude. (2) Full spatial-angular hp -adaptivity is more efficient in solution representation, compared to solely spatial or solely angular -adaptivity. This is most evident when steep gradients are present in both the spatial and angular distribution. These results suggest that adaptive spatial- hp angular-refinement may perform well in large-scale seven-dimensional applications.

Adaptive refinement↗
NASA NTRS2020

Long Term Stability Monitoring of Aqua MODIS Thermal Emissive Bands through Radiative Transfer Modeling

Moderate Resolution Imaging Spectroradiometer (MODIS) on Aqua has been in operation providing continuous global observations for science research and applications since 2002. The long-term stability of thermal emissive bands (TEBs) of Aqua MODIS was monitored through inter-comparisons with measurements by hyperspectral infrared sensors, such as AIRS on Aqua and IASI on MetOp satellite, or through long term vicarious monitoring over cold targets, such as Dome- C and deep convective clouds (DCC). In this paper, a radiative transfer modeling-based simulation model using Community Radiative Transfer Model (CRTM) is developed to perform the long-term monitoring of the stability of TEBs of Aqua MODIS. CRTM is a fast-radiative transfer model for calculations of radiances for satellite infrared radiometers and is able to output infrared radiance and brightness temperature at spectral bands of MODIS. Long term European Centre for Medium-Range Weather Forecasts (ECMWF) global atmospheric reanalysis data, such as temperature and humidity profiles, are used as inputs to the CRTM simulation. By confining the area of interest to be over low to middle latitude ocean, the long-term stabilities of selected Aqua MODIS TEBs are monitored through Observation–Background (O-B) brightness temperature (BT) bias between MODIS measurements and BT retrieval from CRTM simulation. The consistency and relative stability between Aqua MODIS and ECMWF reanalysis data for surface channels of MODIS are evaluated. In addition, the radiative transfer modeling with CRTM enables us to evaluate the impacts of long term variation of global CO2 distribution on the O-B BT biases for CO2 channels of Aqua MODIS through comparison of simulations with constant or long term variable CO2 as inputs. The O-B analysis with RTM show that Aqua MODIS surface channels are all radiometrically stable with yearly BT bias drift less than 0.004K/year for B20, B22, B23 and ~0.01K/year for B31 and B32. The CO2 absorption channels of Aqua MODIS, e.g. B33-B36, are stable with BT bias drift < 0.005K/year.

MODIS↗
DOE OSTI2024

Light Transfers Through a Koch Shape Cloud

Abstract Modeling radiative transfer in a 3D cloudy atmosphere is critical to climate projections. A recently developed fast 3D radiation parameterization scheme gains some success in quantifying horizontal radiative transfer through cloud sides using cloud area fraction. Based on 3D Monte Carlo simulations of radiative transfer through an idealized single‐layer cloud with Koch‐shaped fractal geometry edges, here we show that radiative energy transport through cloud sides correlates more significantly with cloud area fraction than with cloud perimeter length. The results exemplify the importance of accounting for the horizontal radiative energy exchanges between cloud‐free and cloudy regions with cloud area fraction. Results from additional sensitivity simulations show that increased cloud vertical extent often enhances cloud‐side sunlight leak more significantly than cloud‐side sunlight interception. At low sun elevations, cloud‐side sunlight interception is enhanced more than cloud‐side sunlight leak does with the increase of cloud mass.

58 GEOSCIENCES↗
NASA NTRS2011

Computing Temperatures in Optically Thick Protoplanetary Disks

We worked with a Monte Carlo radiative transfer code to simulate the transfer of energy through protoplanetary disks, where planet formation occurs. The code tracks photons from the star into the disk, through scattering, absorption and re-emission, until they escape to infinity. High optical depths in the disk interior dominate the computation time because it takes the photon packet many interactions to get out of the region. High optical depths also receive few photons and therefore do not have well-estimated temperatures. We applied a modified random walk (MRW) approximation for treating high optical depths and to speed up the Monte Carlo calculations. The MRW is implemented by calculating the average number of interactions the photon packet will undergo in diffusing within a single cell of the spatial grid and then updating the packet position, packet frequencies, and local radiation absorption rate appropriately. The MRW approximation was then tested for accuracy and speed compared to the original code. We determined that MRW provides accurate answers to Monte Carlo Radiative transfer simulations. The speed gained from using MRW is shown to be proportional to the disk mass.

radiative transfer↗
NASA NTRS2023

Developing a Radar Signal Simulator for the Community Radiative Transfer Model

Active radar instruments provide vertically resolved clouds and precipitation measurements that cannot be provided by the passive instruments. These active measurements are not conventionally assimilated into the data assimilation systems because of the lack of fast forward radiative transfer models and also difficulties in the error modelling of the measurements. This paper describes the development, evaluation, and sensitivity analysis for a forward radar model implemented in the Community Radiative Transfer Model (CRTM). The scattering properties required by the forward model are provided by the hydrometeor lookup tables that were generated using the discrete dipole approximation. The model is able to calculate both the reflectivity and the attenuated reflectivity for any given radar instrument at any given zenith angles as long as CRTM instrument specific coefficients are available. The evaluation using CloudSat measurements shows a very good agreement between the simulations and measurements as long as the input profiles of hydrometeors are consistent with the measured reflectivity profiles. Major sources contributing to the differences between the measured and simulated reflectivities are input hydrometeor profiles, scattering lookup tables, lack of melting layer in the forward model, CRTM scattering solvers, and attenuation calculations. In addition to the forward model, both Tangent Linear and Adjoint of the model are also implemented and tested within CRTM. These components may be required by some data assimilation systems for the assimilation of radar measurements.

radar↗
DOE OSTI2021

A Broad Grid of 2D Kilonova Emission Models

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M ⊙ and 0.05 to 0.3c, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (Y e ) ejecta with a robust r-process composition and either a spherical or lobed high-Y e ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Parameterization of the Cloud Scattering Polarization Signal Derived From GPM Observations for Microwave Fast Radative Transfer Models

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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A Comparison of Three Neodymium Atomic Data Sets for Kilonova Modeling

We examine the impact of input neodymium (Nd) atomic data on the light curves and spectra of kilonovae (KNe), probing the sensitivity of kilonova observables to the atomic physics of this important lanthanide element. We use the SuperNu Monte Carlo radiative transfer code, simulating a simple semianalytic 1D kilonova (KN) with a pure Nd atmosphere, fixing the radiative transfer method while using input atomic data generated by three different codes: the LANL suite of atomic physics codes, HULLAC, and Autostructure. We see that the choice of atomic data significantly shapes the resulting light curves and spectra. Peak bolometric luminosities differ by a ratio of nearly 1.5 between HULLAC/Autostructure and LANL data sets. Moreover, we observe significant near- to mid-IR differences in the structure of the spectra. We specifically attribute these differences to the choice of atomic data for neutral Nd I. Many of the results here have been adapted from a presentation at “Radiative Transfer and Atomic Physics of Kilonovae” in Stockholm, 2023. We additionally present a LANL data set with energies calibrated to available values in the NIST Atomic Spectra Database, and demonstrate that this calibration also significantly affects IR spectral structure at late time. The substantial differences in KN observables that arise from tuning the atomic data of just one lanthanide element highlight the special attention that must be paid to atomic physics uncertainties when modeling KNe, from AT2017gfo to beyond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2007

General Relativistic Radiative Transfer and General Relativistic MHD Simulations of Accretion and Outflows of Black Holes

We have calculated the emission from relativistic flows in black hole systems using a fully general relativistic radiative transfer, with flow structures obtained by general relativistic magnetohydrodynamic simulations. We consider thermal free-free emission and thermal synchrotron emission. Bright filament-like features are found protruding (visually) from the accretion disk surface, which are enhancements of synchrotron emission when the magnetic field is roughly aligned with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We propose that variations and location drifts of the features are responsible for certain X-ray quasi-periodic oscillations (QPOs) observed in black-hole X-ray binaries.

Fuerst, Steven V.↗
NASA NTRS1974

Reinterpretation of Mariner 9 IRIS data on the basics of a simulation of radiative-conductive transfer in the dust-laden Martian atmosphere

A reinterpretation of the convective heat transport concept is presented for the IRIS measurements as an essential energy transfer mechanism. Thermal convection in a radiatively heated ground atmosphere system is discussed along with CO2 sublimation and diffusion from the Martian polar ice cap. A systematic reappraisal of the existing model is also included.

Pallmann, A. J.↗
NASA NTRS1989

Analysis and use of VAS satellite data

Four interrelated investigations have examined the analysis and use of VAS satellite data. A case study of VAS-derived mesoscale stability parameters suggested that they would have been a useful supplement to conventional data in the forecasting of thunderstorms on the day of interest. A second investigation examined the roles of first guess and VAS radiometric data in producing sounding retrievals. Broad-scale patterns of the first guess, radiances, and retrievals frequently were similar, whereas small-scale retrieval features, especially in the dew points, were often of uncertain origin. Two research tasks considered 6.7 micron middle tropospheric water vapor imagery. The first utilized radiosonde data to examine causes for two areas of warm brightness temperature. Subsidence associated with a translating jet streak was important. The second task involving water vapor imagery investigated simulated imagery created from LAMPS output and a radiative transfer algorithm. Simulated image patterns were found to compare favorably with those actually observed by VAS. Furthermore, the mass/momentum fields from LAMPS were powerful tools for understanding causes for the image configurations.

Fuelberg, Henry E.↗