SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Real and complex valued geometrical optics inverse ray-tracing for inline field calculations

In this work, a 3-D ray based model for computing laser fields in dissipative and amplifying media is presented. The eikonal equation is solved using inverse ray-tracing on a dedicated nonstructured 3-D mesh. Inverse ray-tracing opens the possibility of using Complex Geometrical Optics (CGO), for which we propose a propagation formalism in a finite element mesh. Divergent fields at caustics are corrected using an etalon integral method for fold-type caustics. This method is successfully applied in dissipative media by modifying the ray-ordering and root selection rules, thereby allowing one to reconstruct the field in the entire caustic region. In addition, we demonstrate how caustics in the CGO framework can disappear entirely for sufficiently dissipative media, making the complex ray approach valid in the entire medium. CGO is shown to offer a more precise modeling of laser refraction and absorption in a dissipative medium when compared to Geometrical Optics (GO). In the framework of Inertial Confinement Fusion (ICF), this occurs mostly at intermediate temperatures or at high temperatures close to the critical density. Additionally, GO is invalid at low temperatures if an approximated expression of the permittivity is used. The inverse ray-tracing algorithm for GO and CGO is implemented in the IFRIIT code, in the framework of a dielectric permittivity described in 3-D using a piecewise linear approximation in tetrahedrons. Fields computed using GO and CGO are compared to results from the electromagnetic wave solver LPSE. Excellent agreement is obtained in 1-D linear and nonlinear permittivity profiles. Good agreement is also obtained for ICF-like Gaussian density profiles in 2-D. Finally, we demonstrate how the model reproduces Gaussian beam diffraction using CGO. The IFRIIT code will be interfaced inline to 3-D radiative hydrodynamic codes to describe the nonlinear laser plasma interaction in ICF and high-energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Minidisks in Binary Black Hole Accretion

Newtonian simulations have demonstrated that accretion onto binary black holes produces accretion disks around each black hole (“minidisks”), fed by gas streams flowing through the circumbinary cavity from the surrounding circumbinary disk. We study the dynamics and radiation of an individual black hole minidisk using 2D hydrodynamical simulations performed with a new general relativistic version of the moving-mesh code Disco. We introduce a comoving energy variable that enables highly accurate integration of these high Mach number flows. Tidally induced spiral shock waves are excited in the disk and propagate through the innermost stable circular orbit, providing a Reynolds stress that causes efficient accretion by purely hydrodynamic means and producing a radiative signature brighter in hard X-rays than the Novikov–Thorne model. Disk cooling is provided by a local blackbody prescription that allows the disk to evolve self-consistently to a temperature profile where hydrodynamic heating is balanced by radiative cooling. We find that the spiral shock structure is in agreement with the relativistic dispersion relation for tightly wound linear waves. We measure the shock-induced dissipation and find outward angular momentum transport corresponding to an effective alpha parameter of order 0.01. We perform ray-tracing image calculations from the simulations to produce theoretical minidisk spectra and viewing-angle-dependent images for comparison with observations.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

NASA NTRS1977

The linear Fresnel lens - Solar optical analysis of tracking error effects

Real sun-tracking solar concentrators imperfectly follow the solar disk, operationally sustaining both transverse and axial misalignments. This paper describes an analysis of the solar concentration performance of a line-focusing flat-base Fresnel lens in the presence of small transverse tracking errors. Simple optics and ray-tracing techniques are used to evaluate the lens solar transmittance and focal-plane imaging characteristics. Computer-generated example data for an f/1.0 lens indicate that less than a 1% transmittance degradation occurs for transverse errors up to 2.5 deg. In this range, solar-image profiles shift laterally in the focal plane, the peak concentration ratio drops, and profile asymmetry increases with tracking error. With profile shift as the primary factor, the ninety-percent target-intercept width increases rapidly for small misalignments, e.g., almost threefold for a 1-deg error. The analytical model and computational results provide a design base for tracking and absorber systems for the linear-Fresnel-lens solar concentrator.

Cosby, R. M.↗

NASA NTRS2003

Near-source and Remote Observations of Kilometric Continuum Radiation from Multispacecraft Observations

[I] Kilometric continuum (KC) radiation was first identified from Geotail plasma wave observations. Past authors have shown that this emission has a frequency range that overlaps that of the auroral kilometric radiation (AKR) but is characterized by a fine structure of narrow-bandwidth, linear features that have nearly constant or drifting frequency. This fine structure is distinct from that of AKR. KC also apparently has a distinct source region probably associated with the low-latitude inner magnetosphere, consistent with direction-finding and ray-tracing results. We present new high-resolution electric and magnetic field observations of KC obtained by the Polar plasma wave instrument in the near-source region. These observations show intense electrostatic and less intense electromagnetic emissions near the magnetic equator at the plasmapause. Simultaneously, Geotail, located at 20 to 30 RE in radial distance, observes KC in the same frequency range. These data support a possible mode-conversion source mechanism near a region of high-density gradient. High-resolution data obtained from wideband receivers on board both Polar and Cluster show closely spaced bands of emission near the magnetic equator that may be due to many nearby independent sources of EM emission perhaps associated with density fluctuations or cavities in the plasmasphere.

Menietti, J. D.↗

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 OSTI2015

ALGORITHMS AND PROGRAMS FOR STRONG GRAVITATIONAL LENSING IN KERR SPACE-TIME INCLUDING POLARIZATION

Active galactic nuclei (AGNs) and quasars are important astrophysical objects to understand. Recently, microlensing observations have constrained the size of the quasar X-ray emission region to be of the order of 10 gravitational radii of the central supermassive black hole. For distances within a few gravitational radii, light paths are strongly bent by the strong gravity field of the central black hole. If the central black hole has nonzero angular momentum (spin), then a photon’s polarization plane will be rotated by the gravitational Faraday effect. The observed X-ray flux and polarization will then be influenced significantly by the strong gravity field near the source. Consequently, linear gravitational lensing theory is inadequate for such extreme circumstances. We present simple algorithms computing the strong lensing effects of Kerr black holes, including the effects on polarization. Our algorithms are realized in a program “KERTAP” in two versions: MATLAB and Python. The key ingredients of KERTAP are a graphic user interface, a backward ray-tracing algorithm, a polarization propagator dealing with gravitational Faraday rotation, and algorithms computing observables such as flux magnification and polarization angles. Our algorithms can be easily realized in other programming languages such as FORTRAN, C, and C++. The MATLAB version of KERTAP is parallelized using the MATLAB Parallel Computing Toolbox and the Distributed Computing Server. The Python code was sped up using Cython and supports full implementation of MPI using the “mpi4py” package. As an example, we investigate the inclination angle dependence of the observed polarization and the strong lensing magnification of AGN X-ray emission. We conclude that it is possible to perform complex numerical-relativity related computations using interpreted languages such as MATLAB and Python.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

NASA NTRS1996

Modeling of Jovian Hectometric Radiation Source Locations: Ulysses Observations

The Unified Radio and Plasma Wave (URAP) experiment on Ulysses has provided unique high latitude measurements of Jovian hectometric radiation (HOM) during its encounter with Jupiter in February 1992. URAP was the first radio instrument in the Jovian environment with radio direction-finding capability, which was previously used to determine the HOM source locations in the Jovian magnetosphere. These initial source location determinations were based on several assumptions, including the neglect of refractive effects, which may be tested. We have, for the first time, combined the measured incident ray-direction at the spacecraft with a model magnetosphere to directly trace the rays back to the HOM source. We concentrate on the observations of HOM from high northern latitudes when Ulysses was at distances less than 15 R(sub j). The three- dimensional ray-tracing calculations presented here indicate that the HOM sources probably lie on L shells in the range 3 less than or approximately equal to L less than 7 (tilted dipole magnetic field model) consistent with previous determinations that ignored the effects of refraction. The ray-tracing results, however, indicate that wave refraction due to the Io torus and the magnetic field can significantly influence the precise source location. We show that constraints on the locations imposed by the gyroemission mechanism suggest that the lo torus density may have experienced temporal and/or spatial fluctuations during the Ulysses observations of HOM. Finally, in the cold plasma approximation we demonstrate that even if the emission were nearly linearly polarized near the source region, almost circular polarization will be observed at Ulysses, in agreement with observations.

Menietti, J. D.↗

NASA NTRS1990

Optical modeling for dynamics and control analysis

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

Redding, David C.↗

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 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 OSTI2018

An electron Bernstein wave heating scheme for the Proto-MPEX linear device

Linear plasma devices, such as the Prototype Material Plasma Exposure eXperiment (Proto-MPEX), are an economic method to study plasma-material interactions under high heat and particle fluxes. The Proto-MPEX device at Oak Ridge National Laboratory is a high-density helicon plasma generator with additional resonant electron heating to study plasma-material interactions in ITER-like conditions. Here, we present modeling results of an electron Bernstein wave (EBW) heating scheme for this device using a microwave-based system. A 28 GHz injection system is available to provide power levels up to 100 kW injected into low electron temperature and high electron density plasmas. A modified version of the GENRAY ray-tracing code, GENRAY-C, has been used to determine the EBW and electron cyclotron heating wave accessibility for these overdense plasmas. Modeling has shown that greater than 80% of mode conversion from the ordinary mode to the EBW is possible for 28 GHz injection in Proto-MPEX. Calculations show that by utilizing the effects of Doppler-shifted resonance absorption, power deposition near the core is possible. Finally, significant collisional damping, leading to edge absorption, is expected to occur for neutral pressures >1.0 mTorr; however, experiments can be designed to minimize these effects.

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↗