Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Characterizing Impacts of Dry Coal Feeding in High Pressure Oxy-Coal Combustion Systems

Reaction Engineering International (REI) has managed a team of experts from the University of Utah, Southeast University (SEU) in Nanjing, China, Electric Power Research Institute (EPRI), Corrosion Management Ltd. (C-M), Praxair, and Brigham Young University (BYU) to investigate dry pulverized coal feeding for pressurized oxy-coal combustion. Dry feed firing systems for entrained flow, pressurized, oxy-coal combustors have not been well developed, although related technologies have been used in the Shell Gasification Process and for pressurized fluid bed combustion. DOE-funded research recently completed at REI and the University of Utah focused on characterizing impacts of high temperatures and pressures in oxy-coal combustion systems. For high pressure combustion, that research used a coal slurry feed into a 17 bar pressurized combustor. As a consequence of that research, it was identified that fuel feeding and firing system flexibility are challenges that require attention. Based on that experience, the approach of using a coal slurry feed system leads to challenges in producing consistent atomization of the slurry, which causes burnout problems, especially at high pressures. In addition, slurry atomization processes may be difficult to scale to sizes appropriate for practical commercial use. Dry pressurized coal burner systems, on the other hand, have the potential to yield efficiency gains, improve flexibility and facilitate applications at larger scales. Experimental work was conducted at the University of Utah Industrial Combustion and Gasification Research Facility as well as the 100 kW pressurized oxy-coal combustor (POC) facility at Brigham Young University. Mechanism development and CFD-based combustion and dense-phase flow modeling were performed at REI. Successful completion of the project objectives has resulted in the following key deliverables: 1) Design and prototype of a pressurized pulverized coal feeding and oxy-firing system 2) Data from a 100kW, pressurized (15 bar) entrained flow reactor with a dry feeding delivery and burner system that describes flame characteristics, radiative heat flux profiles, carbon burnout, along with characteristics of ash aerosols, fouling, and slagging. 3) Validated and transportable models that describe the relevant conditions in pressurized oxy-combustion systems and that can be used for scale-up and optimization. 4) Principles to guide design of high pressure, pilot-scale and full-scale coal oxy-firing systems. 5) Assessment of pressurized oxy-combustion impacts on key parameters relevant to oxy-coal fired utility boilers such as coal devolatilization, char oxidation, mineral matter transformation, deposition, and corrosion. The experimental data, pressurized oxy-firing system principles, and process mechanisms provided by this work can be used by electric utilities, boiler OEMs, equipment suppliers, design firms, software vendors, consultants and government agencies to assess the use of high temperature and high pressure oxy-combustion in current research and to guide development of new oxy-coal boiler designs.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Taylor dispersion in osmotically driven laminar flows in phloem

Sucrose is among the main products of photosynthesis that are deemed necessary for plant growth and survival. It is produced in the mesophyll cells of leaves and translocated to different parts of the plant through the phloem. Progress in understanding this transport process remains fraught with experimental difficulties, thereby prompting interest in theoretical approaches and laboratory studies. The Münch pressure and mass flow model is one of the accepted hypotheses describing the physics of sucrose transport in the phloem. It is based on osmosis creating an energy potential difference between the source and the sink. The flow responding to this energy potential is assumed laminar and described by the Hagen–Poiseuille equation. This study revisits such osmotically driven flows in tubes with membrane walls by including the effects of Taylor dispersion on mass transport. This effect has been overlooked in phloem flow studies. Taylor dispersion can increase the effective transport of solutes by increasing the apparent diffusion coefficient. It is shown that, in addition to the conventional diffusive correction derived for impermeable tube walls, a new adjustment to the mean advective terms arises because of osmotic effects. Because the molecular Schmidt number is very large for sucrose in water, the sucrose front speed and travel times have a direct dependence on the Péclet number for different ranges of the Münch number. Fially, this study establishes upper limits on expected Taylor dispersion enhancement of sucrose transport.

42 ENGINEERING↗
DOE OSTI2021

DNS/LES Study of Representative Wall-Bounded Turbulent Flows using SIERRA/Fuego

This report summarizes a series of SIERRA/Fuego validation efforts of turbulent flow models on canonical wall-bounded configurations. In particular, direct numerical simulations (DNS) and large eddy simulations (LES) turbulence models are tested on a periodic channel, a periodic pipe, and an open jet for which results are compared to the velocity profiles obtained theoretically or experimentally. Velocity inlet conditions for channel and pipe flows are developed for application to practical simulations. To show this capability, LES is performed over complex terrain in the form of two natural hills and the results are compared with other flow solvers. The practical purpose of the report is to document the creation of inflow boundary conditions of fully developed turbulent flows for other LES calculations where the role of inflow turbulence is critical.

42 ENGINEERING↗
DOE OSTI2020

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗
DOE OSTI2022

Toward a machine-guided approach to energetic material discovery

In this article, we trained a machine learning (ML) model to connect microstructural details of an energetic material formulation to its performance for the purpose of guiding the discovery of new explosive formulations. Our hypothesis was that the algorithm would robustly learn the training data and produce an accurate surrogate model. Specifically, the algorithm learned the relationship between details of the void size distribution (VSD), initiating shock pressure, and the energetic material performance. We used realistic constraints on the VSD and a range of cases were ingested by a physically informed reactive flow model working within a hydrodynamic solver running on high-performance computing resources. The ML algorithm produced a surrogate model that accurately predicted known test points around the parameter space. In addition to the utility of the model and the process used for its development, we noted interesting comparisons between what we, the authors—subject matter experts, would heuristically conclude from the training data and the surrogate model predictions. We detected nuanced details that were missed by the surrogate model; however, these details are not important to an energetic material formulator. We concluded that the algorithm did indeed robustly learn the training data and produce an accurate surrogate model. We further concluded that the surrogate model is a powerful tool to guide the formulator in the absence of subject matter experts and limited-access computing resources.

42 ENGINEERING↗
DOE OSTI2024

Modeling Radionuclide Inventories in MSR Off-Gas Systems with Radiochemical Transport Analysis

Molten salt reactors (MSR) contain unique characteristics that may require enhancements to modeling tools to accurately predict phenomena. One characteristic that may be advantageous to leverage during normal operation is on-line processing of the circulating fuel salt, such as an off-gas system (OGS) to remove volatile fission products. Therefore, new modeling tools must be developed to integrate spatial resolution and chemistry effects into fuel depletion tools to be able to account for these non-core sources of radioactivity. Such types of radiochemical transport analysis tools were used to estimate the removal rates for 12 elements within a flow model of the Molten Salt Reactor Experiment (MSRE) by optimizing against legacy experimental data of the gas-borne (GB) percentages of 12 nuclides. The removal rates were used in a depletion model to calculate the FP inventory that enters the MSRE OGS. Calculations are in good agreement with the empirical GB percentages reported for the 12 nuclides, which validates the approach and verifies each tool’s treatment of the radiochemical flow effects. The OGS inventory is discussed in terms of the largest nuclide contributors to activity, dose consequence, decay heat, and elemental composition. Finally, insights from the study allow recommendations to be made for future code development activities.

Shahbazi, Shayan↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Can Large Language Models Understand Intermediate Representations?

Intermediate Representations (IRs) are essential in compiler design and program analysis, yet their comprehension by Large Language Models (LLMs) remains underexplored. This paper presents a pioneering empirical study to investigate the capabilities of LLMs, including GPT-4, GPT-3, Gemma 2, LLaMA 3.1, and Code Llama, in understanding IRs. We analyze their performance across four tasks: Control Flow Graph (CFG) reconstruction, decompilation, code summarization, and execution reasoning. Our results indicate that while LLMs demonstrate competence in parsing IR syntax and recognizing high-level structures, they struggle with control flow reasoning, execution semantics, and loop handling. Specifically, they often misinterpret branching instructions, omit critical IR operations, and rely on heuristic-based reasoning, leading to errors in CFG reconstruction, IR decompilation, and execution reasoning. The study underscores the necessity for IR-specific enhancements in LLMs, recommending fine-tuning on structured IR datasets and integration of explicit control flow models to augment their comprehension and handling of IR-related tasks.

Jiang, Hailong↗
DOE OSTI2026

Hydrogen Production Cost from Proton-Conducting Solid Oxide Electrolysis

Rigorous stakeholder-vetted techno-economic analysis (TEA) was conducted to estimate the cost of hydrogen (H 2 ) production using Proton-Conducting Solid Oxide (PSO) electrolysis. The analysis evaluates Current (2025) and Future (2035) technology cases at centralized plant scales of 50 and 500 metric tonnes per day (MTD), assuming electricity, water, and air as the only system inputs. Untaxed, unsubsidized levelized cost of hydrogen (LCOH) is projected to range from 2020 $\$$1.81 to $\$$2.47/kg H 2 at an electricity price of $\$$0.03/kWh and 97% capacity factor under Nth-of-a-kind (NOAK) deployment assumptions. Total installed capital cost was developed using bottom-up Design for Manufacture and Assembly (DFMA) stack cost modeling and detailed balance-of-plant estimates, including mechanical and electrical subsystems, installation, site preparation, engineering, and contingency. Stack performance assumptions include thermoneutral operation, degradation over time, and periodic replacement. LCOH was calculated using the Hydrogen Analysis (H2A) discounted cash flow model in constant 2020 dollars. Results indicate PSO electrolysis has potential for competitive hydrogen production costs under low-cost electricity and mature manufacturing conditions.

08 HYDROGEN↗
DOE OSTI2022

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗
DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗
DOE OSTI2024

Analysis of AP1000 Small-Break Loss-of-Coolant Accident Using Reactor Transient Simulator

The Westinghouse Electric Company’s Advanced Passive Reactor (AP1000) is characterized by the incorporation of passive safety systems (PSSs) designed to ensure core cooling during transient events. The assessment of PSSs requires evaluation of their performance through a combination of experiments and simulations employing various thermal-hydraulic codes. In addition, detailed evaluation of PSSs for a specific reactor system transient analysis such as loss-of-coolant-accident analysis supports understanding representative integral effects test facility development and the further evolution model development and assessment process. Developing a reactor system code is a complex and time-consuming process that requires significant engineering expertise and effort. It can take several months to even years to complete in the early stages of reactor system design and analysis. However, this process can be expedited through the use of transient simulator models for similar reactor systems, which can be used for lesson learning and training purposes. This study uses the Personal Computer Transient Analyzer (PCTRAN) code. The main advantage of PCTRAN is its ease of use and ability to run faster than real time. This study presents the results obtained for a small-break loss-of-coolant accident (SBLOCA) for two breaks using the full version (licensed) of PCTRAN. The purpose of this investigation is to evaluate the overall system behavior during the postulated SBLOCA event as well as assess the capability of the PCTRAN code to reproduce the system response during transient events. The obtained results were compared with the Westinghouse NOTRUMP system code. The PCTRAN code proved to be reliable in predicting the qualitative behavior of the system in both transient cases. As for the system response, it was found that it is contingent on the activation time of the PSSs. The differences in reactor coolant system pressure between the two codes were attributed to the critical flow model and simplification of mass and energy balance. Despite PCTRAN’s limitations, it can still provide a reasonable prediction of various reactor parameters such as pressure, mass flow rate, and void fraction during a SBLOCA scenario. It is worth noting that PCTRAN currently employs a bulk approach similar to that of the Modular Accident Analysis Program (MAAP) and MELCOR codes. However, the upcoming version of PCTRAN will include an artificial intelligence–based detection and accident prevention system, as well as different models for different reactor components. Consequently, PCTRAN has the potential to be upgraded to match the system thermal-hydraulic codes of the U.S. Nuclear Regulatory Commission and become more widely used in cybersecurity to safeguard nuclear power plants from cyberattacks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Normalizing flows for likelihood-free inference with fusion simulations

Fluid-based scrape-off layer transport codes, such as UEDGE, are heavily utilized in tokamak analysis and design, but typically require user-specified anomalous transport coefficients to match experiments. Determining the uniqueness of these parameters and the uncertainties in them to match experiments can provide valuable insights to fusion scientists. Here, we leverage recent work in the area of likelihood-free inference (‘simulation-based inference’) to train a neural network, which enables accurate statistical inference of the anomalous transport coefficients given experimental plasma profile input. UEDGE is treated as a black-box simulator and runs multiple times with anomalous transport coefficients sampled from priors, and the neural network is trained on these simulations to emulate the posterior. The neural network is trained as a normalizing flow model for density estimation, allowing it to accurately represent complicated, high-dimensional distribution functions. With a fixed simulation budget, we compare a single-round procedure to a multi-round approach that guides the training simulations toward a specific target observation. Finally, we discuss the future possibilities for use of amortized models, which train on a wide range of simulations and enable fast statistical inference for results during experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Operation-adversarial scenario generation

This paper proposes a modified conditional generative adversarial network (cGAN) model to generate net load scenarios for power systems that are statistically credible, conditioned by given labels (e.g., seasons), and, at the same time, “stressful” to the system operations and dispatch decisions. The measure of stress used in this paper is based on the operating cost increases due to net load changes. The proposed operation-adversarial cGAN (OA-cGAN) internalizes a DC optimal power flow model and seeks to maximize the operating cost and achieve a worst-case data generation. The training and testing stages employed in the proposed OA-cGAN use historical day-ahead net load forecast errors and has been implemented for the realistic NYISO 11-zone system. In conclusion, our numerical experiments demonstrate that the generated operation-adversarial forecast errors lead to more cost-effective and reliable dispatch decisions.

42 ENGINEERING↗