Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ray-tracing simulation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗
NASA NTRS2020

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

We present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature there by integrating over the thermal seed photon flux from the disk surface. The method either assumes ion and electron tem- peratures are equal (1T), or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to per- form. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. Both 1T and 2T approaches lead to increased radiative efficiency and a larger fraction of total cooling in the corona as compared to the original target-temperature cooling function used by harm3d, especially in the 1T case. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

Brooks E. Kinch↗
NASA NTRS1988

Structure of scintillations in Neptune's occultation shadow

An exceptionally high-quality data set from a Neptune occultation is used here to derive a number of new results about the statistical properties of the fluctuations of the intensity distribution in various parts of Neptune's occultation shadow. An approximate numerical ray-tracing model which successfully accounts for many of the qualitative aspects of the observed intensity fluctuation distribution is introduced. Strong refractive scintillation is simulated by including the effects of 'turbulence' with projected atmospheric properties allowed to vary in both the direction perpendicular and parallel to the limb, and an explicit two-dimensional picture of a typical intensity distribution throughout an occulting planet's shadow is presented. The results confirm the existence of highly anisotropic turbulence.

Hubbard, W. B.↗
NASA NTRS1993

CAD-based stand-alone spacecraft radiation exposure analysis system: An application of the early man-tended Space Station

Major improvements have been completed in the approach to analyses and simulation of spacecraft radiation shielding and exposure. A computer-aided design (CAD)-based system has been developed for determining the amount of shielding provided by a spacecraft and simulating transmission of an incident radiation environment to any point within or external to the vehicle. Shielding analysis is performed using a customized ray-tracing subroutine contained within a standard engineering modeling software package. This improved shielding analysis technique has been used in several vehicle design programs such as a Mars transfer habitat, pressurized lunar rover, and the redesigned international Space Station. Results of analysis performed for the Space Station astronaut exposure assessment are provided to demonastrate the applicability and versatility of the system.

Appleby, M. H.↗
NASA NTRS1993

Electromagnetic field strength prediction in an urban environment: A useful tool for the planning of LMSS

A model for the prediction of the electromagnetic field strength in an urban environment is presented. The ray model, that is based on the Uniform Theory of Diffraction (UTD), includes effects of the non-perfect conductivity of the obstacles and their surface roughness. The urban environment is transformed into a list of standardized obstacles that have various shapes and material properties. The model is capable of accurately predicting the field strength in the urban environment by calculating different types of wave contributions such as reflected, edge and corner diffracted waves, and combinations thereof. Also, antenna weight functions are introduced to simulate the spatial filtering by the mobile antenna. Communication channel parameters such as signal fading, time delay profiles, Doppler shifts and delay-Doppler spectra can be derived from the ray-tracing procedure using post-processing routines. The model has been tested against results from scaled measurements at 50 GHz and proves to be accurate.

Vandooren, G. A. J.↗
NASA NTRS1990

Optical modeling for dynamics and control analysis

This paper presents a coordinate-free ray-trace analysis of optical beam trains consisting of mirrors, lenses and reference surfaces. The analysis leads to optical models that can be directly integrated with standard structure and control models for integrated instrument design, analysis and simulation. This capability is required for a coming generation of spaceborn optical instruments which use controlled optical elements supported by flexible structures. New results include analytic formulas for optical sensitivities as functions of structural and control geometric parameters. An example problem calculates a linear (small-motion) optical model, derives optical gains for small- and large-angle controlled mirrors, and shows the propagation of a wavefront.

Redding, David C.↗
NASA NTRS2000

Ion Transmission to the Focal Plane of the Chandra X-Ray Observatory

The severe degradation in spectrometric performance observed by the AXAF CCD Imaging Spectrometer (ACIS) instrument on-board the Chandra X-ray Observatory (CXO) was associated with a sharp increase in charge transfer inefficiency combined with relatively constant dark current. It was also observed that damage occurred only during passage through the earth's radiation belts, and only when ACIS remained in the focal plane during the passage. Subsequent measurements and analyses support the conjecture that the damaging radiation entered through the Observatory's High-Resolution Mirror Assembly (HRMA) aperture. A mechanism whereby low energy magnetospheric protons and/or heavier ions are scattered through the HRMA and reach the focal plane with just enough energy to stop in the CCD's charge transfer region currently provides a reasonably consistent explanation of all observed phenomena. In this paper, we will describe analyses which support this conclusion. Simulations using standard ion transmission codes were used to generate a bi-directional response function (BDRF) for the mirror surfaces and transmission through the various path elements. The BDRF was convolved with the geometry via ray-trace. Damage estimates using measured proton fluence and ground measurements of ACIS- type CCD damage vs. proton energy will be presented and compared with observed on- orbit damage.

Kolodziejczak, J. J.↗
NASA NTRS2014

Features of Afterbody Radiative Heating for Earth Entry

Radiative heating is identified as a major contributor to afterbody heating for Earth entry capsules at velocities above 10 km/s. Because of rate-limited electron-ion recombination processes, a large fraction of the electronically-excited N and O atoms produced in the high temperature/pressure forebody remain as they expand into the afterbody region, which results in significant afterbody radiation. Large radiative heating sensitivities to electron-impact ionization rates and escape factors are identified. Ablation products from a forebody ablator are shown to increase the afterbody radiation by as much as 40%. The tangent-slab radiation transport approach is shown to over-predict the radiative flux by as much as 40% in the afterbody, therefore making the more computationally expensive ray-tracing approach necessary for accurate radiative flux predictions. For the Stardust entry, the afterbody radiation is predicted to be nearly twice as large as the convective heating during the peak heating phase of the trajectory. Comparisons between simulations and the Stardust Echelle observation measurements, which are shown to be dominated by afterbody emission, indicate agreement within 20% for various N and O lines. Similarly, calorimeter measurements from the Fire II experiment are identified as a source of validation data for afterbody radiation. For the afterbody calorimeter measurement closest to the forebody, which experiences the largest afterbody radiative heating component, the convective heating alone is shown to under-predict the measurement, even for the fullycatalytic assumption. Agreement with the measurements is improved with the addition of afterbody radiation. These comparisons with Stardust and Fire II measurements provide validation that the significant afterbody radiation values proposed in this work are legitimate.

Johnston, Christopher O.↗
NASA NTRS2011

A Deterministic Electron, Photon, Proton and Heavy Ion Radiation Transport Suite for the Study of the Jovian System

A deterministic suite of radiation transport codes, developed at NASA Langley Research Center (LaRC), which describe the transport of electrons, photons, protons, and heavy ions in condensed media is used to simulate exposures from spectral distributions typical of electrons, protons and carbon-oxygen-sulfur (C-O-S) trapped heavy ions in the Jovian radiation environment. The particle transport suite consists of a coupled electron and photon deterministic transport algorithm (CEPTRN) and a coupled light particle and heavy ion deterministic transport algorithm (HZETRN). The primary purpose for the development of the transport suite is to provide a means for the spacecraft design community to rapidly perform numerous repetitive calculations essential for electron, proton and heavy ion radiation exposure assessments in complex space structures. In this paper, the radiation environment of the Galilean satellite Europa is used as a representative boundary condition to show the capabilities of the transport suite. While the transport suite can directly access the output electron spectra of the Jovian environment as generated by the Jet Propulsion Laboratory (JPL) Galileo Interim Radiation Electron (GIRE) model of 2003; for the sake of relevance to the upcoming Europa Jupiter System Mission (EJSM), the 105 days at Europa mission fluence energy spectra provided by JPL is used to produce the corresponding dose-depth curve in silicon behind an aluminum shield of 100 mils ( 0.7 g/sq cm). The transport suite can also accept ray-traced thickness files from a computer-aided design (CAD) package and calculate the total ionizing dose (TID) at a specific target point. In that regard, using a low-fidelity CAD model of the Galileo probe, the transport suite was verified by comparing with Monte Carlo (MC) simulations for orbits JOI--J35 of the Galileo extended mission (1996-2001). For the upcoming EJSM mission with a potential launch date of 2020, the transport suite is used to compute the traditional aluminum-silicon dose-depth calculation as a standard shield-target combination output, as well as the shielding response of high charge (Z) shields such as tantalum (Ta). Finally, a shield optimization algorithm is used to guide the instrument designer with the choice of graded-Z shield analysis.

Norman, Ryan B.↗
NASA NTRS1997

Comet Shoemaker-Levy 9 Fragment Size Estimates: How Big was the Parent Body?

The impact of Comet Shoemaker-Levy 9 on Jupiter in July, 1994 was the largest, most energetic impact event on a planet ever witnessed. Because it broke up during a close encounter with Jupiter in 1992, it was bright enough to be discovered more than a year prior to impact, allowing the scientific community an unprecedented opportunity to assess the effects such an event would have. Many excellent observations were made from Earth-based telescopes, the Hubble Space Telescope (HST), and the Galileo spacecraft en route to Jupiter. In this paper, these observations are used in conjunction with computational simulations performed with the CTH shock-physics hydrocode to determine the sizes of the fifteen fragments that made discernible impact features on the planet. To do this, CTH was equipped with a radiative ablation model and a postprocessing radiative ray-trace capability that enabled light-flux predictions (often called the impact flash) for the viewing geometries of Galileo and ground-based observers. The five events recorded by Galileo were calibrated to give fragment size estimates. Compared against ground-based and HST observations, these estimates were extended using a least-squares analysis to assess the impacts of the remaining ten fragments. Some of the largest impacts (L, G, and K) were greater that 1 km in diameter, but the density of the fragments was low, about 0.25 g/cm(exp 3). The volume of the combined fifteen fragments would make a sphere 1.8 km in diameter. Assuming a prebreakup density of 0.5 g/cm(exp 3), the parent body of Shoemaker-Levy 9 had a probable diameter of 1.4 km. The total kinetic energy of all the impacts was equivalent to the explosive yield of 300 Gigatons of TNT.

Crawford, David A.↗
NASA NTRS2012

GERMcode: A Stochastic Model for Space Radiation Risk Assessment

A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and high charge and energy (HZE) particles that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of HZE particles in tissue and shielding materials is made with a stochastic approach that includes both particle track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections. For NSRL applications, the GERMcode evaluates a set of biophysical properties, such as the Poisson distribution of particles or delta-ray hits for a given cellular area and particle dose, the radial dose on tissue, and the frequency distribution of energy deposition in a DNA volume. By utilizing the ProE/Fishbowl ray-tracing analysis, the GERMcode will be used as a bi-directional radiation transport model for future spacecraft shielding analysis in support of Mars mission risk assessments. Recent radiobiological experiments suggest the need for new approaches to risk assessment that include time-dependent biological events due to the signaling times for activation and relaxation of biological processes in cells and tissue. Thus, the tracking of the temporal and spatial distribution of events in tissue is a major goal of the GERMcode in support of the simulation of biological processes important in GCR risk assessments. In order to validate our approach, basic radiobiological responses such as cell survival curves, mutation, chromosomal aberrations, and representative mouse tumor induction curves are implemented into the GERMcode. Extension of these descriptions to other endpoints related to non-targeted effects and biochemical pathway responses will be discussed.

Kim, Myung-Hee Y.↗
NASA NTRS2006

Optics Program Modified for Multithreaded Parallel Computing

A powerful high-performance computer program for simulating and analyzing adaptive and controlled optical systems has been developed by modifying the serial version of the Modeling and Analysis for Controlled Optical Systems (MACOS) program to impart capabilities for multithreaded parallel processing on computing systems ranging from supercomputers down to Symmetric Multiprocessing (SMP) personal computers. The modifications included the incorporation of OpenMP, a portable and widely supported application interface software, that can be used to explicitly add multithreaded parallelism to an application program under a shared-memory programming model. OpenMP was applied to parallelize ray-tracing calculations, one of the major computing components in MACOS. Multithreading is also used in the diffraction propagation of light in MACOS based on pthreads [POSIX Thread, (where "POSIX" signifies a portable operating system for UNIX)]. In tests of the parallelized version of MACOS, the speedup in ray-tracing calculations was found to be linear, or proportional to the number of processors, while the speedup in diffraction calculations ranged from 50 to 60 percent, depending on the type and number of processors. The parallelized version of MACOS is portable, and, to the user, its interface is basically the same as that of the original serial version of MACOS.

Lou, John↗
NASA NTRS2016

Modeling and Simulation of Radiative Compressible Flows in Aerodynamic Heating Arc-Jet Facility

Numerical simulations of an arc heated flow inside NASA's 20 [MW] Aerodynamics heating facility (AHF) are performed in order to investigate the three-dimensional swirling flow and the current distribution inside the wind tunnel. The plasma is considered in Local Thermodynamics Equilibrium(LTE) and is composed of Air-Argon gas mixture. The governing equations are the Navier-Stokes equations that include source terms corresponding to Joule heating and radiative cooling. The former is obtained by solving an electric potential equation, while the latter is calculated using an innovative massively parallel ray-tracing algorithm. The fully coupled system is closed by the thermodynamics relations and transport properties which are obtained from Chapman-Enskog method. A novel strategy was developed in order to enable the flow solver and the radiation calculation to be preformed independently and simultaneously using a different number of processors. Drastic reduction in the computational cost was achieved using this strategy. Details on the numerical methods used for space discretization, time integration and ray-tracing algorithm will be presented. The effect of the radiative cooling on the dynamics of the flow will be investigated. The complete set of equations were implemented within the COOLFluiD Framework. Fig. 1 shows the geometry of the Anode and part of the constrictor of the Aerodynamics heating facility (AHF). Fig. 2 shows the velocity field distribution along (x-y) plane and the streamline in (z-y) plane.

Pasma flows↗
NASA NTRS2016

Normalizing a Relativistic Model of X-Ray Reflection Definition of the Reflection Fraction and Its Implementation in relxill

Aims. The only relativistic reflection model that implements a parameter relating the intensity incident on an accretion disk to the observed intensity is relxill. The parameter used in earlier versions of this model, referred to as the reflection strength, is unsatisfactory; it has been superseded by a parameter that provides insight into the accretion geometry, namely the reflection fraction. The reflection fraction is defined as the ratio of the coronal intensity illuminating the disk to the coronal intensity that reaches the observer. Methods. The relxill model combines a general relativistic ray-tracing code and a photoionization code to compute the component of radiation reflected from an accretion that is illuminated by an external source. The reflection fraction is a particularly important parameter for relativistic models with well-defined geometry, such as the lamp post model, which is a focus of this paper. Results. Relativistic spectra are compared for three inclinations and for four values of the key parameter of the lamp post model,namely the height above the black hole of the illuminating, on-axis point source. In all cases, the strongest reflection is produced for low source heights and high spin. A low-spin black hole is shown to be incapable of producing enhanced relativistic reflection. Results for the relxill model are compared to those obtained with other models and a Monte Carlo simulation. Conclusions. Fitting data by using the relxill model and the recently implemented reflection fraction, the geometry of a system can be constrained. The reflection-fraction is independent of system parameters such as inclination and black hole spin. The reflection-fraction parameter was implemented with the name reflec_frac all flavours of the relxill model, and the non-relativistic reflection model xillver, in v0.4a (18 January 2016).

accretion↗
NASA NTRS2002

Correlations among the Optical Properties of Cirrus-Cloud Particles: Microphysical Interpretation

Cirrus measurements obtained with a ground-based polarization Raman lidar at 67.9 deg N in January 1997 reveal a strong positive correlation between the particle optical properties, specifically depolarization ratio delta(sub par) and extinction- to-backscatter (lidar) ratio S, for delta(sub par) less than approximately 40%, and an anti-correlation for delta(sub par) greater than approximately 40%. Over the length of the measurements the particle properties vary systematically. Initially, delta (sub par) approximately equals 60% and S approximately equals 10sr are observed. Then, with decreasing delta(sub par), S first increases to approximately 27sr (delta(sub par) approximately equals 40%) before decreasing to values around 10sr again (delta(sub par) approximately equals 20%). The analysis of lidar humidity and radiosonde temperature data shows that the measured optical properties stem from scattering by dry solid ice particles, while scattering by supercooled droplets, or by wetted or subliming ice particles can be excluded. For the microphysical interpretation of the lidar measurements, ray-tracing computations of particle scattering properties have been used. The comparison with the theoretical data suggests that the observed cirrus data can be interpreted in terms of size, shape, and, under the assumption that the lidar measurements of consecutive cloud segments can be mapped on the temporal development of a single cloud parcel moving along its trajectory, growth of the cirrus particles: Near the cloud top in the early stage of cirrus development, light scattering by nearly isometric particles that have the optical characteristics of hexagonal columns (short, column-like particles) is dominant. Over time the ice particles grow, and as the cloud base height extends to lower altitudes characterized by warmer temperatures they become morphologically diverse. For large S and depolarization values of approximately 40%, the scattering contributions of column- and plate-like particles are roughly the same. In the lower ranges of the cirrus clouds, light scattering is predominantly by plate-like ice particles. This interpretation assumes random orientation of the cirrus particles. Simulations with a simple model suggest, however, that the positive correlation between S and delta(sub par) which is observed for depolarization ratios less than 40% mainly at low cloud altitudes, can be alternatively explained by horizontal alignment of a fraction of the cirrus particle population.

Reichardt, J.↗
NASA NTRS2016

Advancements in Afterbody Radiative Heating Simulations for Earth Entry

Four advancements to the simulation of backshell radiative heating for Earth entry are presented. The first of these is the development of a flow field model that treats electronic levels of the dominant backshell radiator, N, as individual species. This is shown to allow improvements in the modeling of electron-ion recombination and two-temperature modeling, which are shown to increase backshell radiative heating by 10 to 40%. By computing the electronic state populations of N within the flow field solver, instead of through the quasi-steady state approximation in the radiation code, the coupling of radiative transition rates to the species continuity equations for the levels of N, including the impact of non-local absorption, becomes feasible. Implementation of this additional level of coupling between the flow field and radiation codes represents the second advancement presented in this work, which is shown to increase the backshell radiation by another 10 to 50%. The impact of radiative transition rates due to non-local absorption indicates the importance of accurate radiation transport in the relatively complex flow geometry of the backshell. This motivates the third advancement, which is the development of a ray-tracing radiation transport approach to compute the radiative transition rates and divergence of the radiative flux at every point for coupling to the flow field, therefore allowing the accuracy of the commonly applied tangent-slab approximation to be assessed for radiative source terms. For the sphere considered at lunar-return conditions, the tangent-slab approximation is shown to provide a sufficient level of accuracy for the radiative source terms, even for backshell cases. This is in contrast to the agreement between the two approaches for computing the radiative flux to the surface, which differ by up to 40%. The final advancement presented is the development of a nonequilibrium model for NO radiation, which provides significant backshell radiation at velocities below 10 km/s. The developed model reduces the nonequilibrium NO radiation by 50% relative to the previous model.

Johnston, Christopher O.↗
NASA NTRS2008

How Very Massive Metal-Free Stars Start Cosmological Reionization

The initial conditions and relevant physics for the formation of the earliest galaxies are well specified in the concordance cosmology. Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include nonequilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the case B approximation using adaptively ray-traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of approx.10(exp 6). These first sources of reionization are highly intermittent and anisotropic and first photoionize the small-scale voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf-sized galaxies, the escape fraction of UV radiation decreases and the H II regions only break out on some sides of the galaxies, making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately five ionizing photons are needed per sustained ionization when star formation in 10(exp 6) stellar Mass halos is dominant in the calculation. As the halos become larger than approx.10(exp 7) Stellar Mass, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Radiative feedback decreases clumping factors by 25% when compared to simulations without star formation and increases the average temperature of ionized gas to values between 3000 and 10,000 K.

Wise, John H.↗
NASA NTRS1990

Linearized ray-trace analysis

A new, coordinate-free version of the exact ray-trace equations for optical systems consisting of conic reflecting, refracting and reference surfaces is presented. These equations are differentiated to obtain closed-form optical sensitivity dyadics. For computation, the sensitivities are evaluated in a single global coordinate frame and combined in linearized ray-trace matrix difference equations that propagate the rays and the sensitivities from element to element. One purpose of this analysis is to create optical models that can be directly integrated with models of the instrument structure and control systems for dynamic simulation.

Redding, David C.↗