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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗
DOE OSTI2017

Linear feature detection algorithm for astronomical surveys – II. Defocusing effects on meteor tracks

Given the current limited knowledge of meteor plasma micro-physics and its interaction with the surrounding atmosphere and ionosphere, meteors are a highly interesting observational target for high-resolution wide-field astronomical surveys. Such surveys are capable of resolving the physical size of meteor plasma heads, but they produce large volumes of images that need to be automatically inspected for possible existence of long linear features produced by meteors. Here in this paper, we show how big aperture sky survey telescopes detect meteors as defocused tracks with a central brightness depression. We derive an analytic expression for a defocused point source meteor track and use it to calculate brightness profiles of meteors modelled as uniform brightness discs. We apply our modelling to meteor images as seen by the Sloan Digital Sky Survey and Large Synoptic Survey Telescope telescopes. The expression is validated by Monte Carlo ray-tracing simulations of photons travelling through the atmosphere and the Large Synoptic Survey Telescope telescope optics. We show that estimates of the meteor distance and size can be extracted from the measured full width at half-maximum and the strength of the central dip in the observed brightness profile. However, this extraction becomes difficult when the defocused meteor track is distorted by the atmospheric seeing or contaminated by a long-lasting glowing meteor trail. The full width at half-maximum of satellite tracks is distinctly narrower than meteor values, which enables removal of a possible confusion between satellites and meteors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Covariance matrices for galaxy cluster weak lensing: from virial regime to uncorrelated large-scale structure

Next-generation optical imaging surveys will revolutionize the observations of weak gravitational lensing by galaxy clusters and provide stringent constraints on growth of structure and cosmic acceleration. In these experiments, accurate modelling of covariance matrices of cluster weak lensing plays the key role in obtaining robust measurements of the mean mass of clusters and cosmological parameters. We use a combination of analytical calculations and high-resolution N -body simulations to derive accurate covariance matrices that span from the virial regime to linear scales of the cluster-matter cross-correlation. We validate this calculation using a public ray-tracing lensing simulation and provide a software package for calculating covariance matrices for a wide range of cluster and source sample choices. We discuss the relative importance of shape noise and density fluctuations, the impact of radial bin size, and the impact of off-diagonal elements. For a weak lensing source density n s = 10 arcmin -2 , shape noise typically dominates the variance on comoving scales $r_{\rm p}\lesssim 5\ h^{-1} \, \rm Mpc$. However, for n s = 60 arcmin -2 , potentially achievable with future weak lensing experiments, density fluctuations typically dominate the variance at $r_{\rm p}\gtrsim 1\ h^{-1} \, \rm Mpc$ and remain comparable to shape noise on smaller scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Optimizing injection for the storage ring proton-EDM experiment

The proposed proton electric dipole moment (pEDM) experiment at Brookhaven National Laboratory (BNL), to be built inside the alternating gradient synchrotron (AGS) tunnel, aims to measure the proton’s electric dipole moment with a sensitivity of 10 −29 𝑒 cm. This paper presents the design of the injection line from the AGS booster to the pEDM storage ring, utilizing portions of the existing booster-to-AGS (BtA) transfer line. Building on the symmetric-hybrid lattice design [, Comprehensive symmetric-hybrid ring design for a proton-EDM experiment at below 10 −29 𝑒 cm, Phys. Rev. D 105, 032001 (2022)], our study emphasizes rigorous optics matching, detailed particle and spin tracking, and systematic error mitigation essential for achieving a target sensitivity of 10 −29 𝑒 cm. This design preserves the proton’s vertical spin orientation within ±20 mrad, a critical requirement for the pEDM measurement. Particle and spin tracking simulations using the ray-tracing code Zgoubi [F. Méot, Zgoubi users’ guide, Technical Report, Brookhaven National Laboratory (BNL), Relativistic Heavy Ion Collider (RHIC), Upton, NY, 2012] validate the design’s performance, demonstrating its feasibility for this precision experiment. The simulation results demonstrate that both clockwise (CW) and counterclockwise (CCW) injection lines meet the stringent beam envelope and polarization requirements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Light atom quantum oscillations in UC and US

High energy vibrational scattering in the binary systems UC and US is measured using time-of-flight inelastic neutron scattering. A clear set of well-defined peaks equally separated in energy is observed in UC, corresponding to harmonic oscillations of the light C atoms in a cage of heavy U atoms. The scattering is much weaker in US and only a few oscillator peaks are visible. We show how the difference between the materials can be understood by considering the neutron scattering lengths and masses of the lighter atoms. Monte Carlo ray tracing is used to simulate the scattering, with near quantitative agreement with the data in UC, and some differences with US. The possibility of observing anharmonicity and anisotropy in the potentials of the light atoms is investigated in UC. Lastly, the observed data is well accounted for by considering each light atom as a single atom isotropic quantum harmonic oscillator.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2014

Impact of spurious shear on cosmological parameter estimates from weak lensing observables

We research, residual errors in shear measurements, after corrections for instrument systematics and atmospheric effects, can impact cosmological parameters derived from weak lensing observations. Here we combine convergence maps from our suite of ray-tracing simulations with random realizations of spurious shear. This allows us to quantify the errors and biases of the triplet (Ω m ,w,σ 8 ) derived from the power spectrum (PS), as well as from three different sets of non-Gaussian statistics of the lensing convergence field: Minkowski functionals (MFs), low-order moments (LMs), and peak counts (PKs). Our main results are as follows: (i) We find an order of magnitude smaller biases from the PS than in previous work. (ii) The PS and LM yield biases much smaller than the morphological statistics (MF, PK). (iii) For strictly Gaussian spurious shear with integrated amplitude as low as its current estimate of σ sys 2 ≈ 10 -7 , biases from the PS and LM would be unimportant even for a survey with the statistical power of Large Synoptic Survey Telescope. However, we find that for surveys larger than ≈ 100 deg 2 , non-Gaussianity in the noise (not included in our analysis) will likely be important and must be quantified to assess the biases. (iv) The morphological statistics (MF, PK) introduce important biases even for Gaussian noise, which must be corrected in large surveys. The biases are in different directions in (Ωm,w,σ8) parameter space, allowing self-calibration by combining multiple statistics. Our results warrant follow-up studies with more extensive lensing simulations and more accurate spurious shear estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg2 CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters Ω m , σ 8 , and w, and replicating the galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of ≤ 5%, and compute the likelihood in the three-dimensional parameter space (Ω m , σ 8 , w) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain w without external data. When the power spectrum and peak counts are combined, the area of the error “banana” in the (Ω m , σ 8 ) plane reduces by a factor of ≈ two, compared to using the power spectrum alone. For a flat Λ cold dark matter model, combining both statistics, we obtain the constraint σ 8 (Ω m /0.27)0.63 = 0.85 +0.03 -0.03 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Sample variance in weak lensing: How many simulations are required?

Constraining cosmology using weak gravitational lensing consists of comparing a measured feature vector of dimension N b with its simulated counterpart. An accurate estimate of the N b × N b feature covariance matrix C is essential to obtain accurate parameter confidence intervals. When C is measured from a set of simulations, an important question is how large this set should be. To answer this question, we construct different ensembles of N r realizations of the shear field, using a common randomization procedure that recycles the outputs from a smaller number N s ≤ N r of independent ray-tracing N-body simulations. We study parameter confidence intervals as a function of (N s , N r ) in the range 1 ≤ N s ≤ 200 and 1 ≤ N r ≲ 105. Previous work [S. Dodelson and M. D. Schneider, Phys. Rev. D 88, 063537 (2013)] has shown that Gaussian noise in the feature vectors (from which the covariance is estimated) lead, at quadratic order, to an O(1/N r ) degradation of the parameter confidence intervals. Using a variety of lensing features measured in our simulations, including shear-shear power spectra and peak counts, we show that cubic and quartic covariance fluctuations lead to additional O(1/N 2 r ) error degradation that is not negligible when N r is only a factor of few larger than N b . We study the large N r limit, and find that a single, 240 Mpc/h sized 512 3 -particle N-body simulation (N s = 1) can be repeatedly recycled to produce as many as N r = few × 10 4 shear maps whose power spectra and high-significance peak counts can be treated as statistically independent. Lastly, a small number of simulations (N s = 1 or 2) is sufficient to forecast parameter confidence intervals at percent accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Geometry and growth contributions to cosmic shear observables

We explore the sensitivity of weak-lensing observables to the expansion history of the Universe and to the growth of cosmic structures, as well as the relative contribution of both effects to constraining cosmological parameters. We utilize ray-tracing dark-matter-only N-body simulations and validate our technique by comparing our results for the convergence power spectrum with analytic results from past studies. We then extend our analysis to non-Gaussian observables which cannot be easily treated analytically. We study the convergence (equilateral) bispectrum and two topological observables, lensing peaks and Minkowski functionals, focusing on their sensitivity to the matter density Ω m and the dark energy equation of state w. We find that a cancellation between the geometry and growth effects is a common feature for all observables and exists at the map level. It weakens the overall sensitivity by factors of up to 3 and 1.5 for w and Ω m , respectively, with the bispectrum worst affected. However, combining geometry and growth information alleviates the degeneracy between Ω m and w from either effect alone. As a result, the magnitudes of marginalized errors remain similar to those obtained from growth-only effects, but with the correlation between the two parameters switching sign. Furthermore, these results shed light on the origin of the cosmology sensitivity of non-Gaussian statistics and should be useful in optimizing combinations of observables.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Lensing reconstruction in post-Born cosmic microwave background weak lensing

The study of the cosmic microwave background (CMB) lensing potential has established itself by now as a robust way of probing the physics of large-scale structure growth. The most common estimators of the lensing potential are derived under the assumption of Gaussianity of the matter distribution and in the Born approximation of the photon diffusion. In this paper we study the performance of quadratic estimators when applied to realistic sky maps extracted from multiple-lens ray tracing techniques in cosmological $N$-body simulations. These are expected to model accurately the effects due to the non-Gaussianity of the matter distribution induced by its nonlinearity and the deviation from the Born approximation. We show that both these effects on their own lead to reconstruction biases, but these tend to partially cancel each other when both these effects are considered together. We forecast the impact of these biases on the estimation of cosmological parameters for future high-sensitivity CMB experiments like CMB-S4. We find that the cold dark matter density, $Ω_{cdm}$, the optical depth to reionization $τ$, the amplitude of primordial inflationary perturbations, $A_s$, and the sum of neutrino masses, $M_ν$, could be biased at the $1–2σ$ level, if no external data set is used. We also observe a reduction of the bias if external data like baryon acoustic oscillations is included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Effect of non-Gaussian lensing deflections on CMB lensing measurements

We investigate the impact of non-Gaussian lensing deflections on measurements of the CMB lensing power spectrum. We find that the false assumption of their Gaussianity significantly biases these measurements in current and future experiments at the percent level. The bias is detected by comparing CMB lensing reconstructions from simulated CMB data lensed with Gaussian deflection fields to reconstructions from simulations lensed with fully non-Gaussian deflection fields. The non-Gaussian deflections are produced by ray-tracing through snapshots of an N-body simulation and capture both the non-Gaussianity induced by non-linear structure formation and by multiple correlated deflections. We find that the amplitude of the measured bias is in agreement with analytical predictions by Bohm et al. 2016. The bias is largest in temperature-based measurements and we do not find evidence for it in measurements from a combination of polarization fields (EB;EB). We argue that the non-Gaussian bias should be even more important for measurements of cross-correlations of CMB lensing with low-redshift tracers of large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive “full-wave” simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. Furthermore, the numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Expansion-Driven Self-Magnetization of High-Energy-Density Plasmas

Understanding plasma self-magnetization is one of the fundamental challenges in both laboratory and astrophysical plasmas. Self-magnetization can modify plasma transport properties, altering the dynamical evolution of plasmas. Multiple high-energy-density (HED) experiments have observed the formation of ion-scale magnetic filaments of megagauss strength, though their origin remains debated. Here, in this study, we conduct 2D collisional particle-in-cell (PIC) simulations with a laser ray-tracing module for a fully self-consistent simulation of the plasma ablation, expansion, and magnetization. The simulations use a planar geometry, effectively suppressing the Biermann magnetic fields, to focus on anisotropy-driven instabilities. The laser intensity is varied between 10 13 and 10 14 W/cm 2 , which is relevant to HED and inertial fusion experiments where collisions must be considered. We find that, above a critical intensity, the plasma rapidly self-magnetizes via an expansion-driven Weibel process, producing a plasma beta of 100 (𝛽 = 8⁢𝜋⁢𝑘 𝐵 ⁢𝑛 𝑒 ⁢𝑇 𝑒 /𝐵 2 ) and Hall parameter 𝜔 ce ⁢𝜏 𝑒 >1 within the first few hundred picoseconds. The magnetic field is sufficiently strong to modify plasma heat transport, and simulations with an artificially suppressed magnetic field show noticeably different temperature profiles.

Lezhnin, K. V. [Princeton Plasma Physics Laborator↗
DOE OSTI2021

Simulating and benchmarking neutron total scattering instrumentation from inception of events to reduced and fitted data

In the design and realization of modern neutron scattering instrumentation, particularly when designing beamline concepts from the ground up, it is desirable to fully benchmark against realistically simulated data. This is especially true for total scattering beamlines, where the future deliverable data is to be analysed in both reciprocal- and real-space representations, and needs must be carefully balanced to ensure sufficient range, resolution and flux of the instrument. An approach to optimize the design of neutron scattering instrumentation via a workflow including ray-tracing simulations, event-based data reduction, heuristic analysis and fitting against realistically simulated spectra is demonstrated here. The case of the DISCOVER beamline concept at the Spallation Neutron Source is used as an example. The results of the calculations are benchmarked through simulation of existing instrumentation and subsequent direct comparison with measured data. Finally, on the basis of the validated models, the ability to explore design characteristics for future beamline concepts or future instrument improvements is demonstrated through the examples of detector tube size and detector layout.

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

A simple optical system delivering a tunable micrometer pink beam that can compensate for heat-induced deformations

The radiation from an undulator reflected from one or more optical elements (usually termed `pink-beam') is used in photon-hungry experiments. The optical elements serve as a high-energy cutoff and for focusing purposes. One of the issues with this configuration is maintaining the focal spot dimension as the energy of the undulator is varied, since this changes the heat load absorbed by the first optical element. Finite-element analyses of the power absorbed by a side water-cooled mirror exposed to the radiation emitted by an undulator at the Advanced Photon Source (APS) and at the APS after the proposed upgrade (APSU) reveals that the mirror deformation is very close to a convex cylinder creating a virtual source closer to the mirror than the undulator source. Here a simple optical system is described based on a Kirkpatrick–Baez pair which keeps the focus size to less than 2 µm (in the APSU case) with a working distance of 350 mm despite the heat-load-induced change in source distance. Detailed ray tracings at several photon energies for both the APS and APSU show that slightly decreasing the angle of incidence on the mirrors corrects the change in the `virtual' position of the source. The system delivers more than 70% of the first undulator harmonic with very low higher-orders contamination for energies between 5 and 10 keV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Mutual optical intensity propagation through non-ideal mirrors

The mutual optical intensity (MOI) model is extended to include the propagation of partially coherent radiation through non-ideal mirrors. The propagation of the MOI from the incident to the exit plane of the mirror is realised by local ray tracing. The effects of figure errors can be expressed as phase shifts obtained by either the phase projection approach or the direct path length method. Using the MOI model, the effects of figure errors are studied for diffraction-limited cases using elliptical cylinder mirrors. Figure errors with low spatial frequencies can vary the intensity distribution, redistribute the local coherence function and distort the wavefront, but have no effect on the global degree of coherence. The MOI model is benchmarked against HYBRID and the multi-electron Synchrotron Radiation Workshop ( SRW ) code. The results show that the MOI model gives accurate results under different coherence conditions of the beam. Other than intensity profiles, the MOI model can also provide the wavefront and the local coherence function at any location along the beamline. The capability of tuning the trade-off between accuracy and efficiency makes the MOI model an ideal tool for beamline design and optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Optimization of a phase-space beam position and size monitor for low-emittance light sources

The recently developed vertical phase-space beam position and size monitor (ps-BPM) system has proven to be able to measure the electron-source position, angle, size and divergence simultaneously in the vertical plane at a single location of a beamline. The optimization of the ps-BPM system is performed by ray-tracing simulation to maximize the instrument sensitivity and resolution. The contribution of each element is studied, including the monochromator, the K -edge filter, the detector and the source-to-detector distance. An optimized system is proposed for diffraction-limited storage rings, such as the Advanced Photon Source Upgrade project. The simulation results show that the ps-BPM system can precisely monitor the source position and angle at high speed. Precise measurements of the source size and divergence will require adequate resolution with relatively longer integration time.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

IRIXS: a resonant inelastic X-ray scattering instrument dedicated to X-rays in the intermediate energy range

A new resonant inelastic X-ray scattering (RIXS) instrument has been constructed at beamline P01 of the PETRA III synchrotron. This instrument has been named IRIXS (intermediate X-ray energy RIXS) and is dedicated to X-rays in the tender-energy regime (2.5–3.5 keV). The range covers the L 2,3 absorption edges of many of the 4delements (Mo, Tc, Ru, Rh, Pd and Ag), offering a unique opportunity to study their low-energy magnetic and charge excitations. The IRIXS instrument is currently operating at the Ru L 3 -edge (2840 eV) but can be extended to the other 4d elements using the existing concept. The incoming photons are monochromated with a four-bounce Si(111) monochromator, while the energy analysis of the outgoing photons is performed by a diced spherical crystal analyzer featuring (10$\bar{2}$ lattice planes of quartz (SiO 2 ). A total resolution of 100 meV (full width at half-maximum) has been achieved at the Ru L 3 -edge, a number that is in excellent agreement with ray-tracing simulations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗