Search NASA⌕ Search

SEARCH · Search NASA

Results for “Linearized 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.

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

New modelling capabilities in IDT

This work concerns the enhanced modelling capabilities of the discrete ordinates transport solver IDT. The novelties introduced allow for modelling unstructured geometries composed by a collection of X/Y segments and circles, and the use of reciprocity and conservation relations reduce the memory imprint as well as the computational cost of the method. IDT decomposes geometries in modular Cartesian patterns, which are the so-called Heterogeneous Cartesian Cells (HCCs), containing a chunk of the original unstructured geometries. Each HCC can be then discretized by superimposing a XY grid to refine locally the HCC. Unlike the most popular MOC, IDT performs the spatial sweeping by directional collision probabilities instead of trajectories. The sources and interface angular fluxes are expanded up to linear order. The accuracy of ray-tracing, the memory imprints together with the novel mesh refinement capabilities have been verified. A first set of preliminary results on PWR lattice problems will be presented in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Modeling and implementation of vertical excursion FFA in the Zgoubi ray-tracing code

Vertical Fixed Field Accelerators (vFFAs) feature complex and highly non-linear magnetic fields, which require simulation codes allowing step-wise particle tracking. Methods to model the 3D magnetic field of scaling vFFAs have been developed in the ray-tracing code Zgoubi. The field modeling and particle tracking methods include the field non-linearities, the fringe fields, and the field superposition of neighboring magnets. The procedure implements the vFFA analytical field expressions, allowing design studies and parameter optimizations using the Zgoubi built-in fit method. The vFFA procedure has been applied to a ten-fold symmetry ring with a triplet focusing structure designed to accelerate protons from 3 MeV to 12 MeV and studied under the ISIS-II proton driver prototype project. Results from particle tracking in externally generated 3D semi-analytical field maps and the developed vFFA analytical model are shown to be in excellent agreement.

43 PARTICLE ACCELERATORS↗
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 OSTI2021

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (ne,lin≈2.0×1019 m−3), simultaneous injection of a 65/35 mix of 2.45/4.6 GHz power achieved an lower hybrid current drive efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker–Planck solver CQL3D. For each discharge, profiles of current and hard x-rays from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Choi, W. (ORCID:0000000268912300)↗
DOE OSTI2021

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (n_e,lin ≈ 2.0 × 10^19m^−3), simultaneous injection of a 65/35 mix of 2.45 GHz/4.6 GHz power achieved an LHCD efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker-Planck solver CQL3D. For each discharge, profiles of current and hard x-ray from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Lower hybrid synergy↗
DOE OSTI2021

Cathode to Target Simulations for Scorpius: I. Simulation Codes and Models

The new Scorpius linear induction electron accelerator is under development for multi-pulse flash radiography of large, explosively-driven hydrodynamic experiments. Beam physics from the cathode to the target was examined with computer simulations, including envelope, ray-trace, and particle-in cell (PIC) codes. Beam instabilities investigated included beam breakup (BBU), image displacement, diocotron, parametric envelope, ion hose, and the resistive wall instability. Beam corkscrew motion and emittance growth from beam mismatch were also studied. The results of these simulations is documented in a series of reports. In this report the computer codes and physical models used for these simulations are described. The conclusion of this study is that Scorpius will produce and accelerate a beam with radiographic quality equivalent to the present accelerators at Los Alamos National Laboratory if the same engineering standards and construction details are upheld.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Techno-economic design of a linear Fresnel reflector for industrial process heat

A techno-economic model of a Concentrating Solar Thermal (CST) system using a Linear Fresnel Reflector (LFR) has been developed. LFRs can deliver process heat suitable for a range of industries, including food and beverage production. This model uses an adaptive algorithm to calculate the optimal secondary reflector shape given the geometry and optical properties of the rest of the system. A ray-tracing program is used to calculate optical efficiency over a wide range of longitudinal and transversal incidence angles and subsequently evaluate the annual efficiency at a given geographical location. A specific LFR design developed by Hyperlight Energy was modelled, and this industrial partner provided a detailed cost breakdown which was used as the basis of an economic model. Combining the technical and economic data facilitates the calculation of the Levelized Cost of Heat (LCOH). The influence of a number of parameters on the annual efficiency and LCOH is explored; notable parameters include the absorber height, the number, width, and spacing of the primary mirrors, the aim point, and the secondary reflector shape and width. By identifying an optimal combination of these parameters, we reduce the LCOH of the industry partner’s system design by 9.2%, from 14.4 $\$/MWh_{th}$ to 13.0 $\$/MWh_{th}$. In comparison, the LCOH of a natural gas boiler delivering the same annual quantity of heat is 29 $\$/MWh_{th}$, which indicates that LFRs can be a competitive heat source for industrial processes.

14 SOLAR ENERGY↗
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 OSTI2021

Initial conditions for simulations of beam physics in linear induction accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the DualAxis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Correct Initial Conditions for Simulations of Beam Physics in Linear Induction Accelerators

Flash radiography of hydrodynamic experiments driven by high explosives is a well-known diagnostic technique in use at many laboratories. At Los Alamos, the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility two linear induction electron accelerators (LIAs) make the bremsstrahlung radiographic source spots for point projection radiographs from orthogonal views. A new LIA, called Scorpius, is presently under development to advance this technology. To better understand electron-beam physics in these LIAs, numerical simulations are frequently performed with the objective of improving the radiography. At Los Alamos we frequently use the TRAK ray-trace and LSP particle-in-cell (PIC) codes to simulate the injector, and the XTR and LAMDA envelope/centroid codes along with LSP to simulate transport of the accelerated beam through the LIAs. The LIA simulations need the injected beam parameters as initial conditions for calculating beam transport and stability. The determination of these initial conditions is the topic of this note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Single-Crystal X-Ray Spectropolarimeter

X-ray spectropolarimetry is a diagnostics technique suitable for measuring intense electron beams and strong magnetic fields in hot, dense plasmas. To increase the accuracy achievable, we investigated the possibility of measuring both orthogonal components of linear polarization simultaneously using a single crystal. We demonstrated experimentally for the first time that crystals with cubic symmetry can be used as polarizing beam splitters. We also found that the polarization-splitting planes in cubic crystals can be accessed in two crystal orientations, a property that can be used to decrease the clutter in polarized spectra. In parallel, we developed a ray-tracing code and are working to add the crystal reflectivity to instrument modeling codes, to help us interpret the spectra and to design spectroscopic and spectropolarimetric instruments. This work is motivated by the need to perform measurements of important quantities in difficult experimental conditions, such as the axial magnetic field at stagnation in magnetized liner inertial fusion loads.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Stochastic Averaging of Radiative Transfer Coefficients for Relativistic Electrons

Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.

Astronomy & Astrophysics↗
DOE OSTI2025

A Scalable Gaussian Process Approach to Shear Mapping with MuyGPs

Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in turn provides us with information about the nature of dark energy. Conventional methods generate shear maps from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser–Squires inversion) for generating these maps, however, are tricky to implement and can introduce bias. Recent alternatives construct a spatial process prior for the lensing potential, which allows for inference of the convergence and shear parameters given lensing shear measurements. Realizing these spatial processes, however, scales cubically in the number of observations—an unacceptable expense as near-term surveys expect billions of correlated measurements. Therefore, we present a linearly scaling shear map construction alternative using a scalable Gaussian process prior called MuyGPs. MuyGPs avoids cubic scaling by conditioning interpolation on only nearest neighbors and fits hyperparameters using batched leave-one-out cross-validation. This work is the first step toward a full, scalable mass mapping method. We work in a simplified regime where we validate our method by interpolating and analyzing maps given noisy point-estimate data from all three shear fields, taken from a suite of N -body ray-tracing simulations. We also show that we can perform these operations at the scale of billions of galaxies on high-performance computing platforms.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗