Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “synthetic grid”

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 307 records · Page 17

Heavy-Duty Truck Electrification and the Impacts of Depot Charging on Electricity Distribution Systems

Major technological advancements and recent policy support are improving the outlook for heavy-duty truck electrification in the United States. In particular, short-haul operations (=200 miles (=322 km)) are prevalent and early candidates for plug-in electric vehicles (EVs) given their short, predictable routes and return-to-base applications, which allows vehicles to recharge when off shift at their depots. Although previous studies investigated the impacts of added electrical loads on distribution systems, which included light-duty EVs, the implications for heavy-duty EV charging are underexplored. Here we summarize the causes, costs and lead times of distribution system upgrades anticipated for depot charging. We also developed synthetic depot charging load profiles for heavy-duty trucks from real-world operating schedules, and found that charging requirements are met at common light-duty EV charging rates (=100 kW per vehicle). Finally, we applied depot charging load profiles to 36 distribution real-world substations, which showed that most can accommodate high levels of heavy-duty EV charging without upgrades.

47 OTHER INSTRUMENTATION↗

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

Multi-area parameter error identification for large power systems

Power grid model parameters may contain errors due to various reasons. Detecting and correcting parameter errors typically requires significant computational effort due to the size and complexity of the parameter database. While the normalized Lagrange multiplier (NLM) method can effectively detect, identify and correct parameter errors, its computational burden could rapidly grow with increasing system size. This paper addresses this issue by proposing a multi-area parameter error identification method. Each area has its own outlier detection tool for detecting the incorrect parameters and measurements within the area. On the other hand, due to the reduced redundancy at area boundaries, parameter errors on branches incident to boundary buses may not be detected. Such errors are subsequently detected by a coordination level estimator completing the system-wide parameter detection procedure. In conclusion, performance of the developed method is demonstrated using the IEEE 118-bus and 2000-bus Texas synthetic systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Life Cycle Cost Analysis of Prestressed Concrete Poles Subjected to Wind, Surges, and Waves

Prestressed concrete (PC) poles are becoming popular choices to support coastal power transmission systems. However, the existing literature does not offer a detailed analysis of the effectiveness of PC poles in terms of long-term vulnerabilities and the direct and indirect costs. This is due to (1) lack of fragility models for PC poles and (2) lack of probabilistic wind, storm surge, and wave models in coastal settings. In this study, we address these gaps through a series of Monte Carlo simulations to estimate fragility of PC poles as a function of age and hazard (wind, surges, and waves) intensity, and the development of a probabilistic hazard model based on 10,000 years of synthetic tropical cyclone data. The probabilistic hazard model is used in conjunction with high-resolution hydrodynamic models to generate realizations of coastal wind, storm surges, and waves for the Louisiana and Mississippi coasts. A comprehensive life cycle cost analysis for a service life of 70 years considering direct and indirect losses is conducted to compare the performance of a transmission line located in Pascagoula, Mississippi, when wood poles are replaced by PC poles. Results showed that aging has a minor effect on the reliability of PC poles, highlighting the advantages of replacing wood poles with PC poles, especially in coastal areas. In addition, PC poles are significantly more cost-effective compared with wood poles over their life cycle, leading to an estimated saving of $11.55 million (68.17% reduction). The results of this study provide key insight to inform decision-making processes to keep the coastal power grids resilient and cost-effective against future storm hazards.

natural disasters↗

Forecasting of Storm-Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm-surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm-surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm-surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and, most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite-element-based codes. It affords a capability for simulating tidal circulation and storm-surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm-surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate. To keep the computational load manageable, this work was conducted using only the wind stress, calculated by using historical data from Hurricane Camille, as the input condition for the model. Hurricane storm-surge simulations were performed on an eight-node Linux computer cluster. Each node contained dual 2-GHz processors, 2GB of memory, and a 40GB hard drive. The digital elevation model (DEM) for this region was specified using a combination of Navy data (over water), NOAA data (for the coastline), and optimized Interferometric Synthetic Aperture Radar data (over land). This high-resolution topographical data of the Mississippi coastal region provided the ADCIRC model with improved input with which to calculate improved storm-surge forecasts.

Valenti, Elizabeth↗

Achieving Accuracy Requirements for Forest Biomass Mapping: A Data Fusion Method for Estimating Forest Biomass and LiDAR Sampling Error with Spaceborne Data

The synergistic use of active and passive remote sensing (i.e., data fusion) demonstrates the ability of spaceborne light detection and ranging (LiDAR), synthetic aperture radar (SAR) and multispectral imagery for achieving the accuracy requirements of a global forest biomass mapping mission. This data fusion approach also provides a means to extend 3D information from discrete spaceborne LiDAR measurements of forest structure across scales much larger than that of the LiDAR footprint. For estimating biomass, these measurements mix a number of errors including those associated with LiDAR footprint sampling over regional - global extents. A general framework for mapping above ground live forest biomass (AGB) with a data fusion approach is presented and verified using data from NASA field campaigns near Howland, ME, USA, to assess AGB and LiDAR sampling errors across a regionally representative landscape. We combined SAR and Landsat-derived optical (passive optical) image data to identify forest patches, and used image and simulated spaceborne LiDAR data to compute AGB and estimate LiDAR sampling error for forest patches and 100m, 250m, 500m, and 1km grid cells. Forest patches were delineated with Landsat-derived data and airborne SAR imagery, and simulated spaceborne LiDAR (SSL) data were derived from orbit and cloud cover simulations and airborne data from NASA's Laser Vegetation Imaging Sensor (L VIS). At both the patch and grid scales, we evaluated differences in AGB estimation and sampling error from the combined use of LiDAR with both SAR and passive optical and with either SAR or passive optical alone. This data fusion approach demonstrates that incorporating forest patches into the AGB mapping framework can provide sub-grid forest information for coarser grid-level AGB reporting, and that combining simulated spaceborne LiDAR with SAR and passive optical data are most useful for estimating AGB when measurements from LiDAR are limited because they minimized forest AGB sampling errors by 15 - 38%. Furthermore, spaceborne global scale accuracy requirements were achieved. At least 80% of the grid cells at 100m, 250m, 500m, and 1km grid levels met AGB density accuracy requirements using a combination of passive optical and SAR along with machine learning methods to predict vegetation structure metrics for forested areas without LiDAR samples. Finally, using either passive optical or SAR, accuracy requirements were met at the 500m and 250m grid level, respectively.

Montesano, P. M.↗

The Essential Role of Nuclear Energy in Achieving Economy-wide Net-Zero Solutions

Governments and private industry around the world have established aggressive goals to achieve net-zero emissions for the power, industrial, and transportation sectors by 2050. These aggressive goals demand immediate action if we are to be successful, and they require us to think more holistically about our clean energy options. Programs within the U.S. Department of Energy (DOE) are addressing these holistic solutions. Traditionally, electricity generation and management and meeting energy demands for industry and transportation are considered independently. As we seek to achieve net-zero, we need to reassess energy demands. When we consider overall energy use, only one-third is in the form of electricity. Additional energy demands are in the form of heat or steam for industrial processes, as well as transportation. These sectors are much harder to abate, and electrification may not be the best option. Reducing environmental emissions at an affordable cost, while maintaining grid reliability and resilience, will require us to use all of the clean energy resources that we have available. That means coordinating the use of nuclear, renewables, and fossil fuels with carbon capture to meet growing demands for electricity, industrial applications, and mobility. The DOE Office of Nuclear Energy (DOE-NE) program on Integrated Energy Systems (IES) is led by researchers at Idaho National Laboratory (INL) and work is conducted in partnership with an array of other DOE laboratories, industry, and academia. The primary focus of IES research is to assess the technical and economic potential of IES to enhance the flexibility and utilization of nuclear reactors working alongside renewable generators to meet an array of energy demands—thereby maximizing the utilization of clean energy resources across all energy sectors. Various energy applications and product streams beyond electricity are being evaluated, ranging from generation of potable water to production of hydrogen, fertilizers, synthetic fuels, and chemicals. The DOE-NE program additionally partners with the Hydrogen and Fuel Cell Technologies Office under the DOE Office of Energy Efficiency and Renewable Energy to jointly fund the development of analysis tools, technologies, and nuclear-integrated hydrogen demonstration projects. This presentation highlights the wide array of R&D being conducted across multiple DOE-funded programs to develop and deploy nuclear-based IES that will be key to achieving our net-zero goals. By working with key collaborators in the nuclear industry, analytical studies are now becoming a reality in demonstration projects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Advancing Nuclear Energy to Support Economy-wide Net-zero Solutions: Challenges and Opportunities

Governments and private industry around the world have established aggressive goals to achieve net-zero emissions for the power, industrial, and transportation sectors by 2050. These aggressive goals demand immediate action if we are to be successful, and they require us to think more holistically about our clean energy options. Programs within the U.S. Department of Energy (DOE) are addressing these holistic solutions. Traditionally, electricity generation and management and meeting energy demands for industry and transportation are considered independently. As we seek to achieve net-zero, we need to reassess energy demands. When we consider overall energy use, only one-third is in the form of electricity. Additional energy demands are in the form of heat or steam for industrial processes, as well as transportation. These sectors are much harder to abate, and electrification may not be the best option. Reducing environmental emissions at an affordable cost, while maintaining grid reliability and resilience, will require us to use all of the clean energy resources that we have available. That means coordinating the use of nuclear, renewables, and fossil fuels with carbon capture to meet growing demands for electricity, industrial applications, and mobility. The DOE Office of Nuclear Energy (DOE-NE) program on Integrated Energy Systems (IES) is led by researchers at Idaho National Laboratory (INL) and work is conducted in partnership with an array of other DOE laboratories, industry, and academia. The primary focus of IES research is to assess the technical and economic potential of IES to enhance the flexibility and utilization of nuclear reactors working alongside renewable generators to meet an array of energy demands—thereby maximizing the utilization of clean energy resources across all energy sectors. Various energy applications and product streams beyond electricity are being evaluated, ranging from generation of potable water to production of hydrogen, fertilizers, synthetic fuels, and chemicals. The DOE-NE program additionally partners with the Hydrogen and Fuel Cell Technologies Office under the DOE Office of Energy Efficiency and Renewable Energy to jointly fund the development of analysis tools, technologies, and nuclear-integrated hydrogen demonstration projects. This presentation highlights the wide array of R&D being conducted across multiple DOE-funded programs to develop and deploy nuclear-based IES that will be key to achieving our net-zero goals. By working with key collaborators in the nuclear industry, analytical studies are now becoming a reality in demonstration projects. The presentation specifically focuses on motivations for the paradigm shift in how nuclear energy is deployed and used, opportunities for nuclear energy to support an array of non-electric application, integration options, computational tools and experimental systems to support evaluation and to accelerate deployment, and industry partnerships.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Formulations and Valid Inequalities for Optimal Black Start Allocation in Power Systems

The restoration of a power system after a blackout starts around units with enhanced technical capabilities, referred to as black start units (BSUs). We examine the planning problem of optimally allocating these units on the grid subject to a budget constraint. We present a mixed integer programming model based on current literature in power systems. Binary variables are associated with the allocation of BSUs and with the energization state of buses, branches, and generators of the power system over a time horizon. We extract a substructure of the feasible region which imposes the requirement that each island that appears during the restoration process must have at least one operational generator. We discuss three equivalent reformulations for this requirement. We introduce a family of exponentially many, polynomially separable, valid inequalities to strengthen the formulation. Under simplifying assumptions, we show that the convex hull of the feasible region is a full-dimensional polyhedron, and prove that some of the constraints we introduced are facet-defining. We perform experiments to examine the difference in strength between the formulations as well as the computational times to solve the problem to near optimality for synthetic instances of the IEEE-39, IEEE-118, Illinois-200, WECC-225, IEEE-300, South Carolina-500, and Texas-2000 power systems. We illustrate a use case of the model. We conclude by suggesting extensions of the current work for future research.

Black Start Allocation↗

Demand Response Analysis for Different Residential Personas in a Comfort-Driven Behavioral Context

Low demand response (DR) participation and high program drop-out rates continue to impede DR goals that could save up to $13 billion in annual grid expansion and electricity demand costs. Yet, the literature lacks a thorough understanding of how different residential customer segments enrolled in DR programs respond to utility signals in view of occupant comfort considerations. The objective of this study is to gain a clear understanding of the effects of four different customer personas on residential DR. Given current data limitations, this work developed an array of hypothetical personas with varied priorities, activity levels, and comfort thresholds based on demographic variables that have been found in previous studies to influence energy consumption. A BEopt DR model for a reference residential single-family building located in Colorado was built to isolate the effect of differences in buildings or climate. The results provide useful evidence on how persona-comfort differences lead to significant deviations in DR goals (especially peak demand reduction), ranging from 0.1% to 20%. This work presents a novel framework representing comfort preferences in DR models. The data generated, albeit synthetic, and the results could inform DR program design considerations of how different people respond to different comfort priorities.

BEopt↗

Comprehensive Life Cycle Analysis of Methanol Production and Methanol-to-Diesel Conversion

Methanol is a strategic chemical and intermediate in the manufacture of synthetic diesel, due to its versatility, diesel’s compatibility with existing infrastructure, and their role in industrial and transport applications. Conventional production methods for methanol, primarily steam methane reforming (SMR), rely on natural gas and are subject to the price variability due to market conditions and geopolitical events. They are also associated with greenhouse gas (GHG) emissions. Methanol and synthetic diesel production could be integrated with nuclear energy to stabilize fuel prices and insulate pricing from outside geopolitical events due to the relative stability of nuclear fuel as compared to natural gas. This could lead to increased transportation fuel security, reliability and resilience. An added benefit is the abatement of emissions when substituting nuclear energy for conventional energy from natural gas. This report presents a comprehensive life cycle analysis(LCA) framework which was developed to evaluate the GHG emissions reduction potential associated with nuclear integrated methanol production, methanol-to-diesel upgrading, and end-use combustion. Gate-to-gate methanol production and cradle-to-grave emissions were evaluated in, starting with a business-as-usual (BAU) SMR-based methanol plant, and then considering stepwise nuclear integration. Methanol-to-diesel (MTD) conversion was evaluated accounting for nuclear energy integration and hydrogen production via high-temperature steam electrolysis (HTSE) using electricity either from the grid or from a dedicated nuclear power system. This multi-step process diverts stable and reliable nuclear energy into the transportation sector by upgrading low energy dense natural gas into liquid fuels fully compatible with existing infrastructure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Supporting ARPA-E Power Grid Optimization (Final Report)

Pacific Northwest National Laboratory (PNNL), Arizona State University (ASU), Georgia Institute of Technology (Georgia Tech), Los Alamos National Laboratory (LANL), National Renewable Energy Laboratory (NREL), Texas A&M University (TAMU), The University of Texas at Austin (UT), and the University of Wisconsin-Madison (UW-M) supported the ARPA-E Grid Optimization (GO) Competition by providing a common problem formulation, data format, datasets, evaluation mechanism, scoring, rules, and results that resulted in the awarding of $\$9.24$ million dollars to teams from academia, industry, and national labs for solving three sets of increasingly difficult non-linear, security- constrained AC Optimal Powerflow (AC-OPF) optimization problems in order to increase the efficiency of the US Electric Grid. It is estimated that a 1% increase in efficiency can save $\$1$ billion. Current industry practices typically use a linear DC model (DC-OPF) in order solve the OPF problem within the time constraints of the operation schedule. The GO Competition challenges the best power engineers, mathematicians, and computer scientists to make possible operational decisions based on accurate physical models. To accomplish this, the GO Competition created a series of Challenges and funded teams to produce the best solver. Challenge 1 was to solve the security constrained Alternating Current Optimal Power Flow (ACOPF) problem. Challenge 2 extended that to by adding adjustable transformer tap ratios, phase shifting transformers, switchable shunts, price-responsive demand, ramp rate constrained generators and loads, and fast-start unit commitment (UC). Furthermore, Challenge 2 was a maximization problem while Challenge 1 was a minimization problem. While Challenge 3 was being developed, the entrants were invited to find better solutions to the Challenge 2 synthetic datasets with no restrictions on time, hardware, or algorithms. The Challenge 2 solutions turned out to be very good. Challenge 3 expanded the Challenge 2 problem further by using multiperiod dynamic markets, including advisory models for extreme weather events, day-ahead markets, and the real-time markets with an extended look-ahead. These problems included active bid-in demand and topology optimization. Together the Challenges used nearly 30 million CPU hours. Since each team was working on the same problem, using the same data, and running on the same hardware, fair comparisons could be drawn as to the best solver. The datasets were varied enough, however, that the best solver for one dataset was not necessarily the best at another, so cumulative scores were used. The process was managed by the PNNL maintained website https://GOCompetition.energy.gov, where Entrants could find information about the problem, the data, the rules, submit their solver for evaluation, and see the scores of all the competing teams on a Leaderboard. Interest was world-wide but only American teams were eligible for prizes. The Competition has produced 34 journal articles 115 papers and been cited over 500 times in the literature, including 12 dissertations (4 from foreign countries; Columbia (2), Germany, and Italy) and 3 from the DOE ExaScale project. Software developed by Pearl Street Technologies for Challenges 1 and 2 is now deployed by Southwest Power Pool (SPP) and Midcontinent Independent Service Operator (MISO). Other teams have received inquiries from venture capitalists. Google DeepMind has thanked the Competition for making the datasets developed for the Competition public. They are using it to train machine learning models. The larger datasets have billions of unknowns to be solved for, but only a small percent matter in the final solution. Knowing what unknowns are important can dramatically speedup the solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Compactly‐Supported Nonstationary Kernels for Computing Exact Gaussian Processes on Big Data

The Gaussian process (GP) is a widely used method for analyzing large-scale data sets, including spatio-temporal measurements of nonlinear processes that are now commonplace in the environmental sciences. Traditional implementations of GPs involve stationary kernels (also termed covariance functions) that limit their flexibility, and exact methods for inference that prevent application to data sets with more than about 10,000 points. Modern approaches to address stationarity assumptions generally fail to accommodate large data sets, while all attempts to address scalability focus on approximating the Gaussian likelihood, which can involve subjectivity and lead to inaccuracies. In this work, we explicitly derive an alternative kernel that can discover and encode both sparsity and nonstationarity. We embed the kernel within a fully Bayesian GP model and leverage high-performance computing resources to enable the analysis of massive data sets. We demonstrate the favorable performance of our novel kernel relative to existing exact and approximate GP methods across a variety of synthetic data examples. Furthermore, we conduct space–time prediction based on more than 1 million measurements of daily maximum temperature and verify that our results outperform state-of-the-art methods in the Earth sciences. More broadly, having access to exact GPs that use ultra-scalable, sparsity-discovering, nonstationary kernels allows GP methods to truly compete with a wide variety of machine learning methods.

Gaussian processes↗

Comparison of source-location algorithms for atmospheric samplers

Numerous algorithms have been developed to determine the source characteristics for an atmospheric release, e.g., (Bieringer et al., 2017). This study compares three models that have been applied to the data collected by the International Monitoring System operated by the Comprehensive Nuclear-Test-Ban Treaty Preparatory Commission to estimate source event parameters. Each model uses a different approach to estimate the parameters. A deterministic model uses a possible source region (PSR) approach (Ringbom et al., 2014) that is based on the correlation between predicted and measured sample values. An unnamed model developed at Pacific Northwest National Laboratory (herein called BAYEST) uses a Bayesian formulation (Eslinger et al., 2019; Eslinger et al., 2020; Eslinger and Schrom, 2016). The FREAR model uses a different Bayesian formulation (De Meutter and Hoffman, 2020; De Meutter et al., 2021; Pieter De Meutter, 2021). The performance of the three source-location models is evaluated with 100 synthetic release cases for the single xenon isotope, 133 Xe. The release cases resulted in detections in a fictitious network with 120 noble gas samplers. All three source-location models use the same sampling data. The two Bayesian models yield more accurate location estimates than the deterministic PSR model, with FREAR having slightly better location performance than BAYEST. Samplers with collection periods of 3, 6, 8, 12, and 24-h were used. Results from BAYEST show that location accuracy improves with each reduction in sample collection length. The BAYEST model is slightly better for estimating the start time of the release. The PSR model has about the same spread in start times as the FREAR model, but the PSR results have a better average start time. The Bayesian source-location algorithms give more accurate results than the PSR approach, and provide release magnitude estimates, while the base PSR model does not estimate the release magnitude. This investigation demonstrates that a reasonably dense sampling grid will sometimes yield poor location and time estimates regardless of the model. The poor estimates generally coincide with cases where there is a much larger distance between the release point and the first detecting sampler than the average sampler spacing.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗

Projected increases in tropical cyclone-induced U.S. electric power outage risk

Abstract While power outages caused by tropical cyclones (TCs) already pose a great threat to coastal communities, how—and why—these risks will change in a warming climate is poorly understood. To address this need, we develop a robust machine learning model to capture TC-induced power outage risk. When applied to 900 000 synthetic TCs downscaled from simulated historical and future climate conditions under a strong warming scenario, we find outage risk in the United States and Puerto Rico is expected to increase broadly by the end of the century, with some states seeing increases of 60% and higher. Further, we discover that rising rainfall rates will play an increasingly important role in TC-induced power outage risk as the climate changes, explaining more than 50% of the projected change in risk in some regions. These insights are important for guiding decision-makers in their future outage risk investment and mitigation plans.

Grid Resilience↗