Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geometrical optics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Geometrical optics without singularities: using the ray time as the coordinate space

Geometrical optics (GO) is widely used for reduced modelling of waves in plasmas, but it fails near reflection points, where it predicts a spurious singularity of the wave amplitude. We show how to avoid this singularity by adopting a different representation of the wave equation. Instead of the physical coordinate 𝑥 and the wavevector 𝑘, we use the ray time 𝜏 as the new canonical coordinate and the ray energy ℎ as the associated canonical momentum. To derive the envelope equation in the 𝜏-representation, we construct the Weyl symbol calculus on the (𝜏,ℎ) space and show that the corresponding Weyl symbols are related to their (𝑥,𝑘) counterparts by the Airy transform. This allows us to express the coefficients in the envelope equation through the known properties of the original dispersion operator. When necessary, solutions of this equation can be mapped to the 𝑥-space using a generalised metaplectic transform. However, the field per se might not even be needed in practice. Instead, knowing the corresponding Wigner function usually suffices for linear and quasilinear calculations. As a Weyl symbol itself, the Wigner function can be mapped analytically, using the aforementioned Airy transform. We show that the standard Airy patterns that form in regions where conventional GO fails are successfully reproduced within metaplectic GO (MGO) simply by remapping the field from the 𝜏-space to the 𝑥-space. An extension to mode-converting waves is also presented. This formulation, which we call generalised MGO, can be particularly useful, for example, for reduced modelling of the O–X conversion in inhomogeneous plasma near the critical density, an effect that is important for fusion applications and also occurs in the ionosphere. Overall, MGO can replace GO for any practical purposes, because it better handles cutoffs and is similar otherwise.

plasma waves↗
DOE OSTI2025

Combining geometric-optical and spectral invariants theories for modeling canopy fluorescence anisotropy

The spectral invariants theory ( p -theory) has received much attention in the field of quantitative remote sensing over the past few decades and has been adopted for modeling of canopy solar-induced chlorophyll fluorescence (SIF). However, the spectral invariant properties (SIP) in simple analytical formulae have not been applied for modeling canopy fluorescence anisotropy primarily because they are parameterized in terms of leaf total scattering, which precludes the differentiation between forward and backward leaf SIF emissions. In this study, we have developed the canopy-SIP SIF model by combining geometric-optical (GO) theory to account for asymmetric leaf SIF forward and backward emissions at the first-order scattering and by modeling multiple scattering based on the p-theory, thus avoiding the dependence on radiative transfer models. The applicability of the model simulations especially over 3D heterogeneous canopies was improved by incorporating canopy structure through multi-angular clumping index, and by modeling single scattering from the four components of the scene in view according to the GO approach. The results show good consistency with both the state-of-the-art SIF models and multi-angular field SIF observations over grass and chickpea canopies. Further, the coefficient of determination (R²) between the simulated SIF and field measurements was 0.75 (red) and 0.74 (far-red) for chickpea, and 0.65 (both red and far-red) for grass. The average relative error was approximately 3% for 1D homogeneous scenes when comparing the canopy-SIP SIF model simulations to the SCOPE model simulations, and around 4% for the 3D heterogeneous scene when comparing to the LESS model simulations. The results indicate that the proposed approach for separating asymmetric leaf SIF emissions is a robust way to keep a balance between satisfactory simulation accuracy and efficiency. Model simulations suggest that neglecting the leaf SIF asymmetry can lead to an underestimation of canopy red SIF by 6.3% to 42.6% for various leaf biochemical and canopy structural parameters. This study presents a simple but efficient analytical approach for canopy fluorescence modeling, with potential for large-scale canopy fluorescence simulations.

3D heterogeneous structure↗
DOE OSTI2026

Geometric optics model of Thomson scattering enhanced by parametric coupling

We formulate a theory of Thomson scattering (TS) from inhomogeneous plasmas using a ray-tracing approach. The scattered radiation is described within the eikonal approximation, including convective gain arising from scattering instabilities driven by the TS probe laser. In the limit of a weak probe and vanishing gain, the theory recovers the standard TS cross section and scattered power due to volumetric particle noise. The ray-tracing formalism is implemented in a TS fitting framework, enabling accurate analysis of inhomogeneous plasmas whose background profiles can be provided by large-scale hydrodynamic simulations or experimental diagnostics. Here, several examples of synthetic TS spectra are presented, and further extensions of the theory are discussed.

Carleton, D. E. [University of Alberta, Edmonton, ↗
DOE OSTI2025

Measuring the Burgers vector of dislocations with dark-field X-ray microscopy

The subsurface dynamics of dislocations are essential to many properties of bulk crystalline materials. However, it is challenging to characterize a bulk crystal by conventional transmission electron microscopy (TEM) due to the limited penetration depth of electrons. A novel X-ray imaging technique – dark-field X-ray microscopy (DFXM) – was developed to image hierarchical dislocation structures in bulk crystals. While today's DFXM can effectively map the line structures of dislocations, it is still challenging to quantify the Burgers vectors, the key characterization governing the dislocation behaviors. Here, we extend the 'invisibility criterion' formalism from the TEM theory to the geometrical-optics model of DFXM and demonstrate the consistency between DFXM and dark-field TEM using multi-diffraction-peak imaging for a single edge dislocation. Due to the practical difficulty of multi-peak DFXM experiments, we further study how the Burgers vector effect is encoded for a single-peak DFXM experiment. Using the geometrical-optics DFXM simulation, we explore the asymmetry of rocking tilt scans at different rolling tilts and develop a new method to characterize the Burgers vector. The conclusions of this study advance our understanding of the use of DFXM in characterizing individual dislocations, enabling the connection from bulk DFXM imaging to dislocation mechanics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Theoretical Benchmark of the Geometric and Optical Properties for 3d Transition Metal Nanoclusters via Density Functional Theory

Understanding structure–property relationships in atomically precise metal nanoclusters is vital in finding selective and tunable catalysts. Here, in this study, density functional theory (DFT) was used to benchmark seven exchange correlation functionals at different basis sets for 17 atomically precise nanoclusters against experimentally determined geometries, band gaps, and optical gaps. The set contains both monometallic and bimetallic clusters that possess at least two types of 3d transition metals (specifically, Cu, Ni, Fe, or Co). The benchmark highlights that PBE0 is a good functional to use regardless of the basis set, and Minnesota functionals do well with respect to specific metals. Further, while long-range corrected functionals overestimate band and optical gaps, they model absorption features better than the other considered functionals. The study additionally looks at the photoinduced hydrogen evolution reaction (HER) and the CO 2 reduction mechanism on nanoclusters reported from the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

A confocal laser-induced fluorescence diagnostic with an annular laser beam

Here, in this work, we report an annular beam confocal laser-induced fluorescence (LIF) configuration, which allows for high spatial resolution measurements of plasma properties in plasma setups and sources with limited optical access. The proposed LIF configuration utilizes the annular laser beam generated by a pair of diffractive axicons. The LIF signal is collected along the main optical axis within the ring region. It is shown experimentally that at a focal distance of 300 mm, a spatial resolution of ~5.3 mm can be achieved. Using geometric optics estimations, we showed that ~1 mm resolution at the same focal distance could potentially be achieved by modifying laser beam parameters. This approaches the localization accuracy of conventional LIF collection methods (with crossing laser beam injection and fluorescence collection optical paths). Measurements of the ion velocity distribution function in an argon plasma using both the confocal LIF with an annular laser beam and conventional LIF demonstrate a satisfactory agreement. The proposed LIF setup has potential applications for diagnostics in various plasma processing equipment and plasma sources, such as hollow cathodes, microplasmas, electric propulsion, etc.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Quasi-optical beam tracing module development for millimeter-wave high-wavenumber collective scattering on the NSTX-U and EAST tokamaks

A Python3-based beam tracing code utilizing Quasi-Optics has been developed to track both incident and receiving beams in high-k collective millimeter wave scattering systems within magnetic fusion plasmas. In contrast to existing ray tracing codes that solely consider refraction, this beam tracing code incorporates diffraction phenomena, providing a more comprehensive calculation. Here, this enhanced capability allows for a more accurate calculation of the scattering volume and spatial resolution in high-k collective scattering systems, crucial for evaluating system performance and facilitating data analysis. Unlike Geometrical Optics, Quasi-Optics employs the complex eikonal method, representing a Gaussian beam as a collection of coupled rays to accurately preserve diffraction characteristics. The developed code is intended for application in NSTX-Upgrade and EAST high-k beam tracing analyses, targeting frequencies of 693 GHz and 270 GHz, respectively. The high-k system's primary objective is the observation of electron-scale instabilities. Employing a symplectic integrator, the code ensures numerical accuracy, assessed through the conservation of the Hamiltonian. With its precision and efficiency, the code facilitates rapid inter-shot analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The role of stereochemistry in combustion processes

Stereochemical effects significantly influence chemical processes, yet it is not well understood if they are a leading source of uncertainty in combustion modeling. Stereochemistry influences a combustion model (i) at the earliest stage of its construction when mapping the reaction network, (ii) in the computation of individual thermochemical and rate parameters, and (iii) in the prediction of combustion observables. The present work reviews the importance of enumerating stereochemical species and reactions at each of these steps. Further, it analyzes the separate influence of several types of stereochemistry, including geometric, optical, and fleeting transition state diastereomers. Three reaction networks serve to examine which stages of low-temperature oxidation are most affected by stereochemistry, including the first and second oxidation of n-butane, the third oxidation of n-pentane, and the early stages of pyrolysis of 1- and 2-pentene. The 149 reactions in the n-butane mechanism are expanded to 183 reactions when accounting for diastereomerism. Each of these 183 reactions is parameterized with ab initio kinetics computations to determine that, for the n-butane mechanism, the median factor of diastereomeric deviation is 3.5 at 360 K for rate constants and as high as 1.6 for mechanism reactivity, in terms of ignition delay times, as opposed to a mechanism without stereochemical expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Device response principles and the impact on energy resolution of epitaxial quantum dot scintillators with monolithic photodetector integration

Abstract Epitaxial quantum dot (QD) scintillator crystals with picosecond-scale timing and high light yield have been created for medical imaging, high energy physics and national security applications. Monolithic photodetector (PD) integration enables the sensing of photons generated within the waveguiding crystal and allows a wide range of scintillator-photodetector coupling geometries. Until recently, these doubly novel devices have suffered from complex, high variance responses to monoenergetic sources which significantly reduces their precision and accuracy. The principles governing the overall device response have now been discerned and embodied by an expression derived within a geometrical optics framework which considers optical properties, surface roughness and photodetector coupling geometry. Response variation due to these factors was sufficiently reduced to obtain material-related energy resolution values of 2.4% with alpha particles. These findings place energy resolution alongside luminescence timescale, photon yield, and radiation hardness as outstanding properties of these engineered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gauge-invariant gravitational waves in matter beyond linearized gravity

Modeling the propagation of gravitational waves (GWs) in media other than vacuum is complicated by the gauge freedom of linearized gravity in that, once nonlinearities are taken into consideration, gauge artifacts can cause spurious acceleration of the matter. To eliminate these artifacts, we propose how to keep the theory of dispersive GWs gauge-invariant beyond the linear approximation and, in particular, obtain an unambiguous gauge-invariant expression for the energy–momentum of a GW in a dispersive medium. Using analytic tools from plasma physics, we propose an exactly gauge-invariant 'quasilinear' theory, in which GWs are governed by linear equations and also affect the background metric on scales large compared to their wavelength. As a corollary, the gauge-invariant geometrical optics of linear dispersive GWs in a general background is formulated. As an example, we show how the well-known properties of vacuum GWs are naturally and concisely yielded by our theory in a manifestly gauge-invariant form. We also show how the gauge invariance can be maintained within a given accuracy to an arbitrary order in the GW amplitude. These results are intended to form a physically meaningful framework for studying dispersive GWs in matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Self-consistent interaction of linear gravitational and electromagnetic waves in non-magnetized plasma

Here, this paper explores the hybridization of linear metric perturbations with linear electromagnetic (EM) perturbations in non-magnetized plasma for a general background metric. The local wave properties are derived from first principles for inhomogeneous plasma, without assuming any symmetries of the background metric. First, we derive the effective (“oscillation-center”) Hamiltonian that governs the average dynamics of plasma particles in a prescribed quasimonochromatic wave that involves metric perturbations and EM fields simultaneously. Then, using this Hamiltonian, we derive the backreaction of plasma particles on the wave itself and obtain gauge-invariant equations that describe the resulting self-consistent gravito-electromagnetic (GEM) waves in a plasma. The transverse tensor modes of gravitational waves are found to have no interaction with the plasma and the EM modes in the geometrical-optics limit. However, for $longitudinal$ GEM modes with large values of the refraction index, the interplay between gravitational and EM interactions in plasma can have a strong effect. In particular, the dispersion relation of the Jeans mode is significantly affected by electrostatic interactions. As a spin-off, our calculation also provides an alternative resolution of the so-called Jeans swindle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Computing virtual dark-field X-ray microscopy images of complex discrete dislocation structures from large-scale molecular dynamics simulations

Dark-field X-ray microscopy (DFXM) is a novel diffraction-based imaging technique that non-destructively maps the local deformation from crystalline defects in bulk materials. While studies have demonstrated that DFXM can spatially map 3D defect geometries, it is still challenging to interpret DFXM images of the high-dislocation-density systems relevant to macroscopic crystal plasticity. This work develops a scalable forward model to calculate virtual DFXM images for complex discrete dislocation structure(s) (DDS) obtained from atomistic simulations. Our new DDS-DFXM model integrates a non-singular formulation for calculating the local strain from the DDS and an efficient geometrical optics algorithm for computing the DFXM image from the strain field. We apply the model to complex DDS obtained from a large-scale mol­ecular dynamics simulation of compressive loading on single-crystal silicon. Simulated DFXM images exhibit prominent contrast for dislocation features between the multiple slip systems, demonstrating the potential of DFXM to resolve features from dislocation multiplication. In conclusion, the integrated DDS-DFXM model provides a toolbox for DFXM experimental design and image interpretation in the context of bulk crystal plasticity for a range of measurements across shock plasticity and the broader materials science community.

X-ray imaging↗
DOE OSTI2023

Optical information transfer through random unknown diffusers using electronic encoding and diffractive decoding

Free-space optical information transfer through diffusive media is critical in many applications, such as biomedical devices and optical communication, but remains challenging due to random, unknown perturbations in the optical path. We demonstrate an optical diffractive decoder with electronic encoding to accurately transfer the optical information of interest, corresponding to, e.g., any arbitrary input object or message, through unknown random phase diffusers along the optical path. This hybrid electronic-optical model, trained using supervised learning, comprises a convolutional neural network-based electronic encoder and successive passive diffractive layers that are jointly optimized. After their joint training using deep learning, our hybrid model can transfer optical information through unknown phase diffusers, demonstrating generalization to new random diffusers never seen before. The resulting electronic-encoder and optical-decoder model was experimentally validated using a 3D-printed diffractive network that axially spans <70λ, where λ = 0.75 mm is the illumination wavelength in the terahertz spectrum, carrying the desired optical information through random unknown diffusers. The presented framework can be physically scaled to operate at different parts of the electromagnetic spectrum, without retraining its components, and would offer low-power and compact solutions for optical information transfer in free space through unknown random diffusive media.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analytical Meridional, Non-Paraxial Ray Tracing

We present an analytical method for non-paraxial, meridional ray tracing through a system of basic geometric surfaces derivable from a generalized ellipse; i.e. the most general ray tracing method for a 2-dimensional system comprised of elliptical, circular, and linear surfaces. This work is meant as a non-small angle generalization of standard ray tracing for geometric optics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2025

Mathematical Foundation for Quantum Computing of Electromagnetic Wave Propagation in Dielectric Media

Can quantum computers effectively simulate the propagation and scattering of electromagnetic waves in a classical plasma? This chapter introduces some of the basic concepts in mathematics and physics essential to answering that question. The numerical simulations of Maxwell equations for wave propagation in dielectrics are constrained by technological limitations of the present-day computers. In contrast, there has been ample fanfare around quantum computers and their potential to far exceed the performance of traditional computers. Whether the enhanced capabilities of a quantum computer can be put to use for simulating topics in classical physics is a source of intrigue and curiosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Celestial gluon amplitudes from the outside in

We show that, given a two-dimensional realization of the celestial OPE in self-dual Yang-Mills, we can find a scalar source around which scattering amplitudes replicate correlation functions computed from the 2D ‘gluon’ operators in a limit where a dynamic massless scalar decouples. We derive conditions on the two-dimensional three-point correlation function so that such a source exists and give two particular examples of this construction, one in which gluons are constructed from vertex operators in the semiclassical limit of Liouville theory and another in which the soft gluons arise from generalized free fields. Finally, we identify a bulk dual to the level of the boundary Kac-Moody algebra and discuss moving beyond the decoupling limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Massive twistor worldline in electromagnetic fields

We study the (ambi-)twistor model for spinning particles interacting via electromagnetic field, as a toy model for studying classical dynamics of gravitating bodies including effects of both spins to all orders. We compute the momentum kick and spin kick up to one-loop order and show precisely how they are encoded in the classical eikonal. The all-orders-in-spin effects are encoded as a dynamical implementation of the Newman-Janis shift, and we find that the expansion in both spins can be resummed to simple expressions in special kinematic configurations, at least up to one-loop order. We confirm that the classical eikonal can be understood as the generator of canonical transformations that map the in-states of a scattering process to the out-states. We also remark that cut contributions for converting worldline propagators from time-symmetric to retarded amount to the iterated action of the leading eikonal at one-loop order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗