Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference models”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Generating Co-expression Networks for Three Cyanobacteria: Synechococcus sp. PCC 7942, Synechococcus sp. PCC 7002, Synechocystis sp. PCC 6803

Cyanobacteria are photosynthetic organisms capable of high growth rate and represent a promising bioplatform for harnessing the sun’s energy to make biofuel. Additionally, the process of photosynthesis absorbs CO2 from the environment. Understanding the metabolic processes involved in photosynthesis could lead to solutions to the recent rise of CO2 concentration in Earth’s atmosphere and the associated climate change. More research on the transcriptional regulation of these cells is needed to learn how to harness the untapped potential of cyanobacteria for these applications. Transcriptional analysis via RNA-seq provides an understanding of how gene expression changes at the mRNA level under diverse growing conditions. I systematically collected and analyzed RNA-Seq data obtained under a variety of conditions and available on the NCBI database for three cyanobacteria model organisms: Synechococcus elongatus sp. PCC 7942, Synechococcus sp. PCC 7002, and Synechocystis sp. PCC 6803. For each organism, the data was mapped to a reference genome to characterize the RNA expression profile. Samples were checked for quality based on the number of reads and the correlation of the expression profile between labeled replicates. All samples were transformed into transcripts per million reads, followed by a log transformation to account for the wide range of sample sizes. Gene co-expression networks were generated and analyzed for each species using Cytoscape. These networks provide a base level of gene expression for each species. The network topology and high-betweeness nodes of these networks need to be analyzed further to provide insight on potential ways to harness cyanobacteria genetics. Additionally, these datasets can be used together to form a core genome network analysis- one that includes only the genes that are homologous between the three species. This project has prepared the way for a more in-depth study on photosynthetic microbes on a genetic level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The first ensemble of kilometer-scale simulations of a hydrological year over the third pole

An accurate understanding of the current and future water cycle over the Third Pole is of great societal importance, given the role this region plays as a water tower for densely populated areas downstream. An emerging and promising approach for skillful climate assessments over regions of complex terrain is kilometer-scale climate modeling. As a foundational step towards such simulations over the Third Pole, we present a multi-model and multi-physics ensemble of kilometer-scale regional simulations for the hydrological year of October 2019 to September 2020. The ensemble consists of 13 simulations performed by an international consortium of 10 research groups, configured with a horizontal grid spacing ranging from 2.2 to 4 km covering all of the Third Pole region. These simulations are driven by ERA5 and are part of a Coordinated Regional Climate Downscaling EXperiment Flagship Pilot Study on Convection-Permitting Third Pole. The simulations are compared against available gridded and in-situ observations and remote-sensing data, to assess the performance and spread of the model ensemble compared to the driving reanalysis during the cold and warm seasons. Although ensemble evaluation is hindered by large differences between the gridded precipitation datasets used as a reference over this region, we show that the ensemble improves on many warm-season precipitation metrics compared with ERA5, including most wet-day and hour statistics, and also adds value in the representation of wet spells in both seasons. As such, the ensemble will provide an invaluable resource for future improvements in the process understanding of the hydroclimate of this remote but important region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Model-Based Approaches to Generate Knowledge from Data in a Plant Reliability Context

One challenge that nuclear power plant system engineers are facing is continuous generation of an extremely large amount of equipment reliability (ER) data. These data elements come in textual (e.g., condition reports) and numeric (e.g., generated by monitoring systems) forms. They provide system engineers with valuable insights and information by discovering anomalous behaviors or degradation trends, identifying possible causes behind such behaviors and trends, and predicting their direct consequences. This paper directly targets the knowledge generation from ER data by putting “data into context.” We employ model-based system engineering (MBSE) of systems and assets to represent and capture their architecture and functional (i.e., cause-effect) relations. ER data elements are processed by first identifying which of the developed MBSE elements they are referring to. This task is harder for textual data since the information contained in issue or maintenance reports needs to be “understood” by a computational tool. We called this process “knowledge extraction” since our methods extract knowledge from textual data. Last, once numeric and textual ER data elements have been processed and “understood,” we discover possible cause-effect relations among them. This is performed by observing whether a logical connection through the MBSE models exists, and if there is a temporal relationship among them. The logic and temporal are the two main ingredients to perform “machine reasoning” from ER data.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Predicting Damages to Remainder Parcels in Right-of-Way Acquisitions for Expanding Transportation Infrastructure: Using a Truncated Finite-Mixture Model

Right-of-way acquisition is a critical component of transportation infrastructure development. Transportation infrastructure projects cannot proceed without proper right-of-way acquisition or may face significant delays. State Departments of Transportation frequently acquire parcels of land for roadway expansion projects. A majority of these acquisitions can be partial takings, referring to a portion of a parcel that is acquired. The remainder of the property usually suffers economic changes due to the partial acquisition, which can be calculated as damage percentages. The damage percentage represents the extent to which the remaining land or property value has been diminished due to the acquisition. It reflects the remaining property value percentage that may have been lost or compromised due to the acquisition. Here, this study aims to provide a robust model to estimate damage percentages to the remainder parcels that may help state Departments of Transportation appraisers make early predictions about the damages in cases involving partial takings. The research uses 509 appraisal reports from the Tennessee Department of Transportation to identify the key parcel attributes that influence the percentage of damages. Three regression models are developed: a linear regression model, a finite-mixture model (FMM), and a truncated FMM with two latent classes. The modeling results show that the truncated FMM with two classes outperforms the other models. To validate the models, actual sales data is collected and analyzed for 59 properties, and the results suggest that the model predictions are fairly accurate. A predictive tool is developed based on the models to help appraisers anticipate right-of-way damages under different scenarios and can provide early predictions about the damages.

42 ENGINEERING↗
DOE OSTI2021

Terry Turbopump Expanded Operating Band Modeling and Simulation Efforts in Fiscal Year 2021 Extended Period of Performance (Final Report)

This report documents the progress made under the Terry Turbine Expanded Operating Band (TTEXOB) program's modeling and simulation (MODSIM) initiative at Sandia National Laboratories (SNL ). It describes the US Federal Fiscal Year 2021 (FY21) extended period-of-performance MODSIM work completed since the closure of FY20 with due reference to the Texas A&M University (TAMU) hybrid milestone 5/6 experimental program. This work, which falls under Milestone 7 of the program, provides a counterpart to the various experiments. The overall TTEXOB program and its milestone-based approach are described in the program's Summary Plan. Details of the individual milestone test plans can be found in the corresponding detailed test plan, e.g. the Milestone 3 and 4 Detailed Test Plan. SNL MODISM is conducted alongside experiments performed at TAMU, and SNL technical staff regularly consults with TAMU on the experimental program. In FY21, MELCOR code models and capabilities were exercised in two different contexts: experimental comparisons to the TAMU ZS-1 and GS-2, and stand-alone analyses of a station black-out (SBO) scenario in a generic boiling water reactor (BWR). Code to experiment comparisons met with fair success when turbine losses were well characterized as for the ZS-1 turbine. Both deterministic and Bayesian calibration processes were used to find a recommended turbine torque multiplier for ZS-1 type turbines. This process could be repeated for GS-2 type turbines if GS-2 losses were better understood. Stand-alone generic BWR SBO calculations revealed that three different modes of self-regulating turbopump behavior may be observed depending on certain modeling parameters and choices having to do with turbine nozzles. Aspects of this predicted behavior may have been observed in TAMU GS-2 experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

ADEPT: A Pedagogical Framework for Integrating Agentic AI with Deterministic Scientific Workflows

The integration of Large Language Models (LLMs) into scientific research promises to accelerate discovery, yet a significant gap remains between the dynamic reasoning of Artificial Intelligence (AI) agents and the static, deterministic nature of canonical scientific workflows. This paper introduces ADEPT (Agentic Discovery and Exploration Platform for Tools), a reference architecture and pedagogical framework explicitly designed to bridge this gap. ADEPT's primary mission is to provide a transparent, "glass-box" environment where researchers and engineers can learn to effectively wrap established scientific software (e.g., BLAST, Nextflow pipelines) and compose it into reliable, agent-driven workflows. We describe its modular, multi-server architecture, which leverages the Model Context Protocol (MCP) for tool serving, LangGraph for robust agentic orchestration, and a secure nsjail-based sandbox for safe code execution. By prioritizing architectural clarity, safety, and modularity, ADEPT serves as an extensible blueprint for building trustworthy AI-augmented systems and fosters the collaborative development necessary to responsibly employ agentic AI for science. We provide practical examples of how to adapt and extend this framework, highlighting its utility in workforce development and AI-readiness capabilities across research and development projects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Comparative Study of Nonlinear Black-Box Modeling for Power Electronics Converters

With the increasing penetration level of renewable sources and power electronics loads in modern power systems, accurate and computationally efficient models are needed. Black-box model (BBM) could be a useful method in such systems. However, not very extensive research efforts have been made for power electronics BBM so far, and existing works mostly focus on steady-state operation, neglecting the important transient behaviors such as load transients, voltage transients, and faults. This paper presents a comparative study of three commonly used nonlinear BBM approaches for transient behaviors of power electronics converters. Comparison methods are proposed, and the evaluations are conducted under different transients using a grid-connected single-phase photovoltaic inverter. The findings of this study provide valuable references for further feasibility investigations on implementing BBMs in large-scale power electronics-rich power systems.

Qiao, Liang↗
DOE OSTI2024

Modeling aerosol transmission spectra from n(λ) and k(λ) infrared optical constants measurements of organic liquids and solids

The effects of light scattering and refraction play significantly different roles for aerosols than for bulk materials, making it challenging to identify aerosolized chemicals using traditional spectral methods or spectral reference libraries. Due to a potentially infinite number of particle morphologies, sizes, and compositions, constructing a database of laboratory-measured aerosol spectra is not a practical solution. Here, as an alternative approach, the measured n / k optical vectors of two example organic materials (diethyl phthalate and D-mannitol) are used in combination with particle absorption / scattering theory (Mie theory and FDTD) and the Beer-Lambert law to generate a series of synthetic infrared transmission / scattered light spectra. The synthetic spectra show significant differences versus simple slab transmission spectra, even for small changes in particle size (e.g., 5 vs. 10 µm) for both single particles and ensembles, potentially serving as useful reference data for aerosol sensing. For spherical single particles with diameters of 1 to 10 µm, FDTD simulations predict changes in the magnitudes of spectral shifts and the shapes of the peaks vs. particle size with only small deviations from Mie theory predictions, yet reliably capture the direction of the shifts. Typical spectral peak shifts in the longwave infrared correspond to Δλ ∼0.20 µm (∼34 cm -1 ) when compared to corresponding slab transmission spectra. Additionally, synthetic spectra generated from the n / k values derived using two different methods (KBr pellet transmission and single-angle reflectance) are compared using the Mie theory model.

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

Hybrid RANS-LES of the Atmospheric Boundary Layer for Wind Farm Simulations: Preprint

Wind farm simulations often do not accurately represent wake-atmospheric boundary layer (ABL) interactions, blade boundary layer (BL) dynamics, and turbine-turbine interactions. In this work, we use Active Model Split (AMS), a new hybrid Reynolds-Averaged Navier Stokes (RANS)-large eddy simulation (LES) model, which is well suited to capture these effects because the model can (i) accurately simulate the ABL with the Coriolis effect, (ii) is accurate in adverse pressure gradients such as those near wind turbine blades, and (iii) has sufficiently low computational cost to simulate multiple turbines while resolving the blade BL. For simplicity and consistency we develop AMS to be used throughout the domain rather than in a zonal method. We implement our work in the massively parallel flow solver, Nalu-Wind, so that our model can access the compute resources needed for blade-resolved simulations of multiple wind turbines. To accomplish these aims, we modify the baseline AMS by changing the RANS contribution to SST k - omega with a length scale limiter, adding the Coriolis effect, and developing an appropriate wall treatment. We show that AMS of the ABL with the Coriolis effect matches LES reference results better than those obtained with RANS. We describe our plans to add buoyancy effects and wind turbines to our AMS simulations.

atmospheric boundary layer↗
DOE OSTI2022

Modeling the non-Schmid crystallographic slip in MAX phases

We present a crystal plasticity constitutive relation for the description of experimentally observed non-Schmid crystallographic slip in a class of ternary carbides and nitrides commonly referred to as MAX phases. In the constitutive relation, we assume that the evolution of the slip system strength in MAX phases has two components – a classical component that depends on the Taylor cumulative shear strain and a non-Schmid component that depends on the stress normal to the slip plane. The non-Schmid crystal plasticity constitutive relation is then used to carry out finite element simulations of micropillar compression of single crystals of two MAX phases, Ti 2 AlC and Ti 3 AlC 2 . The finite element simulations not only quantitatively predict the stress – strain response of a wide range of crystallographic orientations of the micropillars but also rationalize the non-uniform deformation and the deformed shape of the micropillars observed in the experiments for the two materials. As a result, parametric studies are also carried out to quantify the role of the non-Schmid effect and understand the effects of key experimental parameters on the stress – strain response of the micropillars of the two MAX phases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Energy Requirements for Integration of Nuclear Reactors with Iron and Steel Plants

This report identifies energy needs of heavy energy users within the domestic iron and steel industry and suggests solutions for integrating nuclear energy. The iron and steel industry, composed of several types of plants which perform different processes with varied energy demands and vectors, is a heavy consumer of electric power and fossil fuels including coke and natural gas. Almost all major process temperatures exceed the temperatures of direct heat available from advanced reactors, and so electricity and hydrogen were considered instead. Reference units were adopted and estimated energy demands computed for the blast furnace (BF), direct reduced iron (DRI) unit, electric arc furnace (EAF), and reheat furnaces. By utilizing production capacity data from industry reports, the ranges of power demands were estimated, including for hydrogen production by high temperature steam electrolysis (HTSE). For the EAF and DRI unit, more detailed integration studies with thermodynamic modeling were also conducted and determined a possible solution with a specific reactor design and number of modules. Furthermore, because many unit processes are co-located, entire plants were considered by adding the energy demands of the unit processes to form four hypothetical reference plants. The range of power needs for the reference plants is compatible with multi-unit banks of microreactors at the low end, and would create a need for multiple larger-capacity SMRs at the high end (1 GWe plus 0.14 GWt). Although the overall power need at the high end is well-matched with one present-day large reactor offering (1.1 GWe), redundancy considerations may require a minimum of two reactors, potentially eliminating the single large reactor from consideration. Auxiliary or house loads would increase the reference plant estimates. Finally, the report provides total estimated energy needs under integration of all U.S. units of each process (BF, DRI, EAF, and reheat furnaces), representing a national potential for nuclear energy in the industry. U.S. iron and steel plants may be candidates for integration with nuclear reactors via electricity and hydrogen, and many sites have energy requirements that correspond well to the capacities of several advanced nuclear power designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A review on the sizing and selection of control valves for thermal hydraulics for reactor system applications

This study presents an overview regarding how to model/size, choose, and maintain control valves (CVs) for a reactor system thermal hydraulics experimentation—specifically simulating pipe break(s) and a loss of coolant accident (LOCA) analysis. In a thermal hydraulics test loop, CVs modulate fluid flow according to the control signal generated by the controller to control the process parameter. The selection of a CV in industrial applications involves many factors, including process diversity, safety, and reliability. Further, the key issues and challenges for performing proper mathematical modeling, installation, servicing, and calibration are also discussed for the benefit of operators and end-users. Additionally, this research includes case studies illustrating major plant-level accidents associated with the malfunctioning of valves. Therefore, the study will serve as a knowledge base or reference guide for young professionals and operators who wish to better understand process control applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗
DOE OSTI2022

A pair of early- and late-forming galaxy cluster samples: A novel way of studying halo assembly bias assisted by a constrained simulation

The halo assembly bias, a phenomenon referring to dependencies of the large-scale bias of a dark matter halo other than its mass, is a fundamental property of the standard cosmological model. First discovered in 2005 from the Millennium Run simulation, it has been proven very difficult to be detected observationally, with only a few convincing claims of detection so far. The main obstacle lies in finding an accurate proxy of the halo formation time. In this study, by utilizing a constrained simulation that can faithfully reproduce the observed structures larger than 2 Mpc in the local universe, for a sample of 634 massive clusters at z ≤ 0.12, we found their counterpart halos in the simulation and used the mass growth history of the matched halos to estimate the formation time of the observed clusters. This allowed us to construct a pair of early- and late-forming clusters, with a similar mass as measured via weak gravitational lensing, and large-scale biases differing at the ≈3σ level, suggestive of the signature of assembly bias, which is further corroborated by the properties of cluster galaxies, including the brightest cluster galaxy and the spatial distribution and number of member galaxies. Our study paves a way to further detect assembly bias based on cluster samples constructed purely on observed quantities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Coordinate transformation methodology for simulating quasistatic elastoplastic solids

Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypo-elastoplastic mechanical model, and develop a projection method to enforce quasi-static equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A Distributionally Robust Resilience Enhancement Strategy for Distribution Networks Considering Decision-Dependent Contingencies

When performing the resilience enhancement for distribution networks, there are two obstacles to reliably model the uncertain contingencies: 1) decision-dependent uncertainty (DDU) due to various line hardening decisions, and 2) distributional ambiguity due to limited outage information during extreme weather events (EWEs). Here, to address these two challenges, this paper develops scenario-wise decision-dependent ambiguity sets (SWDD-ASs), where the DDU and distributional ambiguity inherent in EWE-induced contingencies are simultaneously captured for each possible EWE scenario. Then, a two-stage tri-level decision-dependent distributionally robust resilient enhancement (DD-DRRE) model is formulated, whose outputs include the optimal line hardening, distributed generation (DG) allocation, and proactive network reconfiguration strategy under the worst-case distributions in SWDD-ASs. Subsequently, the DD-DRRE model is equivalently recast to a mixed-integer linear programming (MILP)-based master problem and multiple scenario-wise subproblems, facilitating the adoption of a customized column-and-constraint generation (C&CG) algorithm. Finally, case studies demonstrate a remarkable improvement in the out-of-sample performance of our model, compared to its prevailing stochastic and robust counterparts. Moreover, the potential values of incorporating the ambiguity and distributional information are quantitatively estimated, providing a useful reference for planners with different budgets and risk-aversion levels.

decision-dependent uncertainty↗