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 595 records · Page 33

Faster spectral density calculation using energy moments

Abstract Accurate predictions of inclusive scattering cross sections in the linear response regime require efficient and controllable methods to calculate the spectral density in a strongly-correlated many-body system. In this work we reformulate the recently proposed Gaussian Integral Transform technique in terms of Fourier moments of the system Hamiltonian which can be computed efficiently on a quantum computer. One of the main advantages of this framework is that it allows for an important reduction of the computational cost by exploiting previous knowledge about the energy moments of the spectral density. For a simple model of medium mass nucleus like $$^{40}$$ 40 Ca and target energy resolution of 1 MeV we find an expected speed-up of $$\approx 125$$ ≈ 125 times for the calculation of the giant dipole response and of $$\approx 50$$ ≈ 50 times for the simulation of quasi-elastic electron scattering at typical momentum transfers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A laboratory investigation of the variability of cloud reflected radiance fields

A method to determine the radiative properties of complex cloud fields was developed. A Cloud field optical simulator (CFOS) was constructed to simulate the interaction of cloud fields with visible radiation. The CFOS was verified by comparing experimental results from it with calculations performed with a Monte Carlo radiative transfer model. A software library was developed to process, reduce, and display CFOS data. The CFSOS was utilized to study the reflected radiane patterns from simulated cloud fields.

Mckee, T. B.↗

Ablation and Heating During Atmospheric Entry and Its Effect on Airburst Risk

Large meteoroids and asteroids entering the atmosphere endure tremendous heating from the shock heated air, and thereby lose a significant fraction of their mass during atmospheric entry a process known as ablation. The predicted evolution of the asteroids mass as it passes through the atmosphere can affect both the predicted energy deposition profile relevant to an airburst event, or the residual mass that strikes the ground in the case of an impact event. This presentation is divided roughly into two parts. In the first part, an overview of traditional models for heat transfer and ablation that are historically used in the meteor physics community is presented, and the validity in the asteroid entry regime discussed. Sensitivity analyses performed using the recently developed Fragment-Cloud Model (FCM) will be presented which show illustrate the range of sizes and entry parameters for which the predicted asteroid threat is most sensitive to the models for ablation and heat transfer. The second part of the presentation shall focus on recent work done under NASAs Asteroid Threat Assessment Project (ATAP) to develop new models for heat transfer and ablation using high-fidelity numerical simulation in concert with state-of-the-art experiments. Coupled computational fluid dynamics (CFD)radiation transport simulations preformed using the state-of-the-art entry modeling tools at NASA show that, for large meteoroids and asteroids, there can significant attenuation of the heat transfer to the surface (95 in some cases) by the products of ablation. In addition to the heat transfer, new models for the material response and ablation of asteroidal material have been developed [cite]. In the current work, we present finding from recent novel experiments performed in the arc jet facility at NASA Ames, which allows us to, in part, simulate the extreme environment experienced by the asteroid during entry. Briefly, the experimental set-up was comprised of a 1.5 conical article of machined H5 chondrite, which was exposed to a high-enthalpy flow resulting in approximately 4 kWcm2 of heating to the surface. A still frame capture from high-speed video taken during this experiment can be seen in Figure 1. In this figure, we can observe some of the major mechanisms for meteoroid ablation, such as melt flow, spallation (mechanical removal of material), and vaporization. Major findings from this, and other experiments will be discussed, as well progress on utilizing the data from the experiments to inform and develop improved models for ablation.

Asteroids↗

Laboratory evaluation and application of microwave absorption properties under simulated conditions for planetary atmospheres

Laboratory measurements of microwave and millimeter wave properties of the simulated atmosphere of the outer planets and their satellites has continued. One of the focuses is on the development of a radiative transfer model of the Jovian atmosphere at wavelengths from 1 mm to 10 cm. This modeling effort led to laboratory measurements of the millimeter wave opacity of hydrogen sulfide (H2S) under simulated Jovian conditions. Descriptions of the modeling effort, the Laboratory experiment, and the observations are presented. Correlative studies of measurements with Pioneer-Venus radio occultation measurements with longer wavelength emission measurements have provided new ways for characterizing temporal and spatial variations in the abundance of both gases H2SO4 and SO2, and for modeling their roles in the subcloud atmosphere. Laboratory measurements were conducted on 1.35 cm (and 13 cm) opacity of gaseous SO2 and absorptivity of gaseous SO2 at the 3.2 mm wavelength under simulated Venus conditions. Laboratory measurements were completed on millimeter wave dielectric properties of liquid H2SO4, in order to model the effects of the opacity of the clouds of Venus onto millimeter wave emission spectrum.

Steffes, Paul G.↗

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗

CRTM Support to GMAO, Validation and Coefficient Generation

Radiative transfer (RT) models play a very critical role in assimilating satellite radiances into NWP models. The RT models are used as forward operator to simulate satellite radiances from atmopspheric control variables such as pressure, temperature, water vapor, and ozone. However because line-by-line RT models are computationally very expensive, fast RT models have been developed and advanced especially in past two decades to overcome these limitations. Community Radiative Transfer Model (CRTM) developed by Joint Center for Satellite Data Assimilation is widely used in the U.S. as the forward operator for the assimilation of microwave and infrared satellite radiances. This abstract summarizes the GMAO activities in the support of CRTM including generating training coefficients for new instruments as well as developments for assimilating satellite radiances from shortwave infrared channels.

Moradi, Isaac↗

Improving a Spectral Bin Microphysical Scheme Using TRMM Satellite Observations

Comparisons between cloud model simulations and observations are crucial in validating model performance and improving physical processes represented in the mod Tel.hese modeled physical processes are idealized representations and almost always have large rooms for improvements. In this study, we use data from two different sensors onboard TRMM (Tropical Rainfall Measurement Mission) satellite to improve the microphysical scheme in the Goddard Cumulus Ensemble (GCE) model. TRMM observed mature-stage squall lines during late spring, early summer in central US over a 9-year period are compiled and compared with a case simulation by GCE model. A unique aspect of the GCE model is that it has a state-of-the-art spectral bin microphysical scheme, which uses 33 different bins to represent particle size distribution of each of the seven hydrometeor species. A forward radiative transfer model calculates TRMM Precipitation Radar (PR) reflectivity and TRMM Microwave Imager (TMI) 85 GHz brightness temperatures from simulated particle size distributions. Comparisons between model outputs and observations reveal that the model overestimates sizes of snow/aggregates in the stratiform region of the squall line. After adjusting temperature-dependent collection coefficients among ice-phase particles, PR comparisons become good while TMI comparisons worsen. Further investigations show that the partitioning between graupel (a high-density form of aggregate), and snow (a low-density form of aggregate) needs to be adjusted in order to have good comparisons in both PR reflectivity and TMI brightness temperature. This study shows that long-term satellite observations, especially those with multiple sensors, can be very useful in constraining model microphysics. It is also the first study in validating and improving a sophisticated spectral bin microphysical scheme according to long-term satellite observations.

Li, Xiaowen↗

Deriving Soil Moisture with the Combined L-Band Radar and Radiometer Measurements

In this study, we develop a combined active/passive technique to estimate surface soil moisture with the focus on the short vegetated surfaces. We first simulated a database for both active and passive signals under SMAP s sensor configurations using the radiative transfer model with a wide range of conditions for surface soil moisture, roughness and vegetation properties that we considered as the random orientated disks and cylinders. Using this database, we developed 1) the techniques to estimate surface backscattering and emission components and 2) the technique to estimate soil moisture with the estimated surface backscattering and emission components. We will demonstrate these techniques with the model simulated data and its validation with the airborne PALS image data from the soil moisture SGP 99 and SMEX 02 experiments.

Shi, Jiancheng↗

Computational analysis of Variable Thrust Engine (VTE) performance

The Variable Thrust Engine (VTE) of the Orbital Maneuvering Vehicle (OMV) uses a hypergolic propellant combination of Monomethyl Hydrazine (MMH) and Nitrogen Tetroxide (NTO) as fuel and oxidizer, respectively. The performance of the VTE depends on a number of complex interacting phenomena such as atomization, spray dynamics, vaporization, turbulent mixing, convective/radiative heat transfer, and hypergolic combustion. This study involved the development of a comprehensive numerical methodology to facilitate detailed analysis of the VTE. An existing Computational Fluid Dynamics (CFD) code was extensively modified to include the following models: a two-liquid, two-phase Eulerian-Lagrangian spray model; a chemical equilibrium model; and a discrete ordinate radiation heat transfer model. The modified code was used to conduct a series of simulations to assess the effects of various physical phenomena and boundary conditions on the VTE performance. The details of the models and the results of the simulations are presented.

Giridharan, M. G.↗

Constraining the Properties of the Eta Carinae System via 3-D SPH Models of Space-Based Observations: The Absolute Orientation of the Binary Orbit

The extremely massive (> 90 Solar Mass) and luminous (= 5 x 10(exp 6) Solar Luminosity) star Eta Carinae, with its spectacular bipolar "Homunculus" nebula, comprises one of the most remarkable and intensely observed stellar systems in the galaxy. However, many of its underlying physical parameters remain a mystery. Multiwavelength variations observed to occur every 5.54 years are interpreted as being due to the collision of a massive wind from the primary star with the fast, less dense wind of a hot companion star in a highly elliptical (e approx. 0.9) orbit. Using three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the binary wind-wind collision in Eta Car, together with radiative transfer codes, we compute synthetic spectral images of [Fe III] emission line structures and compare them to existing Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) observations. We are thus able, for the first time, to constrain the absolute orientation of the binary orbit on the sky. An orbit with an inclination of i approx. 40deg, an argument of periapsis omega approx. 255deg, and a projected orbital axis with a position angle of approx. 312deg east of north provides the best fit to the observations, implying that the orbital axis is closely aligned in 3-1) space with the Homunculus symmetry axis, and that the companion star orbits clockwise on the sky relative to the primary.

Madura, Thomas I.↗

Radiation Coupling with the FUN3D Unstructured-Grid CFD Code

The HARA radiation code is fully-coupled to the FUN3D unstructured-grid CFD code for the purpose of simulating high-energy hypersonic flows. The radiation energy source terms and surface heat transfer, under the tangent slab approximation, are included within the fluid dynamic ow solver. The Fire II flight test, at the Mach-31 1643-second trajectory point, is used as a demonstration case. Comparisons are made with an existing structured-grid capability, the LAURA/HARA coupling. The radiative surface heat transfer rates from the present approach match the benchmark values within 6%. Although radiation coupling is the focus of the present work, convective surface heat transfer rates are also reported, and are seen to vary depending upon the choice of mesh connectivity and FUN3D ux reconstruction algorithm. On a tetrahedral-element mesh the convective heating matches the benchmark at the stagnation point, but under-predicts by 15% on the Fire II shoulder. Conversely, on a mixed-element mesh the convective heating over-predicts at the stagnation point by 20%, but matches the benchmark away from the stagnation region.

Wood, William A.↗

Influence of Ice Particle Surface Roughening on the Global Cloud Radiative Effect

Ice clouds influence the climate system by changing the radiation budget and large-scale circulation. Therefore, climate models need to have an accurate representation of ice clouds and their radiative effects. In this paper, new broadband parameterizations for ice cloud bulk scattering properties are developed for severely roughened ice particles. The parameterizations are based on a general habit mixture that includes nine habits (droxtals, hollow/solid columns, plates, solid/hollow bullet rosettes, aggregate of solid columns, and small/large aggregates of plates). The scattering properties for these individual habits incorporate recent advances in light-scattering computations. The influence of ice particle surface roughness on the ice cloud radiative effect is determined through simulations with the Fu-Liou and the GCM version of the Rapid Radiative Transfer Model (RRTMG) codes and the National Center for Atmospheric Research Community Atmosphere Model (CAM, version 5.1). The differences in shortwave (SW) and longwave (LW) radiative effect at both the top of the atmosphere and the surface are determined for smooth and severely roughened ice particles. While the influence of particle roughening on the single-scattering properties is negligible in the LW, the results indicate that ice crystal roughness can change the SW forcing locally by more than 10 W m(exp −2) over a range of effective diameters. The global-averaged SW cloud radiative effect due to ice particle surface roughness is estimated to be roughly 1-2 W m(exp −2). The CAM results indicate that ice particle roughening can result in a large regional SW radiative effect and a small but nonnegligible increase in the global LW cloud radiative effect.

clouds↗

3D Radiative Hydrodynamics Modeling of Convection of Stars to Probe Their Interiors and Photospheric Properties

The dramatic flow of data from the Kepler and K2 missions opens the opportunity to significantly improve our knowledge of stellar interiors, surface dynamics, and structure. However, interpretation of these observations is a challenging task because it depends on tiny effects that can be studied only with advanced first-principles modeling. We present results of 3D time-dependent radiative hydrodynamic simulations of stellar outer convection zones and atmospheres taking into account chemical composition, radiative transfer, turbulence effects, and a realistic equation of state for main sequence stars. We will discuss properties of convective structure and dynamics, convective overshoot, effects of magnetic fields and rotation, as well as the potential influence of turbulent surface dynamics on high-precision RV measurements.

Photospheric↗

Assessment of Retrieved GMI Emissivity Over Land, Snow and Sea Ice in the GEOS System

Measurements from microwave sounders and imagers provide a valuable source of information including atmospheric temperature and water vapor in Numerical Weather Prediction (NWP) systems that assimilate these observations directly over water surfaces (oceans and other large water bodies). In a recent decadal survey, targeted observables in the Planetary Boundary Layer (PBL) were cited as a key need for future observations (NASEM, 2018). Microwave observations which sense in the PBL are currently available, however, utilizing surface-sensitive microwave observations for atmospheric data assimilation remains a challenge over land, snow and sea ice. This is in part due to the inability of surface emissivity models used by NWP data assimilation systems to simulate observations with sufficient accuracy. The GEOS-ADAS (Todling and el Akkraoui, 2018) which utilizes the Community Radiative Transfer Model (CRTM) (Han, 2006; Chen 2009) is no exception. The ECMWF system has retrieved instantaneous surface emissivity from surface-sensitive channels for SSMI/S and MHS radiance observations, and apply these estimates to the closest channels higher in frequency (Baordo and Geer 2016) in the calculation of simulated radiances. This approach currently is also being tested in the GEOS-ADAS for AMSU-A and ATMS radiances (Zhu et al. 2021). No or minimal emissivity spectral variability has been assumed in the above-mentioned studies. Recently, work by Munchak et al., 2020 (hereby referred to as M2020) provided a new database for emissivity over land, snow and sea ice retrieved from the NASA Global Precipitation Mission (GPM). Compared with Tool to Estimate Land Surface Emissivities at Microwave (TELSEM2; Wang et al., 2017), M2020 provides emissivities for more frequencies(i.e., 10.7 GHz V/H). Moreover, this database is unique in that it utilizes both active and passive data to retrieve surface emissivity and normalized radar cross section. While the emissivity values may be useful for other sensors, they are most applicable to the GPM Microwave Imager (GMI). In this work the GEOS-ADAS is modified to utilize emissivity values from Munchak et al, 2020 in place of values used by CRTM. Presently, only GMI radiances over ocean are used in the operational GEOS-ADAS. This study will focus on the GMI radiances over land, snow, and ice, as a first attempt to evaluate GMI radiances over these non-water surface types. Two cases are then presented, one with one week of observation minus background departures using the modified GEOS-ADAS, and one utilizing the original GEOS-ADAS. It should be noted that the surface emissivity models in CRTM are not state of the art and are scheduled to be replaced by the Community Surface Emissivity Module (CSEM; Chen and Weng, 2016). Simulations using default CRTM emissivity values are used merely as reference comparing against M2020, and is not a thorough comparison against other more state of the art modules such as CSEM.

GMI↗

Space-Environment Simulation at the Lewis Research Center

A knowledge of the environmental conditions in space and an understanding of the possible effect on objects in that environment is necessary in order to make a full-scale, systematic exploration of outer space. Once these conditions have been determined, the understanding of the behavior of objects in space may be gained by simulating the environment in the laboratory and conducting experiments in the simulated environment. Among the conditions in space which have been determined and may be simulated are the pressure, the temperature, and the radiant energy. Quantitatively, simulating these conditions would require a pressure as low as 10 (exp -16) mm Hg, a sink temperature of approximately 30 to 40 K, and radiation similar to that from a 6,000 K black body. The features of the space environment directly affect the thermal balance of the spacecraft. The temperature of the object is that temperature which provides equilibrium between the energy absorbed by the object and the energy emitted by the object. Hence, a knowledge of the absorption and emission characteristics for materials under space environmental conditions must be obtained. High-speed particles may impact and erode the material surfaces and in that way alter the radiation properties of the materials and finally change the equilibrium temperature. Also, since most of the radiant energy (that from the sun) is from one direction, the-heat-transfer within the body will provide surface-temperature variations which depend upon the shape and orientation of the object and the properties of the material. These variations may be effectively studied in the simulated environment. Finally, because of the vacuum the heat transferred to and from the object must be accomplished solely by radiation. Thus, if heat-cycle power-generation equipment is to be used, the necessary rejection of heat must be accomplished by radiation from waste heat radiators whose design must be investigated in order to determine optimum configurations.

Nichols, Lester D.↗

Measurement of early time outer laser beam reflection inside a cylindrical hohlraum at the National Ignition Facility

In a National Ignition Facility (NIF) cylindrical hohlraum, any specular reflection (“glint”) off the wall from the outer cone (incident angle of 50° and 44° relative to the surface of the wall) laser beam will irradiate the capsule poles. If the glint power is sufficiently large during the picket (early time) of the laser pulse, it may seed high-mode perturbations on the capsule surface that can grow during the implosion. To quantify the glint power on the capsule during the picket by the outer beams, we performed dedicated experiments on NIF using a flat witness foil as a surrogate for the capsule in a half-hohlraum target. We found that the measured glint power is approximately comparable to simulations using a nominal electron conduction flux limiter of f = 0.15, and over an order of magnitude lower than that predicted using f = 0.03 in the wall. Based on our current understanding, we conclude that the glint from the outer beams plays an insignificant role in capsule drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiative Heat Transport and Optical Characterization of High Temperature Molten Salts

The goal of this project was to advance the field of molten fluoride radiative heat transfer (RHT), through optical property measurements, heat transfer loop construction and operation, and computational fluid dynamics simulations. The following outcomes have been accomplished. A comprehensive review has been conducted to compile all available RHT-relevant optical absorption data for molten fluorides. The various measurement methods have also been compiled and critiqued. Thermal emission measurement methods have been refined and documented. A unique thermal emissivity measurement facility has been constructed and used to measure the emissivity of flowing FLiBe-exposed SS316 samples. A slight increase has been observed in the emissivity compared to unexposed SS316. A dual reflectance-emission setup was partially demonstrated for the purpose of molten fluoride absorption measurement. Preliminary infrared edge data for FLiNaK was obtained. A visible near-infrared (NIR) measurement setup has been constructed and will be used to measure absorption of corrosion products in FLiBe. A new multi-band thermal radiation method using blackbody-weighted spectral bands for the absorption-coefficient was deployed to investigate the thermal performance of the UW Natural Circulation FLiBe Loop under various power levels and operating conditions. The results showed an increasing and decreasing trend in the overall heat transfer with optical thickness with a maximum at optical thickness τ D ≈Ο(1). A forced convection FLiBe loop was constructed for the purpose of studying corrosion and thermal hydraulic phenomena. A vertical cantilever pump is used to achieve up to 5 gpm and 5 m/s flow velocity in hot and cold test sections. A methodology has been established for modeling radiative heat transfer (RHT) in molten salt coolants for the purpose of quantifying the contribution of thermal radiation in conditions relevant to FHR or MSR systems. The comparison between the results obtained for forced and mixed convection showed that in the studied ranges of optical thickness, the maximum effect of RHT on overall heat transfer was 60% improvement observed at z/D = 250 for forced convection; whereas for mixed convection, that value was higher at 170% observed at z/D = 30. The highest effect of RHT on the overall energy distribution was observed at τ D ≈3.4. Negligible effect (<20%) was observed above τ D =18 and below τ D ≈0.002. A CFD code-to-code comparison was conducted using COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT, for a forced and mixed convection laminar flow for a grey medium in a vertical heated using both COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT. The relative error was the largest at τ D =40 for the radiative heat flux, bulk temperature, and Nusselt number. On the other hand, the relative error was the least at τ D =0.002 for the conductive heat flux, bulk temperature, and Nusselt number.

Anderson, Mark↗

A numerical study of a marine subtropical stratus cloud layer and its stability

A model for numerical simulation of stratus cloud layers is constructed by combining a second-order closure, turbulent transfer model with a thermal radiative transfer model. The turbulent transfer model allows water vapor saturation. The combined turbulence-radiation model is applied to both a horizontally uniform one-dimensional case and a horizontally nonuniform two-dimensional case. In the latter, the dynamics of mesoscale circulations are also incorporated. Results of the two-dimensional simulation show that the layer cloud instability occurs where the sea surface temperature is high and the large-scale subsidence is weak. The simulated instability is analyzed in view of an instability criterion, the eddy kinetic energy budget, and evaporative cooling near the cloud top.

Moeng, C.-H.↗