Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverse 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Uncertainty Quantification in GADRAS Inverse Modeling

The Gamma Detector Response and Analysis Software (GADRAS) package includes an inverse modeling tool that is helpful in identifying characteristics of unknown radioactive materials. Traditionally, uncertainties in this analysis were derived solely from measurement data quality and the fit of synthetic spectra. This paper aims to rigorously quantify additional sources of uncertainty, focusing on uncertainties arising from measurements being analyzed, Detector Response Function (DRF) characterization, and DRF extrapolation. Applying these findings to the BeRPBall benchmark data set, we demonstrated the impact of these uncertainties on plutonium and polyethylene estimates. The results underscore the importance of incorporating diverse uncertainty sources to enhance the accuracy and reliability of GADRAS’s inverse modeling capabilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗
DOE OSTI2024

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Development of the TREAT M2 experiment as a transient benchmark

A transient benchmark based on the Transient Reactor Test Facility (TREAT) M2 Calibration experiment (M2-CAL) is under development. TREAT, at Idaho National Laboratory, is a graphite moderated air-cooled research reactor which has been used extensively for fuel material testing under extreme and accident conditions. Accurate benchmark models are a beneficial component in the operations, experimental planning, and development of TREAT. In this work, we present a transient benchmark model for the M2-CAL experiment core loading. The benchmark model incorporates the coupling between the Monte Carlo Code SERPENT and the computational fluid dynamics code OpenFOAM to capture the temperature feedback mechanism. The M2-CAL transient 2580 was simulated in this work. A pre-transient analysis was performed to determine the optimum core composition and conditions before the beginning of the transient. The analysis was validated against historic TREAT kinetic measurements and the worth of the transient rod T-2. The axial power distribution in the flux wire for the M2-CAL experiment was determined and contrasted with the experiment. The model was then used to simulate the M2-CAL transient 2580 experiment based on the reported pre-transient and transient conditions. Several transient observables were calculated and compared to the experiment. The model shows good agreement with the experimental power traces, period, and average increase of power at the power ramp (time > 7.8 s). Inverse point kinetics analysis was performed during the period of the power ramp for the model and the experiment based on that, the temperature feedback component was isolated and contrasted with the experimental feedback. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Turbulence and Cloudiness in Cumulus-Topped Marine Boundary Layers

Turbulence and cloudiness in cumulus-topped marine boundary layers are studied using data from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site. The analysis includes eight periods of nonprecipitating shallow cumulus clouds, spanning 141 h and encompassing 603 individual clouds. On average, the cumulus had bases at 558 m, were 99 m thick, had a chord length of 500 m, and exhibited an hourly base-layer cloudiness of 12%. Changes in cloud fraction were primarily driven by variations in cloud number not chord length. High-resolution Doppler lidar and cloud radar observations were used to estimate updraft and downdraft mass fluxes in the cumulus base layer under both clear and cloudy regions. Cloudy updrafts contributed only ∼23% of the total updraft mass flux at cloud base, indicating that most upward mass flux originated from clear-air eddies. Updraft strength, rather than updraft fraction, was found to predominantly control both clear and cloudy updraft mass fluxes. Cloud-base cloudiness strongly correlated with cloudy updraft mass flux but showed negligible correlation with clear-air updraft mass flux. Clear-air downdraft mass flux exhibited a strong relationship with the ratio of surface buoyancy to inversion strength. Mesoscale analysis revealed that moist patches had greater cloudiness, updraft mass flux, and vertical velocity variance compared to dry patches. Additionally, mesoscale base-layer cloud fraction was highly and significantly correlated with cloudy updraft mass flux in both dry and moist environments. Results presented herein have implications for cumulus parameterization development along with climatological and model evaluation studies conducted at the ENA site.

cumulus clouds↗
DOE OSTI2023

On unifying randomized methods for inverse problems

This work unifies the analysis of various randomized methods for solving linear and nonlinear inverse problems with Gaussian priors by framing the problem in a stochastic optimization setting. By doing so, we show that many randomized methods are variants of a sample average approximation (SAA). More importantly, we are able to prove a single theoretical result that guarantees the asymptotic convergence for a variety of randomized methods. Additionally, viewing randomized methods as an SAA enables us to prove, for the first time, a single non-asymptotic error result that holds for randomized methods under consideration. Another important consequence of our unified framework is that it allows us to discover new randomization methods. Here, we present various numerical results for linear, nonlinear, algebraic, and PDE-constrained inverse problems that verify the theoretical convergence results and provide a discussion on the apparently different convergence rates and the behavior for various randomized methods.

42 ENGINEERING↗
DOE OSTI2020

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Electrical Resistivity Changes During Heating Experiments Unravel Heterogeneous Thermal‐Hydrological‐Mechanical Processes in Salt Formations

Abstract Rock salt is considered a suitable medium for the permanent disposal of heat‐generating radioactive waste due to its isolation properties. However, excavation damage and heating induce complex and heterogeneous thermal‐hydrological‐mechanical (THM) processes across different zones. Quantifying this heterogeneity is crucial for accurate long‐term performance assessment models, but traditional methods lack the necessary resolution. This study employs 4D electrical resistivity tomography (ERT) monitoring during controlled heating experiments in a salt formation to unravel the spatiotemporal dynamics of THM processes. Advanced time‐lapse inversion and clustering analysis quantify subsurface properties and map the heterogeneity of THM dynamics. The ERT results can estimate subsurface properties and delineate the damaged and intact zones, enabling appropriate parameterization and representation of processes for long‐term modeling. This approach may be used in further improving the predictive models and ensuring the safe long‐term disposal of radioactive waste in rock salt.

58 GEOSCIENCES↗
DOE OSTI2021

Deformation of binary and boron-doped Ni 3 Al alloys at high pressures studied with synchrotron x-ray diffraction

In situ x-ray synchrotron diffraction experiments were carried out on nickel-based high-strength superalloys under pressure to understand their deformation mechanism using a diamond anvil cell (DAC). Radial x-ray diffraction determines the room-temperature equations of state and yield strengths of binary Ni 3 Al alloy and 500 ppm boron-doped Ni 3 Al to pressures of 20 and 46 GPa, respectively. Crystallographic preferred orientations observed in these superalloys due to anisotropic stress field in DAC indicate the onset of plastic deformation. Inverse pole figure analysis reveals that the underlying deformation mechanisms change from an octahedral slip to a simultaneous activation of octahedral and cube slips upon doping with boron. The yield-strength values were found to increase with pressure and are comparable to those determined from axial diffraction experiments. The results indicate that the yield strength of Ni 3 Al:B is about 0.5 GPa higher (at pressures below 20 GPa) due to grain boundary strengthening by boron. It is shown that due to high elastic anisotropy of Ni 3 Al alloy, the yield-strength estimations from diffraction experiments strongly depend on the micromechanical model used to convert the measured elastic strains into stresses.

42 ENGINEERING↗
DOE OSTI2023

Application of Electromagnetic Methods for Reservoir Monitoring with Emphasis on Carbon Capture, Utilization, and Storage

The Controlled-Source ElectroMagnetic (CSEM) method provides crucial information about reservoir fluids and their spatial distribution. Carbon dioxide (CO2) storage, enhanced oil recovery (EOR), geothermal exploration, and lithium exploration are ideal applications for the CSEM method. The versatility of CSEM permits its customization to specific reservoir objectives by selecting the appropriate components of a multi-component system. To effectively tailor the CSEM approach, it is essential to determine whether the primary target reservoir is resistive or conductive. This task is relatively straightforward in CO2 monitoring, where the injected fluid is resistive. However, for scenarios involving brine-saturated (water-wet) or oil-wet (carbon capture, utilization, and storage—CCUS) reservoirs, consideration must also be given to conductive reservoir components. The optimization of data acquisition before the survey involves analyzing target parameters and the sensitivity of multi-component CSEM. This optimization process typically includes on-site noise measurements and 3D anisotropic modeling. Based on our experience, subsequent surveys tend to proceed smoothly, yielding robust measurements that align with scientific objectives. Other critical aspects to be considered are using magnetotelluric (MT) measurements to define the overall background resistivities and integrating real-time quality assurance during data acquisition with 3D modeling. This integration allows the fine tuning of acquisition parameters such as acquisition time and necessary repeats. As a result, data can be examined in real-time to assess subsurface information content while the acquisition is ongoing. Consequently, high-quality data sets are usually obtained for subsequent processing and initial interpretation with minimal user intervention. The implementation of sensitivity analysis during the inversion process plays a pivotal role in ensuring that the acquired data accurately respond to the target reservoirs’ expected depth range. To elucidate these concepts, we present an illustrative example from a CO2 storage site in North Dakota, USA, wherein the long-offset transient electromagnetic method (LOTEM), a variation of the CSEM method, and the MT method were utilized. This example showcases how surface measurements attain appropriately upscaled log-scale sensitivity. Furthermore, the sensitivity of the CSEM and MT methods was examined in other case histories, where the target reservoirs exhibited conductive properties, such as those encountered in enhanced oil recovery (EOR), geothermal, and lithium exploration applications. The same equipment specifications were utilized for CSEM and MT surveys across all case studies.

Barajas-Olalde, César↗
DOE OSTI2020

Calculation of Critical Heat Flux Using an Inverse Heat Transfer Method to Support TREAT Experiment Analysis

Heat transfer between cladding and coolant during transient scenarios remains a critical area of uncertainty in understanding nuclear reactor safety. To advance the understanding of transient and accident scenarios involving critical heat flux (CHF), an in-pile experiment for the Transient Reactor Test (TREAT) facility at Idaho National Laboratory was developed. The experiment named CHF-Static Environment Rodlet Transient Test Apparatus (CHF-SERTTA), consists of a hollow borated stainless steel heater rod submerged in a static water pool which is heated via the (n, a) reaction in boron-10. This paper presents a novel inverse heat transfer method to determine CHF by using the optimization and uncertainty software DAKOTA to calibrate a RELAP5-3D model of CHF-SERTTA to temperature measurements obtained from a thermocouple welded to the surface of the rod.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantitative texture analysis at the WAND 2 and HIDRA diffractometers

Data collection and analysis strategies have been developed for efficient and reliable crystallographic texture measurements at two recently upgraded neutron diffractometers: the Wide Angle Neutron Diffractometer Squared (WAND 2 ) and the High Intensity Diffractometer for Residual Stress Analysis (HIDRA) at the High Flux Isotope Reactor located at Oak Ridge National Laboratory. These methods are demonstrated using measurements on a variety of textured samples, including multi-phase steel composites and polycrystalline calcite (CaCO 3 ). Reference measurements were also made at VULCAN, the engineering diffractometer located at the Spallation Neutron Source. The texture data obtained on the different instruments are in agreement, and WAND 2 is more time efficient than HIDRA. Two analysis methods were investigated, single-peak fitting to obtain individual pole figures for inversion and Rietveld texture analysis using MAUD. Here, the impact of the differences between the various textures obtained was evaluated through the calculation of diffraction elastic constants, which is one application of the texture data collected. Both instruments were found to provide texture data that are suitable for complementing other analyses, such as residual stress mapping.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗
DOE OSTI2024

Multiwavelength investigations of PKS 2300–18: S-shaped radio quasar with precessing jets and double-peaked broad emission-line spectrum

ABSTRACT S-shaped radio galaxy jets are prime sources for investigating the dynamic interplay between the central active galactic nucleus (AGN), the jets, and the ambient intergalactic medium. These sources are excellent candidates for studying jet precession, as their S-shaped inversion symmetry strongly indicates underlying precession. We present a multiwavelength analysis of the giant inversion-symmetric S-shaped radio galaxy PKS 2300$-$18, which spans 0.76 Mpc. The host is a quasar at a redshift of 0.128, displaying disturbed optical morphology due to an ongoing merger with a companion galaxy. We conducted a broad-band radio spectral study using multifrequency data ranging from 183 MHz to 6 GHz, incorporating dedicated observations with the upgraded Giant Metrewave Radio Telescope (uGMRT) and Karl G. Jansky Very Large Array (JVLA) alongside archival radio data. A particle injection model was fitted to the spectra of different regions of the source to perform ageing analysis, which was supplemented with a kinematic jet precession model. The ageing analysis revealed a maximum plasma age of $\sim$ 40 Myr, while the jet precession model indicated a precession period of $\sim$ 12 Myr. ROentgen SATellite (ROSAT) data revealed an X-ray halo of Mpc size, and from Chandra the AGN X-ray spectrum was modelled using thermal and power-law components. The optical spectrum displaying double-peaked broad emission lines was modelled, indicating complex broad-line region kinematics at the centre with the possibility of a binary SMBH. We present the results of our multiwavelength analysis of the source, spanning scales from a few light-days to a few Mpc, and discuss its potential evolutionary path.

Misra, Arpita↗
DOE OSTI2021

Towards robust autonomous impedance spectroscopy analysis: A calibrated hierarchical Bayesian approach for electrochemical impedance spectroscopy (EIS) inversion

Distribution-based analyses, such as the distribution of relaxation times (DRT) and the distribution of diffusion times (DDT), present model-free alternatives to equivalent circuit modeling for analysis of electrochemical impedance spectroscopy (EIS) data. However, reconstructing such distributions from noisy impedance data is an ill-posed problem that must be solved with specialized inversion algorithms, requiring careful control and tuning. Furthermore, most inversion algorithms developed to date can only solve problems of limited complexity. Herein, we present a new hierarchical Bayesian method for EIS inversion, leveraging efficient algorithms for optimization and Hamiltonian Monte Carlo (HMC) sampling to solve models of arbitrary complexity. We overcome the challenge of ad-hoc parameter tuning by encoding intrinsic characteristics of the DRT and DDT into flexible prior distributions and “pre-calibrating” the model to simulated data. This approach is versatile, highly robust to noise, and provides quantitative estimates of both the error structure of the data and the uncertainty in the recovered distributions. The model is validated with simulated data to demonstrate accurate recovery of the DRT and the DDT. The method also shows promise for simultaneous recovery of multiple distributions, raising the intriguing possibility of semi-autonomous EIS analysis and ad-hoc model construction. Finally, the practical utility of the method is illustrated with experimental data. Throughout, we draw comparisons to several recently published EIS inversion methodologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Simulations and analysis tools for charge-exchange (d, 2 He) reactions in inverse kinematics with the AT-TPC

Charge-exchange (d, 2 He) reactions in inverse kinematics at intermediate energies are a very promising method to investigate the Gamow–Teller transition strength in unstable nuclei. A simulation and analysis software based on the attpcroot package was developed to study this type of reactions with the active-target time projection chamber (AT-TPC). The simulation routines provide a realistic detector response that can be used to understand and benchmark experimental data. Analysis tools and correction routines can be developed and tested from simulations in ATTPCROOT , because they are processed in the same way as the real data. In particular, we study the feasibility of using coincidences with beam-like particles to unambiguously identify the (d, 2 He) reaction channel, and to develop a kinematic fitting routine for future applications. More technically, the impact of space-charge effects in the track reconstruction, and a possible correction method are investigated in detail. Finally, this analysis and simulation package constitutes an essential part of the software development for the fast-beams program with the AT-TPC.

(d,2He)↗
DOE OSTI2023

Emulating the Deutsch-Josza algorithm with an inverse-designed terahertz gradient-index lens

An all-dielectric photonic metastructure is investigated for application as a quantum algorithm emulator (QAE) in the terahertz frequency regime; specifically, we show implementation of the Deustsh-Josza algorithm. The design for the QAE consists of a gradient-index (GRIN) lens as the Fourier transform subblock and patterned silicon as the oracle subblock. First, we detail optimization of the GRIN lens through numerical analysis. Then, we employed inverse design through a machine learning approach to further optimize the structural geometry. Through this optimization, we enhance the interaction of the incident light with the metamaterial via spectral improvements of the outgoing wave.

Blackwell, Ashley N.↗