Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid Adaptation”

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 1,513 records · Page 84

Subsurface Stress Fields in FCC Single Crystal Anisotropic Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent high cycle fatigue (HCF) failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and non-crystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is presented for evaluating the subsurface stresses in the elastic half-space, based on the adaptation of a stress function method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis (FEA). Effects of crystal orientation on stress response and fatigue life are examined. Obtaining accurate subsurface stress results for anisotropic single crystal contact problems require extremely refined three-dimensional (3-D) finite element grids, especially in the edge of contact region. Obtaining resolved shear stresses (RSS) on the principal slip planes also involves considerable post-processing work. For these reasons it is very advantageous to develop analytical solution schemes for subsurface stresses, whenever possible.

Arakere, Nagaraj K.↗

QuakeSim and the Solid Earth Research Virtual Observatory

We are developing simulation and analysis tools in order to develop a solid Earth science framework for understanding and studying active tectonic and earthquake processes. The goal of QuakeSim and its extension, the Solid Earth Research Virtual Observatory (SERVO), is to study the physics of earthquakes using state-of-the-art modeling, data manipulation, and pattern recognition technologies. We are developing clearly defined accessible data formats and code protocols as inputs to simulations, which are adapted to high-performance computers. The solid Earth system is extremely complex and nonlinear resulting in computationally intensive problems with millions of unknowns. With these tools it will be possible to construct the more complex models and simulations necessary to develop hazard assessment systems critical for reducing future losses from major earthquakes. We are using Web (Grid) service technology to demonstrate the assimilation of multiple distributed data sources (a typical data grid problem) into a major parallel high-performance computing earthquake forecasting code. Such a linkage of Geoinformatics with Geocomplexity demonstrates the value of the Solid Earth Research Virtual Observatory (SERVO) Grid concept, and advances Grid technology by building the first real-time large-scale data assimilation grid.

virtual observatory↗

An Automatic Cloud Mask Algorithm Based on Time Series of MODIS Measurements

Quality of aerosol retrievals and atmospheric correction depends strongly on accuracy of the cloud mask (CM) algorithm. The heritage CM algorithms developed for AVHRR and MODIS use the latest sensor measurements of spectral reflectance and brightness temperature and perform processing at the pixel level. The algorithms are threshold-based and empirically tuned. They don't explicitly address the classical problem of cloud search, wherein the baseline clear-skies scene is defined for comparison. Here, we report on a new CM algorithm which explicitly builds and maintains a reference clear-skies image of the surface (refcm) using a time series of MODIS measurements. The new algorithm, developed as part of the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm for MODIS, relies on fact that clear-skies images of the same surface area have a common textural pattern, defined by the surface topography, boundaries of rivers and lakes, distribution of soils and vegetation etc. This pattern changes slowly given the daily rate of global Earth observations, whereas clouds introduce high-frequency random disturbances. Under clear skies, consecutive gridded images of the same surface area have a high covariance, whereas in presence of clouds covariance is usually low. This idea is central to initialization of refcm which is used to derive cloud mask in combination with spectral and brightness temperature tests. The refcm is continuously updated with the latest clear-skies MODIS measurements, thus adapting to seasonal and rapid surface changes. The algorithm is enhanced by an internal dynamic land-water-snow classification coupled with a surface change mask. An initial comparison shows that the new algorithm offers the potential to perform better than the MODIS MOD35 cloud mask in situations where the land surface is changing rapidly, and over Earth regions covered by snow and ice.

Lyapustin, Alexei↗

Remote Off-grid Microgrid Design Support Tool (ROMDST) - An Optimal Design Support Tool for Remote, Resilient, and Reliable Microgrids Phase II - Final Report

The objective of the project is to develop an advanced optimization-based design support tool for AC or DC microgrids in remote locations, where utility grids may not be accessible. The mathematical model and the interface are being developed such that multiple design objectives and criteria/constraints can be easily enabled or disabled, to deliver a flexible tool relevant to a large, diverse user base, and to facilitate future feature developments. The new tool delivered by this project leverages the team’s extensive expertise in the development, testing, deployment, and commercialization of the state-of-the-art Distributed Energy Resources Customer Adoption Model (DER-CAM) – which is the foundation for the new tool. Phase I of the project leverages the microgrid design and decision support tool DER-CAM. Key Phase I developments include the following adaption to the DER-CAM model: a) AC or DC microgrid architectures, b) multiple economic objectives (e.g. single customer vs community vs utility perspective), c) constraints on fuel availability (e.g. total fuel available during outage periods), d) active and reactive power flow constraints for normal and N-1 contingency cases, e) component part load efficiencies (e.g. power electronic devices), and f) interactive data visualization capabilities (e.g. design topology elements). A component library including components’ default technical and economical characteristics was also developed. In Phase II, the project is focused on tool testing and validation, use for real-world microgrid designs, and outreach. Fully developing this tool will have the following impacts: 1) optimal off-grid microgrid designs that replace existing back of the envelope or non-optimal calculations, reducing capital costs and risk of microgrid deployment; 2) removing barriers to microgrid assessments by lowering microgrid soft costs, as the tool is freely usable; 3) reliable and resilient designs that reduce the cost of critical load shedding due to component outages.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Remote Off-grid Microgrid Design Support Tool (ROMDST) - An Optimal Design Support Tool for Remote, Resilient, and Reliable Microgrids (Phase II, Final Report)

The objective of the project is to develop an advanced optimization-based design support tool for AC or DC microgrids in remote locations, where utility grids may not be accessible. The mathematical model and the interface are being developed such that multiple design objectives and criteria/constraints can be easily enabled or disabled, to deliver a flexible tool relevant to a large, diverse user base, and to facilitate future feature developments. The new tool delivered by this project leverages the team’s extensive expertise in the development, testing, deployment, and commercialization of the state-of-the-art Distributed Energy Resources Customer Adoption Model (DER-CAM) – which is the foundation for the new tool. Phase I of the project leverages the microgrid design and decision support tool DER-CAM. Key Phase I developments include the following adaption to the DER-CAM model: a) AC or DC microgrid architectures, b) multiple economic objectives (e.g. single customer vs community vs utility perspective), c) constraints on fuel availability (e.g. total fuel available during outage periods), d) active and reactive power flow constraints for normal and N-1 contingency cases, e) component part load efficiencies (e.g. power electronic devices), and f) interactive data visualization capabilities (e.g. design topology elements). A component library including components’ default technical and economical characteristics was also developed. In Phase II, the project is focused on tool testing and validation, use for real-world microgrid designs, and outreach. Fully developing this tool will have the following impacts: 1) optimal off-grid microgrid designs that replace existing back of the envelope or non-optimal calculations, reducing capital costs and risk of microgrid deployment; 2) removing barriers to microgrid assessments by lowering microgrid soft costs, as the tool is freely usable; 3) reliable and resilient designs that reduce the cost of critical load shedding due to component outages.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DownScaleBench for developing and applying a deep learning based urban climate downscaling- first results for high-resolution urban precipitation climatology over Austin, Texas

Abstract Cities need climate information to develop resilient infrastructure and for adaptation decisions. The information desired is at the order of magnitudes finer scales relative to what is typically available from climate analysis and future projections. Urban downscaling refers to developing such climate information at the city (order of 1 – 10 km) and neighborhood (order of 0.1 – 1 km) resolutions from coarser climate products. Developing these higher resolution (finer grid spacing) data needed for assessments typically covering multiyear climatology of past data and future projections is complex and computationally expensive for traditional physics-based dynamical models. In this study, we develop and adopt a novel approach for urban downscaling by generating a general-purpose operator using deep learning. This ‘DownScaleBench’ tool can aid the process of downscaling to any location. The DownScaleBench has been generalized for both in situ (ground- based) and satellite or reanalysis gridded data. The algorithm employs an iterative super-resolution convolutional neural network (Iterative SRCNN) over the city. We apply this for the development of a high-resolution gridded precipitation product (300 m) from a relatively coarse (10 km) satellite-based product (JAXA GsMAP). The high-resolution gridded precipitation datasets is compared against insitu observations for past heavy rain events over Austin, Texas, and shows marked improvement relative to the coarser datasets relative to cubic interpolation as a baseline. The creation of this Downscaling Bench has implications for generating high-resolution gridded urban meteorological datasets and aiding the planning process for climate-ready cities.

Singh, Manmeet (ORCID:0000000233747149)↗

Reducing Demand Charges and Onsite Generation Variability Using Behind-the-Meter Energy Storage

Electric utilities in the United States are increasingly employing demand charges and/or real-time pricing. This directive is bringing potential opportunities in deploying behindthe-meter energy storage (BMES) systems for various grid functionalities. This study quantifies techno-economic benefits of BMES in reducing demand charge and smoothing load/generation intermittencies, and determines how those benefits vary with onsite distributed photovoltaic. We proposed a two-stage control algorithm, whereby the first stage proactively determines costoptimal BMES configuration for reducing peak-demands and demand charges, and the second stage adaptively compensates intermittent generations and short load spikes that may otherwise increase the demand charges. The performance of the proposed algorithm is evaluated through a 24 hours time sweep simulation performed using data from smart microgrid testbed at Idaho National Laboratory (INL). The simulation results demonstrated that this research provides a simple but effective solution for peak shaving, demand charge reductions, and smoothing onsite PV variability.

25 ENERGY STORAGE↗

An adaptive multigrid model for hurricane track prediction

This paper describes a simple numerical model for hurricane track prediction which uses a multigrid method to adapt the model resolution as the vortex moves. The model is based on the modified barotropic vorticity equation, discretized in space by conservative finite differences and in time by a Runge-Kutta scheme. A multigrid method is used to solve an elliptic problem for the streamfunction at each time step. Nonuniform resolution is obtained by superimposing uniform grids of different spatial extent; these grids move with the vortex as it moves. Preliminary numerical results indicate that the local mesh refinement allows accurate prediction of the hurricane track with substantially less computer time than required on a single uniform grid.

Fulton, Scott R.↗

High-spatial-resolution passive microwave sounding systems

The principal contributions of this combined theoretical and experimental effort were to advance and demonstrate new and more accurate techniques for sounding atmospheric temperature, humidity, and precipitation profiles at millimeter wavelengths, and to improve the scientific basis for such soundings. Some of these techniques are being incorporated in both research and operational systems. Specific results include: (1) development of the MIT Microwave Temperature Sounder (MTS), a 118-GHz eight-channel imaging spectrometer plus a switched-frequency spectrometer near 53 GHz, for use on the NASA ER-2 high-altitude aircraft, (2) conduct of ER-2 MTS missions in multiple seasons and locations in combination with other instruments, mapping with unprecedented approximately 2-km lateral resolution atmospheric temperature and precipitation profiles, atmospheric transmittances (at both zenith and nadir), frontal systems, and hurricanes, (3) ground based 118-GHz 3-D spectral images of wavelike structure within clouds passing overhead, (4) development and analysis of approaches to ground- and space-based 5-mm wavelength sounding of the upper stratosphere and mesosphere, which supported the planning of improvements to operational weather satellites, (5) development of improved multidimensional and adaptive retrieval methods for atmospheric temperature and humidity profiles, (6) development of combined nonlinear and statistical retrieval techniques for 183-GHz humidity profile retrievals, (7) development of nonlinear statistical retrieval techniques for precipitation cell-top altitudes, and (8) numerical analyses of the impact of remote sensing data on the accuracy of numerical weather predictions; a 68-km gridded model was used to study the spectral properties of error growth.

Staelin, D. H.↗

WRF-UCM CDD (cooling degree days) data

The MATLAB (.mat) files (spatial grids of daily CDD data) are generated by a suit of WRF-UCM simulations over California at 4.5 km resolution (1.5 km resolution over cities) for 5 scenarios

54 ENVIRONMENTAL SCIENCES↗

Coordination of Wind Farm and Pumped-Storage Hydro for a Self-Healing Power Grid

The increasing penetration of wind energy poses great challenges to the operation of power systems in normal and emergency states. However, energy storage technologies can help accommodate wind power uncertainty and variability due to their flexible characteristics. This paper focuses on the restoration phase, and provides a novel coordination strategy of wind and pumped-storage hydro (PSH) units for a faster and reliable self-healing process. The wind-PSH assisted power system restoration is formulated as a two-stage adaptive robust optimization problem. The first-stage problem determines the start-up sequence of generators and the energization times of transmission paths; and the second-stage problem decides load pickup sequences, wind power dispatch levels, and PSH units' operating modes. The column-and-constraint generation decomposition algorithm is applied to solve the two-stage adaptive robust optimization problem, which has a mixed-integer optimization in the inner-level problem. The developed coordination strategy is tested on the modified IEEE 39-bus system. Numerical results demonstrate that the coordinated wind and PSH units can increase the total energy served, enhance wind power dispatchability, and reduce wind power curtailment.

17 WIND ENERGY↗

SU-C-202-05: Pilot Study of Online Treatment Evaluation and Adaptive Re-Planning for Laryngeal SBRT

Purpose: We have instigated a phase I trial of 5-fraction stereotactic body radiotherapy (SBRT) for advanced-stage laryngeal cancer. We conducted this pilot dosimetric study to confirm the potential utility of online adaptive re-planning to preserve treatment quality. Methods: Ten cases of larynx cancer were evaluated. Baseline and daily SBRT treatment plans were generated per trial protocol. Daily volumetric images were acquired prior to every fraction of treatment. Reference simulation CT images were deformably registered to daily volumetric images using Eclipse. Planning contours were then deformably propagated to daily images. Reference SBRT plans were directly copied to calculate delivered dose distributions on deformed reference CT images. In-house software platform has been developed to calculate cumulative dose over a course of treatment in four steps: 1) deforming delivered dose grid to reference CT images using deformation information exported from Eclipse; 2) generating tetrahedrons using deformed dose grid as vertices; 3) resampling dose to a high resolution within every tetrahedron; 4) calculating dose-volume histograms. Our inhouse software was benchmarked with a commercial software, Mirada. Results: In all ten cases including 49 fractions of treatments, delivered daily doses were completely evaluated and treatment could be re-planned within 10 minutes. Prescription dose coverage of PTV was less than intended in 53% of fractions of treatment (mean: 94%, range: 84%–98%) while minimum coverage of CTV and GTV was 94% and 97%, respectively. Maximum bystander point dose limits to arytenoids, parotids, and spinal cord remained respected in all cases, although variances in carotid artery doses were observed in a minority of cases. Conclusion: Although GTV and CTV coverage is preserved by in-room 3D image guidance of larynx SBRT, PTV coverage can vary significantly from intended plans. Online adaptive treatment evaluation and re-planning is potentially necessary and our procedure is clinically applicable to fully preserve treatment quality. This project is supported by CPRIT Individual Investigator Research Award RP150386.

60 APPLIED LIFE SCIENCES↗

Technoeconomic Design Optimization for Fast Reactors. Part II: Impact of Technoeconomic Constraints on Optimal Design

There is a current drive toward optimizing reactors, particularly small/micro reactors to minimize cost and maximize performance. Previous work has investigated the development of technoeconomic workflows for the design optimization of pool-type fast reactors that aim to deploy into district energy grids. Initial scoping studies verified that the workflow was capable of capturing design trends throughout a variety of design configurations and problem formulations while remaining sufficiently flexible. In this paper, this methodology is applied to understand how cost functions and technoeconomic constraints can drive optimal reactor design. Specifically, the UPu10Zr-fueled fast reactor model from Part I is adapted to include changes in the fissile content limits, control rod worth limits, control rod drive cost, and assumed fuel form. In the case of constraint relaxation at fixed power (fissile content and control rod worth limits), cost sensitivities of 5% to 10% were uncovered. Multi-objective optimization at varying reactor power levels with individualized control rod drives for each assembly (as opposed to one operational and one safety drive) increased cost by $\$10$ to $\$25$ million and substantially altered the optimal core geometry, favoring geometries with substantially fewer control rod placements relative to baseline optimization. Finally, a multi-objective optimization was performed at varying power levels with the fuel form overhauled to metallic, high-assay low-enriched uranium–based U10Zr with more refined fuel cost models. In the case of uranium fueling, the costs increased by at least $50 million relative to the baseline case. Furthermore, economic fuel zoning and lower reactivity swing cores were recovered. Each case serves to demonstrate the value of applying technoeconomic workflows to initial reactor design scoping studies to better understand the trade-off for a proposed concept between different design options.

Argonne Reactor Computation (ARC) codes↗

Near-surface ice on Mercury and the Moon: A topographic thermal model

A thermal model that can be easily adapted to craters of arbitrary shape is developed and applied to high-latitude impact craters on Mercury and the Moon, Chao Meng Fu crater at -87.5 deg L on Mercury, an unnamed bowl-shaped crater at 86.7 deg L on Mercury, and Peary crater at 88.6 deg L on the Moon. For an assumed input topography and grid of surface elements, the model computes for each element the irradiation from direct insolation and reflected and emitted radiation from other elements, taking into account shadowing by walls of the crater, partial obscuration of the solar disk near the poles and the diurnal, orbital, and seasonal cycles. Temperatures are computed over the surface grid as functions of depth and time from the surface to a specified depth and over the pertinent astronomical cycles, including the effects of direct and indirect surface irradiation, infrared radiation, heat conduction, and interior heating. Vapor fluxes and ice recession times are computed as functions of ice depth over the surface grid. Temperatures profiles, vapor fluxes, and ice recession times were computed for flat surfaces not associated with craters near the poles of Mercury and the Moon. It was found that water ice could have existed throughout geologic time within the maximum radar detection depth of recent observation of Mercury (J. K. Harmon and M. A. Slade, 1992, Science 258, 640-643) poleward of approximately 87 - 88 deg L on Mercury and poleward of approximately 73 deg L on the Moon. For Chao Meng Fu crater it was found that approximately 40% of the crater floor is permanently shadowed from direct solar insolation, while the remainder of the crater floor is periodically illuminated by a partially obscured Sun. Temperatures at the upper levels of the south wall can slightly exceed 550 K. Surface temperatures in the permanently shadowed region of the crater floor are under approximately 130 K, which could have allowed water ice to exist throughout geologic time within the radar detection depth of recent observation of Mercury. For small bowl-shaped crater on Mercury, it was found that most of the crater is permanently shadowed from direct solar radiation, except for a narrow semicircular band bordering the north rim. However, temperatures in the permanently shadowed region periodically reach a maximum near approximately 315 K due to efficient heating of the small crater by thermal emission and reflection from the small sunlit region, which periodically reaches temperatures exceeding 630 K. Water ice could not have existed throughout geologic time anywhere in this crater within the radar detection depth. For Peary crater on the Moon, the entire crater floor is permanently shadowed from direct solar insolation with maximum temperature under 120 K. The upper level of the north wall periodically reaches a maximum temperature near 310 K. The low temperatures on the crater floor would have allowed water ice to exist near the surface throughout geologic time, provided that the Moon's obliquity was always as low as it is at present.

Salvail, James R.↗

Annual Research Briefs - 1996, Center for Turbulence Research

Partial contents include: (1) Simulation of jet impingement heat transfer with the k-e-V2 model; (2) A new methodology for turbulence modelers using Direct Numerical Simulation (DNS) database analysis; (3) A particle representation model for the deformation of homogeneous turbulence; (4) Modeling of inhomogeneous compressible turbulence using a two-scale statistical theory; (5) Lewis number effects on partially premixed flames; (6) Conditional moment closure of mixing and reaction in turbulent non-premixed combustion; (7) Dynamic models for LES of turbulent front propagation with a spectral method; (8) Jet noise using large eddy simulation; (9) Towards numerical simulations of trailing-edge aeroacoustics; (10) A model for the onset of vortex breakdown; (11) Compressing an elliptic vortex: transition to turbulence by tumble breakdown; (12) A new method for the adaptive control of vortex-wall interactions; (13) Symmetries in turbulent boundary layer flows; (14) Near-wall models in large eddy simulations of flow behind a backward-facing step; (15) Assessment of the subgrid-scale models at low and high Reynolds numbers; (16) Large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids; (17) Progress on LES of flow past a circular cylinder; (18) Experimental investigation of flow through an asymmetric plane diffuser; (19) Progress in the large-eddy simulation of an asymmetric plane diffuser; (20) A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation; (21) On the generation of turbulent inflow conditions for boundary layer simulations, and (22) A new class of finite difference schemes.

COMPUTATIONAL FLUID DYNAMICS↗

Evaluation of a CFD Method for Aerodynamic Database Development using the Hyper-X Stack Configuration

A computational fluid dynamic (CFD) study is performed on the Hyper-X (X-43A) Launch Vehicle stack configuration in support of the aerodynamic database generation in the transonic to hypersonic flow regime. The main aim of the study is the evaluation of a CFD method that can be used to support aerodynamic database development for similar future configurations. The CFD method uses the NASA Langley Research Center developed TetrUSS software, which is based on tetrahedral, unstructured grids. The Navier-Stokes computational method is first evaluated against a set of wind tunnel test data to gain confidence in the code s application to hypersonic Mach number flows. The evaluation includes comparison of the longitudinal stability derivatives on the complete stack configuration (which includes the X-43A/Hyper-X Research Vehicle, the launch vehicle and an adapter connecting the two), detailed surface pressure distributions at selected locations on the stack body and component (rudder, elevons) forces and moments. The CFD method is further used to predict the stack aerodynamic performance at flow conditions where no experimental data is available as well as for component loads for mechanical design and aero-elastic analyses. An excellent match between the computed and the test data over a range of flow conditions provides a computational tool that may be used for future similar hypersonic configurations with confidence.

Parikh, Paresh↗

Transparent Electrodes for Efficient Optoelectronics

With the development of new generations of optoelectronic devices that combine high performance and novel functionalities (e.g., flexibility/bendability, adaptability, semi or full transparency), several classes of transparent electrodes have been developed in recent years. These range from optimized transparent conductive oxides (TCOs), which are historically the most commonly used transparent electrodes, to new electrodes made from nano‐ and 2D materials (e.g., metal nanowire networks and graphene), and to hybrid electrodes that integrate TCOs or dielectrics with nanowires, metal grids, or ultrathin metal films. Here, the most relevant transparent electrodes developed to date are introduced, their fundamental properties are described, and their materials are classified according to specific application requirements in high efficiency solar cells and flexible organic light‐emitting diodes (OLEDs). This information serves as a guideline for selecting and developing appropriate transparent electrodes according to intended application requirements and functionality.

Morales‐Masis, Monica↗

Space‐Time Causal Discovery in Earth System Science: A Local Stencil Learning Approach

Causal discovery tools enable scientists to infer meaningful relationships from observational data, spurring advances in fields as diverse as biology, economics, and climate science. Despite these successes, the application of causal discovery to space-time systems remains immensely challenging due to the high-dimensional nature of the data. For example, in climate sciences, modern observational temperature records over the past few decades regularly measure thousands of locations around the globe. To address these challenges, we introduce Causal Space-Time Stencil Learning (CaStLe), a novel meta-algorithm for discovering causal structures in complex space-time systems. CaStLe leverages regularities in local space-time dependencies to learn governing global dynamics. This local perspective eliminates spurious confounding and drastically reduces sample complexity, making space-time causal discovery practical and effective. For causal discovery, CaStLe flexibly accepts any appropriately adapted time series causal discovery algorithm to recover local causal structures. These advances enable causal discovery of geophysical phenomena that were previously unapproachable, including non-periodic, transient phenomena such as volcanic eruption plumes. Regularities in local space-time dependencies are transformed into informative spatial replicates, which actually improve CaStLe's performance when applied to ever-larger spatial grids. We successfully apply CaStLe to discover the atmospheric dynamics governing the climate response to the 1991 Mount Pinatubo volcanic eruption. We provide validation experiments to demonstrate the effectiveness of CaStLe over existing causal-discovery frameworks on a range of geophysics-inspired benchmarks while identifying the method's limitations and domains where its assumptions may not hold.

Nichol, J. Jake [Univ. of New Mexico, Albuquerque,↗