Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ray trace”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

DWARF GALAXIES WITH IONIZING RADIATION FEEDBACK. I. ESCAPE OF IONIZING PHOTONS

We describe a new method for simulating ionizing radiation and supernova feedback in the analogs of low-redshift galactic disks. In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self-regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. By simulating a galactic disk in a halo of 2.3 × 10{sup 11} M{sub ☉}, we find that the average escape fraction from all radiating sources on the spiral arms (excluding the central 2.5 kpc) fluctuates between 0.08% and 5.9% during a ∼20 Myr period with a mean value of 1.1%. The flux of escaped photons from these sources is not strongly beamed, but manifests a large opening angle of more than 60° from the galactic pole. Further, we investigate the escape fraction per SFMC particle, f{sub esc}(i), and how it evolves as the particle ages. We discover that the average escape fraction f{sub esc} is dominated by a small number of SFMC particles with high f{sub esc}(i). On average, the escape fraction from an SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle lifetime, 6 Myr. This is because SFMC particles drift away from the dense gas clumps in which they were born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova feedback. The framework established in this study brings deeper insight into the physics of photon escape fraction from an individual star-forming clump and from a galactic disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Diffuse ionizing radiation within HH jets

We present numerical hydrodynamical simulations of a time-dependent ejection velocity precessing jet. The parameters used in our models correspond to a high excitation Herbig-Haro object, such as HH 80/81. We have included the transfer of ionizing radiation produced within the shocked regions of the jet. The radiative transfer is computed with a ray-tracing scheme from all the cells with an emissivity above a certain threshold. We show the development of a radiative precursor, and compare the morphology with a model without the diffuse radiation. Our simulations show that the morphology of the Hα emission is affected considerably if the diffuse ionizing radiation is accounted for. The predicted Hα position-velocity diagram (i.e., spatially resolved emission line profiles) from a model with the transfer of ionizing radiation has a relatively strong component at zero velocity, corresponding to the radiative precursor. Qualitatively similar 'zero velocity components' are observed in HH 80/81 and in the jet from Sanduleak's star in the Large Magellanic Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Improved reflection models of black hole accretion disks: Treating the angular distribution of X-rays

X-ray reflection models are used to constrain the properties of the accretion disk, such as the degree of ionization of the gas and the elemental abundances. In combination with general relativistic ray tracing codes, additional parameters like the spin of the black hole and the inclination to the system can be determined. However, current reflection models used for such studies only provide angle-averaged solutions for the flux reflected at the surface of the disk. Moreover, the emission angle of the photons changes over the disk due to relativistic light bending. To overcome this simplification, we have constructed an angle-dependent reflection model with the XILLVER code and self-consistently connected it with the relativistic blurring code RELLINE. The new model, relxill, calculates the proper emission angle of the radiation at each point on the accretion disk and then takes the corresponding reflection spectrum into account. We show that the reflected spectra from illuminated disks follow a limb-brightening law highly dependent on the ionization of disk and yet different from the commonly assumed form I∝ln (1 + 1/μ). A detailed comparison with the angle-averaged model is carried out in order to determine the bias in the parameters obtained by fitting a typical relativistic reflection spectrum. These simulations reveal that although the spin and inclination are mildly affected, the Fe abundance can be overestimated by up to a factor of two when derived from angle-averaged models. The fit of the new model to the Suzaku observation of the Seyfert galaxy Ark 120 clearly shows a significant improvement in the constraint of the physical parameters, in particular by enhancing the accuracy in the inclination angle and the spin determinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A three-dimensional model of Tangential YORP

Tangential YORP, or TYORP, has recently been demonstrated to be an important factor in the evolution of an asteroid's rotation state. It is complementary to normal YORP, or NYORP, which used to be considered previously. While NYORP is produced by non-symmetry in the large-scale geometry of an asteroid, TYORP is due to heat conductivity in stones on the surface of the asteroid. To date, TYORP has been studied only in a simplified one-dimensional model, substituting stones with high long walls. This article for the first time considers TYORP in a realistic three-dimensional model, also including shadowing and self-illumination effects via ray tracing. TYORP is simulated for spherical stones lying on regolith. The model includes only five free parameters and the dependence of the TYORP on each of them is studied. The TYORP torque appears to be smaller than previous estimates from the one-dimensional model, but is still comparable to the NYORP torques. These results can be used to estimate TYORP of different asteroids and also as a basis for more sophisticated models of TYORP.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The direct collapse of a massive black hole seed under the influence of an anisotropic Lyman-Werner source

The direct collapse model of supermassive black hole seed formation requires that the gas cools predominantly via atomic hydrogen. To this end we simulate the effect of an anisotropic radiation source on the collapse of a halo at high redshift. The radiation source is placed at a distance of 3 kpc (physical) from the collapsing object and is set to emit monochromatically in the center of the Lyman-Werner (LW) band. The LW radiation emitted from the high redshift source is followed self-consistently using ray tracing techniques. Due to self-shielding, a small amount of H{sub 2} is able to form at the very center of the collapsing halo even under very strong LW radiation. Furthermore, we find that a radiation source, emitting >10{sup 54} (∼ 10{sup 3} J{sub 21}) photons s{sup –1}, is required to cause the collapse of a clump of M ∼ 10{sup 5} M {sub ☉}. The resulting accretion rate onto the collapsing object is ∼0.25 M {sub ☉} yr{sup –1}. Our results display significant differences, compared to the isotropic radiation field case, in terms of the H{sub 2} fraction at an equivalent radius. These differences will significantly affect the dynamics of the collapse. With the inclusion of a strong anisotropic radiation source, the final mass of the collapsing object is found to be M ∼ 10{sup 5} M {sub ☉}. This is consistent with predictions for the formation of a supermassive star or quasi-star leading to a supermassive black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THE POWER OF IMAGING: CONSTRAINING THE PLASMA PROPERTIES OF GRMHD SIMULATIONS USING EHT OBSERVATIONS OF Sgr A*

Recent advances in general relativistic magnetohydrodynamic simulations have expanded and improved our understanding of the dynamics of black-hole accretion disks. However, current simulations do not capture the thermodynamics of electrons in the low density accreting plasma. This poses a significant challenge in predicting accretion flow images and spectra from first principles. Because of this, simplified emission models have often been used, with widely different configurations (e.g., disk- versus jet-dominated emission), and were able to account for the observed spectral properties of accreting black holes. Exploring the large parameter space introduced by such models, however, requires significant computational power that exceeds conventional computational facilities. In this paper, we use GRay, a fast graphics processing unit (GPU) based ray-tracing algorithm, on the GPU cluster El Gato, to compute images and spectra for a set of six general relativistic magnetohydrodynamic simulations with different magnetic field configurations and black-hole spins. We also employ two different parametric models for the plasma thermodynamics in each of the simulations. We show that, if only the spectral properties of Sgr A* are used, all 12 models tested here can fit the spectra equally well. However, when combined with the measurement of the image size of the emission using the Event Horizon Telescope, current observations rule out all models with strong funnel emission, because the funnels are typically very extended. Our study shows that images of accretion flows with horizon-scale resolution offer a powerful tool in understanding accretion flows around black holes and their thermodynamic properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Impact of baryonic processes on weak-lensing cosmology: Power spectrum, nonlocal statistics, and parameter bias

We study the impact of baryonic physics on cosmological parameter estimation with weak-lensing surveys. We run a set of cosmological hydrodynamics simulations with different galaxy formation models. We then perform ray-tracing simulations through the total matter density field to generate 100 independent convergence maps with a field of view of 25 deg{sup 2}, and we use them to examine the ability of the following three lensing statistics as cosmological probes: power spectrum (PS), peak counts, and Minkowski functionals (MFs). For the upcoming wide-field observations, such as the Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of 1400 deg{sup 2}, these three statistics provide tight constraints on the matter density, density fluctuation amplitude, and dark energy equation of state, but parameter bias is induced by baryonic processes such as gas cooling and stellar feedback. When we use PS, peak counts, and MFs, the magnitude of relative bias in the dark energy equation of state parameter w is at a level of, respectively, δw∼0.017, 0.061, and 0.0011. For the HSC survey, these values are smaller than the statistical errors estimated from Fisher analysis. The bias could be significant when the statistical errors become small in future observations with a much larger survey area. We find that the bias is induced in different directions in the parameter space depending on the statistics employed. While the two-point statistic, i.e., PS, yields robust results against baryonic effects, the overall constraining power is weak compared with peak counts and MFs. On the other hand, using one of peak counts or MFs, or combined analysis with multiple statistics, results in a biased parameter estimate. The bias can be as large as 1σ for the HSC survey and will be more significant for upcoming wider-area surveys. We suggest to use an optimized combination so that the baryonic effects on parameter estimation are mitigated. Such a “calibrated” combination can place stringent and robust constraints on cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Reapproaching the spin estimate of GX 339-4

We systematically reanalyze two previous observations of the black hole (BH) GX 339-4 in the very high and intermediate state taken with XMM-Newton and Suzaku. We utilize up-to-date data reduction procedures and implement the recently developed, self-consistent model for X-ray reflection and relativistic ray tracing, relxill. In the very high and intermediate state, the rate of accretion is high and thus the disk remains close to the innermost stable circular orbit. We require a common spin parameter and inclination when fitting the two observations since these parameters should remain constant across all states. This allows for the most accurate determination of the spin parameter of this galactic BH binary from fitting the Fe Kα emission line and provides a chance to test previous estimates. We find GX 339-4 to be consistent with a near maximally spinning BH with a spin parameter a{sub ∗}>0.97 with an inclination of 36° ± 4°. This spin value is consistent with previous high estimates for this object. Further, if the inner disk is aligned with the binary inclination, this modest inclination returns a high BH mass, but they need not be aligned. Additionally, we explore how the spin is correlated with the power of the jet emitted but find no correlation between the two.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

TIME-DEPENDENT MULTI-GROUP MULTI-DIMENSIONAL RELATIVISTIC RADIATIVE TRANSFER CODE BASED ON SPHERICAL HARMONIC DISCRETE ORDINATE METHOD

We develop a time-dependent, multi-group, multi-dimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids that are involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) which evaluates a source function including anisotropic scattering in spherical harmonics and implicitly solves the static radiative transfer equation with ray tracing in discrete ordinates. We implement treatments of time dependence, multi-frequency bins, Lorentz transformation, and elastic Thomson and inelastic Compton scattering to the publicly available SHDOM code. Our code adopts a mixed-frame approach; the source function is evaluated in the comoving frame, whereas the radiative transfer equation is solved in the laboratory frame. This implementation is validated using various test problems and comparisons with the results from a relativistic Monte Carlo code. These validations confirm that the code correctly calculates the intensity and its evolution in the computational domain. The code enables us to obtain an Eddington tensor that relates the first and third moments of intensity (energy density and radiation pressure) and is frequently used as a closure relation in radiation hydrodynamics calculations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Modelling of the EAST lower-hybrid current drive experiment using GENRAY/CQL3D and TORLH/CQL3D

The coupled GENRAY-CQL3D code has been used to do systematic ray-tracing and Fokker–Planck analysis for EAST Lower Hybrid wave Current Drive (LHCD) experiments. Despite being in the weak absorption regime, the experimental level of LH current drive is successfully simulated, by taking into account the variations in the parallel wavenumber due to the toroidal effect. The effect of radial transport of the fast LH electrons in EAST has also been studied, which shows that a modest amount of radial transport diffusion can redistribute the fast LH current significantly. Taking advantage of the new capability in GENRAY, the actual Scrape Off Layer (SOL) model with magnetic field, density, temperature, and geometry is included in the simulation for both the lower and the higher density cases, so that the collisional losses of Lower Hybrid Wave (LHW) power in the SOL has been accounted for, which together with fast electron losses can reproduce the LHCD experimental observations in different discharges of EAST. We have also analyzed EAST discharges where there is a significant ohmic contribution to the total current, and good agreement with experiment in terms of total current has been obtained. Also, the full-wave code TORLH has been used for the simulation of the LH physics in the EAST, including full-wave effects such as diffraction and focusing which may also play an important role in bridging the spectral gap. The comparisons between the GENRAY and the TORLH codes are done for both the Maxwellian and the quasi-linear electron Landau damping cases. In conclusion, these simulations represent an important addition to the validation studies of the GENRAY-CQL3D and TORLH models being used in weak absorption scenarios of tokamaks with large aspect ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A novel approach to neutron scattering instrumentation for probing multiscale dynamics in soft and biological matter

We present a concept and ray-tracing simulation results of a mechanical device that will enable inelastic neutron scattering measurements where the data at energy transfers from a few eV to several hundred meV can be collected in a single, gapless spectrum. Besides covering 5 orders of magnitude on the energy (time) scale, the device provides data over 2 orders of magnitude on the scattering momentum (length) scale in a single measurement. Such capabilities are geared primarily toward soft and biological matter, where the broad dynamical features of relaxation origin largely overlap with vibration features, thus necessitating gapless spectral coverage over several orders of magnitude in time and space. Furthermore, neutron scattering experiments with such a device are performed with a fixed neutron final energy, which enables measurements, with neutron energy loss in the sample, at arbitrarily low temperatures over the same broad spectral range. Lastly, this capability is also invaluable in biological and soft matter research, as the variable temperature dependence of different relaxation components allows their separation in the scattering spectra as a function of temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Development of a tomographic reconstruction method for axisymmetric D α emission profiles in the ITER plasma boundary

A tomographic method for reconstructing the axially symmetric D α emission profile in the ITER scrape-off layer and divertor from D α intensity measurements in the fields of view (FoV) of Vis/IR TV and Divertor Impurity Monitor diagnostics is under development. The method takes into account the strong background signal due to stray light from divertor produced by reflections from the metal wall of the vacuum chamber. The method allows to filter out the reflections in D α signals measured in the main chamber. Here, the accuracy of this method is tested within the framework of a synthetic diagnostic, which uses the Raysect and Cherab numerical codes for ray-tracing simulations. The D α emission profiles, simulated with the SOLPS and OEDGE codes, are used as synthetic experimental data. Uncertainty in light reflection properties of the first wall (FW) is the major factor affecting the accuracy of the D α emission reconstruction. The dependence on the uncertainty in the FW surface roughness and on the ratio of diffusive to specular reflection of the error in recovering useful D α signals (the direct light without reflections) in the FoV of the ITER main chamber H-alpha and Visible Spectroscopy Diagnostic is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Power transport efficiency during O-X-B 2nd harmonic electron cyclotron heating in a helicon linear plasma device

The principal objective of this work is to report on the power coupled to a tungsten target in the prototype-material plasma exposure experiment device during oblique injection of a microwave beam (<70 kW at 28 GHz) into a high-power (~100 kW at 13.56 MHz) over-dense (${n_e} > 1 \times {10^{19}}{\text{ }}{{\text{m}}^{ - 3}}$) deuterium helicon plasma column. The experimental setup, electron heating system, electron heating scheme, and IR thermographic diagnostic for quantifying the power transport is described in detail. It is demonstrated that the power transported to the target can be effectively controlled by adjusting the magnetic field profile. Using this method, heat fluxes up to 22 MW m -2 and power transport efficiencies in the range of 17%–20% have been achieved using 70 kW of microwave power. It is observed that most of the heat flux is confined to a narrow region at the plasma periphery. Ray-tracing calculations are presented which indicate that the power is coupled to the plasma electrons via an O-X-B mode conversion process. Calculations indicate that the microwave power is absorbed in a single pass at the plasma periphery via collisions and in the over-dense region via 2nd harmonic cyclotron resonance of the electron Bernstein wave. The impact of our results is discussed in the context of MPEX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Optimization of the N ∥ Upshift in the DIII-D high field side lower hybrid current drive experiment

High field side lower hybrid current drive (LHCD) is one potential candidate for efficient non-inductive current drive in tokamak power plants, and the first test of this technology will occur on the DIII-D tokamak during the 2024 campaign. Previous LFS launch experiments operated in the multi-pass regime and relied on scrape-off layer interactions to close the spectral gap. In the DIII-D experiment, single-pass damping is achievable via an upshift in the parallel refractive index N ∥ caused by mode converting twice (slow → fast → slow). This mode conversion affects the ray trajectories and can lead to enhanced N ∥ upshift depending on where mode conversion occurs. Compared to multi-pass absorption experiments, the optimization of launched N ∥ and plasma parameters can be counter-intuitive: increased density may increase efficiency and smaller N ∥,launch tend to damp closer to the separatrix. A hard x-ray camera installed to measure the bremsstrahlung (50–250 keV) radiation from LHCD-generated fast electrons is capable of verifying the trends reporting in this paper through comparison to the ray-tracing/Fokker–Planck codes GENRAY/CQL3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

TorbeamNN: machine learning-based steering of ECH mirrors on KSTAR

We have developed TorbeamNN: a machine learning surrogate model for the TORBEAM ray tracing code to predict electron cyclotron heating (ECH) and current drive locations in tokamak plasmas. TorbeamNN provides more than a 100 times speed-up compared to the highly optimized and simplified real-time implementation of TORBEAM without any reduction in accuracy compared to the offline, full fidelity TORBEAM code. The model was trained using KSTAR ECH mirror geometries and works for both O-mode and X-mode absorption. The TorbeamNN predictions have been validated both offline and real-time in experiment. TorbeamNN has been utilized to track an ECH absorption vertical position target in dynamic KSTAR plasmas as well as under varying toroidal mirror angles and with a minimal average tracking error of 0.5 cm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Parametric raytracing modeling for NSTX-U scenario development with high harmonic fast waves and neutral beam injection

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker–Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k φ phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Imaging the flow of holes from a collimating contact in graphene

A beam of holes formed in graphene by a collimating contact is imaged using a liquid-He cooled scanning probe microscope (SPM). The mean free path of holes is greater than the device dimensions. A zigzag shaped pattern on both sides of the collimating contact absorbs holes that enter at large angles. The image charge beneath the SPM tip defects holes, and the pattern of flow is imaged by displaying the change in conductance between contacts on opposite sides, as the tip is raster scanned across the sample. Collimation is confirmed by bending hole trajectories away from the receiving contact with an applied magnetic field. The SPM images agree well with ray-tracing simulations.

36 MATERIALS SCIENCE↗