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 829 records · Page 46

Characterization and simulation of soft gamma-ray mirrors for their use with spent fuel rods at reprocessing facilities

The use of a grazing incidence optic to selectively reflect K-shell fluorescence emission and isotope-specific lines from special nuclear materials is a highly desirable nondestructive analysis method for use in reprocessing fuel environments. Preliminary measurements have been performed, and a simulation suite has been developed to give insight into the design of the x ray optics system as a function of the source emission, multilayer coating characteristics, and general experimental configurations. As a result, the experimental results are compared to the predictions from our simulation toolkit to illustrate the ray-tracing capability and explore the effect of modified optics in future measurement campaigns.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design of a broadband, high-power, off-axis spiral polarization reflector

Spirally polarized beams are known to produce a strong longitudinal field under high numerical aperture focusing. Here, we report on the design of a reflective optical element with an off-axis geometry that is capable of transforming a uniformly polarized beam into a predominantly spirally polarized beam. By utilizing the retardance present in Fresnel reflection, our device can operate over a broad range of wavelengths. The output polarization and intensity structure of a beam passing through the reflector are calculated using polarization and scalar ray tracing. Vector diffraction is used to calculate the polarization and intensity structure under $f/0.656$ focusing. The off-axis reflective design is well-suited for high-power laser applications and could enable current and future laser-based particle accelerators to achieve strong longitudinal fields at focus.

Black, A. Nicholas [Univ. of Rochester, NY (United↗

Altering apparent optical properties with an array of semitransparent mesoscale structures

The ability to control and optimize interactions between light and matter has much utility in engineering design. A well-researched way to achieve optical property modulation is via the use of optical metamaterials, which feature sub-wavelength scale surface structures. In this work, an alternative approach for modulating optical properties is presented using a composite surface modified with a periodic array of semitransparent hemispherical shell mesoscale structures which are larger than the incident light wavelength. A ray-tracing simulation approach is used to predict the optical behavior for an arrayed surface. At oblique angles of incidence, significant increases and decreases in apparent absorptance are achieved via the use of optically thick and thin shells, respectively. Additionally, a potential application to solar cells is described with optimal spectral behavior achieved via the use of semitransparent external structures.

Kucuktas, Onur A.↗

Modal phase measuring deflectometry

Here in this work, a model based method is applied to phase measuring deflectometry, which is named as modal phase measuring deflectometry. The height and slopes of the surface under test are represented by mathematical models and updated by optimizing the model coefficients to minimize the discrepancy between the reprojection in ray tracing and the actual measurement. The pose of the screen relative to the camera is pre-calibrated and further optimized together with the shape coefficients of the surface under test. Simulations and experiments are conducted to demonstrate the feasibility of the proposed approach.

42 ENGINEERING↗

Reflective all-sky thermal infrared cloud imager

A reflective all-sky imaging system has been built using a long-wave infrared microbolometer camera and a reflective metal sphere. This compact system was developed for measuring spatial and temporal patterns of clouds and their optical depth in support of applications including Earth-space optical communications. The camera is mounted to the side of the reflective sphere to leave the zenith sky unobstructed. The resulting geometric distortion is removed through an angular map derived from a combination of checkerboard-target imaging, geometric ray tracing, and sun-location-based alignment. A tape of high-emissivity material on the side of the reflector acts as a reference that is used to estimate and remove thermal emission from the metal sphere. In conclusion, once a bias that is under continuing study was removed, sky radiance measurements from the all-sky imager in the 8-14 μm wavelength range agreed to within 0.91 W/(m 2 sr) of measurements from a previously calibrated, lens-based infrared cloud imager over its 110° field of view.

47 OTHER INSTRUMENTATION↗

Effects of a nuclear-disturbed environment on electromagnetic wave propagation through the atmosphere

This paper investigates the effects of a nuclear-disturbed environment on the transmission of electromagnetic (EM) waves through the atmosphere. An atmospheric nuclear detonation can produce heightened free electron densities in the surrounding atmosphere that can disrupt EM waves that propagate through the disturbed region. Radiation transport models simulated the ionization and free electron densities created in the atmosphere from a 1 MT detonation at heights of burst of 5 km, 25 km, and 75 km. Recombination rates for the free electrons in the atmosphere were applied, from previous work in the literature, to determine the nuclear-induced electron densities as a function of time and space after the detonation. A ray-tracing algorithm was applied to determine the refraction and reflection of waves propagating in the different nuclear-disturbed environments. The simulation results show that the free electron plasma created from an atmospheric nuclear detonation depend on the height of burst of the weapon, the weapon yield, and the time after detonation. Detonations at higher altitudes produce higher free electron densities for greater durations and over larger ranges. The larger the free electron densities, the greater the impact on EM wavelengths in regards to refraction, reflection, and absorption in the atmosphere. An analysis of modern infrastructure and the effects of nuclear-disturbed atmospheres on different signal wavelengths and systems is discussed.

42 ENGINEERING↗

Manufacturability-based optical design optimization for advanced Kirkpatrick–Baez X-ray focusing mirrors

The advanced Kirkpatrick–Baez (AKB) mirror setup is an effective and compelling solution to provide stable X-ray nano-focusing for synchrotron radiation or free-electron laser beamlines. We propose an AKB mirror design optimization approach to mitigate the difficulties associated with mirror fabrication by minimizing the total slope ranges of the four curved mirrors while achieving the expected focusing performance. In the optimization, we have considered geometry constraints to ensure the beam acceptance with the required clear aperture, the diffraction-limited focal size with the adequate numerical aperture, and the desired mirror gaps for adjustment and the necessary working distance for the sample stage. Additionally, practical constraints linked to mirror metrology and fabrication, such as mirror length limits and curvature uncertainty in measurement, are taken into account. Furthermore, progressive objective optimization eliminates the need for any initial guess, fully automating the AKB optimization process. This approach facilitates the development of an elegant Wolter-I or Wolter-III type AKB design solution that satisfies these multiple constraints. In cases where constraints cannot be simultaneously satisfied, the optimization results provide valuable insights into areas where trade-offs need to be considered. Simulations with ray tracing and wavefront propagation validate the optimized AKB design showing high tolerance to the beam incident angle.

36 MATERIALS SCIENCE↗

Offner radial group delay compensator for ultra-broadband laser beam transport

Conventional lens-based image relays introduce radial group delay that significantly reduces focal-spot intensity of ultra-broadband short pulses. A direct compensation scheme that can be used with conventional singlet image relays is proposed. The compensation scheme is based on using two negative lenses embedded within an Offner imaging system. The setup allows a single-pass configuration and eliminates the need for polarization switching optics. The performance of this setup is calculated using ray tracing. The compensation improves the spatiotemporal Strehl ratio from 0.02 to 0.97.

Bahk, Seung-Whan↗

CHESS 2025: Orthorectified airborne RGB imagery from NEON AOP surveys

This dataset provides Level 1 (L1) and Level 3 (L3) orthorectified Red-Green-Blue (RGB) imagery collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). This high-resolution imagery is a photographic record of red, green, and blue visible light from sunlight reflected off of the Earth’s surface. The data comprise full-color images of the ground surface and are primarily intended to provide context to imaging spectroscopy and light detection and ranging (LiDAR) data. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. RGB images were acquired using the PhaseOne IXM-RS150F high-resolution digital camera onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). The package data include both an L1 product comprising one camera frame per file and an L3 mosaic aligned to the Universal Transverse Mercator (UTM) Zone 13N grid and the World Geodetic System (WGS) 84 projection. Both products are provided in geotif (.tif) format at 0.1 m ground resolution. The bulk of the imagery was collected during the main CHESS field campaign from June 13 to July 15, 2025. Additional images of a portion of the Upper Taylor (UPTA) domain were collected on September 18, 2025, to fill gaps in imagery identified after the main campaign was complete. RGB camera imagery is not radiometrically calibrated, and therefore pixel values should not be exploited for scientific analysis. Pixel values have undergone a manual adjustment to enhance feature identification. The imagery is rigorously geolocated which does allow for reliable geometric information to be retrieved. To generate the orthorectified imagery, the NEON AOP camera captured visible spectrum in red, green, and blue bands. The raw images were then processed using NEON’s camera orthorectification workflow. A boresight calibration flight was made to build a complete camera, distortion, and alignment model. Color balance/white balance and exposure correction were applied to the raw RGB images. The corrected images were orthorectified by ray-tracing image pixels to a lidar-derived digital surface model (DSM) mesh using the refined camera model, outputting orthorectified raster pixels on a regular grid. Flightline-level data were mosaicked by selecting per-pixel contributions from overlapping orthorectified images using line-of-sight (LOS) zenith angle minimization to reduce edge distortions. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Investigation of structural heterogeneity at the SPE site using combined P –wave travel times and Rg phase velocities

Here, we present analyses of the 2D seismic structure beneath Source Physics Experiments (SPE) geophone lines that extended radially at 100 m spacing from 100 to 2000 m from the source borehole. With seismic sources at only one end of the geophone lines, standard refraction profiling methods cannot resolve seismic velocity structures unambiguously. In previous work, we demonstrated overall agreement between body-wave refraction modeling and Rg dispersion curves for the least complex of the five lines. A more detailed inspection supports a 2D reinterpretation of the structure. We obtained Rg phase velocity measurements in both the time and frequency domains, then used iterative adjustment of the initial 1D body-wave model to predict Rg dispersion curves to fit the observed values. Our method applied to the most topographically severe of the geophone lines is supplemented with a 2D ray-tracing approach, whose application to P-wave arrivals supports the Rg analysis. In addition, midline sources will allow us to refine our characterization in future work.

58 GEOSCIENCES↗

SKS Polarization Anomalies Due to the Coriolis Force

The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

58 GEOSCIENCES↗

Simulated NMIS Imaging Data for an Unknown Object

This report presents simulated Nuclear Materials Identification System (NMIS) imaging data for an unclassified object, whose characteristics are initially unknown to the reader. This data will be used to test various analysis capabilities and was created with a simple deterministic ray-tracing algorithm. NMIS is a time-dependent coincidence counting system that is used to characterize both fissile and non-fissile materials undergoing nondestructive assay. NMIS characterizes materials by interrogating them with neutrons, either from an associated-particle deuterium-tritium (DT) neutron generator, which produces a time and directionally tagged monoenergetic beam of 14.1 MeV neutrons, or a time-tagged spontaneous fission source in an ionization chamber.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Zgoubi-ing AGS : spin motion with snakes and jump-quads, Gγ = 43.5 → Gγ = 46.5 and beyond

This Note reports on the first, and successful, simulations of particle and spin dynamics in the AGS in presence of the two helical snakes and of the tune-jump quadrupoles, using the ray-tracing code Zgoubi. It includes DA tracking in the absence or in the presence of the two helical snakes, simulation of particle and spin motion in the snakes using their magnetic field maps, spin flipping at integer resonances in the 36+Qy depolarizing resonance region, with and without tune-jump quadrupole gymnastics. It also includes details on the setting-up of Zgoubi input data files and on the various numerical methods of concern in and available from Zgoubi.

43 PARTICLE ACCELERATORS↗

SciDAC Cooperative Agreement: Center for Wave Interactions with Magnetohydrodynamics (Final Technical Report)

Final technical report for research performed by Dr. Thomas G. Jenkins in collaboration with Professor Dalton D. Schnack on SciDAC Cooperative Agreement: Center for Wave Interactions with Magnetohydrodyanics, DE-FC02-06ER54899, for the period of 8/15/06 - 8/14/11. This report centers on the Slow MHD physics campaign work performed by Dr. Jenkins while at UW-Madison and then at Tech-X Corporation. To make progress on the problem of RF induced currents affect magnetic island evolution in toroidal plasmas, a set of research approaches are outlined. Three approaches can be addressed in parallel. These are: (1) Analytically prescribed additional term in Ohm's law to model the effect of localized ECCD current drive; (2) Introduce an additional evolution equation for the Ohm's law source term. Establish a RF source 'box' where information from the RF code couples to the fluid evolution; and (3) Carry out a more rigorous analytic calculation treating the additional RF terms in a closure problem. These approaches rely on the necessity of reinvigorating the computation modeling efforts of resistive and neoclassical tearing modes with present day versions of the numerical tools. For the RF community, the relevant action item is - RF ray tracing codes need to be modified so that general three-dimensional spatial information can be obtained. Further, interface efforts between the two codes require work as well as an assessment as to the numerical stability properties of the procedures to be used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development Status of PROTEUS-MOC

The PROTEUS-MOC solver is intended to resolve the outstanding computational challenges that the PROTEUS-SN is limited by and bridge the gap between PROTEUS-SN and existing homogenized structured grid solvers. The primary push behind all of this work is to carry out coupled thermal-mechanical, fluid dynamics, and neutronics physics problems. The PROTEUS-MOC development has taken advantage of the combined experience with PROTEUS-SN and DeCART with the goal of producing a tool which combines the best features of both tools. Taking the DeCART tool as the model, it merges a subgroup cross section interpolation technique with a 2D method of characteristics approach and a 1D SN formulation. DeCART has been specifically demonstrated to be accurate on PWR and VHTR problems, but it has displayed several issues with regard to accuracy on thin axial mesh sizes which are attributable to the 2D1D scheme it is built upon. There are also notable problems with performance and scalability all of which we hope to resolve in PROTEUS-MOC. PROTEUS-MOC is at its core, based upon the success achieved with MOCFE which demonstrated excellent parallelization with respect to angle and good parallelization with respect to space. More recent efforts have indicated parallelization with respect to energy is also feasible. The work completed yielded a novel methodology where the axial basis functions are treated separately from the radial ones, but the geometry is extruded. This eliminates the need to carry out and store 3D ray tracing data and simultaneously allows the coarse axial node sizes typical with DeCART. To rapidly build a prototype, a standard power iteration algorithm was constructed using an unpreconditioned GMRES solution algorithm for solving the flux vector in each fission source iteration. The new PROTEUS-MOC solver was deployed on two test problems: Takeda and a heterogeneous cylindrical heavy water system. The Takeda benchmark problem is a standard homogeneous test problem which was previously found to be too difficult for the 3D MOCFE code to solve. Using the PROTEUS-MOC we were able to show space-angle convergence although full convergence requires further improvements in performance. The second benchmark problem was constructed more for future work than this year, but the results generated are promising indicating that it will be very useful in highlighting the weaknesses and strengths of the PROTEUS-MOC tool. While we are far from a production tool, the results included in this report give us confidence that we are on the right path to constructing an interceding code between the plethora of structured grid homogenized assembly solvers and the supercomputer consuming PROTEUS-SN code. With time, the creation of a suitable preconditioner along with typical optimization should reduce the performance shortcomings that are observable today.

Marin-Lafleche, A.↗

ECRH/EBWH SYSTEM FOR NSTX-

The National Spherical Torus Experiment Upgrade (NSTX-U) will operate at an axial toroidal field of up to 1 T, about twice the field available on NSTX. A 28 GHz electron cylotron resonance heating (ECRH) system is currently being planned for NSTX-U. A 1 MW 28 GHz gyrotron will be employed. Intially the system will use short, 10-50 ms, 1 MW pulses for ECRH-assisted discharge start-up. Later the pulse length will be extended to 1-5 s to study electron Bernstein wave heating (EBWH) during the plasma current flat top. A mirror launcher will be used to couple microwave power to the plasma via O-mode to the slow X-mode to EBW (O-X-B) double mode conversion. This paper presents a pre-conceptual design for the ECRH/EBWH system proposed for NSTX-U and includes ray tracing and Fokker-Planck modeling results for 28 GHz ECRH during plasma start-up and EBW heating and current drive during the plasma current flattop of a NSTX-U advanced H-mode plasma scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Report for "Tech-X Corporation work for the SciDAC Center for Simulation of RF Wave Interactions with Magnetohydrodynamics (SWIM)"

Work carried out by Tech-X Corporation for the DoE SciDAC Center for Simulation of RF Wave Interactions with Magnetohydrodynamics (SWIM; U.S. DoE Office of Science Award Number DE-FC02-06ER54899) is summarized and is shown to fulfil the project objectives. The Tech-X portion of the SWIM work focused on the development of analytic and computational approaches to study neoclassical tearing modes and their interaction with injected electron cyclotron current drive. Using formalism developed by Hegna, Callen, and Ramos [Phys. Plasmas 16, 112501 (2009); Phys. Plasmas 17, 082502 (2010); Phys. Plasmas 18, 102506 (2011)], analytic approximations for the RF interaction were derived and the numerical methods needed to implement these interactions in the NIMROD extended MHD code were developed. Using the SWIM IPS framework, NIMROD has successfully coupled to GENRAY, an RF ray tracing code; additionally, a numerical control system to trigger the RF injection, adjustment, and shutdown in response to tearing mode activity has been developed. We discuss these accomplishments, as well as prospects for ongoing future research that this work has enabled (which continue in a limited fashion under the SciDAC Center for Extended Magnetohydrodynamic Modeling (CEMM) project and under a baseline theory grant). Associated conference presentations, published articles, and publications in progress are also listed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Role of Gravity Waves in the Formation and Organization of Clouds During TWPICE

All convective clouds emit gravity waves. While it is certain that convectively-generated waves play important parts in determining the climate, their precise roles remain uncertain and their effects are not (generally) represented in climate models. The work described here focuses mostly on observations and modeling of convectively-generated gravity waves, using the intensive observations from the DoE-sponsored Tropical Warm Pool International Cloud Experiment (TWP-ICE), which took place in Darwin, from 17 January to 13 February 2006. Among other things, the research has implications the part played by convectively-generated gravity waves in the formation of cirrus, in the initiation and organization of further convection, and in the subgrid-scale momentum transport and associated large-scale stresses imposed on the troposphere and stratosphere. The analysis shows two groups of inertia-gravity waves are detected: group L in the middle stratosphere during the suppressed monsoon period, and group S in the lower stratosphere during the monsoon break period. Waves belonging to group L propagate to the south-east with a mean intrinsic period of 35 h, and have vertical and horizontal wavelengths of about 5-6 km and 3000-6000 km, respectively. Ray tracing calculations indicate that these waves originate from a deep convective region near Indonesia. Waves belonging to group S propagate to the south-south-east with an intrinsic period, vertical wavelength and horizontal wavelength of about 45 h, 2 km and 2000-4000 km, respectively. These waves are shown to be associated with shallow convection in the oceanic area within about 1000 km of Darwin. The intrinsic periods of high-frequency waves are estimated to be between 20-40 minutes. The high-frequency wave activity in the stratosphere, defined by mass-weighted variance of the vertical motion of the sonde, has a maximum following the afternoon local convection indicating that these waves are generated by local convection. The wave activity is strongest in the lower stratosphere below 22 km and, during the suppressed monsoon period, is modulated with a 3-4-day period. The concentration of the wave activity in the lower stratosphere is consistent with the properties of the environment in which these waves propagate, whereas its 3-4-day modulation is explained by the variation of the convection activity in the TWP-ICE domain. At low rainfall intensity the wave activity increases as rainfall intensity increases. At high values of rainfall intensity, however, the wave activity associated with deep convective clouds is independent of the rainfall intensity. The convection and gravity waves observed during TWP-ICE are simulated with the Weather Research and Forecasting (WRF) Model. These simulations are compared with radiosonde observations described above and are used to determine some of the properties of convectively generated gravity waves. The gravity waves appear to be well simulated by the model. The model is used to explore the relationships between the convection, the gravity waves and cirrus.

54 ENVIRONMENTAL SCIENCES↗