Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid impact analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

Cross-Border Electricity Trade between India and Sri Lanka: Impact on Power System Operations

Recent developments in South Asia have the potential to increase the incentives for cross-border energy trade (CBET) between India and Sri Lanka. India's surplus power, including from its rapid expansion of wind and solar generation, provides an incentive for it to increase investment returns through exports. Sri Lanka is considering options to increase its power supply, and CBET could alleviate the need for investments in new thermal generation capacity, while giving it access to India's large and diverse power system. Finally, CBET can facilitate renewable energy (RE) integration by allowing the two countries to more efficiently harness balancing resources and smooth out wind and solar variability over larger geographic areas. This study builds on several prior analyses of a potential high-voltage direct current transmission link between India and Sri Lanka. Uniquely, the study makes use of a detailed one-year, 30-minute resolution unit commitment and economic dispatch (UC&ED) model of both countries to examine the technical and economic impacts of CBET between India and Sri Lanka in the year 2025. The study finds that CBET, enabled by a 500-MW high voltage direct current transmission link, could generate annual production cost savings of USD 180 million and improve power system operations. Through better utilization of the most efficient generators, generation from relatively expensive natural gas plants in Sri Lanka is displaced by lower cost imports from India. In addition, increased flexibility through exports reduces RE curtailment in India's Southern region by 400 GWh annually. While the study captures many previously unconsidered intertemporal operational constraints, it does not provide a complete characterization of system operations with a high-voltage direct current interconnection. Additional work should include a dynamic analysis of system stability following interruptions to the direct current interconnection, an AC power flow analysis and consideration of any non-technical barriers to trade which might prevent the optimization of CBET between India and Sri Lanka at 30-minute timescales.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A multi-scale calibration approach for process-oriented aggregated building energy demand models

Long-term energy planning relies largely on projections of future energy demand and hourly load profiles. Aggregate building models are increasingly being utilized to characterize the sensitivity of current and future building stocks to changes in climate, population, and building technology on city-to-regional-to-national scales. Due to challenges in the availability of data, those analyses have been limited to projection of energy demand at a single scale, usually state or country, while long-term planning power system models and production cost models might operate at different spatial scales such as energy regions which are focused on ensuring adequate generation infrastructure. In this work, we propose and evaluate a novel method to calibrate an aggregate building energy demand model (PNNL’s BEND model) against the best available data at the spatial scale of balancing authorities. This approach extends previous work on aggregated building energy demand by facilitating analysis of building energy demand across scales, in particular policy and operational decision-making scales. We show that the bias-corrected model estimates building electric loads reasonably well compared with estimates from a statistical model, but has the additional feature of flexibility across spatial scales. While the calibration approach is presently U.S.-centric and associated with U.S. energy regions, it can be extrapolated to other worldwide regions with similar scale challenges between policy and operational implementation decision making. We discuss the significant challenges involved in formulating and calibrating a complex physical model based on simulations of roughly 100,000 individual buildings against available aggregate regional electric load data and highlight areas for potential future work and improved data collection.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Biomass Burning Research Using DOE ARM Single-Particle Soot Photometer (SP2) Field Campaign Report

The focus of this laboratory study was to investigate the chemical and optical properties, and the detection efficiencies, of tar balls generated in the laboratory using the same instruments deployed on the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Climate Research Facility Gulfstream-1 (G-1) aircraft during the 2013 Biomass Burning Observation Project (BBOP) field study, during which tar balls were observed in wildland biomass burning particulate emissions. Key goals of this laboratory study were: (a) measuring the chemical composition of tar balls to provide insights into the atmospheric processes that form (evaporation/oxidation) and modify them in biomass burning plumes, (b) identifying whether tar balls contain refractory black carbon, (c) determining the collection efficiencies of tar balls impacting on the 600oC heated tungsten vaporizer in the Aerodyne Soot Particle Aerosol Mass Spectrometer (SP-AMS) (i.e., given the observed low volatilities, AMS measurements might underestimate organic biomass burning plume loadings), and (d) measuring the wavelength-dependent, mass-specific absorption cross-sections of brown carbon components of tar balls. This project was funded primarily by the DOE Atmospheric System Research (ASR) program, and the ARM Facility made their single-particle soot photometer (SP2) available for September 1-September 31, 2016 in the Aerodyne laboratories. The ARM mentor (Dr. Sedlacek) requested no funds for mentorship or data reduction. All ARM SP2 data collected as part of this project are archived in the ARM Data Archive in accordance with established protocols. The main objectives of the ARM Biomass Burning Observation Period (BBOP, July-October, 2013) field campaign were to (1) assess the impact of wildland fires in the Pacific Northwest on climate, through near-field and regional intensive measurement campaigns, and (2) investigate agricultural burns to determine how those biomass burn plumes differ from those from wildland fires. During BBOP, tar balls, small solid particles of organic substances, were observed downwind from wildland fires (at plume ages of 0-3 hours), but not agricultural burns. Observations of the tar balls on transmission electron microscope (TEM) grids suggest that they formed during atmospheric transport, likely due to the same atmospheric processes that increased the oxidation levels of the organic aerosol. Preliminary analyses suggest that tar balls may account for almost 50% of the total particle number, and 30% of the total organic particle mass, of the aerosol emitted from the burning events. These BBOP observations are described in detail in a manuscript in preparation (Sedlacek et al., 2017). The current laboratory study lasted four weeks and was conducted in the aerosol laboratories located at Aerodyne Research, Inc. in Billerica, Massachusetts. Tar balls were generated from several different biomass fuels, including samples from BBOP-related field sites, following literature procedures (Hoffer, Tóth, Nyirö-Kósa, Pósfai, and Gelencsér, 2016; Tóth, Hoffer, Nyirö-Kósa, Pósfai, and Gelencsér, 2014), and they were characterized using the same equipment used during the 2013 BBOP study, specifically the SP-AMS, SP2 and TEM. This study determined that laboratory-generated tar balls (1) are refractory with respect to TEM analysis in a similar manner to those collected during BBOP from wildland fires, (2) are composed of organic material with some refractory carbon components, (3) can be measured quantitatively by the SP-AMS, strengthening observations during BBOP, (4) absorb visible light, and (4) are dominated by unsaturated hydrocarbons that may be responsible for their light-absorbing properties. The results from this project are already being incorporated into our analysis of the formation processes and emission rates of tar balls as a function of fuel and combustion conditions from wildland fires.

Aerodyne Soot Particle Aerosol Mass Spectrometer↗

Modeling and Analysis for Spent Nuclear Fuel Seismic Testing

Note: This is the final public release version of PNNL-31671 Draft, which was previously released to the sponsor for review. The intent is the public release version will be PNNL-31671. There are no significant changes from the earlier draft version. This unlimited distribution report is the deliverable for M3SF-21PN010202014: Modeling and Analysis for Spent Nuclear Fuel Seismic Testing. This report summarizes the modeling, analysis, and test plan support completed by Pacific Northwest National Laboratory for the spent nuclear fuel dry storage system seismic test plan through May of 2021. Test plan preparation is planned to continue until the seismic test is completed in July of 2022. This report covers preliminary structural dynamic model development and computer-aided design of test hardware. Based on preliminary modeling, the strongest earthquakes under consideration in this test program provide mechanical loading on the spent nuclear fuel that is comparable to the 30 cm cask drop scenario, although the loads do not appear to be strong enough to cause significant permanent deformation of the fuel assembly spacer grids. The potential for grid deformation during the test will be assessed when the proposed shake table motion becomes available. The weakest earthquakes considered in this test are expected to be comparable to the mechanical loads witnessed in the multimodal transportation test of 2017. The horizontal canister system is predicted to provide nearly-uniform loading condition on the fuel assemblies it contains, while the vertical cask system is predicted to cause non-uniform dynamic loads on the fuel assemblies. In the horizontal case, gravity keeps the fuel assemblies in contact with one basket wall surface unless the loads are strong enough to cause a separation. In the vertical case, the fuel assemblies are relatively long and slender, and seismic motion in the anticipated test range is predicted to cause the assemblies to lean, tilt, and impact the basket walls throughout the seismic event. These gap closures are a nonlinear force transmission condition, so the vertical system is expected to have more variation and variability than the horizontal system. While the horizontal system is expected to have a relatively more linear response than the vertical system, there is still the potential for nonlinear behavior in the horizontal system because the fuel assemblies are free to slide, bounce, and impact the basket walls if the seismic loads are strong enough. One major conclusion of this study is that the use of mixed fuel assemblies in the canister will be acceptable. There are differences in overall system response when the mass, center of gravity, or gaps are changed, but the changes in response are within the bounds of a system that contains completely homogeneous fuel assemblies. One important observation is that the loading for each individual fuel assembly within the vertical canister is expected to be different from that of the others because of nonlinearities in the system. The horizontal canister system is expected to have a more uniform and more predictable response than the vertical canister, making variations in fuel assembly characteristics easier to account for. One important recommendation that comes from this modeling work is to repeat some of the strongest shake tests at different angles, particularly for the vertical cask system. Modeling predicts that strong seismic motion will cause the fuel assemblies to close initial gaps and impact the fuel basket walls when the canister is in the vertical cask configuration. Each shake test will impose a pre-defined three-dimensional time history on the cask system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Demonstration of Wake Steering Through Yaw Control in a Wind Plant Field Experiment: Cooperative Research and Development Final Report, CRADA Number CRD-16-00629

Over the last few decades, wind energy has evolved into a large international industry involving major players in the manufacturing, construction, and utility sectors. Coinciding with the industry’s growth, significant innovation in the technology has resulted in larger turbines with lower associated costs of energy and more complex designs in all subsystems. However, as the deployment of the technology has grown and its role within the electricity sector become more prominent, so have the expectations of the technology in terms of performance, reliability, and cost. The industry currently partitions its efforts into separate paths for turbine design, plant design and development, finance, grid interaction and operation, mitigation of adverse community and environmental impacts, and other areas. One prominent area where this partition is evident is in wind turbine control. Traditionally, each wind turbine in a wind plant has been controlled separately – via its own internal controller using only its own sensors. However, wind turbines in a plant interact with each other through the plant-level fluid dynamics. Wake losses (due to upstream turbines extracting energy from the winds and “waking” downstream turbines) can be up to 10% or even 20% of the gross energy production (if each turbine experienced the free stream wind inflow to the plant). A series of studies and experiments have demonstrated that there is potential for improving energy output at existing plants through plant control methods which seek to optimize total wind plant energy production over the current “greedy” approach where each turbine maximizes its own production. Wake steering induced by yaw offsets (turning the turbine to be out of the plane perpendicular to wind inflow) for upstream turbines has shown significant promise in simulations and wind tunnel experiments. In simulation studies, annual energy production has been shown to increase by 2% or more depending on the particular aspects of the wind plant (turbine spacing, meteorological conditions, etc). This project seeks to demonstrate the potential of plant-level controls via wake steering at a commercial wind plant. This is an important step towards commercialization and industry adoption of this plant-level modeling and analysis capability.

17 WIND ENERGY↗

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

Advanced Measurements for Resilient Integration of Inverter-Based Resources: PROGRESS MATRIX Final Report

As nearly every aspect of the electric power grid undergoes rapid change, measurement technologies that support grid operation and planning must evolve as well. The rapid large-scale deployment of inverter-based resources (IBRs) vital to achieving the nation’s clean energy goals has in some cases led to negative impacts on the reliability and security of the bulk power system (BPS). Advanced power system measurements, including synchronized phasor and waveform measurements, are key to making IBR integration secure and reliable. To this end, the Department of Energy (DOE) initiated the PROGRESS MATRIX project to develop advanced measurement capabilities and analytics that will accelerate adoption of IBRs while improving the reliability and resilience of the BPS. This report discusses the outcomes of the project, which was a joint effort between the Pacific Northwest National Laboratory (PNNL), Oak Ridge National Laboratory (ORNL), the National Renewable Energy Laboratory (NREL), and Lawrence Berkeley National Laboratory (LBNL). In the project’s first year, PNNL, NREL, and ORNL partnered with the Bonneville Power Administration (BPA), the Western Area Power Administration (WAPA), and Kauai Island Utility Cooperative (KIUC) to understand their existing measurement capabilities and the gaps limiting deployment of IBR-focused measurement systems and analytics. The other primary activity in the first year was deployment of GridSweep instruments, which provide unprecedented precision in waveform measurement while probing distribution systems. The instruments were deployed at Dominion Energy and the University of California, Riverside. In the project’s second year, the input from partner utilities and collected measurements were used to advance measurement capabilities. Twelve analytical methods spanning disturbance analysis, power plant evaluation, feeder evaluation, and modeling were developed. Two software tools were developed, one to analyze GridSweep measurements and another to automatically evaluate the control performance of power plants connected to the BPS. Testbeds at ORNL and NREL were augmented to better enable studies of IBR integration. The project culminated in demonstrations of these analytical methods, software tools, and testbeds, both in the field and in the laboratory. This report discusses these various accomplishments and documents the significant progress in developing advanced measurement capabilities to support the secure, reliable, and accelerated adoption of IBRs in the BPS.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Multilevel in Space and Energy Solver for 3-D Multigroup Diffusion and Coarse-Mesh Finite Difference Eigenvalue Problems

The Multilevel in Space and Energy Diffusion (MSED) method accelerates the iterative convergence of multigroup diffusion eigenvalue problems by performing work on lower-order equations with only one group and/or coarser spatial grids. It consists of two primary components: (1) a grey (one-group) diffusion eigenvalue problem that is solved via Wielandt-shifted power iteration (PI) and (2) a multigrid-in-space linear solver. In previous work, the efficiency of MSED was verified using Fourier analysis and numerical results from a one-dimensional multigroup diffusion code. Since that work, MSED has been implemented as a solver for the coarse-mesh finite difference (CMFD) system in the three-dimensional Michigan Parallel Characteristics Transport (MPACT) code. In this paper, the results from the implementation of MSED in MPACT are introduced, and the changes needed to make MSED more suitable for MPACT are described. For problems without feedback, the conclusions in this paper show that MSED can reduce the CMFD run time by an order of magnitude and the overall run time by a factor of 2 to 3 compared to the default CMFD solver in MPACT [PI with the generalized minimal residual (GMRES) method]. For problems with feedback, the convergence of the outer Picard iteration scheme is worsened by the well-converged CMFD solutions produced by the standard MSED method. To overcome this unintuitive deficiency, MSED may be run with looser convergence criteria (a modified version of the MSED method called MSED-L) to circumvent the issue until the multiphysics iteration in MPACT is improved. Results show that MSED-L can reduce the CMFD run time in MPACT by an order of magnitude, without negatively impacting the outer Picard iteration scheme.

3-D↗

Electromagnetic Transient Modeling of Large Data Centers for Grid-Level Studies: Beta Release

The magnitude and complexity of electricity usage patterns from large data centers are having significant impacts on the operation and dynamics of the power grid; grid operators and planners require a range of specialized data center models to properly evaluate these impacts and specify technical solutions as needed. Towards addressing this need, Pacific Northwest National Laboratory (PNNL) has developed a library of electromagnetic transient (EMT) models for grid-level studies of data centers called the data center model library (DML). This report describes how the DML was created and how it may properly be used. This report details the DML’s beta release, completed in July 2026. This is a revision and expansion of the alpha release, which was made available in January 2026 The models present in the DML are generic models; subject matter expertise and additional technical data are needed to modify these models before they can represent any real data center. However, they will significantly reduce the level of effort required to develop site-specific models and can serve as a common starting point to guide industry towards a more refined consensus. Most of the models within DML are dedicated to representing the power electronics interfaces commonly used in modern data centers, such as double-conversion uninterruptible power supplies and single-phase power factor correction converters. These models are intended for use in grid-level studies and are a simplified aggregation of many small components. That said, background material on the physical and electrical design of large data centers is provided as companion material so that users can be aware of many of the details which have been omitted or streamlined as a matter of practical necessity. Additionally, guidance on the application of EMT analysis for data center interconnection studies is provided, which aids users in identifying when the DML is necessary and what sort of additional model development may be necessary for conducting real-world studies.

electromagnetic transients↗

Feasibility and strategic implications of deploying nuclear power reactors in Africa

This report assesses the feasibility and strategic implications of deploying nuclear power reactors, including large-scale plants, advanced small modular reactors (SMRs), and microreactors, in African countries. Case studies focus on South Africa, Egypt, Kenya, Ghana, and Nigeria, examining nuclear energy’s role in Africa’s rapidly evolving energy landscape, marked by fast-growing demand, significant electricity access gaps, increasing renewable penetration, and strong policy commitments to industrialization and energy security. Several U.S. reactor technologies and designs are considered based on their development status and readiness for deployment. The analysis finds that nuclear power can provide reliable, clean baseload and flexible generation, as well as high-temperature process heat for desalination, hydrogen production, and industrial applications. However, suitability is highly country-specific, depending on grid size and stability, transmission capacity, cooling water availability, regulatory readiness, and fuel supply chains. Near-term deployment opportunities are strongest for light-water reactors (such as NuScale, BWRX-300, AP300, and SMR-300) that use low-enriched uranium and build on proven technology. More advanced concepts, including gas-cooled, sodium-cooled, molten-salt cooled reactors, and microreactors, will likely be relevant for African deployment in the 2030s or later, contingent on demonstration projects, high-assay low-enriched uranium (HALEU) fuel availability, and mature international licensing frameworks. Economic analysis shows that SMRs are capital-intensive, with projected overnight costs for 300 MWe units in 2025 ranging from approximately 1.4 to 2.6 billion USD per module. The levelized cost of electricity (LCOE) is highly sensitive to the weighted average cost of capital (WACC). Given typically higher financing costs and utility balance-sheet weaknesses in many African countries, bankable project structures will require sovereign guarantees, robust offtake arrangements, and layered financing from export credit agencies, development finance institutions, and vendor nations. Comparisons with recent large nuclear projects in the United Arab Emirates (UAE) and Egypt underscore the central role of state-backed loans, long tenors, and concessional terms. Country case studies illustrate a spectrum of readiness and opportunity. South Africa operates two 920 MWe pressurized light water reactors (totaling 1,840 MWe) at Koeberg and has the most mature regulatory and industrial base, positioning it as a prime candidate for both large reactors and SMRs to replace coal, support desalination, and anchor industrial hubs. Egypt is constructing four VVER-1200 units at El Dabaa with strong state leadership and could later complement this fleet with SMRs for coastal and industrial applications. Kenya and Ghana are advancing through IAEA Milestones with growing institutional capacity and clear interest in SMRs that match their smaller grids and industrialization plans. Nigeria has the largest demand potential but faces acute constraints in grid reliability, project bankability, and regulatory capacity; targeted deployments of large reactors and SMRs near coastal or industrial sites could have high impact if accompanied by major grid upgrades and institutional reforms. The report identifies cross-cutting challenges such as financing, political continuity, public acceptance, nonproliferation and security, waste and back-end management, regulatory capacity, grid adequacy, and long deployment timelines for first-of-a-kind designs, and ANL/NSE-26/3 ii proposes broad directions for resolution. These include stronger multifaceted financing for nuclear, long-term national energy strategies that transcend electoral cycles, proactive stakeholder engagement, strengthened regional and national regulators, and systematic workforce development through centers of excellence and expanded training. The United States should develop partnerships with African countries and offer end-to-end nuclear package similar to those used effectively by competitors: coordinated project development, state-backed financing, long-term fuel services, and durable in-country support through regional offices and sustained workforce/regulatory training. With timely planning, sustained political commitment, and appropriate financing and institutional support, nuclear energy, both large reactors and advanced SMRs, can become a meaningful, though not dominant, pillar of Africa’s future power mix, enhancing energy security, enabling industrial growth, and supporting climate goals.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Water-Energy Considerations in California’s Agricultural Sector and Opportunities to Provide Flexibility to California’s Grid

Water and energy use on farms are deeply linked, and trends changing the way water is used on farms have implications for the electrical grid, and vice versa. In 2014, California’s agricultural sector accounted for approximately 80% of the state’s water consumption and 5% of its electricity use. The majority of the electricity was for surface conveyance and groundwater pumping. Diminishing surface water availability and increasing reliance on groundwater pumping are likely to increase future energy demand. These trends are coincident with farms embracing new precision irrigation technologies that increase water productivity and crop yield, but may have the potential of other collateral effects, such as increasing energy intensity of irrigation and total electricity use, and thus demand on the grid. The impacts of these technologies and their loads on the electrical grid have been explored only to a limited degree, and the potential grid benefits afforded by these technologies have been studied even less. For example, studies show that California’s increasing adoption of renewable electricity may result in excess generation during certain periods of the day, and that responsive load management (i.e. demand response) may be a cost-effective method to shift loads to utilize this excess generation. The degree to which farms in California do or can provide responsive loads has not been explored. In this paper, we explore whether farms that are adopting technologies and automation for irrigation management can provide co-benefits to the grower and the electrical grid. To do so, we define the characteristics of future grid needs specifically related to responsive load management, explore existing market mechanisms and policies for proving those resources, estimate the temporal distribution of energy use for irrigation in California, and develop scenarios of how irrigation loads might meet these needs. We conclude by exploring potential ancillary service benefits and associated revenues that California’s agricultural industry may provide. Ancillary services include using loads to dynamically adjust demand on the system to alleviate short-run ramps and disturbances at timescales ranging from seconds up to an hour. With higher penetrations of renewable resources, larger quantities of ancillary services will be necessary and may become more valuable to the grid than traditional demand response. Ancillary services from agricultural loads is a nascent topic, and no such analysis has been conducted to quantify its benefits to our knowledge.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Four-Dimensional Oceanic and Atmosperic Data Assimilation with Tropical Rainfall Measuring Mission Data

An oceanic data assimilation system which allows to utilize the forthcoming Tropical Rainfall Measuring Mission (TRMM) data has been developed and applied to the Pacific Ocean to produce the velocity field. The assimilated data will be indispensable to examine the effects of rainfall and its variability on the structure and circulation of the tropical oceans and to assess the impact of global warming due to the increase of carbon dioxide on the ocean circulation system and the marine pollution caused by oil spill and ocean damping of radionuclide. The data will also provide the verification for the oceanic and ocean-atmosphere coupled General Circulation Models (GCM's). The system consists of oceanic GCM, analysis scheme and data. In the system the flow field has been determined to be physically consistent with the observed density field and the sea surface winds derived from the Special Sensor Microwave Imagery (SSM/I) data which drive the ocean current. The time integration has been performed for five years until the flow field near the surface attained the steady state starting from the rest ocean with observed temperature and salinity fields, and the SSM/I surface wind velocity. The resultant flow field showed high producibility of the system. Especially the flow near the ocean surface agreed well with available observed data. The system, for the first time, succeeded to produce the eastward subtropical current which has been discovered in the joint investigation on Kuroshio current (CSK) in the 1960s. To verify the quality of the flow field a trajectory analysis has been carried out and compared with the Algos buoy data. BRIEF DESCRIPTION OF THE DATA ASSIMILATION SYSTEM ## Oceanic GCM and analysis scheme--The basic equations are much the same as used for the GCM's, except for the Newtonian damping terms introduced into the prediction equations for the potential temperature and salinity to maintain these fields as observed. The C grid of 2'lat. by 2'long. in horizontal and the 11 vertical levels are applied to the entire Pacific Ocean. At the east and west ocean boundaries the periodic boundary conditions are applied creating fictitious ocean there. The SMAC Method is used to increase the accuracy of mass conservation. * Data--The JODC temperature and salinity data obtained from 1906 to 1988 are used in the system between Long.100'E. and 60'W. The surface wind data are derived from the SSM/I data by Dr-R. Atlas of NASA/GSFC. The data set contains every 6 hours data from July 1987 to June 1989 on the grid of 2'lat. by 2.5'long. The averaged for the whole period and then interpolated into the 2'lat. by 2'long. grid data are used to force the system. The sea bottom topography data was based on the General Bathymetric Chart of the Ocean (GEBCO) supplied by the Canadian Hydrographic Service under contract with the International Hydrographic Organization and International Oceanographic Commission of UNESCO.

Takano, Kenji↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

The path to clean energy: direct coupling of nuclear and renewable technologies for thermal and electrical applications

The U.S. Department of Energy (DOE) recognizes the need to transform the energy infrastructure of the U.S. and elsewhere to systems that can significantly reduce environmental impacts in an efficient and economically viable manner while utilizing both clean energy generation sources and hydrocarbon resources. Thus, DOE is supporting research and development that could lead to more efficient utilization of clean nuclear and renewable energy generation sources. A concept being advanced by the DOE Offices of Nuclear Energy (NE) and Energy Efficiency and Renewable Energy (EERE) is tighter coupling of nuclear and renewable energy sources in a manner that better optimizes energy use for the combined electricity, industrial manufacturing, and the transportation sectors. This integration concept has been referred to as a “hybrid system” that is capable of providing energy (thermal or electrical) where it is needed, when it is needed. For the purposes of this work, the hybrid system would integrate two or more energy resources to generate two or more products, one of which must be an energy commodity, such as electricity or transportation fuel. This definition requires coupling of subsystems ‘‘behind’’ the electrical transmission bus, where energy flows are dynamically apportioned as necessary to meet demand and the system has a single connection to the grid that provides dispatchable electricity as required while capital intensive generation assets operate at full capacity. Development of integrated energy systems for an “energy park” must carefully consider the intended location and the associated regional resources, traditional industrial processes, energy delivery infrastructure, and markets to identify viable region-specific system configurations. This paper will provide an overview of the current status of regional hybrid energy system design, development and application of dynamic analysis tools to assess technical and economic performance, and roadmap development to identify and prioritize component, subsystem and system testing that will lead to prototype demonstration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Precipitation characteristics of CAM5 physics at mesoscale resolution during MC3E and the impact of convective timescale choice

Abstract The physics suite of the Community Atmosphere Model version 5 (CAM5) has recently been implemented in the Weather Research and Forecasting (WRF) model to explore the behavior of the parameterization suite at high resolution and within the controlled setting of a limited area model. The initial paper documenting this capability characterized the behavior for a northern high‐latitude region. This paper characterizes the precipitation characteristics for continental, midlatitude, springtime conditions during the Midlatitude Continental Convective Clouds Experiment (MC3E) over the central United States. This period exhibited a range of convective conditions from those driven strongly by large‐scale synoptic regimes to more locally driven convection. The study focuses on the precipitation behavior at 32 km grid spacing to better anticipate how the physics will behave in a global model when used at similar grid spacing in the coming years. Importantly, one change to the Zhang‐McFarlane deep convective parameterization when implemented in WRF was to make the convective timescale parameter an explicit function of grid spacing. This study examines the sensitivity of the precipitation to the default value of the convective timescale in WRF, which is 600 s for 32 km grid spacing, to the value of 3600 s used for 2° grid spacing in CAM5. For comparison, a 1200 s and an infinite convective timescale are also used. The results show that the 600 s timescale gives the most accurate precipitation amount over the central United States. However, this setting has the worst precipitation diurnal cycle, with the convection too tightly linked to the daytime surface heating. Longer timescales greatly improve the diurnal cycle but result in less precipitation and produce a low bias. An analysis of rain rates shows the accurate precipitation amount with the shorter timescale is assembled from an over abundance of drizzle combined with too few heavy rain events. With longer timescales, one can improve the frequency distribution, particularly for the extreme rain rates. Ultimately, without changing other aspects of the physics, one must decide between accurate diurnal timing and rain amount when choosing an appropriate convective timescale.

54 ENVIRONMENTAL SCIENCES↗

Memory forensic analysis of a programmable logic controller in industrial control systems

In industrial control systems (ICS), programmable logic controllers (PLCs) are used to automate physical processes such as nuclear plants and power grid stations, and are often subject to cyber attacks. As in conventional IT domain, the memory analysis of the PLCs can help answer important forensic questions about the attack, such as the presence of malicious firmware, injection of modified control logic (the program running on the PLC), and manipulation of I/O devices (e.g., sensors and actuators). Unlike conventional IT domain, PLCs have heterogeneous hardware architecture, proprietary firmware and control software, making it challenging to employ a unified framework for their memory forensics. For merely extracting artifacts of forensic importance, reverse-engineering the firmware is a tedious task, and the effort needs to be repeated for every PLC model. As a community, a step-wise approach to tackle this challenge is to analyze the memory of specific PLCs, and subsequently find a generic framework applicable to all PLCs. Our work is a step forward in this direction. By following a methodology that focuses on the functional layer of PLCs instead of reverse engineering the firmware, we analyze the digital forensic artifacts available in a common PLC, Allen-Bradley ControlLogix 1756-L61. Before diving into the memory dump, we analyze the PLC control software to create a list of important artifacts that are sure to exist in the PLC memory dump. The approach employs a setup where PLC control software RSLogix-5000 is connected to the PLC, and the memory dump can be obtained as and when needed. We create test cases that sequentially highlight each category of artifacts, followed by an examination of the resultant impact on memory. After attaining the listed artifacts, we employ conventional string and known data searches to extract interesting information present in this PLC's memory. The memory analysis profile, presented as a Python library and shared with the community, can help a forensic investigator to readily extract forensic artifacts from the same model's controller. The adopted approach may help researchers in creating memory profile of other PLCs, and ultimately formulating a generic PLC memory analysis framework.

Rais, Muhammad Haris↗

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗

The National Climate Data Base (NCDB): A Bias-Corrected High-Resolution Climate Dataset

Assessing renewable energy resources under future climate scenarios has been highlighted in recent years to analyze and understand potential impacts of future change in renewable generation on the power sector. Solar energy is well-known as the most plentiful among various renewable resources and usually converted to electricity using photovoltaics (PV) technologies, and the global deployment of PV technology has increased rapidly in recent decades. In this study, we develop a statistical technique to downscale the future projection of solar irradiance for PV energy-related applications. A set of Regional Climate Model (RCM)-based projections obtained from the North American Coordinated Regional Climate Downscaling Experiment (NA-CORDEX) are used as inputs to statistical methods to generate high-resolution global horizontal irradiance (GHI) over the contiguous United States (CONUS). The main steps of the statistical downscaling method include (1) regridding RCM output (0.22 degree and daily resolutions) to handle the modeled-observed data sets on a common grid, (2) correcting bias of RCM GHI using satellite-derived observation, and (3) implementing temporal and spatial downscaling to generate GHI at 8-km and hourly resolution. Basically, complex physical processes and interactions between solar radiation and various atmospheric constituents lead solar irradiance to be highly variable and uncertain. Underrepresentation of clouds from the RCM parameterizations is the main source of error and uncertainty in modeling solar irradiance. Thus, we adapt and use the high-quality satellite-derived data from the National Solar Radiation Database (NSRDB) to analyze the bias and error of RCM GHI as well as estimate the statistical parameters for spatial and temporal downscaling. This presentation will summarize the comprehensive analysis conducted to produce and assess the results under two climate scenarios (RCP4.5 and RCP8.5). We will also present a detailed validation demonstrating the strengths of the proposed downscaling method and future extension of this research.

climate data↗