Search NASA⌕ Search

SEARCH · Search NASA

Results for “GRID MODERNIZATION”

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

Reinforcement Learning-based Output Structured Feedback for Distributed Multi-Area Power System Frequency Control

Load frequency control (LFC) is a key factor to maintain the stable frequency in multi-area power systems. As the modern power systems evolve from centralized to decentralized paradigm, LFC needs to consider the decentralized scheme that considers limited information from the information-exchange graph for the generator control of each interconnected area. This paper aims to solve a data-driven constrained LQR problem with mean-variance risk constraints and output structured feedback, and applies this framework to solve the LFC problem in multi-area power systems. By reformulating the constrained optimization problem into a minimax problem, the stochastic gradient descent max-oracle (SGDmax) algorithm with zero-order policy gradient (ZOPG) is adopted to find the optimal feedback gain from the learning, while guaranteeing the convergence. In addition, to improve the adaptation of the proposed learning method to new or varying models, we construct an emulator grid that approximates the dynamics of a physical grid and performs training based on this model. Once the feedback gain is obtained from the emulator grid, it is applied to the physical grid with a robustness test to check whether the controller from the approximated emulator applies to the actual system. Numerical tests show that the obtained feedback controller can successfully control the frequency of each area, while mitigating the uncertainty from the loads, with reliable robustness that ensures the adaptability of the obtained feedback gain to the actual physical grid.

Kwon, Kyung-bin↗
DOE OSTI2023

PNNL DSSE Docker Image

This is the docker image for Pacific Northwest National Laboratory's DSSE (distribution system state estimator) used for the demo of OEDI-SI platform. To support the operation of modern distribution systems, operators require real-time visibility into system states. Due to a lack of measurements and unbalanced operation, the state estimation in distribution systems is challenging. This submission is related to an OEDI-SI use case which demonstrates the application of an extended Kalman filter-based state estimator on an IEEE-123 bus system. The state estimator uses the measurements to generate voltage estimates for the system.

Array↗
NASA NTRS2006

Computational Analysis of a Chevron Nozzle Uniquely Tailored for Propulsion Airframe Aeroacoustics

A computational flow field and predicted jet noise source analysis is presented for asymmetrical fan chevrons on a modern separate flow nozzle at take off conditions. The propulsion airframe aeroacoustic asymmetric fan nozzle is designed with an azimuthally varying chevron pattern with longer chevrons close to the pylon. A baseline round nozzle without chevrons and a reference nozzle with azimuthally uniform chevrons are also studied. The intent of the asymmetric fan chevron nozzle was to improve the noise reduction potential by creating a favorable propulsion airframe aeroacoustic interaction effect between the pylon and chevron nozzle. This favorable interaction and improved noise reduction was observed in model scale tests and flight test data and has been reported in other studies. The goal of this study was to identify the fundamental flow and noise source mechanisms. The flow simulation uses the asymptotically steady, compressible Reynolds averaged Navier-Stokes equations on a structured grid. Flow computations are performed using the parallel, multi-block, structured grid code PAB3D. Local noise sources were mapped and integrated computationally using the Jet3D code based upon the Lighthill Acoustic Analogy with anisotropic Reynolds stress modeling. In this study, trends of noise reduction were correctly predicted. Jet3D was also utilized to produce noise source maps that were then correlated to local flow features. The flow studies show that asymmetry of the longer fan chevrons near the pylon work to reduce the strength of the secondary flow induced by the pylon itself, such that the asymmetric merging of the fan and core shear layers is significantly delayed. The effect is to reduce the peak turbulence kinetic energy and shift it downstream, reducing overall noise production. This combined flow and noise prediction approach has yielded considerable understanding of the physics of a fan chevron nozzle designed to include propulsion airframe aeroacoustic interaction effects.

Massey, Steven J.↗
DOE OSTI2022

A low-cost centralized HVAC control system solution for energy savings, load shedding, and improved maintenance

University campuses rely on centralized controls for managing and optimizing complex HVAC systems in larger buildings. However, most campuses also have many smaller buildings with packaged HVAC systems controlled by a stand-alone thermostat. Even when these distributed and often overlooked systems have modern programmable thermostats, they cannot be centrally monitored or controlled, and they are typically not programmed adequately. This paper describes the implementation of a low-cost centralized control solution for these systems serving smaller campus buildings, mostly under 5,000 sf and representative of light commercial spaces. Thanks to advances in technology spurred by residential and commercial IoT developments, simple networked thermostat solutions exist that can easily replace original thermostats, and, connect these systems to a web-based portal for monitoring and control. We show that, with small customizations, these platforms can be integrated into facility management workflows. Beyond the energy savings potential from improved scheduling and closer management of these systems, there are significant advantages for maintenance crews since these systems can now be monitored on smart phones or tablets. A grid-responsive load-shedding program has also been implemented for additional cost savings. The networked thermostats can also be connected to additional systems such as economizer controls for improved ventilation management and energy savings. With data from these systems integrated centrally, it can also be used for improved analytics and fault detection. A toolkit has been developed to share the program with other campuses, whether for energy savings, improved management of ventilation, or a more proactive maintenance approach.

Fauchier-Magnan, Nicolas↗
NASA NTRS2008

NASA's Modern Era Retrospective-Analysis for Research and Applications (MERRA): Early Results and Future Directions

This talk will review the status and progress of the NASA/Global Modeling and Assimilation Office (GMAO) atmospheric global reanalysis project called the Modern Era Retrospective-Analysis for Research and Applications (MERRA). An overview of NASA's emerging capabilities for assimilating a variety of other Earth Science observations of the land, ocean, and atmospheric constituents will also be presented. MERRA supports NASA Earth science by synthesizing the current suite of research satellite observations in a climate data context (covering the period 1979-present), and by providing the science and applications communities with of a broad range of weather and climate data with an emphasis on improved estimates of the hydrological cycle. MERRA is based on a major new version of the Goddard Earth Observing System Data Assimilation System (GEOS-5), that includes the Earth System Modeling Framework (ESMF)-based GEOS-5 atmospheric general circulation model and the new NOAA National Centers for Environmental Prediction (NCEP) unified grid-point statistical interpolation (GST) analysis scheme developed as a collaborative effort between NCEP and the GMAO. In addition to MERRA, the GMAO is developing new capabilities in aerosol and constituent assimilation, ocean, ocean biology, and land surface assimilation. This includes the development of an assimilation capability for tropospheric air quality monitoring and prediction, the development of a carbon-cycle modeling and assimilation system, and an ocean data assimilation system for use in coupled short-term climate forecasting.

Schubert, Siegfried↗
NASA NTRS1990

Unsteady transonic Navier-Stokes computations for an oscillating wing using single and multiple zones

Modern jet transports and maneuvering tactical fighters operating in the transonic regime often give rise to time-dependent fluid physics that interact with flexible structural components, e.g., vortical flow, shocks, and separation. Efficient computational fluid dynamic (CFD) methods are required to study such computationally intensive problems. In this work a numerical method is presented to address this problem. Time-dependent, compressible, Navier-Stokes equations are used to simulate unsteady transonic flow about a three-dimensional rigid wing undergoing a forced periodic motion in angle of attack. An efficient, implicit, diagonal algorithm is utilized because of its low operation count per time step compared to other methods that solve systems of block matrix equations. The formal time accuracy is addressed theoretically and demonstrated numerically by comparison of computational results with experimental data. A zonal grid approach, capable of treating complex geometries, is presented and its time accuracy is demonstrated by comparing a zonal computation with a single grid computation and experimental data.

Chaderjian, Neal M.↗
NASA NTRS2019

Assimilation of Satellite-Based Snow Cover and Freeze/Thaw Observations Over the High Mountain Asia

Toward qualifying hydrologic changes in the High Mountain Asia (HMA) region, this study explores the use of a hyper-resolution (1 km) land data assimilation (DA)framework developed within the NASA Land Information System using the Noah Multi-parameterization Land Surface Model (Noah-MP) forced by the meteorological boundary conditions from Modern-Era Retrospective analysis for Research and Applications, Version 2 data. Two different sets of DA experiments are conducted:(1) the assimilation of a satellite-derived snow cover map (MOD10A1) and (2) the assimilation of the NASA MEaSUREs landscape freeze/thaw product from 2007 to 2008. The performance of the snow cover assimilation is evaluated via comparisons with available remote sensing-based snow water equivalent product and ground-based snow depth measurements. For example, in the comparison against ground-based snow depth measurements, the majority of the stations (13 of 14) show slightly improved goodness-of-fit statistics as a result of the snow DA, but only four are statistically significant. In addition, comparisons to the satellite-based land surface temperature products (MOD11A1 and MYD11A1) show that freeze/thaw DA yields improvements (at certain grid cells) of up to 0.58 K in the root-mean-square error (RMSE) and 0.77K in the absolute bias (relative to model-only simulations). In the comparison against three ground-based soil temperature measurements along the Himalayas, the bias and the RMSE in the 0-10 cm soil temperature are reduced (on average) by 10 and 7%,respectively. The improvements in the top layer of soil estimates also propagate through the deeper soil layers, where the bias and the RMSE in the 10-40 cm soil temperature are reduced (on average) by 9 and 6%, respectively. However, no statistically significant skill differences are observed for the freeze/thaw DA system in the comparisons against ground-based surface temperature measurements at mid-to-low altitude. Therefore, the two proposed DA schemes show the potential of improving the predictability of snow mass, surface temperature, and soil temperature states across HMA, but more ground-based measurements are still required, especially at high-altitudes, in order to document a more statistically significant improvement as a result of the two DA schemes.

High Mountain Asia↗
DOE OSTI2019

Laboratory-Based X-ray Absorption Spectroscopy on a Working Pouch Cell Battery at Industrially-Relevant Charging Rates

Li-ion battery (LIB) research has continuing importance for the entire range of applications from consumer products to vehicle electrification and grid stabilization. In many cases, standard electrochemical methods only provide an overall voltage or specific capacity, giving an inadequate description of parallel redox processes or chemical gradients at the particle and pack level. X-ray absorption fine structure (XAFS) is frequently used to augment bulk electrochemical data, as it provides element-specific changes in oxidation state and local atomic structure. Such microscopic descriptors are crucial for elucidating charge transfer and structural changes associated with bonding or site mixing, two key factors in evaluating state of charge and modes of cell failure. However, the impact of XAFS on LIB research has been significantly constrained by a logistical barrier: contemporary XAFS work is performed almost exclusively at synchrotron X-ray light sources, where beamtime is infrequent and experiment time-frames are limited. In this report we show that modern laboratory-based XAFS cannot only be applied to, e.g., characterization of ex situ LIB electrode materials, but can also be used for operando studies at industrially-relevant charging rates in a standard pouch cell preparation. Such capability enables accelerated discovery of new materials and improved operation modes for LIBs.

25 ENERGY STORAGE↗
DOE OSTI2026

The 2025 Evaluation of Experimental Thermonuclear Reaction Rates (ETR25)

This work describes the formalism for estimating thermonuclear reaction rates for astrophysical applications, emphasizing modern statistical approaches such as Monte Carlo sampling and Bayesian models. We discuss related topics including the calculation of resonance energies from nuclear Q values, indirect estimates of particle partial widths, and matching of reaction rates at elevated temperatures to statistical model results. We have evaluated available experimental data on cross sections, resonance energies and strengths, partial widths, lifetimes, spin-parities, and spectroscopic factors. Based on these results, we have estimated numerical values of 78 experimental charged-particle thermonuclear reaction rates for target nuclei in the A = 2–40 mass region, for temperatures ranging from 1 MK to 10 GK. For each reaction, three rate values are provided: low, median, and high, corresponding to the 16th, 50th, and 84th percentiles, respectively, of the cumulative reaction rate probability density distribution. Additionally, we present the factor uncertainty of each rate at each temperature grid point. These results enable users to sample the reaction rate probability density in nucleosynthesis calculations, facilitating uncertainty estimates of nuclidic abundances. The rates presented here refer to their laboratory values. For use in stellar model simulations, these values need to be corrected for the effects of thermal excitations of the interacting nuclei. For each reaction, we include graphs that illustrate the fractional contributions to the overall reaction rate along with the associated uncertainty. These visuals are designed to assist both stellar modelers and nuclear experimentalists by identifying the primary sources of rate uncertain=^texttx);ty at specific stellar temperatures. A graphical comparison with earlier Monte Carlo rates is also provided.

Nuclear astrophysics↗
NASA NTRS2014

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗
DOE OSTI2025

Hawaii Threat Brief

The Digital Energy Transformation is redefining how power systems operate, integrating physical infrastructure with digital technologies to enhance efficiency, visibility, and resilience. For islanded and digitally modernizing systems like Hawai‘i’s, this shift presents both new opportunities and evolving challenges in cybersecurity, supply chain assurance, and environmental resilience. This brief provides an overview of critical infrastructure dependencies and systemic risks associated with increasing digital integration. It summarizes recent energy-sector threat activity and case studies that illustrate adversary tactics and supply chain vulnerabilities, and identifies pathways to strengthen resilience.

25 - ENERGY STORAGE↗
NASA NTRS2016

Extremely Severe Space Weather and Geomagnetically Induced Currents in Regions with Locally Heterogeneous Ground Resistivity

Large geomagnetically induced currents (GICs) triggered by extreme space weather events are now regarded as one of the serious natural threats to the modern electrified society. The risk is described in detail in High-Impact, Low-Frequency Event Risk, A Jointly-Commissioned Summary Report of the North American Electric Reliability Corporation and the US Department of Energy's November 2009 Workshop, June 2010. For example, the March 13-14,1989 storm caused a large-scale blackout affecting about 6 million people in Quebec, Canada, and resulting in substantial economic losses in Canada and the USA (Bolduc 2002). Therefore, European and North American nations have invested in GIC research such as the Solar Shield project in the USA (Pulkkinen et al. 2009, 2015a). In 2015, the Japanese government (Ministry of Economy, Trade and Industry, METI) acknowledged the importance of GIC research in Japan. After reviewing the serious damages caused by the 2011 Tohoku-Oki earthquake, METI recognized the potential risk to the electric power grid posed by extreme space weather. During extreme events, GICs can be concerning even in mid- and low-latitude countries and have become a global issue.

Fujita, Shigeru↗
DOE OSTI2025

Spacetime pq theory for AC and DC electric power systems

The 50/60 Hz alternating current (AC) electric power has been the standard and most flexible energy source powering our modern societies for one and a half centuries since the war of the currents: AC versus direct current (DC). A reactive power concept that was introduced at the beginning of the AC power was very useful for circuit/system analysis, design, control, optimization, and ultimately for more efficient and stable generation, transmission, distribution, and consumption. The initial reactive power theory was based on single-phase sinusoidal AC power to capture inductive and capacitive power that yields to net-zero average power over one fundamental cycle. Soon it was expanded to non-sinusoidal AC power and finally to instantaneous three-phase AC power. However, these reactive power theories remain separate and limited to special cases and have never been consolidated and made valid to all cases. Today, more widespread adoption of power electronics and renewable energy is bringing back DC power into the electric grids. The reactive power concept has never been applied to DC power systems. There is no reactive power in DC power systems according to the existing reactive power theories. Do DC power systems really have no reactive power? Capacitors and inductors are widely used in DC just like in AC power systems. Are they not reactive power components? Why are they different from their AC counterparts? Furthermore, are batteries active or reactive power components? What about active devices like power converters (or inverters) with AC (or DC) on one side and DC (or AC) on the other? Do they generate or consume reactive power? Finally, what about AC and DC hybrid power systems? How to define reactive power in such a complex power system that has a multitude of loads, buses, and sources? Is there reactive power between any two loads, any two buses, or any two sources in a power system and what is the total reactive power in such a complex power system as a whole? As the motivation and goal of this paper to answer the above basic questions, to unify the existing AC reactive power theories and to ultimately provide theoretical and insightful guidance for system analysis, design, control, efficiency, optimization, and operation of complex power systems, a concept of spacetime (both spatial and temporal) active and reactive power (pq) theory—the spatiotemporal aspect of active and reactive power—is developed for both AC and DC power systems. The theoretical definitions and physical meanings of the spacetime reactive power will be developed, and real applications and thought experiments/cases/exercises will be explored and discussed. The developed mathematics to define the active (or real) and reactive (or imaginary) power— p and q respectively by dot (scalar) and cross (vector) products of multi-dimension spacetime vectors and time-space mapping principle/law can have some fundamental implications as well.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

SU-E-T-629: Feasibility Study of Treating Multiple Brain Tumors with Large Number of Noncoplanar IMRT Beams

Purpose: To study the feasibility of treating multiple brain tumors withlarge number of noncoplanar IMRT beams. Methods: Thirty beams are selected from 390 deliverable beams separated by six degree in 4pi space. Beam selection optimization is based on a column generation algorithm. MLC leaf size is 2 mm. Dose matrices are calculated with collapsed cone convolution and superposition method in a 2 mm by 2mm by 2 mm grid. Twelve brain tumors of various shapes, sizes and locations are used to generate four plans treating 3, 6, 9 and 12 tumors. The radiation dose was 20 Gy prescribed to the 100% isodose line. Dose Volume Histograms for tumor and brain were compared. Results: All results are based on a 2 mm by 2 mm by 2 mm CT grid. For 3, 6, 9 and 12 tumor plans, minimum tumor doses are all 20 Gy. Mean tumor dose are 20.0, 20.1, 20.1 and 20.1 Gy. Maximum tumor dose are 23.3, 23.6, 25.4 and 25.4 Gy. Mean ventricles dose are 0.7, 1.7, 2.4 and 3.1 Gy.Mean subventricular zone dose are 0.8, 1.3, 2.2 and 3.2 Gy. Average Equivalent uniform dose (gEUD) values for tumor are 20.1, 20.1, 20.2 and 20.2 Gy. The conformity index (CI) values are close to 1 for all 4 plans. The gradient index (GI) values are 2.50, 2.05, 2.09 and 2.19. Conclusion: Compared with published Gamma Knife treatment studies, noncoplanar IMRT treatment plan is superior in terms of dose conformity. Due to maximum limit of beams per plan, Gamma knife has to treat multiple tumors separately in different plans. Noncoplanar IMRT plans theoretically can be delivered in a single plan on any modern linac with an automated couch and image guidance. This warrants further study of using noncoplanar IMRT as a viable treatment solution for multiple brain tumors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Opacity Table Generation

Radiation opacity is a key property of material systems which describes the absorption and scattering of photons. Such processes are an important part of energy transport, diagnostics, and astrophysical phenomena. Recent experiments in iron opacity and continuum lowering have sparked renewed interest in and careful vetting of opacity models. Even within the local thermodynamic equilibrium approximation, opacity varies over a large parameter space: material, temperature, density, and energy. Opacity tables that span wide regimes of this parameter space are of prime importance to radiation-hydrodynamics codes. Accurate opacity calculations cannot be achieved inline during a radiation hydrodynamics simulation since they involve accounting for all possible electron transitions of all important atomic states of a thermal ensemble. Livermore’s current framework utilizes precalculated tables at many densities, temperatures, and energies. Last year, the opacity theory team completed an L2 milestone concerning the development of a new code, Opus, which served the dual purpose of creating a modern code infrastructure and version control as well as the basis for training new members of that team, a critical need as key members near or are passed retirement. For the current L2 milestone, we exercised and validated Opus on three elements: boron, carbon, and nitrogen, which have few enough electrons to allow a careful convergence study over many input parameters. Here, we report on the production of new tables, show their convergence with respect to important parameters, state some of their limitations, and document a mostly automated production process. In this section, we address the milestone completion criteria, found on the title page and the subsection headings. In section 2, we validate Opus against Tycho and compare with TabOp and Atomic. Section 3 details some of the physics implemented in Opus and relevant to this report. We then cover details of the density-temperature grid in section 4 and do a spectral comparison between different parameter settings and different codes in section 5. After summarizing these results in section 6, we give instructions how to run the opacity table scripts in appendix A and document the scripts and input/output files in appendix B.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Use and Limitations of a Climate-Quality Data Record to Study Temperature Trends on the Greenland Ice Sheet

Enhanced melting of the Greenland Ice Sheet has been documented in recent literature along with surface-temperature increases measured using infrared satellite data since 1981. Using a recently-developed climate-quality data record, 11- and 12-year trends in the clear-sky ice-surface temperature (IST) of the Greenland Ice Sheet have been studied using the Moderate-Resolution Imaging Spectroradiometer (MODIS) IST product. Daily and monthly MODIS ISTs of the Greenland Ice Sheet beginning on 1 March 2000 and continuing through 31 December 2010 are now available at 6.25-km spatial resolution on a polar stereographic grid as described in Hall et al. (submitted). This record will be elevated in status to a climate-data record (CDR) when more years of data become available either from the MODIS on the Terra or Aqua satellites, or from the Visible Infrared Imager Radiometer Suite (VIIRS) to be launched in October 2011. Maps showing the maximum extent of melt for the entire ice sheet and for the six major drainage basins have been developed from the MODIS IST dataset. Twelve-year trends of the duration of the melt season on the ice sheet vary in different drainage basins with some basins melting progressively earlier over the course of the study period. Some (but not all) of the basins also show a progressively-longer duration of melt. IST 12-year trends are compared with in-situ data, and climate data from the Modern Era Retrospective-Analysis for Research and Applications (MERRA) Reanalysis.

Hall, Dorothy K.↗
DOE OSTI2025

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,↗
DOE OSTI2020

Discretization limits of lattice‐Boltzmann methods for studying immiscible two‐phase flow in porous media

Summary Digital images of porous media often include features approaching the image resolution length scale. The behavior of numerical methods at low resolution is therefore important even for well‐resolved systems. We study the behavior of the Shan‐Chen (SC) and Rothman‐Keller (RK) multicomponent lattice‐Boltzmann models in situations where the fluid‐fluid interfacial radius of curvature and/or the feature size of the medium approaches the discrete unit size of the computational grid. Various simple, small‐scale test geometries are considered, and a drainage test is also performed in a Bentheimer sandstone sample. We find that both RK and SC models show very high ultimate limits: in ideal conditions the models can simulate static fluid configuration with acceptable accuracy in tubes as small as three lattice units across for RK model (six lattice units for SC model) and with an interfacial radius of curvature of two lattice units for RK and SC models. However, the stability of the models is affected when operating in these extreme discrete limits: in certain circumstances the models exhibit behaviors ranging from loss of accuracy to numerical instability. We discuss the circumstances where these behaviors occur and the ramifications for larger‐scale fluid displacement simulations in porous media, along with strategies to mitigate the most severe effects. Overall we find that the RK model, with modern enhancements, exhibits fewer instabilities and is more suitable for systems of low fluid‐fluid miscibility. The shortcomings of the SC model seem to arise predominantly from the high, strongly pressure‐dependent miscibility of the two fluid components.

Li, Zhe↗