Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed model predictive control”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

WE-DE-202-03: Modeling of Biological Processes - What Happens After Early Molecular Damage?

Radiation therapy for the treatment of cancer has been established as a highly precise and effective way to eradicate a localized region of diseased tissue. To achieve further significant gains in the therapeutic ratio, we need to move towards biologically optimized treatment planning. To achieve this goal, we need to understand how the radiation-type dependent patterns of induced energy depositions within the cell (physics) connect via molecular, cellular and tissue reactions to treatment outcome such as tumor control and undesirable effects on normal tissue. Several computational biology approaches have been developed connecting physics to biology. Monte Carlo simulations are the most accurate method to calculate physical dose distributions at the nanometer scale, however simulations at the DNA scale are slow and repair processes are generally not simulated. Alternative models that rely on the random formation of individual DNA lesions within one or two turns of the DNA have been shown to reproduce the clusters of DNA lesions, including single strand breaks (SSBs), double strand breaks (DSBs) without the need for detailed track structure simulations. Efficient computational simulations of initial DNA damage induction facilitate computational modeling of DNA repair and other molecular and cellular processes. Mechanistic, multiscale models provide a useful conceptual framework to test biological hypotheses and help connect fundamental information about track structure and dosimetry at the sub-cellular level to dose-response effects on larger scales. In this symposium we will learn about the current state of the art of computational approaches estimating radiation damage at the cellular and sub-cellular scale. How can understanding the physics interactions at the DNA level be used to predict biological outcome? We will discuss if and how such calculations are relevant to advance our understanding of radiation damage and its repair, or, if the underlying biological processes are too complex for a mechanistic approach. Can computer simulations be used to guide future biological research? We will debate the feasibility of explaining biology from a physicists’ perspective. Learning Objectives: Understand the potential applications and limitations of computational methods for dose-response modeling at the molecular, cellular and tissue levels Learn about mechanism of action underlying the induction, repair and biological processing of damage to DNA and other constituents Understand how effects and processes at one biological scale impact on biological processes and outcomes on other scales J. Schuemann, NCI/NIH grantsS. McMahon, Funding: European Commission FP7 (grant EC FP7 MC-IOF-623630)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Design of an Out-Of-Pile Experimental Facility to Demonstrate the Feasibility of In Situ Thermal Conductivity Measurements of Nuclear Fuels Under Irradiation

There is substantial merit in quantifying nuclear fuel performance under irradiation. At Oak Ridge National Laboratory (ORNL), the MiniFuel irradiation platform has become the primary test vehicle for conducting separate-effects fuel performance irradiation experiments. The MiniFuel experiment is a passively controlled capsule design deployed in the High Flux Isotope Reactor (HFIR) through which fuel performance data is collected post-irradiation. Separate effects fuels irradiation capabilities are being expanded at ORNL by developing instrumented capsule designs that aim to capture fuel performance phenomena in-situ. One such capsule will specifically target fuel specimen thermal conductivity changes as a function of fuel burnup. Due to the complexity of making this measurement on nuclear fuel in-pile, this paper describes the necessary out-of-pile testing conducted on the thermal conductivity capsule (TCC) design. The measurement is ascertained via a thermopile system with heat transferred unidirectionally through a surrogate fuel specimen sandwiched between two conductive materials. The capsules investigated in this study are representative of the in-pile design, with the primary departure from irradiation conditions being the distribution of heat generation within the capsule. In the out-of-pile experiment, an external heater was used to drive heat through the conductive slug materials and into the specimen. This paper expounds the design of the out-of-pile experimental system and the thermal conductivity measurement technique. Predictive models used to determine the sensitivity of the measurement to variables governing thermal contact conductance between the specimen and slug materials and to predict experimental results are also described. Data from the out-of-pile experiment will be used to validate the readiness of the design for insertion into HFIR for irradiation.

Parker, Trevor [ORNL]↗
DOE OSTI2022

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Elucidating the Radiation-Induced Redox Chemistry of Plutonium Under Used Nuclear Fuel Reprocessing Conditions

Plutonium plays a critical role in the development of sustainable nuclear fuel cycles, and yet, our fundamental understanding of this element’s inherent radiation-induced redox chemistry and associated impacts on nuclear fuel cycle technologies is limited. Unanticipated changes in oxidation state distribution can influence the speciation and transport of plutonium in a given process. Control of these parameters is especially important for used nuclear fuel reprocessing technologies, wherein the separation and recovery of plutonium is typically achieved by the selective formation, maintenance, and complexation of specific oxidation states. Furthermore, plutonium’s inherent radiation-induced redox chemistry has the capacity to influence the radiolytic behavior of its complexes, the longevity of which are critical in the design of efficient and cost-effective advanced reprocessing technologies. These radiation-induced processes are unavoidable under fuel cycle conditions owing to the inherency of ionizing radiation fields to the decay of plutonium’s isotopes and to the various other radioisotopes generated by nuclear fission and neutron-capture process and the subsequent radioactive decay of their products. As such, mechanistically understanding the response of plutonium’s multiple oxidation states to multi-component ionizing radiation fields is essential for predicting the behavior of this critical element under used nuclear fuel reprocessing conditions. Here, through a combination of time-resolved (electron pulse) and steady-state (alpha and gamma) irradiation experiments complemented by quantitative, multiscale modeling calculations, we present advances in our understanding of radiation-induced plutonium redox chemistry!

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Integrated system modeling analysis of a cryogenic multi-cell deflecting-mode cavity resonator

A deflecting mode cavity is the integral element for six-dimensional phase-space beam control in bunch compressors and emittance transformers at high energy beam test facilities. RF performance of a high-Q device is, however, highly sensitive to operational conditions, in particular in a cryo-cooling environment. Using analytic calculations and RF simulations, we examined cavity parameters and deflecting characteristics of TM{sub 110,π} mode of a 5 cell resonator in a liquid nitrogen cryostat, which has long been used at the Fermilab A0 Photoinjector (A0PI). The sensitivity analysis indicated that the cavity could lose 30%–40% of deflecting force due to defective input power coupling accompanying non-uniform field distribution across the cells with 40 ∼ 50 MeV electron beam and 70–80 kW klystron power. Vacuum-cryomodules of the 5 cell cavity are planned to be installed at the Fermilab Advanced Superconducting Test Accelerator facility. Comprehensive modeling analysis integrated with multi-physics simulation tools showed that RF loading of 1 ms can cause a ∼5 K maximum temperature increase, corresponding to a ∼4.3 μm/ms deformation and a 1.32 MHz/K maximum frequency shift. The integrated system modeling analysis will improve design process of a high-Q cavity with more accurate prediction of cryogenic RF performance under a high power pulse operation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Quantifications of growth rate and particle size of magnesium silicate hydrate cement using small‐angle X‐ray scattering

The focus of recent research on cement binders has shifted to alternatives such as magnesium silicate hydrate (M-S-H) to mitigate the environmental impact of conventional concrete production. However, a detailed understanding of the morphology and growth rate of M-S-H, which informs its applicability as a cementing agent, is lacking. These properties are important as they influence cement hydration and consequently the development of macroscopic properties. Here, in this study, the particle size distribution and growth rate of M-S-H synthesized by mixing aqueous stock solutions were analyzed using synchrotron-based small-angle X-ray scattering. The reaction solutions contained dissolved magnesium and silicon with concentrations ranging from 2 to 500 mM, corresponding to saturation indices with respect to M-S-H endmembers, M 0.75 -S-H and M 1.50 -S-H, between 5.8 and 10.3. The precipitation rates of M-S-H align with affinity-based kinetic models and can be described by the following equations: For M 0.75 -S-H: 𝑟 M 0.75 −S−H =10 −9.54±0.70 ⁢(Ω$^{0.59±0.08}_{M_{0.75}⁢\textrm{SH}}$−1), and for M 1.50 -S-H: 𝑟 M 1.50 −S−H =10 −7.56±0.37 ⁢(Ω$^{0.34±0.04}_{M_{1.50}⁢\textrm{SH}}$−1), where 𝑟 is the rate in mol/L/s, and Ω represents the saturation ratio with respect to M-S-H. Moreover, we observed that changing the [Mg]/[Si] ratio affects particle size, a trend that cannot be fully explained by changes in the saturation index alone. These insights are critical for improving the predictability and control of cementation in concrete that utilizes M-S-H-based binders.

M-S-H↗
DOE OSTI2014

Comparative genomics and evolution of regulons of the LacI-family transcription factors

DNA-binding transcription factors (TFs) are essential components of transcriptional regulatory networks in bacteria. LacI-family TFs (LacI-TFs) are broadly distributed among certain lineages of bacteria. The majority of characterized LacI-TFs sense sugar effectors and regulate carbohydrate utilization genes. The comparative genomics approaches enable in silico identification of TF-binding sites and regulon reconstruction. To study the function and evolution of LacI-TFs, we performed genomics-based reconstruction and comparative analysis of their regulons. For over 1300 LacI-TFs from over 270 bacterial genomes, we predicted their cognate DNA-binding motifs and identified target genes. Using the genome context and metabolic subsystem analyses of reconstructed regulons, we tentatively assigned functional roles and predicted candidate effectors for 78 and 67% of the analyzed LacI-TFs, respectively. Nearly 90% of the studied LacI-TFs are local regulators of sugar utilization pathways, whereas the remaining 125 global regulators control large and diverse sets of metabolic genes. The global LacI-TFs include the previously known regulators CcpA in Firmicutes, FruR in Enterobacteria, and PurR in Gammaproteobacteria, as well as the three novel regulators—GluR, GapR, and PckR—that are predicted to control the central carbohydrate metabolism in three lineages of Alphaproteobacteria. Phylogenetic analysis of regulators combined with the reconstructed regulons provides a model of evolutionary diversification of the LacI protein family. The obtained genomic collection of in silico reconstructed LacI-TF regulons in bacteria is available in the RegPrecise database (http://regprecise.lbl.gov). It provides a framework for future structural and functional classification of the LacI protein family and identification of molecular determinants of the DNA and ligand specificity. The inferred regulons can be also used for functional gene annotation and reconstruction of sugar catabolic networks in diverse bacterial lineages.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Separate-Effects Validation Experiments for Models of Fracture in Ceramic Nuclear Fuel

A number of advances have recently been made in modeling fracture initiation and propagation in ceramic light water reactor (LWR) fuel. These include the use of techniques such as peridynamics, the extended finite element method (XFEM), the discrete element method (DEM), and phase field fracture. Improved predictions of fracture in LWR fuel can benefit fuel performance simulations in a number of ways, with important implications for radial relocation, pellet-cladding mechanical interactions, and fragmentation during accident conditions. Although these modeling approaches appear promising, to-date most of their applications have been to fresh fuel. A major limitation in applying these approaches with confidence to a wider variety of conditions, such as fragmentation of high-burnup fuel, is that there is very limited data that can be directly used to validate the results of these fracture models. While fracture patterns are available from cross-sections of numerous irradiated fuel rods, the fuel rods have typically been subjected to complex irradiation histories. To validate models of fracture propagation, focused experiments that give snapshots of the initiation and propagation of fracture under well characterized prototypical reactor conditions are needed. Under a multi-institutional U.S. Department of Energy Integrated Research Project (IRP), three separate-effects validation experiments are underway to provide information on the propagation of cracks in LWR fuel under controlled conditions. The primary challenge in developing such an experiment is in providing conditions that result in a radial temperature distribution that is prototypical of LWR in-reactor conditions. Two out-of-reactor experiments are being planned that employ very different approaches to achieve such a temperature distribution. The first of these achieves this by employing a combination of inductive and resistive electrical heating to generate volumetric heating. The second experiment involves slowly heating a pellet in a metal rod to a high temperature, and then quenching that rod in a cold bath. The quenching process generates a temperature profile approximating that in an LWR fuel rod for a moment, but has a very different temperature history from the LWR. Finally, an experiment is being planned for the Idaho National Laboratory TREAT reactor that will use nuclear fission as a heat source, and which will employ a metal heat sink to replicate the thermal conditions provided by coolant in an LWR. These experiments all have their advantages and drawbacks, with varying levels of ability to reproduce actual LWR conditions and varying ability to use instruments to monitor quantities of interest over the course of the experiment that can be used to validate fracture models. For that reason, these different experimental approaches are complimentary to each other, and the combined results are expected to significantly improve confidence in the ability of computational models to simulate fracture in ceramic nuclear fuel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Improving Trustworthiness of Data-Driven Power Grid Contingency Analysis With Bayesian Residual Graph Neural Networks

The evolving energy landscape requires novel tools to efficiently perform contingency analysis and reliability assessment of power grids, potentially in real-time. The high computational cost of traditional power flow solvers limits their applicability in practice. Machine learning (ML) surrogates such as deep neural networks (NNs) accelerate power flow solvers computations, enabling high-order contingency analysis and real-time decision-making by learning highly nonlinear functions and integrating grid topology via graph architectures. However, (graph) NNs lack predictive power away from training data and do not provide predictive confidence estimates. Here, we present a Bayesian residual graph NN that integrates knowledge from low-fidelity data via residual training and embeds granular quantification of uncertainties, improving trustworthiness critical for high-consequence decision-making. Applying Bayesian concepts to NNs is challenging due to the high-dimensionality of both the parameter space, complicating derivation of a meaningful prior, and the output space in large grid systems, requiring enhanced techniques to assess the predicted high-dimensional uncertainties. Our contributions include: (1) Deriving a prior for fully connected and graph NNs that leverages low-fidelity data to guide mean predictions and appropriately control prior predictive uncertainty. (2) Integrating this prior within an ensembling with anchoring scheme for efficient approximate posterior inference. (3) Deriving enhanced metrics to assess accuracy of both the mean and uncertainty predictions in high dimensions, appropriately accounting for correlations propagated through graph layers. The resulting Bayesian residual graph NN is tested on a contingency analysis task for 14-bus and 118-bus grids.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Predictive Large-Eddy Simulation of Supercritical-Pressure Reactive Flows in the Cold Ignition Regime

This report describes a one-year study performed under DOE sponsorship, continuing the investigation of high-pressure turbulent reactive flows. The interest is in the effect of the chemical species distribution in high-pressure turbulent flows in the presence of strong temperature gradients as would occur during reactions in realistic flows. The prime example of such flows are boundary layers in which the wall is at a lower temperature than that of the fluid, as would be the case in Diesel engines. Previous DOE BES work in the program further highlighted the importance of the boundary layer as a configuration for fundamental studies: soot formation in boundary layers is still a problem poorly understood and depends on the availability of particular chemical species at that location, and when boundary layer Large Eddy Simulations results were compared with measurements, the agreement was unfavorable, showing that this important ‘unit’ problem is not well understood. To understand this unit problem, modeling and Direct Numerical Simulations of this unit problem were conducted for the simplest possible multispecies system, that is a binary-species system. The results discovered a new phenomenon, that is, Soret effect induced uphill diffusion. The far-reaching implication is that through an imposed wall/fluid temperature difference it is possible to control the distribution of the species in the boundary layer. Because the results have been documented in one paper published in the refereed literature, and also in conference papers, this final report is in the form of an Executive Summary succinctly describing the results and putting them in perspective with respect to existing information. The refereed and conference papers published are individually listed as Appendices and attached to this report. One manuscript is still in progress and is thus not listed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Low- and intermediate-temperature oxidation of dimethyl ether up to 100 atm in a supercritical pressure jet-stirred reactor

Understanding the low- and intermediate-temperature oxidation chemistry of oxygenated fuels like dimethyl ether (DME) at high pressure is paramount to the development of advanced engines with low carbon emissions. The supercritical pressure jet-stirred reactor (SP-JSR) recently developed at Princeton provides a new platform for conducting kinetic studies at low and intermediate temperatures at extremely high pressures with a uniform temperature distribution and a short flow residence time. This paper uses the SP-JSR to investigate DME oxidation at equivalence ratios of 0.175, 1.0, and 1.72, for pressures of 10 and 100 atm, and temperatures ranging from 400 to 900 K. The results demonstrate weakened NTC behavior at 100 atm relative to 10 atm due to increased flux through QOOH + O 2 = O 2 QOOH relative to QOOH = 2 CH 2 O + OH at 100 atm. Furthermore, the intermediate-temperature oxidation window is shifted to lower temperatures at 100 atm. The experimental data are compared with several chemical kinetic models from the literature. The existing models are seen to agree quite well with the experimental data at 10 atm. However, the models fail to properly capture the NTC behavior at 100 atm. Reaction pathway analyses indicate that both the low- and intermediate-temperature chemistries are controlled by RO 2 consumption pathways. The reaction rates for several of the important reactions, such as DME + OH = CH 3 OCH 2 + H 2 O, H 2 O 2 (+M) = 2 OH (+M), and 2 HO 2 = 2 OH + O 2 are updated in this work. Furthermore, the updated model improves the predictability for all key species compared to the original model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Short-Term Forecasting Across a Network for the Autonomous Wind Farm: Preprint

In an autonomous wind farm, turbines will use information from nearby turbines to achieve wind farm-level objectives such as optimizing the overall performance of a wind farm, ensuring resiliency when other sensors fail, and adapting to changing local conditions. In this paper, the wind farm can be modeled as a network within which turbines (nodes) share information across designated communication channels, with a focus on turbines at the outside of the wind farm capturing local effects and sharing that information with downstream turbines. Understanding of varied inflow conditions can be especially important in complex terrain. This information can be used to monitor turbines, self-organize turbines into groups, and predict the power performance of a wind farm. In particular, this paper describes an autonomous wind farm that incorporates information from local sensors in real time to predict wind speed and wind direction at each turbine over a short-term horizon. Results indicate that the estimate of wind direction can be used to improve the knowledge of the wind speed and direction over the persistence method on a 10- 15-minute time horizon. These short-term forecasts can also be used to facilitate advanced control methods such as feedforward control within a wind farm.

17 WIND ENERGY↗
DOE OSTI2019

Short-Term Forecasting Across a Network for the Autonomous Wind Farm

In an autonomous wind farm, turbines will use information from nearby turbines to achieve wind farm-level objectives such as optimizing the overall performance of a wind farm, ensuring resiliency when other sensors fail, and adapting to changing local conditions. In this paper, the wind farm can be modeled as a network within which turbines (nodes) share information across designated communication channels, with a focus on turbines at the outside of the wind farm capturing local effects and sharing that information with downstream turbines. Understanding of varied inflow conditions can be especially important in complex terrain. This information can be used to monitor turbines, self-organize turbines into groups, and predict the power performance of a wind farm. In particular, this paper describes an autonomous wind farm that incorporates information from local sensors in real time to predict wind speed and wind direction at each turbine over a short-term horizon. Results indicate that the estimate of wind direction can be used to improve the knowledge of the wind speed and direction over the persistence method on a 10-15-minute time horizon. These short-term forecasts can also be used to facilitate advanced control methods such as feedforward control within a wind farm.

17 WIND ENERGY↗
DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗
DOE OSTI2017

Plasma-Surface Interaction Studies on DIII-D and Their Implications for Next-Step Fusion Experiments

Unique diagnostic and access features of the DIII-D tokamak, including a sample exposure system, have been used to carry out controlled and well-diagnosed plasma-surface interactions (PSI) experiments. An important contribution of the experiments has been the ability to link a given plasma exposure condition to a measured response of the plasma-facing surface and to thus understand the interaction. This has allowed for benchmarking certain aspects of erosion models, particularly near-surface particle transport. DIII-D has empirically quantified some of the PSI effects that will limit the operation availability and lifetime of future fusion devices, namely, net erosion limiting divertor plate lifetime and hydrogenic fuel retention in deposit layers. Cold divertor plasmas obtained with detachment can suppress net carbon divertor erosion, but many low-temperature divertor PSI phenomena remain poorly understood: nondivertor erosion sources, long-range particle transport, global erosion/deposition patterns, the enhancement of carbon erosion with neon impurity seeding, the sputtered carbon velocity distribution, and the apparent suppression of carbon chemical erosion in detachment. Long-term particle and energy fluences have reduced the chemical erosion yield of lower-divertor tiles. Finally, plasma-caused modification of a material’s erosion properties, including material mixing, will occur quickly and be important in long-pulse fusion devices, making prediction of PSI difficult in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

NRAP-Open-IAM: NRAP Open Source Integrated Assessment Model

Note: This is the last version (a2.6.1) of NRAP-Open-IAM released during NRAP Phase II in 2022. The latest version of NRAP-Open-IAM is available here: https://edx.netl.doe.gov/dataset/phase-iii-nrap-open-iam NRAP-Open-IAM is an open-source software product that enables quantification of containment effectiveness and leakage risk at storage sites in the context of system uncertainties and variability. NRAP-Open-IAM represents the next-generation in a line of systems-based computational models developed for quantitative geological carbon storage (GCS) risk assessment. The model comprises a set of reduced-order and analytical models of various components of the GCS system, potential leakage pathways, receptors of concern including impact to groundwater resources and the atmosphere, a framework to support stochastic simulation, time stepping, uncertainty quantification, other analytical functionality for scenario and risk-performance evaluation, and a basic graphical user interface to support scenario development, data input simulation definition, and basic post-processing and results display. As the NRAP Open-IAM functionality continues to evolve, we continue to add to its capability to develop quantitative, probabilistic, and time-dependent profiles of the evolution of risk at a GCS site and evaluate the influence of uncertain parameters on uncertainty in predicted risk. It can be used to quantify the dynamics of reservoir saturation plume and pressure-affected area, for evaluation of the area of potential groundwater impact (i.e., Area of Review) and monitoring requirements to support cost and regulatory analysis, and for consideration of different post-injection site care and closure scenarios. This submission contains the current version of NRAP-Open-IAM available for evaluation and testing. To use the NRAP-Open-IAM, download the source code (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/4c24a3da-3b40-4ffe-9892-c807ae9f8760) then open the NRAP-Open-IAM user's guide (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669) to read more about the tool. Installation instructions for Windows, Mac, and Linux can be found in the "installers" folder of the extracted NRAP-Open-IAM folder and describe setup of environment (e.g., Python libraries) needed for proper work of the tool. Test of installation can be done by running "python openiam_setup_tests.py" in the "setup" folder. The installation test also runs a test suite to see if the NRAP-Open-IAM has been installed correctly. To run the test suite separately, run "python iam_test.py" in the "test" folder. User's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669 Developer's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/3bc6ee7d-609d-4eb6-80ba-fa6130ee0313 Reservoir simulation data used in some examples distributed with NRAP-Open-IAM: - Kimberlina: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/eb62cece-61b2-4037-9b6d-32407dde2ab8 - Kimberlina (compartmentalized): https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/366f9530-3b32-4b84-affe-ab2df1d9a8b5 - FutureGen 2.0: https://edx.netl.doe.gov/dataset/futuregen-2-0-1008-simulation-reservoir-lookup-table NRAP-Open-IAM GitLab repository: https://gitlab.com/NRAP/OpenIAM Related publications: - Bacon, D., Yonkofski, C., Brown, C., Demirkanli, D. and Whiting, J., 2019. Risk-based post injection site care and monitoring for commercial-scale carbon storage: Reevaluation of the FutureGen 2.0 site using NRAP-Open-IAM and DREAM. International Journal of Greenhouse Gas Control 90: 102784. - Bacon, D. Demirkanli, D., and White, S., 2020. Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site. International Journal of Greenhouse Gas Control 102: 103153. - Harp, D., Oldenburg, C., and Pawar, R., 2019. A metric for evaluating conformance robustness during geologic CO2 sequestration operations. International Journal of Greenhouse Gas Control 85: 100-108. - Lackey, G., Vasylkivska, V., Huerta, N., King, S., and Dilmore, R., 2019. Managing well leakage risks at a geologic carbon storage site with many wells, International Journal of Greenhouse Gas Control, 88 :182-194. - Vasylkivska, V., Dilmore, R., Lackey, G., Zhang, Y., King, S., Bacon, D., Chen, B., Mansoor, K., and Harp, D., 2021. NRAP-Open-IAM: A flexible open-source integrated assessment model for geologic carbon storage risk assessment and management, Environmental Modelling & Software, 143: 105114. Presentations: - Chen, B., Harp, D., and Pawar, R., A data assimilation approach (ES-MDA) coupling with NRAP-Open-IAM for quantifying uncertainty reduction in geological CO2 sequestration. AGUFM 2019: T44A-02. - Chen, B., and Harp, D., Improving risk analysis precision for geologic CO2 sequestration by quantifying the uncertainty reduction before and after acquiring monitoring data. 14th Greenhouse Gas Control Technologies Conference, Melbourne, Australia, 2018, pp. 21-26. - Harp, D., National Risk Assessment Partnership Task 2: Containment Assurance. No. LA-UR-19-28654, Los Alamos National Laboratory (LANL), Los Alamos, NM (United States), 2019. - Vasylkivska, V., King, S., Bacon, D., Harp, D., Chen, B., Mansoor, K., Onishi, T., Yang, Y., Zhang, Y., and Keating, E., NRAP-Open-IAM: An open-source integrated assessment model, poster, Mastering the Subsurface Through Technology Innovation, Partnerships and Collaboration: Carbon Storage and Oil and Natural Gas Technologies Review Meeting, Pittsburgh, PA, August 13-16, 2018. - Vasylkivska, V., Lackey, G., King, S., Wentworth, A., Huerta, N., Creason, C., DiGiulio, J., Yang, Y., and Dilmore, R., Long-term risk analysis of a geologic CO2 storage project during the post-injection period, SIAM Conference on Computational Science and Engineering, Spokane, WA, February 25-March 1, 2019. - Vasylkivska, V., Overview of the NRAP-Open-IAM tool for carbon storage (beta release), 2019 Annual NRAP Tool Users Meeting, Pittsburgh, PA, August 27, 2019. - Vasylkivska, V., Bacon, D., Chen, B., Dilmore, R., Harp, D., King, S., Lackey, G., Lindner, E., Liu, G., Mansoor, K. and Zhang, Y., NRAP-Open-IAM: A new, open-source code for integrated assessment of geologic carbon storage containment effectiveness and leakage risk, poster, American Geophysical Union Fall Meeting 2020 (virtual meeting), December 2020. - Vasylkivska, V., NRAP open-source integrated assessment model and relevant application, oral presentation, NRAP workshop "NRAP Tools for Geologic Carbon Storage Risk-Based Decision Making" held in conjunction with Groundwater Protection Council (GWPC) 2021 Annual Forum (virtual meeting), Salt Lake City, UT, September 2021. - Vasylkivska, V., NRAP-Open-IAM: open-source integrated assessment model, digital poster/demonstration, software demonstration session, 2022 Carbon Management Project Review Meeting, August 16, 2022

AoR↗
DOE OSTI2024

Suppression of diffraction in deep-inelastic scattering on nuclei and dynamical mechanism of leading twist nuclear shadowing

Abstract Using the leading twist approach (LTA) to nuclear shadowing, we calculate the ratios of diffractive and usual parton distributions for a heavy nucleus (Pb) and the proton,$$ {R}_{A/p}=\left({f}_{i/A}^{D(3)}/{f}_{i/A}\right)/\left({f}_{i/p}^{D(3)}/{f}_{i/p}\right) $$ R A / p = f i / A D 3 / f i / A / f i / p D 3 / f i / p , for coherent and summed (coherent plus quasi-elastic) nuclear deep-inelastic scattering. We find thatR A/p ≈ 0.5 − 1 for quarks as well as for the ratio of the diffractive and total cross sections$$ {\left[\left({d\sigma}_{\textrm{diff}}/{d M}_X^2\right)/{\sigma}_{\textrm{tot}}\right]}_{eA}/{\left[\left({d\sigma}_{\textrm{diff}}/{d M}_X^2\right)/{\sigma}_{\textrm{tot}}\right]}_{ep} $$ dσ diff / dM X 2 / σ tot eA / dσ diff / dM X 2 / σ tot ep andR A/p ≈ 0.5 − 1.3 for gluons in a broad range ofx, including the kinematics of the Electron-Ion Collider, which reaffirms the difference from the nuclear enhancement ofR A/p predicted in the gluon saturation framework. We demonstrate that the magnitude ofR A/p is controlled by the cross section of the interaction of hadronic fluctuations of the virtual photon with target nucleons, which explains an enhancement ofR A/p in the color dipole model and its suppression in LTA. We argue that the black disk limit corresponds toR A/p = 1 and$$ {R}_{A/p}^{\textrm{coh}} $$ R A / p coh = 0.86 for the summed and coherent scattering, respectively. Relying on an intuitive definition of the saturation scale, we show that the ratio of the saturation scales of a heavy nucleus and proton$$ {Q}_{sA}^2(b)/{Q}_{sp}^2(b)\approx 1 $$ Q sA 2 b / Q sp 2 b ≈ 1 at small impact parametersbdue to the strong leading twist nuclear shadowing and diluteness of the nuclear density.

Physics↗
DOE OSTI2024

Model-Based Detection of Coordinated Attacks (DCA) in Distribution Systems

The fast-paced growth in digitization of smart grid components enhances system observability and remote-control capabilities through efficient communication. However, enhanced connectivity results in heightened system vulnerability towards cybersecurity risks in the cyber-physical power system. Coordinated cyber-attacks (CCA), when undetected, lead to system-wide impact in terms of large disturbances or widespread outages. Detecting CCA in the cyber layer is critical to thwart cyber-attacks in real-time before the attack impacts the physical system. The challenge of locating CCA stems from the complex grid dynamics, making it difficult to distinguish between normal operational variations and cyber-attack impact. CCA often employs multiple attack vectors targeting geographically distributed components, further complicating CCA identification. Existing research in intrusion detection is primarily focused on the transmission network and limited to detecting individual attacks. In this paper, a novel proactive DCA strategy is proposed for early detection of CCA by establishing correlations among distinct attack events through model-based reinforcement learning that utilizes abductive reasoning to conclude the attacker goal. The solution includes understanding the system model, learning the system dynamics, and correlating individual cyber-attacks to extract the attacker’s objective. The developed learning algorithm identifies the most probable attack path to reach the attacker’s objective by predicting the next attack steps. A DNP3-based cyber-physical co-simulation testbed is developed to test the proposed algorithm using the IEEE 13-node test feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗