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

Maps of CMB lensing deflection from N-body simulations in Coupled Dark Energy Cosmologies

We produce lensing potential and deflection-angle maps in order to simulate the weak gravitational lensing of the Cosmic Microwave Background (CMB) via ray-tracing through the COupled Dark Energy Cosmological Simulations (CoDECS), the largest suite of N-body simulations to date for interacting Dark Energy cosmologies. The constructed maps faithfully reflect the N-body cosmic structures on a range of scales going from the arcminute to the degree scale, limited only by the resolution and extension of the simulations. We investigate the variation of the lensing pattern due to the underlying Dark Energy (DE) dynamics, characterised by different background and perturbation behaviours as a consequence of the interaction between the DE field and Cold Dark Matter (CDM). In particular, we study in detail the results from three cosmological models differing in the background and perturbations evolution at the epoch in which the lensing cross section is most effective, corresponding to a redshift of ∼ 1, with the purpose to isolate their imprints in the lensing observables, regardless of the compatibility of these models with present constraints. The scenarios investigated here include a reference ΛCDM cosmology, a standard coupled DE (cDE) scenario, and a ''bouncing'' cDE scenario. For the standard cDE scenario, we find that typical differences in the lensing potential result from two effects: the enhanced growth of linear CDM density fluctuations with respect to the ΛCDM case, and the modified nonlinear dynamics of collapsed structures induced by the DE-CDM interaction. As a consequence, CMB lensing highlights the DE impact in the cosmological expansion, even in the degenerate case where the amplitude of the linear matter density perturbations, parametrised through σ{sub 8}, is the same in both the standard cDE and ΛCDM cosmologies. For the ''bouncing'' scenario, we find that the two opposite behaviours of the lens density contrast and of the matter abundance lead to a counter-intuitive effect, making the power of the lensing signal in this model lower by 10% than in the ΛCDM scenario. Moreover, we compare the behaviour of CDM and baryons in CoDECS separately, in order to isolate effects coming from the coupling with the DE component. We find that, in the bouncing scenario, baryons show an opposite trend with respect to CDM, due to the coupling of the latter with the DE component. These results confirm the relevance of CMB lensing as a probe for DE at the early stages of cosmic acceleration, and demonstrate the reliability of N-body based large scale CMB lensing simulations in the context of DE studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

N-body lensed CMB maps: lensing extraction and characterization

We reconstruct shear maps and angular power spectra from simulated weakly lensed total intensity (TT) and polarised (EB) maps of the Cosmic Microwave Background (CMB) anisotropies, obtained using Born approximated ray-tracing through the N-body simulated Cold Dark Matter (CDM) structures in the Millennium Simulations (MS). We compare the recovered signal with the ΛCDM prediction, on the whole interval of angular scales which is allowed by the finite box size, extending from the degree scale to the arcminute, by applying a quadratic estimator in the flat sky limit; we consider PRISM-like instrumental specification for future generation CMB satellites, corresponding to arcminute angular resolution of 3.2' and sensitivity of 2.43 μK-arcmin. The noise contribution in the simulations closely follows the estimator prediction, becoming dominated by limits in the angular resolution for the EB signal, at ℓ ≅ 1500. The recovered signal shows no visible departure from predictions of the weak lensing power within uncertainties, when considering TT and EB data singularly. In particular, the reconstruction precision reaches the level of a few percent in bins with Δℓ ≅ 100 in the angular multiple interval 1000∼<ℓ∼<2000 for T, and about 10% for EB. Within the adopted specifications, polarisation data do represent a significant contribution to the lensing shear, which appear to faithfully trace the underlying N-body structure down to the smallest angular scales achievable with the present setup, validating at the same time the latter with respect to semi-analytical predictions from ΛCDM cosmology at the level of CMB lensing statistics. This work demonstrates the feasibility of CMB lensing studies based on large scale simulations of cosmological structure formation in the context of the current and future high resolution and sensitivity CMB experiment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Relativistic weak lensing from a fully non-linear cosmological density field

In this paper we examine cosmological weak lensing on non-linear scales and show that there are Newtonian and relativistic contributions and that the latter can also be extracted from standard Newtonian simulations. We use the post-Friedmann formalism, a post-Newtonian type framework for cosmology, to derive the full weak-lensing deflection angle valid on non-linear scales for any metric theory of gravity. We show that the only contributing term that is quadratic in the first order deflection is the expected Born correction and lens-lens coupling term. We use this deflection angle to analyse the vector and tensor contributions to the E- and B- mode cosmic shear power spectra. In our approach, once the gravitational theory has been specified, the metric components are related to the matter content in a well-defined manner. Specifying General Relativity, we write down a complete set of equations for a GR+ΛCDM universe for computing all of the possible lensing terms from Newtonian N-body simulations. We illustrate this with the vector potential and show that, in a GR+ΛCDM universe, its contribution to the E-mode is negligible with respect to that of the conventional Newtonian scalar potential, even on non-linear scales. Thus, under the standard assumption that Newtonian N-body simulations give a good approximation of the matter dynamics, we show that the standard ray tracing approach gives a good description for a ΛCDM cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The theory of stochastic cosmological lensing

On the scale of the light beams subtended by small sources, e.g. supernovae, matter cannot be accurately described as a fluid, which questions the applicability of standard cosmic lensing to those cases. In this article, we propose a new formalism to deal with small-scale lensing as a diffusion process: the Sachs and Jacobi equations governing the propagation of narrow light beams are treated as Langevin equations. We derive the associated Fokker-Planck-Kolmogorov equations, and use them to deduce general analytical results on the mean and dispersion of the angular distance. This formalism is applied to random Einstein-Straus Swiss-cheese models, allowing us to: (1) show an explicit example of the involved calculations; (2) check the validity of the method against both ray-tracing simulations and direct numerical integration of the Langevin equation. As a byproduct, we obtain a post-Kantowski-Dyer-Roeder approximation, accounting for the effect of tidal distortions on the angular distance, in excellent agreement with numerical results. Besides, the dispersion of the angular distance is correctly reproduced in some regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

DEMNUni: ISW, Rees-Sciama, and weak-lensing in the presence of massive neutrinos

We present, for the first time in the literature, a full reconstruction of the total (linear and non-linear) ISW/Rees-Sciama effect in the presence of massive neutrinos, together with its cross-correlations with CMB-lensing and weak-lensing signals. The present analyses make use of all-sky maps extracted via ray-tracing across the gravitational potential distribution provided by the ''Dark Energy and Massive Neutrino Universe'' (DEMNUni) project, a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as separate collisionless particles. We correctly recover, at 1–2% accuracy, the linear predictions from CAMB. Concerning the CMB-lensing and weak-lensing signals, we also recover, with similar accuracy, the signal predicted by Boltzmann codes, once non-linear neutrino corrections to HALOFIT are accounted for. Interestingly, in the ISW/Rees-Sciama signal, and its cross correlation with lensing, we find an excess of power with respect to the massless case, due to free streaming neutrinos, roughly at the transition scale between the linear and non-linear regimes. The excess is ∼ 5 – 10% at l ∼ 100 for the ISW/Rees-Sciama auto power spectrum, depending on the total neutrino mass M {sub ν}, and becomes a factor of ∼ 4 for M {sub ν} = 0.3 eV, at l ∼ 600, for the ISW/Rees-Sciama cross power with CMB-lensing. This effect should be taken into account for the correct estimation of the CMB temperature bispectrum in the presence of massive neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Light propagation through black-hole lattices

The apparent properties of distant objects encode information about the way the light they emit propagates to an observer, and therefore about the curvature of the underlying spacetime. Measuring the relationship between the redshift z and the luminosity distance D {sub L} of a standard candle, for example, yields information on the Universe's matter content. In practice, however, in order to decode this information the observer needs to make an assumption about the functional form of the D {sub L}( z ) relation; in other words, a cosmological model needs to be assumed. In this work, we use numerical-relativity simulations, equipped with a new ray-tracing module, to numerically obtain this relation for a few black-hole-lattice cosmologies and compare it to the well-known Friedmann-Lema(ȋtre-Robertson-Walker case, as well as to other relevant cosmologies and to the Empty-Beam Approximation. We find that the latter provides the best estimate of the luminosity distance and formulate a simple argument to account for this agreement. We also find that a Friedmann-Lema(ȋtre-Robertson-Walker model can reproduce this observable exactly, as long as a time-dependent cosmological constant is included in the fit. Finally, the dependence of these results on the lattice mass-to-spacing ratio μ is discussed: we discover that, unlike the expansion rate, the D {sub L}( z ) relation in a black-hole lattice does not tend to that measured in the corresponding continuum spacetime as 0μ → .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Plasma heating by electron cyclotron wave and the temperature effects on lower hybrid current drive on EAST

This paper presents the progress in the long—pulse operation of the electron cyclotron (EC) system and the achievements in high—electron temperature plasmas by the combined EC and lower hybrid (LH) waves heating since the EC system was built in 2015. An electron temperature of up to 12 keV with a duration over 100 s was realized by the simultaneous heating of EC and LH waves at the line-averaged density ${\bar n_{\text{e}}}$ ~1.8 × 10 19 m –3 . The plasma heating effect strongly depends on the location of EC power deposition. H-mode plasmas with solely EC wave auxiliary heating have been obtained on EAST for the first time. These H-mode discharges show an enhanced confinement factor H 98(y,2) around 1.0, which is higher than the previous H-mode using LH power alone (Xu et al 2011 Nucl. Fusion 51 072001). The total heating power is very close to the threshold value for L–H transition according to the international tokamak scaling. In addition to the temperature effects inside the separatrix, higher electron temperature produced by EC wave is found to reduce the LH power loss in the scrape-off layer due to collisional absorption, which is beneficial to further increase the LH current drive efficiency. Ray-tracing/Fokker–Planck modeling results indicate that higher electron temperature can shorten the LH wave propagation on EAST in a multiple-pass regime, thus decreasing the collisional dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First high-power helicon results from DIII-D

Abstract More than 0.6 MW of rf power at 476 MHz has been coupled to DIII-D plasmas by launching helicon (whistler) waves with a traveling-wave antenna (comb-line) in the fast-wave polarization (Van Compernolle et al 2021 Nucl. Fusion 61 116034) which resulted in the observation of electron heating of the core plasma with single-pass absorption based on ray-tracing in L-mode discharges. The coupling performance of the 1.5 m wide 30-element comb-line traveling-wave antenna has been consistent with expectations based on the 2015–2016 experiments on DIII-D with a low-power 12-element prototype (Pinsker et al 2018 Nucl. Fusion 58 106007). The conditioning process that was necessary to carry out high-power experiments is discussed; rf-specific impurities have not been observed. Parametric decay instabilities have been observed and are being investigated as a potential edge absorption mechanism (Porkolab et al 2023 AIP Conf. Proc. 2984 070004).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design workflow of a symmetric traveling wave antenna for fast ion production on DD tokamaks

Initial computational plasma physics scoping and a finite element method antenna modeling design workflow for a symmetric center-fed high-field side high harmonic fast wave traveling wave array (TWA) antenna are reported here. The TWA is designed to generate a test population of fast deuterium ions in an existing D–D tokamak by heating neutral beam deuterium ions, accelerating them from 80 keV to several hundred keV. The resulting fast particles are tailored to mimic key reactor energetic particle parameters with regards to exciting Alfven eigenmode instabilities, allowing for a D–D tokamak like DIII-D or ASDEX-U to replicate reactor-relevant conditions experimentally. Initial scenario scoping for high single-pass absorption as well as good preferential fast ion damping relative to electron damping was completed using the ray-tracing/Fokker–Planck codes GENRAY and CQL3D. Python RF network analysis packages were used to create a custom TWA optimization tool to inform a COMSOL flat antenna design, and Petra-M was used to study cold plasma effects. The TWA produced by this workflow has several novel features when compared to previous TWA studies, including symmetric center feeding, and passive end straps for image current cancellation for reduced impurity production. We show here that the antenna design workflow can readily produce TWA antennas optimized for reflection coefficient, image current cancellation, and launched power spectrum shape; and that a population of fast ions can be generated in the correct region of parameter space, warranting future more detailed studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Numerical investigation of liquid wall ablation in inertial fusion energy chambers

This paper presents a novel approach for modeling liquid wall ablation in liquid wall-protected inertial fusion energy (IFE) chambers. These systems are promising candidates for the implementation of fusion technology, yet significant gaps remain in understanding the underlying physical processes and their implications for design. Following target ignition, a portion of the fusion energy is released as x-rays, which deposit their energy into an array of liquid jets, leading to partial vaporization. Accurately modeling this heat deposition and vaporization process remains challenging due to the complex geometries typical of (pre-conceptual) IFE chamber designs. Furthermore, the subsequent expansion of vaporized material into the chamber’s vacuum environment poses difficulties for conventional CFD methods based on continuum assumptions, which can lead to significant inaccuracies. To address some aspects of these challenges, this work introduces a ray-tracing-based methodology to map the spatial distribution of ablated material in liquid wall-protected systems. In addition, a vacuum-tracking scheme is developed to extend the applicability of an OpenFOAM-based solver to gas dynamics in rarefied environments. The proposed approach has been verified through numerical benchmarks and applied to a practical case involving the HYLIFE-II (High Yield Lithium Injection Fusion Energy) chamber. The methodology advances the modeling capabilities for liquid wall-protected IFE systems and provides valuable tools to support their design and optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Parametrizations of coupled betatron motion for strongly coupled lattices

The coupling of transverse motion is a natural occurrence in particle accelerators, either in the form of a residual coupling arising from imperfections or originating by design from strong systematic coupling fields. While the first can be treated perturbatively, the latter requires a robust approach adapted to strongly coupled optics, and a parametrization of the linear optics must be performed to explore beam dynamics in such peculiar lattices. This work highlights the key physical interpretations of the main parametrization formalisms to describe linear coupled optics, along with explicit links and comparisons of these parametrizations. Concepts rarely illustrated in other works, such as forced mode flips and local coupling, are explored in detail, clarifying some anomalies that can arise in lattice functions. The analytical methods have been implemented in a reference Python package and connected with ray-tracing and integration codes to explore examples of strongly coupled lattices, which are discussed in detail to highlight the key physical interpretations of the parametrizations and characteristics of the lattices.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

The pathfinder X-Ray axion telescope for the International Axion Observatory (IAXO) and axion searches at the CERN Axion Solar Telescope (CAST)

The axion, a hypothetical pseudo-scalar particle proposed by Peccei and Quinn to resolve the strong CP problem in quantum chromodynamics, remains one of the most compelling dark matter candidates. Axions and Axion-Like Particles (ALPs) are characterized by a broad and largely unconstrained parameter space in mass and coupling strength, motivating extensive experimental searches. The International Axion Observatory (IAXO) is a next-generation axion helioscope designed to achieve over an order of magnitude improvement in sensitivity to the axion-photon coupling constant gaγ relative to previous experiments such as the CERN Axion Solar Telescope (CAST). IAXO will employ large-scale superconducting magnets, precision x-ray optics, and ultra-low-background detectors to search for solar axions produced primarily via the Primakoff effect. This work presents the design and performance of the IAXO pathfinder x-ray optic, a technological demonstrator that validated the optical concepts for IAXO through axion searches at CAST. Comprehensive ray-tracing simulations and experimental measurements were conducted to assess the optic's effective area, detection efficiency, and background suppression capabilities. Deployed at CAST in conjunction with Micromegas and GridPix detectors, the Pathfinder demonstrated record-setting sensitivity to both gaγ and the axion-electron coupling gae, establishing the most stringent laboratory constraints on solar axions to date. These results confirm the viability of the IAXO optical design and detector concept, providing a critical technological benchmark that informs the ongoing development of the BabyIAXO and full-scale IAXO telescopes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

Enhanced absorption of thin-film photovoltaic cells using an optical cavity

We show here via numerical simulations that the absorption and solar energy conversion efficiency of a thin-film photovoltaic (PV) cell can be significantly enhanced by embedding it into an optical cavity. A reflective hemi-ellipsoid with an aperture for sunlight placed over a tilted PV cell reflects unabsorbed photons back to the cell, allowing for multiple opportunities for absorption. Ray tracing simulations predict that with the proposed cavity a textured thin-film silicon cell can exceed the Yablonovitch (Lambertian) limit for absorption across a broad wavelength range, while the performance of the cavity-embedded planar PV cell approaches that of the cell with the surface texturing.

14 SOLAR ENERGY↗
DOE OSTI2019

Recent developments of MCViNE and its applications at SNS

Abstract MCViNE is an open source, object-oriented Monte Carlo neutron ray-tracing simulation software package. Its design allows for flexible, hierarchical representations of sophisticated instrument components such as detector systems, and samples with a variety of shapes and scattering kernels. Recently this flexible design has enabled several applications of MCViNE simulations at the Spallation Neutron Source (SNS) at Oak Ridge National Lab, including assisting design of neutron instruments at the second target station and design of novel sample environments, as well as studying effects of instrument resolution and multiple scattering. Here we provide an overview of the recent developments and new features of MCViNE since its initial introduction (Jiao et al 2016 Nucl. Instrum. Methods Phys. Res., Sect. A 810 , 86–99), and some example applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗
DOE OSTI2018

Characteristics of infrasound signals from North Korean underground nuclear explosions on 2016 January 6 and September 9

This study investigates infrasound signals from North Korean underground nuclear explosions on 2016 January 6 at 01:30:01 UTC (UNEJ16) and on 2016 September 9 at 00:30:01 UTC (UNES16). Infrasound observations from eight, seismo-acoustic arrays at distances of 300–600 km from the explosions and aligned approximately east–west across South Korea are used in the analysis. Progressive multichannel correlation was used followed by analyst review to identify the infrasound arrivals with estimates of azimuth, phase velocity, celerity, F-statistic, signal-to-noise ratio, correlation and signal duration. These observations are compared with model predictions of travel path, ray turning height and relative amplitudes, using Ground-to-Space atmospheric specifications at the time of each explosion. Stratospheric conditions at the time of the UNEJ16 are favourable to the Korean arrays when all arrays detected infrasound signals, while the few detections for UNES16 are consistent with unfavourable stratospheric winds and reduced amplitude thermospheric arrivals. Infrasound locations of the explosions were estimated using the Bayesian Infrasonic Source Location method (BISL), with wind-corrected backazimuth estimates based on the ray tracing calculations. Adding wind corrections to BISL resulted in infrasound locations up to 50 per cent closer to the seismic epicentre over uncorrected locations. The UNES16 infrasound location estimate using atmospheric corrections had larger differences (~67 km difference from seismic location) than the UNEJ16 location (~11 km difference from seismic location) as a result of the reduced number of observations and poor azimuth coverage for UNES16. Based on stratospheric phases recorded at six arrays and the empirical yield-scaling relations, wind corrected amplitudes from UNEJ16 lead to equivalent infrasound source energy estimates between 0.9 and 16.1 tons of TNT detonated in the atmosphere with an average of 6.4 tons and a standard deviation of 4.6 tons.

58 GEOSCIENCES↗
DOE OSTI2021

Time delay lens modelling challenge

ABSTRACT In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant H0. However, published state-of-the-art analyses require of order 1 yr of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge, we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated data sets. The results in Rungs 1 and 2 show that methods that use only the point source positions tend to have lower precision ($10\!-\!20{{\ \rm per\ cent}}$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic data sets can recover H0 within the target accuracy (|A| < 2 per cent) and precision (<6 per cent per system), even in the presence of a poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

Ding, X.↗
DOE OSTI2023

Pulse profile modelling of thermonuclear burst oscillations – II. Handling variability

ABSTRACT Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii, and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in significant biases in the inferred mass and radius, particularly for the high count rates that are nominally required to obtain meaningful constraints. In this follow-on study, we show that the bias can be mitigated by slicing the bursts into shorter segments where variability can be neglected, and jointly fitting the segments. Using this approach, the systematic uncertainties on the mass and radius are brought within the range of the statistical uncertainty. With about 106 source counts, this yields uncertainties of approximately 10 per cent for both the mass and radius. However, this modelling strategy requires substantial computational resources. We also confirm that the posterior distributions of the mass and radius obtained from multiple bursts of the same source can be merged to produce outcomes comparable to that of a single burst with an equivalent total number of counts.

Astronomy & Astrophysics↗