Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial analysis”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT

Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.

58 GEOSCIENCES↗
DOE OSTI2025

Characterization of Coal Refuse Sites in West Virginia: Stream Loss, Volume Estimation, and Land Cover Analysis for Restoration Feasibility

Coal refuse disposal across Appalachia has resulted in widespread environmental degradation, including stream burial, landscape destabilization, and long-term hydrologic impacts. This study provides a GIS-based spatial characterization of 623 coal refuse sites in West Virginia, evaluating their potential for ecological and economic restoration. Stream loss was quantified through digitization of historic topographic maps, revealing over 194,500 meters of historically present streams lost, with 105,000 meters directly buried beneath refuse piles. Volume estimation using empirical Bayesian kriging and modern LiDAR surfaces indicated approximately 1.96 billion cubic meters of deposited refuse. Land cover analysis based on NAIP-derived classification showed that 48% of the total site area was forested, while 15% remained barren. These results highlight not only the environmental footprint of legacy refuse deposits but also their potential for restoration through stream daylighting and beneficial reuse. Emerging opportunities, including rare earth element recovery and mitigation banking, may offer economic incentives to facilitate reclamation. This work establishes a spatial framework to prioritize refuse sites for restoration, balancing ecological uplift with feasible material management strategies.

58 GEOSCIENCES↗
DOE OSTI2022

Unwrapping Radiograph Target Disassembly Using Adaptive Machine Learning [Slides]

Advanced analysis methods were developed to determine the electron temperature and density evolution of targets disassembly by intense electron beams at DARHT based on shearing interferometer and shadowgraph measurements. New methods have been developed to spatially and temporally resolve electron velocity and temperature, expanding on recently published first results of Aluminum target disassembly. A surrogate model has been created for use in developing Adaptive Machine Learning analysis techniques. The results will also be used to guide the development of the system and experimental approach to future measurements on DARHT Axis-I and a new system to be deployed on DARHT Axis-II.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Helium Ion Microscopy with Secondary Ion Mass Spectrometry for Nanoscale Chemical Imaging and Analysis of Polyolefins

To date, chemical imaging of polymers and polymer blends has been primarily accomplished using time-of-flight secondary ion mass spectrometry (ToF-SIMS) to directly visualize the distribution of components in a complex material with spatial resolution ranging from 100 nm to 5 μm. However, in many cases, this resolution falls far short of visualizing interfaces directly. To overcome these limitations, recent work has focused on developing a SIMS detection system based on a helium ion microscope (HIM) enabling chemical imaging with a demonstrated ~14 nm spatial resolution. Here, we utilize HIM-SIMS for differentiation between the olefin-based polymers of polyethylene (PE) and polypropylene (PP). We illustrate both analyses for separating PE and PP using specific mass fragment ratios as well as demonstrate spatially resolved imaging of phase-separated domains within PE. Altogether, we demonstrate the abilities of HIM-SIMS as a multimodal chemical technique for imaging and quantification of polyolefin interfaces, which could be more broadly applied to the analysis of more complex polymeric systems.

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

Cosmological constraints from higher redshift gamma-ray burst, H ii starburst galaxy, and quasar (and other) data

ABSTRACT We use higher redshift gamma-ray burst (GRB), H ii starburst galaxy (H iiG), and quasar angular size (QSO-AS) measurements to constrain six spatially flat and non-flat cosmological models. These three sets of cosmological constraints are mutually consistent. Cosmological constraints from a joint analysis of these data sets are largely consistent with currently accelerating cosmological expansion and with cosmological constraints derived from a combined analysis of Hubble parameter (H(z)) and baryon acoustic oscillation (BAO, with Planck-determined baryonic matter density) measurements. A joint analysis of the H(z) + BAO + QSO-AS + H iiG + GRB data provides fairly model-independent determinations of the non-relativistic matter density parameter $\Omega _{\rm m_0}=0.313\pm 0.013$ and the Hubble constant $H_0=69.3\pm 1.2\, \rm {km \, s^{-1} \, Mpc^{-1}}$. These data are consistent with the dark energy being a cosmological constant and with spatial hypersurfaces being flat, but they do not rule out mild dark energy dynamics or a little spatial curvature. We also investigate the effect of including quasar flux measurements in the mix and find no novel conclusions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Interlaboratory Study of Digital Volume Correlation Error Due to X-Ray Computed Tomography Equipment and Scan Parameters: an Update from the DVC Challenge

Background: The quality of Digital Volume Correlation (DVC) full-field displacement measurements depends directly on the characteristics of the X-ray Computed Tomography (XCT) equipment, and scan procedures used to acquire the tomographic images. Objective: In this work, we seek to experimentally study the effects of XCT equipment and tomographic scan procedures on the quality of these images for DVC analysis, and to survey the level of DVC error that may be achieved using standard XCT operating procedures. Methods: Six participants in an interlaboratory study acquired high-quality XCT scans of a syntactic foam before and after rigid body motion. The resulting images were correlated using commercial DVC software to quantify error sources due to random image noise, reconstruction artifacts, as well as systematic spatial or temporal distortion. Results: In the absence of rigid body motion, the standard deviation of the displacement measurements ranged from 0.012 to 0.043 voxels using a moderate subvolume size, indicating that subvoxel measurement resolution could readily be achieved with a variety of XCT equipment and scan recipes. Comparison of consecutive scans without rigid body motion showed transient dilatational displacement gradients due to self-heating of the X-ray source and/or thermal expansion of the foam. Evaluation of the scans after rigid body motion showed significant, machine-specific spatial distortion in the displacement fields of up to 0.5 voxels; new approaches to remove this error need to be developed. Conclusions: Analysis of the scan protocols used in the interlaboratory study, as well as a complementary parametric sensitivity study, showed that the DVC error was strongly influenced by the XCT equipment, but could be mitigated by adjusting the total scan duration.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Multivariate discrimination in quantum target detection

In this work, we describe a simple multivariate technique of likelihood ratios for improved discrimination of signal and background in multi-dimensional quantum target detection. Furthermore, the technique combines two independent variables, time difference and summed energy, of a photon pair from the spontaneous parametric down-conversion source into an optimal discriminant. The discriminant performance was studied in experimental data and in Monte-Carlo modelling with clear improvement shown compared to previous techniques. As novel detectors become available, we expect this type of multivariate analysis to become increasingly important in multi-dimensional quantum optics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Daily, 30 m Resolution NDSI Data for the East River Watershed, CO for 2000-2020

This dataset contains daily Normalized Difference Snow Index (NDSI) values at 30 m spatial resolution for the East River watershed in Colorado, USA. The temporal range of these data includes water years 2001-2020. These data were created using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). This model fuses low spatial and high temporal resolution data from MODIS (500 m, daily) with high spatial and low temporal resolution data from Landsat (30 m, 16 days) to create a 30m synthetic daily snow product. This product allows for the analysis of historical snow covered area trends in the East River Watershed at fine spatiotemporal resolutions where it was not available previously. This research was performed as a part of the Department of Energy’s Subsurface Biogeochemical Research Program with the primary intent of better understanding the timing and spatial patterns of water delivery to the Critical Zone in mountain watersheds. Each .zip file contains one "water year" of data (October 1 - September 30; i.e., water year 2010 starts October 1, 2010 and ends September 30, 2011). Each zip file contains the following: STARFM daily Normalized Difference Snow Index (NDSI) fusion data files in GeoTiff format with one layer for each day between Landsat data acquisition dates (i.e., for dates of Landsat acquisition, the Landsat image is included for that date). The study area is located in an area of Landsat path overlap, so Landsat dates acquisitions are every 7-9 days. Landsat NDSI files containing the high spatial (30m), low temporal (7-9 days due to Landsat path overlap) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Dates for which no Landsat data were obtained are included as NoData layers. MODIS NDSI files containing the high temporal (daily), low spatial (500m) resolution data used as input to STARFM in GeoTiff format with one layer for each day. Please note the MODIS data were resampled to 30m pixels for input into the STARFM model. The data have a scale factor of 10,000 and a no data value of -32767. The projection of all datasets is WGS 84 (EPSG: 4326), which has a latitude/longitude based degree resolution of 0.0002694946 X 0.0002694946, and approximates to the 30 m spatial resolution mentioned above. The Layer Index files in .csv format. They contain information for each layer in the above GeoTiff files regarding the corresponding date for each layer, the fraction of pixels in the image that contain valid data (missing data is due to either cloud cover or poor data quality; these values are not percent snow cover). Dates of Landsat overpass are indicated in these files. If no Landsat data were able to be obtained due to cloud cover or lack of Landsat Tier 1 data available on Google Earth Engine, this is also noted.

EARTH SCIENCE > CRYOSPHERE > SNOW/ICE↗
DOE OSTI2024

Spatially Resolved Raman Spectroscopic Investigation of Uranyl Fluoride: A Case Study in the Importance of Instrument Optimization

Raman spectroscopy is an emerging technique for rapid and nondestructive analysis of nuclear materials for forensic and nonproliferation applications as it is a powerful tool for distinguishing multiple chemical forms of materials with similar stoichiometries. Recent developments in spectroscopic software have enabled rapid data collection with high-speed Raman spectroscopic mapping capabilities. However, some uranium-rich materials are susceptible to degradation in humid air and/or laser-induced phase transformations. To mitigate environmental or measurement-related sample degradation of potential samples of interest, we have taken a systematic approach to define optimized data collection parameters for high-throughput measurements of uranyl fluoride (UO 2 F 2 ), which is an important intermediate material in the nuclear fuel cycle. First, we systematically describe the influence of optical magnification (5× to 100×), laser power, and exposure time on obtained signal for identical particles of UO 2 F 2 and find that at low laser power and exposure times, comparable signal is obtained regardless of optical magnification. Second, we ensure sample integrity during data collection, and third, collect spectroscopic maps that employ optimized parameters to reduce the time required to obtain spatially resolved spectroscopic information. Reductions of 90% and 99% in measurement times are discussed as they relate to differences in resolving spectroscopic features of particles in identical mapping areas. Finally, during this work, we found that additional data processing options were needed and thus developed a customized Python script for importing, processing, analyzing, and visualizing Raman spectroscopic map data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Modeling Radionuclide Inventories in MSR Off-Gas Systems with Radiochemical Transport Analysis

Molten salt reactors (MSR) contain unique characteristics that may require enhancements to modeling tools to accurately predict phenomena. One characteristic that may be advantageous to leverage during normal operation is on-line processing of the circulating fuel salt, such as an off-gas system (OGS) to remove volatile fission products. Therefore, new modeling tools must be developed to integrate spatial resolution and chemistry effects into fuel depletion tools to be able to account for these non-core sources of radioactivity. Such types of radiochemical transport analysis tools were used to estimate the removal rates for 12 elements within a flow model of the Molten Salt Reactor Experiment (MSRE) by optimizing against legacy experimental data of the gas-borne (GB) percentages of 12 nuclides. The removal rates were used in a depletion model to calculate the FP inventory that enters the MSRE OGS. Calculations are in good agreement with the empirical GB percentages reported for the 12 nuclides, which validates the approach and verifies each tool’s treatment of the radiochemical flow effects. The OGS inventory is discussed in terms of the largest nuclide contributors to activity, dose consequence, decay heat, and elemental composition. Finally, insights from the study allow recommendations to be made for future code development activities.

Shahbazi, Shayan↗
DOE OSTI2020

Robust framework and software implementation for fast speciation mapping

One of the greatest benefits of synchrotron radiation is the ability to perform chemical speciation analysis through X-ray absorption spectroscopies (XAS). XAS imaging of large sample areas can be performed with either full-field or raster-scanning modalities. A common practice to reduce acquisition time while decreasing dose and/or increasing spatial resolution is to compare X-ray fluorescence images collected at a few diagnostic energies. In this work, several authors have used different multivariate data processing strategies to establish speciation maps. Furthermore, the theoretical aspects and assumptions that are often made in the analysis of these datasets are focused on. A robust framework is developed to perform speciation mapping in large bulk samples at high spatial resolution by comparison with known references. Two fully operational software implementations are provided: a user-friendly implementation within the MicroAnalysis Toolkit software, and a dedicated script developed under the R environment. The procedure is exemplified through the study of a cross section of a typical fossil specimen. Additionally, the algorithm provides accurate speciation and concentration mapping while decreasing the data collection time by typically two or three orders of magnitude compared with the collection of whole spectra at each pixel. Whereas acquisition of spectral datacubes on large areas leads to very high irradiation times and doses, which can considerably lengthen experiments and generate significant alteration of radiation-sensitive materials, this sparse excitation energy procedure brings the total irradiation dose greatly below radiation damage thresholds identified in previous studies. This approach is particularly adapted to the chemical study of heterogeneous radiation-sensitive samples encountered in environmental, material, and life sciences.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

A robust approach to Gaussian process implementation

Abstract. Gaussian process (GP) regression is a flexible modeling technique used to predict outputs and to capture uncertainty in the predictions. However, the GP regression process becomes computationally intensive when the training spatial dataset has a large number of observations. To address this challenge, we introduce a scalable GP algorithm, termed MuyGPs, which incorporates nearest-neighbor and leave-one-out cross-validation during training. This approach enables the evaluation of large spatial datasets with state-of-the-art accuracy and speed in certain spatial problems. Despite these advantages, conventional quadratic loss functions used in the MuyGPs optimization, such as root mean squared error (RMSE), are highly influenced by outliers. We explore the behavior of MuyGPs in cases involving outlying observations and, subsequently, develop a robust approach to handle and mitigate their impact. Specifically, we introduce a novel leave-one-out loss function based on the pseudo-Huber function (LOOPH) that effectively accounts for outliers in large spatial datasets within the MuyGPs framework. Our simulation study shows that the LOOPH loss method maintains accuracy despite outlying observations, establishing MuyGPs as a powerful tool for mitigating unusual observation impacts in the large data regime. In the analysis of US ozone data, MuyGPs provides accurate predictions and uncertainty quantification, demonstrating its utility in managing data anomalies. Through these efforts, we advance the understanding of GP regression in spatial contexts.

Mukangango, Juliette↗
DOE OSTI2021

A field–particle correlation analysis of a perpendicular magnetized collisionless shock

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Spatially Resolved Near Field Spectroscopy of Vibrational Polaritons at the Small N Limit

Vibrational polaritons, which have been primarily studied in Fabry–Pérot cavities with a large number of molecules (N ~ 10 6 –10 10 ) coupled to the resonator mode, exhibit various experimentally observed effects on chemical reactions. However, the exact mechanism is elusively understood from the theoretical side, as the large number of molecules involved in an experimental strong coupling condition cannot be represented completely in simulations. This discrepancy between theory and experiment arises from computational descriptions of polariton systems typically being limited to only a few molecules, thus failing to represent the experimental conditions adequately. To address this mismatch, we used surface phonon polariton (SPhP) resonators as an alternative platform for vibrational strong coupling. SPhPs exhibit strong electromagnetic confinement on the surface and thus allow for coupling to a small number of molecules. As a result, this platform can enhance nonlinearity and slow down relaxation to the dark modes. In this study, we fabricated a pillar-shaped quartz resonator and then coated it with a thin layer of cobalt phthalocyanine (CoPc). By employing scattering-type scanning near-field optical microscopy (s-SNOM), we spatially investigated the dependency of vibrational strong coupling on the spatially varying electromagnetic field strength and demonstrated strong coupling with 38,000 molecules only–reaching to the small N limit. Through s-SNOM analysis, we found that strong coupling was observed primarily on the edge of the quartz pillar and the apex of the s-SNOM tip, where the maximum field enhancement occurs. In contrast, a weak resonance signal and lack of coupling were observed closer to the center of the pillar. This work demonstrates the importance of spatially resolved polariton systems in nanophotonic platforms and lays a foundation to explore polariton chemistry and chemical dynamics at the small N limit–one step closer to reconcile with high-level quantum calculations.

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

Focused‐beam X‐ray fluorescence and diffraction microtomographies for mineralogical and chemical characterization of unsectioned extraterrestrial samples

Abstract This study describes the application of new synchrotron X‐ray fluorescence (XRF) and diffraction (XRD) microtomographies for the 3‐D visualization of chemical and mineralogical variations in unsectioned extraterrestrial samples. These improved methods have been applied to three compositionally diverse chondritic meteorite samples that were between 300 and 400 μm in diameter, including samples prepared from fragments of the CR2 chondrite LaPaz Icefield (LAP) 02342, H5 chondrite MacAlpine Hills (MAC) 88203, and the CM2 chondrite Murchison. The synchrotron‐based XRF and XRD tomographies used are focused‐beam techniques that measure the intensities of fluorescent and diffracted X‐rays in a sample simultaneously during irradiation by a high‐energy microfocused incident X‐ray beam. Measured sinograms of the emitted and diffracted intensities were then tomographically reconstructed to generate 2‐D slices of XRF and XRD intensity through the sample, with reconstructed pixel resolution of 1–2 μm, defined by the resolution of the focused incident X‐ray beam. For sample LAP 02342, primary mineral phases that were visualized in reconstructed slices using these techniques included isolated grains of α‐Fe, orthopyroxene, and olivine. For our sample of MAC 88203, XRF/XRD tomography allowed visualization of forsteritic olivine as a primary mineral phase, a vitrified fusion crust at the sample surface, identification of localized Cr‐rich spinels at spatial resolutions of several micrometers, and imaging of a plagioclase‐rich glassy matrix. In the sample of Murchison, major identifiable phases include clinoenstatite‐ and olivine‐rich chondrules, variable serpentine matrix minerals and small Cr‐rich spinels. Most notable in the tomographic analysis of Murchison is the ability to quantitatively distinguish and visualize the complex mixture of serpentine‐group minerals and associated tochilinite–cronstedtite intergrowths. These methods provide new opportunities for spatially resolved characterization of sample texture, mineralogy, crystal structure, and chemical state in unsectioned samples. This provides researchers an ability to characterize such samples internally with minimal disruption of sample micro‐structures and chemistry, possibly without the need for sample extraction from some types of sampling and capture media.

Geochemistry & Geophysics↗
DOE OSTI2026

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗
DOE OSTI2024

Health and air pollutant emission impacts of net zero CO2 by 2050 scenarios from the energy modeling forum 37 study

Carbon dioxide and non-greenhouse gas air pollutants are emitted from many of the same sources. Decarbonization actions thus typically yield air pollutant emission reductions, resulting in significant air quality benefits. Although several studies have highlighted this connection, including in the context of net zero carbon emission targets, substantial uncertainty remains regarding how alternative technological pathways to this goal will affect the spatial distribution and magnitude of air pollutants. Comprehensive multi-model and multi-scenario analyzes are needed to explore the relative impacts of alternative pathways. Here, our study begins to address this gap by leveraging the results from the recent Energy Modeling Forum 37 inter-model comparison exercise on U.S. decarbonization pathways. Comparing the results of the six teams who submitted air pollutant emissions suggests that strategies that target net zero U.S. carbon emissions would yield significant reductions in many air pollutants, and that this finding is generally robust across pathways. However, some energy sources, such as biomass and fossil fuels with carbon capture, will emit air pollutants and can potentially influence the magnitude, spatial distribution, and even sign of localized air pollutant emission changes. In the second part of this analysis, a simplified air quality and health impacts screening model is used to evaluate the air quality impacts in 2035 of sectoral emission changes from the three models that provided sectoral detail. Relative to a reference scenario, a net zero pathway is estimated to reduce fine particulate matter concentrations across the contiguous U.S., with health benefits from reduced mortality ranging from $\$$65 billion to $\$$250 billion in 2035 alone (2023$\$$s). These benefits would be expected to grow over time as the net zero trajectory becomes more stringent. Both the magnitude of potential benefits and the substantial variation of the projections across models underscore the need for an EMF-like inter-model comparison exercise focused on air quality.

Air pollutants↗
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↗