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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Acute effects of a large bolide impact simulated by a global atmospheric circulation model

The goal is to use a global three-dimensional atmospheric circulation model developed for studies of atmospheric effects of nuclear war to examine the time evolution of atmospheric effects from a large bolide impact. The model allows for dust and NOx injection, atmospheric transport by winds, removal by precipitation, radiative transfer effects, stratospheric ozone chemistry, and nitric acid formation and deposition on a simulated Earth having realistic geography. Researchers assume a modest 2 km-diameter impactor of the type that could have formed the 32 km-diameter impact structure found near Manson, Iowa and dated at roughly 66 Ma. Such an impact would have created on the order of 5 x 10 to the 10th power metric tons of atmospheric dust (about 0.01 g cm(-2) if spread globally) and 1 x 10 to the 37th power molecules of NO, or two orders of magnitude more stratospheric NO than might be produced in a large nuclear war. Researchers ignore potential injections of CO2 and wildfire smoke, and assume the direct heating of the atmosphere by impact ejecta on a regional scale is not large compared to absorption of solar energy by dust. Researchers assume an impact site at 45 N in the interior of present day North America.

Thompson, Starley L.↗
NASA NTRS1978

Carbonaceous materials subjected to extreme heating - A comparison of numerical simulation and experiments

The ablation of carbonaceous materials in a hydrogen-helium stream has been simulated using a charring materials ablator computer code. These results are compared with the first ablation data to be obtained from the Ames-NASA Giant Planet Pilot Facility. Test stream diagnostics and ablation effects on convective and radiative heat transfer are discussed since these parameters constitute important input data to the numerical simulation. Graphite ablation was predicted to within 10 to 20%, and carbon-phenolic somewhat less accurately.

Davy, W. C.↗
NASA NTRS2014

Simulating Wake Vortex Detection with the Sensivu Doppler Wind Lidar Simulator

In support of NASA's Atmospheric Environment Safety Technologies NRA research topic on Wake Vortex Hazard Investigation, Aerospace Innovations (AI) investigated a set of techniques for detecting wake vortex hazards from arbitrary viewing angles, including axial perspectives. This technical report describes an approach to this problem and presents results from its implementation in a virtual lidar simulator developed at AI. Threedimensional data volumes from NASA's Terminal Area Simulation System (TASS) containing strong turbulent vortices were used as the atmospheric domain for these studies, in addition to an analytical vortex model in 3-D space. By incorporating a third-party radiative transfer code (BACKSCAT 4), user-defined aerosol layers can be incorporated into atmospheric models, simulating attenuation and backscatter in different environmental conditions and altitudes. A hazard detection algorithm is described that uses a twocomponent spectral model to identify vortex signatures observable from arbitrary angles.

Ramsey, Dan↗
DOE OSTI2025

RSPECT: A PROSPECT-based model incorporating the real structure of rice leaves

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.

hyperspectral↗
NASA NTRS2012

Validation of Aquarius Measurements Using Radiative Transfer Models at L-Band

Aquarius/SAC-D was launched in June 2011 by NASA and CONAE (Argentine space agency). Aquarius includes three L-band (1.4 GHz) radiometers dedicated to measuring sea surface salinity. We report detailed comparisons of Aquarius measurements with radiative transfer model predictions. These comparisons were used as part ofthe initial assessment of Aquarius data. In particular, they were used successfully to estimate the radiometer calibration bias and stability. Further comparisons are being performed to assess the performance of models in the retrieval algorithm for correcting the effect of sources of geophysical "noise" (e.g. the galactic background, atmospheric attenuation and reflected signal from the Sun). Such corrections are critical in bringing the error in retrieved salinity down to the required 0.2 practical salinity unit (psu) on monthly global maps at 150 km by 150 km resolution. The forward models making up the Aquarius simulator have been very useful for preparatory studies in the years leading to Aquarius' launch. The simulator includes various components to compute effects ofthe following processes on the measured signal: 1) emission from Earth surfaces (ocean, land, ice), 2) atmospheric emission and absorption, 3) emission from the Sun, Moon and celestial Sky (directly through the antenna sidelobes or after reflection/scattering at the Earth surface), 4) Faraday rotation, and 5) convolution of the scene by the antenna gain patterns. Since the Aquarius radiometers tum-on in late July 2011, the simulator has been used to perform a first order validation of the data. This included checking the order of magnitude ofthe signal over ocean, land and ice surfaces, checking the relative amplitude of signal at different polarizations, and checking the variation with incidence angle. The comparisons were also used to assess calibration bias and monitor instruments calibration drift. The simulator is also being used in the salinity retrieval. For example, initial assessments of the salinity retrieved from Aquarius data showed degradation in accuracy at locations where glint from the galactic sky background was important. This was traced to an inaccurate correction for the Sky glint. We present comparisons of the simulator prediction to the Aquarius data in order to assess the performances of the models of various physical processes impacting the measurements, such as the effect of sea surface roughness, the impact of the celestial Sky and the Sun emission scattered at the rough ocean surface. We discuss what components of the simulator appear reliable and which ones need improvements. Improved knowledge on the radiative transfer models at L-band will not only lead to better salinity retrieved from Aquarius data, it will also allow be beneficial for SMOS or the upcoming SMAP mission.

Dinnat, E.↗
DOE OSTI2024

Making Uintah Performance Portable for Department of Energy Exascale Testbeds

To help ease ports to forthcoming Department of Energy (DOE) exascale systems, testbeds have been made available to select users. These testbeds are helpful for preparing codes to run on the same hardware and similar software as in their respective exascale systems. This paper describes how the Uintah Computational Framework, an open-source asynchronous many-task (AMT) runtime system, has been modified to be performance portable across the DOE Crusher, DOE Polaris, and DOE Sunspot testbeds in preparation for portable simulations across the exascale DOE Frontier and DOE Aurora systems. The Crusher, Polaris, and Sunspot testbeds feature the AMD MI250X, NVIDIA A100, and Intel PVC GPUs, respectively. This performance portability has been made possible by extending Uintah’s intermediate portability layer [18] to additionally support the Kokkos::HIP, Kokkos::OpenMPTarget, and Kokkos::SYCL back-ends. This paper also describes notable updates to Uintah’s support for Kokkos, which were required to make this extension possible. Results are shown for a challenging radiative heat transfer calculation, central to the University of Utah’s predictive boiler simulations. These results demonstrate single-source portability across AMD-, NVIDIA-, and Intel-based GPUs using various Kokkos back-ends.

Holmen, John↗
NASA NTRS2016

The Impact of ARM on Climate Modeling

Climate models are among humanity's most ambitious and elaborate creations. They are designed to simulate the interactions of the atmosphere, ocean, land surface, and cryosphere on time scales far beyond the limits of deterministic predictability, and including the effects of time-dependent external forcings. The processes involved include radiative transfer, fluid dynamics, microphysics, and some aspects of geochemistry, biology, and ecology. The models explicitly simulate processes on spatial scales ranging from the circumference of the Earth down to one hundred kilometers or smaller, and implicitly include the effects of processes on even smaller scales down to a micron or so. The atmospheric component of a climate model can be called an atmospheric global circulation model (AGCM). In an AGCM, calculations are done on a three-dimensional grid, which in some of today's climate models consists of several million grid cells. For each grid cell, about a dozen variables are time-stepped as the model integrates forward from its initial conditions. These so-called prognostic variables have special importance because they are the only things that a model remembers from one time step to the next; everything else is recreated on each time step by starting from the prognostic variables and the boundary conditions. The prognostic variables typically include information about the mass of dry air, the temperature, the wind components, water vapor, various condensed-water species, and at least a few chemical species such as ozone. A good way to understand how climate models work is to consider the lengthy and complex process used to develop one. Lets imagine that a new AGCM is to be created, starting from a blank piece of paper. The model may be intended for a particular class of applications, e.g., high-resolution simulations on time scales of a few decades. Before a single line of code is written, the conceptual foundation of the model must be designed through a creative envisioning that starts from the intended application and is based on current understanding of how the atmosphere works and the inventory of mathematical methods available.

atmospheric circulation↗
NASA NTRS2014

Introducing Multisensor Satellite Radiance-Based Evaluation for Regional Earth System Modeling

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.

Planetary Boundary Layer↗
DOE OSTI2020

A Cold Moderator For Sub-Thermal Neutron Flux Enhancement at the RPI-LINAC

The importance of high-quality thermal neutron scattering kernels coupled with a lack of accurate sub-thermal total cross section data prompted researchers at Rensselaer Polytechnic Institute (RPI) to design a cold moderation system to enhance sub-thermal neutron flux (defined as flux below 5 meV), for nuclear data measurements at the RPI Gaerttner LINAC facility. Using a series of neutronic and heat transfer simulations obtained with MCNP and COMSOL, an optimized final design was produced. The final design exceeded current target sub-thermal neutron producing capabilities by a factor of 8 in the sub-thermal region and is easily configurable with existing neutron producing targets.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2014

Simulation for Prediction of Entry Article Demise (SPEAD): an Analysis Tool for Spacecraft Safety Analysis and Ascent/Reentry Risk Assessment

For the purpose of performing safety analysis and risk assessment for a probable offnominal suborbital/orbital atmospheric reentry resulting in vehicle breakup, a synthesis of trajectory propagation coupled with thermal analysis and the evaluation of node failure is required to predict the sequence of events, the timeline, and the progressive demise of spacecraft components. To provide this capability, the Simulation for Prediction of Entry Article Demise (SPEAD) analysis tool was developed. This report discusses the capabilities, modeling, and validation of the SPEAD analysis tool. SPEAD is applicable for Earth or Mars, with the option for 3 or 6 degrees-of-freedom (DOF) trajectory propagation. The atmosphere and aerodynamics data are supplied in tables, for linear interpolation of up to 4 independent variables. The gravitation model can include up to 20 zonal harmonic coefficients. The modeling of a single motor is available and can be adapted to multiple motors. For thermal analysis, the aerodynamic radiative and free-molecular/continuum convective heating, black-body radiative cooling, conductive heat transfer between adjacent nodes, and node ablation are modeled. In a 6- DOF simulation, the local convective heating on a node is a function of Mach, angle-ofattack, and sideslip angle, and is dependent on 1) the location of the node in the spacecraft and its orientation to the flow modeled by an exposure factor, and 2) the geometries of the spacecraft and the node modeled by a heating factor and convective area. Node failure is evaluated using criteria based on melting temperature, reference heat load, g-load, or a combination of the above. The failure of a liquid propellant tank is evaluated based on burnout flux from nucleate boiling or excess internal pressure. Following a component failure, updates are made as needed to the spacecraft mass and aerodynamic properties, nodal exposure and heating factors, and nodal convective and conductive areas. This allows the trajectory to be propagated seamlessly in a single run, inclusive of the trajectories of components that have separated from the spacecraft. The node ablation simulates the decreasing mass and convective/reference areas, and variable heating factor. A built-in database provides the thermo-mechanical properties of For the purpose of performing safety analysis and risk assessment for a probable offnominal suborbital/orbital atmospheric reentry resulting in vehicle breakup, a synthesis of trajectory propagation coupled with thermal analysis and the evaluation of node failure is required to predict the sequence of events, the timeline, and the progressive demise of spacecraft components. To provide this capability, the Simulation for Prediction of Entry Article Demise (SPEAD) analysis tool was developed. This report discusses the capabilities, modeling, and validation of the SPEAD analysis tool. SPEAD is applicable for Earth or Mars, with the option for 3 or 6 degrees-of-freedom (DOF) trajectory propagation. The atmosphere and aerodynamics data are supplied in tables, for linear interpolation of up to 4 independent variables. The gravitation model can include up to 20 zonal harmonic coefficients. The modeling of a single motor is available and can be adapted to multiple motors. For thermal analysis, the aerodynamic radiative and free-molecular/continuum convective heating, black-body radiative cooling, conductive heat transfer between adjacent nodes, and node ablation are modeled. In a 6- DOF simulation, the local convective heating on a node is a function of Mach, angle-ofattack, and sideslip angle, and is dependent on 1) the location of the node in the spacecraft and its orientation to the flow modeled by an exposure factor, and 2) the geometries of the spacecraft and the node modeled by a heating factor and convective area. Node failure is evaluated using criteria based on melting temperature, reference heat load, g-load, or a combination of the above. The failure of a liquid propellant tank is evaluated based on burnout flux from nucleate boiling or excess internal pressure. Following a component failure, updates are made as needed to the spacecraft mass and aerodynamic properties, nodal exposure and heating factors, and nodal convective and conductive areas. This allows the trajectory to be propagated seamlessly in a single run, inclusive of the trajectories of components that have separated from the spacecraft. The node ablation simulates the decreasing mass and convective/reference areas, and variable heating factor. A built-in database provides the thermo-mechanical properties of

Ling, Lisa↗
NASA NTRS2016

Simulation of Satellite, Airborne and Terrestrial LiDAR with DART (I):Waveform Simulation with Quasi-Monte Carlo Ray Tracing

Light Detection And Ranging (LiDAR) provides unique data on the 3-D structure of atmosphere constituents and the Earth's surface. Simulating LiDAR returns for different laser technologies and Earth scenes is fundamental for evaluating and interpreting signal and noise in LiDAR data. Different types of models are capable of simulating LiDAR waveforms of Earth surfaces. Semi-empirical and geometric models can be imprecise because they rely on simplified simulations of Earth surfaces and light interaction mechanisms. On the other hand, Monte Carlo ray tracing (MCRT) models are potentially accurate but require long computational time. Here, we present a new LiDAR waveform simulation tool that is based on the introduction of a quasi-Monte Carlo ray tracing approach in the Discrete Anisotropic Radiative Transfer (DART) model. Two new approaches, the so-called "box method" and "Ray Carlo method", are implemented to provide robust and accurate simulations of LiDAR waveforms for any landscape, atmosphere and LiDAR sensor configuration (view direction, footprint size, pulse characteristics, etc.). The box method accelerates the selection of the scattering direction of a photon in the presence of scatterers with non-invertible phase function. The Ray Carlo method brings traditional ray-tracking into MCRT simulation, which makes computational time independent of LiDAR field of view (FOV) and reception solid angle. Both methods are fast enough for simulating multi-pulse acquisition. Sensitivity studies with various landscapes and atmosphere constituents are presented, and the simulated LiDAR signals compare favorably with their associated reflectance images and Laser Vegetation Imaging Sensor (LVIS) waveforms. The LiDAR module is fully integrated into DART, enabling more detailed simulations of LiDAR sensitivity to specific scene elements (e.g., atmospheric aerosols, leaf area, branches, or topography) and sensor configuration for airborne or satellite LiDAR sensors.

Gastellu-Etchegorry, Jean-Philippe↗
NASA NTRS1991

Modeling of Jupiter's millimeter wave emission utilizing laboratory measurements of ammonia (NH3) opacity

The techniques used to make laboratory measurements of the millimeter wave opacity from gaseous ammonia under simulated Jovian conditions are described. The results are applied to a radiative transfer model, which is used to compute Jupiter's millimeter wavelength emission. The absorptivity of gaseous NH3 is measured to reduce one of the largest uncertainties in modeling Jupiter's millimeter wave emission. Several other millimeter wave opacity sources are examined. New expressions are given for computing the absorptivity of NH3, H2O, cloud condensates, and pressure-induced absorption. Jupiter's reliable millimeter wavelength observations are compared with synthetic emission spectra which utilize these new absorption expressions.

Joiner, Joanna↗
NASA NTRS1995

Effects of energy release on near field flow structure of gas jets

The primary objective is to understand how buoyancy affects the structure of the shear layer, the development of fluid dynamic instabilities, and formation of the coherent structures in the near-nozzle regions of gas jets. The secondary objectives are to study the role of buoyancy in lifting and reattachment process of diffusion flames, to evaluate the scaling behavior of diffusion flames, and to aid development and/or validation of theoretical models by providing quantitative data in the absence of buoyancy. Fast reacting hydrogen or hydrogen-inert fuels are used to isolate the effects of buoyancy on fluid dynamics without masking the flame behavior by soot and radiative heat transfer. This choice of fuel also permits an evaluation of simulating low gravity in low pressure ground experiments because the similarity constraints are relaxed for the fast reacting, nonsooting diffusion flames. The diagnostics consists primarily of a color schlieren system coupled with computer generated rainbow filters, video recording, and image analysis. The project involves (1) drop tower experiments, (2) ground experiments, and (3) theoretical analysis.

Agrawal, Ajay K.↗
NASA NTRS2001

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Agudo, Ivan↗
NASA NTRS2011

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present new global calculations of X-ray spectra from fully relativistic magneto-hydrodynamic (MHO) simulations of black hole (BH) accretion disks. With a self consistent radiative transfer code including Compton scattering and returning radiation, we can reproduce the predominant spectral features seen in decades of X-ray observations of stellar-mass BHs: a broad thermal peak around 1 keV, power-law continuum up to >100 keV, and a relativistically broadened iron fluorescent line. By varying the mass accretion rate, different spectral states naturally emerge: thermal-dominant, steep power-law, and low/hard. In addition to the spectral features, we briefly discuss applications to X-ray timing and polarization.

Schnittman, Jeremy D.↗
NASA NTRS2022

Numerical Results for Polarized Light Scattering in a Spherical Atmosphere

We report numerical results for polarized light reflection from the top of a Rayleigh scattering spherical atmosphere with height-dependent single scattering albedo over a dark surface. Michael Mishchenko considered this scenario back in the 1990’s, for a plane-parallel atmosphere of unit optical thickness (OT), for which radiance errors arising from neglecting polarization reaches their highest values. To further extend Mishchenko’s results, we consider a value of OT = 0.25, for which the effect of atmospheric curvature is pronounced. New results are generated using three state-of-the art radiative transfer (RT) codes. These are: the MYSTIC and MCSSA models, which simulate light scattering in a true- spherical atmosphere using Monte Carlo methods; and the discrete ordinate code VLIDORT, operating with a new multiple-scatter spherical correction designed to deliver reasonable approximations to spherical-medium scattering. In this work, we report results for both single and multiple scattering; this will help to support the validation of existing and future polarized spherical RT codes, especially those using approximative methods to deal with sphericity.

S Korkin↗
DOE OSTI2024

Modeling Molecular Hydrogen in Low-metallicity Galaxies

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which includes radiative transfer, a nonequilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H 2 abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the H I–H 2 transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the ISM drops drastically at Z ≲ 0.2 Z ⊙ as the formation time of H 2 becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected ƒ H 2 measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy, the overall molecular gas mass is reproduced to better than a factor of ≲1.5 across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous ƒ H 2 models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes.

79 ASTRONOMY AND ASTROPHYSICS↗