SEARCH · Engineering Papers

Results for “optimal input design”

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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Mixed Ionic Electronic Conducting Quaternary Perovskites: Materials by Design for Solar Thermochemical Hydrogen

The innovative research conducted by Arizona State University and Princeton University in the project "Mixed Ionic-Electronic Conducting Quaternary Perovskites: Materials by Design for Solar Thermochemical Hydrogen" marks a significant stride forward in thermochemical water splitting. Through an intricate blend of computational design and experimental validation, the project delved into the promising potential of Mixed Ionic Electronic Conducting (MIEC) perovskites. These complex materials, characterized by their unique redox-active nature and adaptability in stoichiometry, present a promising frontier for efficient solar thermochemical hydrogen production. Firstly, the research enhanced the science by utilizing state-of-the-art computational methodologies to unravel the nuanced chemical potentials of MIEC perovskites. By simulating various off-stoichiometric scenarios and redox conditions, the team was able to predict material behaviors under diverse environmental conditions, a feat unachievable through conventional experimental methodologies alone. This approach not only fast-tracks the material screening process, significantly reducing the time from laboratory re-search to practical application, but also uncovers trends and correlations that are pivotal for future materials innovation. Regarding technical effectiveness, the project stands out in its economic feasibility. Traditional methods of materials discovery are often marred by high costs and extensive timeframes, owing to the iterative nature of experimental processes. However, by employing theoretical computations and validating these findings with targeted experiments, the project introduced a cost-effective paradigm for materials discovery and the first ever prediction, synthesis, and preliminary validation of a material solely from computational and theoretical considerations. This synergy between computation and experimentation expedites the discovery of optimal materials conducive to high-efficiency solar-to-hydrogen conversion processes. Furthermore, the public stands to benefit substantially from this research. The success of MIEC perovskites in solar thermochemical applications heralds a shift towards lower cost and lower electricity input for clean hydrogen production, hence potentially impacting climate and energy resilience. By improving the efficiency of solar-to-hydrogen conversions, the research paves the way for reduced dependency on fossil fuels, addressing the urgent global need for accessible and renewable energy sources. Moreover, the project's advancements contribute to scientific literacy in renewable energy technologies, empowering society through knowledge and spurring future innovations. In essence, this research project demonstrates significant progress in the realm of advanced water splitting through solar thermochemistry. Through its groundbreaking approaches in computational materials science and its implications for real-world applications, it holds the promise of a cleaner, more energy-resilient future.

08 HYDROGEN↗

DOE OSTI2023

Rewiring cell-free metabolic flux in E. coli lysates using a block—push—pull approach

Abstract Cell-free systems can expedite the design and implementation of biomanufacturing processes by bypassing troublesome requirements associated with the use of live cells. In particular, the lack of survival objectives and the open nature of cell-free reactions afford engineering approaches that allow purposeful direction of metabolic flux. The use of lysate-based systems to produce desired small molecules can result in competitive titers and productivities when compared to their cell-based counterparts. However, pathway crosstalk within endogenous lysate metabolism can compromise conversion yields by diverting carbon flow away from desired products. Here, the ‘block—push—pull’ concept of conventional cell-based metabolic engineering was adapted to develop a cell-free approach that efficiently directs carbon flow in lysates from glucose and toward endogenous ethanol synthesis. The approach is readily adaptable, is relatively rapid and allows for the manipulation of central metabolism in cell extracts. In implementing this approach, a block strategy is first optimized, enabling selective enzyme removal from the lysate to the point of eliminating by-product-forming activity while channeling flux through the target pathway. This is complemented with cell-free metabolic engineering methods that manipulate the lysate proteome and reaction environment to push through bottlenecks and pull flux toward ethanol. The approach incorporating these block, push and pull strategies maximized the glucose-to-ethanol conversion in an Escherichia coli lysate that initially had low ethanologenic potential. A 10-fold improvement in the percent yield is demonstrated. To our knowledge, this is the first report of successfully rewiring lysate carbon flux without source strain optimization and completely transforming the consumed input substrate to a desired output product in a lysate-based, cell-free system.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Application of Sequential Design of Experiments (SDoE) to Large Pilot-Scale Solvent-Based CO2 Capture Process at Technology Centre Mongstad (TCM)

The United States Department of Energy’s Carbon Capture Simulation for Industry Impact (CCSI2) program has developed a framework for sequential design of experiments (SDoE) that aims to maximize knowledge gained from budget- and schedule-limited pilot scale testing. SDoE was applied to the planning and execution of campaigns for testing CO2 capture systems at pilot-scale in order to optimally allocate resources available for the testing. In this methodology, a stochastic process model is developed by quantifying the parametric uncertainty in submodels of interest; for a solvent-based CO2 capture system, these may include physical properties and equipment performance submodels (e.g., mass transfer, interfacial area). This uncertainty is propagated through the full process model, over variable operating conditions, for estimating the resulting uncertainty in key model outputs (e.g., percentage of CO2 capture, solvent regeneration energy requirement). In developing a data collection plan, the predicted output uncertainty is incorporated into an algorithm that seeks simultaneously to select process operating conditions for which the predicted uncertainty is relatively high and to ensure that the entire space of operation is well represented. This test plan is then used to guide operation of the pilot plant at varying steady-state conditions, with resulting process data incorporated into the existing model using Bayesian inference to refine parameter distributions. The updated stochastic model, with reduced parametric uncertainty from data collected, is then used to guide additional data collection, thus the sequential nature of the experimental design. The SDoE process was implemented at the pilot test unit (12 MWe in scale) at Norway’s Technology Centre Mongstad (TCM) in a summer 2018 test campaign with aqueous monoethanolamine (MEA). During the test campaign, the varied operating conditions included the flowrates of circulated solvent, flue gas, and reboiler steam and the CO2 concentration in the flue gas. The process data were used to update probability distributions of mass transfer and interfacial area parameters of a stochastic process model developed by the CCSI2 team. Two iterations of the SDoE process were executed, resulting in the uncertainty in model predicted CO2 capture percentage decreasing by an average of 58.0 ± 4.7% over the full input space of interest. This work demonstrates the potential of the SDoE process for model refinement through reduction in process model parametric uncertainty, and ultimately risk in scale-up, in CO2 capture technology performance.

carbon capture↗

DOE OSTI2023

Generating MCNP Input Files for Unstructured Mesh Geometries

The Los Alamos National Laboratory’s (LANL) Monte Carlo N-Particle (MCNP)1 transport code version 6.3 (also known as MCNP6.3) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is time-consuming and error-prone as the complexities of geometries increase. A UM geometry model is a collection of finite elements representing a solid geometry. The first step of the MCNP UM calculation is using other software packages to create a finite element mesh representation of a solid 3D geometry because the MCNP code cannot be used to generate a UM model. Computer-aided design (CAD) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. Some mesh generation software packages may also be used to create solid geometries and thus CAD files are not needed. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software suite. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with a 6.3 version, the MCNP code can process HDF5 mesh input files. We only focus on the UM models formatted as Abaqus input files in this report since currently no external software can be used to generate HDF5 mesh input files for MCNP UM calculations. The MCNP code version 6.3 can be used to convert the Abaqus mesh input files into the HDF5 mesh input files, but this option is typically used by the MCNP code development team to test the HDF5 mesh input file feature. Several software packages (such as Abaqus, Attila4MC, or Cubit) can be used to create the Abaqus input files for MCNP UM calculations. An MCNP UM calculation using an Abaqus model requires two input file types: MCNP and Abaqus input files. The Abaqus input files needed for MCNP UM calcu lations must have the correct Abaqus syntax and meet the additional requirements by the MCNP code. The MCNP code can process only Abaqus input files that make use of part and assembly definitions, where elements in each part must be grouped into one or more element sets (i.e., elset) using *Elset keyword lines with specified naming formats. The MCNP and Abaqus input files required for MCNP UM simulations must be related; pseudo-cells in an MCNP input file must be constructed from mesh model data from an Abaqus input file. For large complex UM models, it is tedious to manually create MCNP UM input files. The um pre op (unstructured mesh pre operations) program with the -m option can be used to create a skeleton MCNP input file from an Abaqus input file [6]. Since the um pre op program was written in Fortran and was not written for optimized performance, this program is a deprecated feature in the MCNP code version 6.3 and may be removed in the next release of the code. To improve calculation flow of multiphysics calculations, a Python3 code called write mcnp um input has been developed to generate an MCNP input file instead of using the um_pre_op -m option. This Python code was initially released to the public in 2020. We have updated this Python code for MCNP6.3 and it was used to generate the MCNP input files used to verify the MCNP6.3 code. The write_mcnp_um_input code is included with the MCNP6.3 code package which will be released to the public through the Radiation Safety Information Computational Center (RSICC) at Oak Ridge National Laboratory. This report is a revision of LA-UR-20-27139 report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

High-Level Synthesis of Parallel Specifications Coupling Static and Dynamic Controllers

The increased need for efficient ways to implement domain-specific accelerators is driving design methodologies towards the use of abstractions higher than the Register Transfer Level (RTL). In this scenario, High Level Synthesis (HLS) plays a significant role by enabling the automatic generation of custom hardware accelerators starting from high level descriptions (e.g., C code). Conventional HLS tools exploit parallelism mostly at the Instruction Level (ILP). They statically schedule the input specifications, and build centralized Finite State Machine (FSM) controllers. However, aggressive exploitation of ILP in many applications has diminishing returns and, usually, centralized approaches do not efficiently exploit coarser parallelism because FSMs are inherently serial. In this paper we present a HLS framework able to synthesize applications that, beside ILP, also expose Task Level Parallelism (TLP). An application can expose TLP through annotations that identify the parallel functions (i.e., tasks). To generate accelerators that efficiently execute concur- rent tasks, we need to solve several issues: devise a mechanism to support concurrent execution flows, exploit memory parallelism, and manage synchronization. To support concurrent execution flows, we introduce a novel adaptive controller. The adaptive controller is composed of a set of interacting control elements that independently manage the execution of a single operation or function call. These control elements check dependencies and resource constraints at runtime, enabling as soon as possible execution. To support parallel access to shared memories and synchronization, we introduce a novel Hierarchical Memory Interface (HMI). With respect to previous solutions, the proposed interface supports multi-ported memories and atomic memory operations, which commonly occur in parallel programming. Our framework can generate the hardware implementation of C functions by employing two different approaches, depending on its characteristics. If a function exposes TLP, then the framework generates hardware implementations based on the adaptive controller. Otherwise, the framework implements the function by exploiting a more conventional FSM approach, which is optimized for ILP exploitation. We evaluate our framework on a set of parallel applications, and show substantial performance improvements (average speedup of 4.7) with limited area over- heads (average area increase of 5.48 times).

Castellana, Vito G.↗

DOE OSTI2021

High-Level Synthesis of Parallel Specifications Coupling Static and Dynamic Controllers

Castellana, Vito G.↗

DOE OSTI2022

Fault Characterization and Diagnostics Supporting Condition-Based Operation and Maintenance of Gas Turbine Engines

Condition-Based Operation and Maintenance (CBOM) is the state-of-the-art in maintenance approaches for gas turbine engines. CBOM applies engine sensor information to the optimization of future operation and maintenance procedures; this technique reduces costs for engine operators by minimizing unplanned outages and catastrophic engine degradation. However, due to the complexities of gas turbine operation and the lack of available engine sensors, further development is required to fully realize the benefits of CBOM. In the turbine section, rotating components interact with high temperature flows, which creates intense thermal and mechanical stresses. As a result, there are numerous mechanisms of component degradation in the turbine section. Furthermore, turbine components – like stator vanes and rotor blades – are among the most expensive in the engine because they are complex to design and manufacture. For these reasons, this dissertation addresses two main questions: (i) which parameters or faults within the turbine section are most important to monitor, and (ii) how can these parameters or faults be monitored in an engine-relevant environment? Although many turbine parameters and faults have been investigated in the open literature, there are some faults that are still not well understood. Rotor-casing eccentricity, which causes a non-constant blade tip clearance around the annulus, has not been investigated in terms of its effects on turbine efficiency. Therefore, the first study in this dissertation quantifies overall and local turbine efficiency for varying levels of rotor-casing eccentricity. Results showed negligible variations to overall turbine efficiency, meaning rotor-casing eccentricity only becomes relevant to CBOM when its severity causes rotordynamic issues. Purge flow is critical to turbine hardware longevity because it prevents ingestion of hot main gas path (MGP) flow into the under-platform regions. Despite its importance, there are currently no methods for monitoring purge flow performance in an engine environment. Therefore, the second study in this dissertation develops a predictive model for sealing effectiveness using inputs from two fast-response pressure sensors. Results exhibited low prediction errors across a full range of purge flow rates, which supports the viability of the modelling approach for CBOM. The final two studies in this dissertation address blade coolant flow monitoring. This cooling flow is responsible for protecting the turbine blades from the MGP flow, which exits the combustor at temperatures greater than the blade melting point. These studies showed that temperature measurements on the blade surface can be used to accurately predict blade coolant flow rate, and that defining the candidate features relative to the coolant trajectory is important for maintaining accuracy as coolant flow rate degradation occurs. This work enables blade coolant flow monitoring, which is currently not possible through existing condition monitoring techniques.

condition-based, gas turbines, diagnostics,↗

DOE OSTI2022

FORCE-DISPATCHES Integration - Initial Demonstration

Integrated energy systems (IES) combine, in mutually beneficial ways, power from variable renewable energy sources and nuclear power plants (NPP) to improve economic viability under uncertain market and weather conditions. The open-source Framework for Optimization of Resources and Economics (FORCE) tool suite, developed at Idaho National Laboratory (INL), has enabled comprehensive modeling and simulation of IES. The capabilities within FORCE include grid portfolio optimization through the Holistic Energy Resource Optimization Network (HERON) and the transient process model analysis library HYBRID, among others. Continuous efforts and investments from the IES programs have been made to expand and improve the versatility of the FORCE toolset in fiscal year 2022. Code-coupling and cross-tool communication have been important methods for improving this versatility. This report focuses on an additional workflow in the HERON tool for capacity and dispatch stochastic optimization through integration with the external tool Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems (DISPATCHES). DISPATCHES was primarily developed by the National Energy Technology Laboratory, in collaboration with other national laboratories, which included INL, universities, and industry partners. It is coupled to a library of algebraic models for specific plant components, and to a framework for stochastic optimization different from that provided in the current Risk Analysis Virtual Environment (RAVEN)-running-RAVEN algorithm in HERON. HERON currently conducts stochastic optimization via an outer-inner loop: it optimizes over variable capacity on the outer loop, and at each step within the capacity parameter space, conducts an inner optimization over scenarios (of market signals, demand, and/or weather patterns) and hourly dispatch throughout a user-specified number of years. On the other hand, DISPATCHES conducts stochastic optimization via an “all-at-once” strategy in which capacity variables are optimized at the same level as dispatch variables, as all scenarios are considered at once. The latter method works especially well for projects of limited size and project length, as the necessary computational power and memory increases with the number of variables and scenarios. The new capability to use the DISPATCHES workflow in HERON enhances standalone simulations by leveraging FORCE tools—namely, the economic metrics from the Tool for Economic Analysis (TEAL) and reduced-order model (ROM) sampling from RAVEN. The initial demonstration of the DISPATCHES workflow simulates an existing nuclear-case flowsheet within the DISPATCHES repository—this models a NPP with a secondary revenue stream for hydrogen production. Electrical output from the plant is converted to hydrogen via a proton-exchange membrane (PEM) electrolyzer, hydrogen tanks are used for storage, and an additional turbine is added for hydrogen combustion. Continued work regarding this FORCE-DISPATCHES integration will include automatic generation of DISPATCHES models from HERON inputs, offering analysts the option of using either the RAVEN-runsRAVEN or DISPATCHES workflow to solve technoeconomic optimization problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Technology Demonstration of a High-Pressure Swirl Oxy-Coal Combustor

This technical report presents the exploration of the design and prototyping of a High-Pressure Swirl Oxy-coal Combustor. Pressurized oxy-coal combustion systems have the potential to improve efficiency along with an increased carbon capture rate. Reduction of flue gas at higher pressure, smaller system size, and capital cost reductions render high-pressure oxy-coal systems particularly attractive as next-generation energy-producing systems. High-pressure oxy-coal combustion systems are a recent concept, and thus operability issues of combustor designs for such systems are not fully understood. Significant challenges exist to maintain oxy-coal combustion stability at elevated pressure and a high CO 2 diluent environment. Although a body of knowledge exists for high-pressure oxygen combustion in rocket engines (or similar applications), it is yet to be strategized how these fundamental concepts can be translated to low-temperature CO 2 diluent combustion regimes. The realization of the pressurized oxy-coal based systems requires combustor components to be designed and demonstrated for an operating pressure over 10 bar. However, pressurized oxy-coal combustor design information at this pressure range and scale relevant to validate those proposed systems is currently limited. Experimental data from MWth scale oxy-coal combustors are needed to identify the optimal trade-off between net efficiency and systems size. The proposed effort is aimed at demonstrating a 1 MWth down-fired swirl Oxy-Coal combustor and investigate the interrelation between combustor operating conditions (pressure; flame stability; flue gas recirculation ratio) and conversion efficiencies to minimize oxygen requirements. One of the key challenges is to configure burner design (i.e., swirl number and injector) and operating conditions for high-pressure oxy-coal combustion systems. These experiments differ from current systems partly due to the high theoretical flame temperature and related burner operability issues associated with oxy-combustion. An ASPEN PLUS® model study for 550 MWe TIPS and ENEL pressurized oxy-coal systems with CO 2 recirculation was performed to evaluate system design, subsystems sizing, and operating condition determination. The system analysis effort included TRL and technology gap determination of subsystems and critical components. This information was scaled to develop design requirements (design pressure and flue gas recirculation: RR Flue Gas = $\frac{m_{flue}}{m_{total}}$) for the 1 MWth combustor. The effects of a wide range of carbon dioxide recirculation ratios on the thermal efficiency of ENEL and TIPS cycles are studied. The pressure of 10 bar and 80 bar are used for ENEL and TIPS cycles, respectively. The thermal efficiency of ENEL is significantly higher than the efficiency of TIPS at a pressure of less than 10 bar. The insights from system analysis were then used to design a 1 MWth swirl oxy-coal combustor. Flame temperature analysis and material strength analysis was performed to determine the combustor thickness. The structural integrity of the combustor was validated by finite element analysis using Abacus® and Hypermesh®. Feasibility of igniters and secondary burners are investigated in successful high-pressure oxy-methane combustion. The secondary burners are designed in such a way that it can operate between 100 to 500 kW firing input. Three generations of the pintle injector were designed based on swirl numbers (S=0, 0.9, and 1.2). Key pintle injector parameters such as pintle size, pintle orifice size, spray pattern were investigated by cold flow tests. Information from these tests was used to modify injector design for smooth and successful operation. A 5 mm pintle orifice size was decided upon as the optimum size for oxy-coal operation for the combustor. Shadow sizing experiments were performed to identify the atomization rate of each injector. Different coal water slurry mixtures (30 – 50% coal by wt% in the mixture) at various total momentum ratios (TMR) were investigated for this purpose. These experiments provided decisive information to choose the best design of the injector. The injector with 1.2 swirl provided higher atomization in all cases than other designs. The mean equivalent droplet size of the jet was similar at different TMR and mixture ratios using this injector, thus making it suitable for use in most cases. Therefore, the 1.2 swirl-pintle injector was chosen for the shakedown test. The combustor and other sub-systems, including feed systems and control and data acquisition, have been manufactured, assembled, and integrated. The total system integration and installation began on July 1, 2020. The shake-down tests and initial operational capability demonstration are expected to be completed by September 30, 2020.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

Virtual Engineering Software Framework for Integrated Biomass Conversion Modeling

This presentation covers the design and implementation of a software tool to systematically connect computational models of unit operations to simulate an integrated process of low-temperature conversion of biomass to fuel. This virtual engineering (VE) software was designed with the overarching goal of connecting unit models written in various programming languages and requiring different computational resources within a single, flexible framework. The models and features currently considered for the VE library include mechanistic models for pretreatment, enzymatic hydrolysis, and aerobic bioreaction; high-fidelity computational fluid dynamics (CFD) simulations for enzymatic hydrolysis and aerobic bioreaction; and the capability to perform techno-economic analyses (TEA) using Aspen Plus, a commercial software package. The CFD models require access to high-performance computing (HPC) resources, so in addition to handling multiple programming languages and interfaces, the VE software must also be capable of interacting with an HPC scheduler to submit, run, and post-process jobs. Using the Python programming language, a new VE software package has been developed that contains functionality to manage the input-output communication between various unit models, schedule simulations to run on NREL's HPC and analyze those results, and interface with existing TEA software workflows. A Jupyter-notebook GUI was also created to solicit user input and provide documentation. In cases where multiple models for a particular unit-operation exist, selection between models is accomplished through a simple checkbox, with the appropriate inputs and outputs being parsed and converted seamlessly in the background. Each operation makes use of a different programming language, but the flow of information from pretreatment to enzymatic hydrolysis to bioreaction is managed with an intuitive, centralized file-communication strategy. In this talk, the programming approach and implementation details of the notebook are presented for multiple possibilities of the conversion process, including a demonstration of the ability to manage HPC resources. Additionally, an example of a sensitivity study of treatment parameters governing the overall conversion outcome is shown which highlights the ease of defining new problems using the VE Notebook workflow and leads into a discussion of ongoing work to enable outer-loop optimization studies.

biofuel↗

DOE OSTI2022

Integrated Optimization and Control of a Hybrid Gas Turbine/sCO 2 Power System

During phase-I, the project team led by Echogen Power Systems (EPS) had two primary objectives based on investigating the application of gas turbines with supercritical carbon dioxide (sCO 2 ) power cycles. The first objective was to improve the overall efficiency and performance of a hybrid gas turbine/sCO 2 power system through a joint optimization of the two subsystems (gas turbine and sCO 2 power cycle) using non-linear optimization techniques that simultaneously evaluate thermal performance of the combined cycle. The hybrid power system included several points of interaction, including (but not limited to) gas turbine exhaust, fuel heating, inlet chilling and turbine cooling. The second objective was to establish a baseline transient response model of the hybrid power system and a notional microgrid and begin steps to integrate the control systems of the three major elements (gas turbine, sCO 2 cycle and grid controller). The project team established a baseline performance for a combined cycle power plant using a production gas turbine and scaled sCO 2 power cycle only utilizing exhaust heat recovery. Echogen’s non-linear techno-economic optimization code was extended by adding gas turbine component models derived from a in-house developed gas turbine design code. With the two cycles coupled by the gas turbine exhaust, design parameters of both cycles were allowed to vary simultaneously to determine performance opportunity versus isolated designs. Returning to the baseline gas turbine/sCO 2 power cycle transient models: Echogen had in-house developed sCO 2 cycle transient model in GT-Suite system simulation software, and had partnered with Siemens Finspång for gas turbine transient model, and Siemens PTI group to provide micro-grid load profile as well as hybrid power cycle generated load (power and frequency) analysis. The transient model for the SGT-750 Siemens gas turbine was a “black-box” functional mock-up interface (FMI) model developed by Siemens Industrial Turbomachinery in Finspång, Sweden. The SGT-750 is a twin-shaft gas turbine that produces 40 MW electricity with an efficiency of about 40% at ISO conditions. At 100% gas turbine throttle (load), the SGT-750 has average exhaust conditions of 114.6 kg/s and 469.8°C. The transient model for sCO 2 power cycle was developed by Echogen in GT-SUITE 1D system simulation software platform. The basic CO 2 flow circuit has single-shaft turbomachinery with net 11.5 MW electrical power output at design conditions. The power turbine has a double-ended shaft with one end connected to synchronous generator through a fixed-ratio gearbox. The other end of power turbine is connected to the compressor through a continuously variable transmission. The major components of the sCO 2 power cycle modeled include air cooled condenser/cooler, CO 2 compressor, recuperator, two waste heat exchanger coils, power turbine, continuous variable transmission, gearbox and generator. Integration of SGT-750 transient model and sCO 2 power cycle transient model was done in Matlab Simulink. In the integrated model, the gas turbine and sCO 2 power cycle interacted at two points, first one being the gas turbine exhaust gas flow rate and temperature, which were inputs to sCO 2 power cycle model. The second point was the distribution of micro-grid load demand signal between the SGT-750 generator and sCO 2 cycle generator. For a given combined-cycle load demand, the gas turbine load demand was equal to the total demand minus the sCO 2 cycle power generated. In the present study the integrated model was simulated for two cases of grid load demand: (i) for a step change, both positive-step and negative-step, in grid load demand (ii) for a micro-grid load demand curve provided by Siemens PTI group. Finally, the time series plots representing load demand versus integrated system response were presented including the sCO 2 power cycle control system performance plots. The actual generated power and frequency of both the generators, gas turbine and sCO 2 power cycle, was supplied to Siemens PTI group for dynamic grid assessment, results of which are provided in appendices.

03 NATURAL GAS↗

DOE OSTI2021

Bamboo Bio-composite Truck/Trailer Decking

Bamboo is one of the fastest growing plants in the world. Massively productive, bamboo will maintain that productivity with limited inputs and minimal management resulting in predictable volume and operating margins and with low overhead. The goal of this project was to develop a bamboo bio-composite trailer decking product that replaces apitong, is lighter weight, stronger and luminescent. The collaboration partners, Fontaine, Resource Fiber and IACMI, worked closely to develop the bamboo bio-composite decking to Fontaine Trailer’s (end-user) specifications. The various technological objectives were (a) To conduct a comprehensive Design of Experiments study for bamboo composites to understand the structure-property relationships for different resin composites and bamboo forms (strips, woven, bulk etc); (b) To establish the nail pull out characteristics of trailer deck geometry bamboo composite form. The nail pull out is a critical test in trailer decks; (c) To design, process and prototype select number of trailer decking planks for Fontaine evaluation; and (d) To conduct life cycle analysis to conduct energy calculations from the various conversion steps of the bamboo from crop to product. The project meets DOE/IACMI metrics of reduced embodied energy, lightweighting and lowering the cost of the end product. It also offers a green solution to a value-added application, i.e. trailer decking. The development goals and benefits of bamboo composite decking included: A. Decreased weight when compared to Apitong and Gen 1, thereby reducing operating costs and petroleum usage due to increased miles per gallon. B. Flexibility to allow camber design into trailers. C. Luminescence to increase safety during low light and dark conditions. D. Composite material to maximize use of bamboo bio-composites as substitute for petroleum-based composites. E. Embedded layer(s) of conditioned bamboo for added strength and stability. F. Improved safety over Apitong by temporary cargo indentation in planks for better stability. G. Cradle-to-cradle design so end-of-life becomes beginning-of-life for other products. Resource Fiber’s bamboo biocomposite trailer decking passed nail pullout tests as compared to apitong (incumbent), was lighter weight than apitong, and used less embodied energy particularly when bamboo is sourced domestically. The planks were optimized in the lab setting to the extent possible and successfully installed on the trailer at Fontaine for field testing. Despite the process not being fully optimized due to through-heat and tooling limitations within the budget constraints, the prototype decking withstood 400 cycles of reverse fatigue loading of a 107,000 lb Caterpillar 349F which was extreme conditions of field testing by Fontaine Trailer. Resource Fiber plans to outsource production of a next round of planks with a commercial pultrusion processor, then to do a re-test with Fontaine Trailers. Tooling specific to the part is required. Long-term commercial plans are to continue outsourcing production while supplying bamboo fiber and mats. Commercial markets include heavy haul trailers, military trailers, and decking for marine, industrial and residential use.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Co-synthesis of Hydrogen and High-Value Carbon Products from Methane Pyrolysis

The ARPA-E Methane Pyrolysis Project successfully developed a scalable technology for hydrogen production with low-CO x emission through methane pyrolysis, co-producing high-value carbon nanotubes (CNTs). The project focused on optimizing reactor design, enhancing catalyst performance, and assessing techno-economic feasibility to create a commercially viable and environmentally sustainable process. The fluidized bed reactor achieved over 90% methane (CH 4 ) conversion by using a 5% CO 2 co-feed, which stabilizes carbon yields and minimizes catalyst deactivation. This setup allowed for continuous operation across ten cycles, each consisting of a 14-minute pyrolysis phase followed by a 10-minute dislodging phase to remove a fraction of the accumulated carbon, resulting in stable performance and high-quality CNT production. In parallel, monolith reactors coated with Fe demonstrated a sustained methane conversion of 73% while producing CNTs with high crystallinity. Although promising for continuous operation, monolith reactors face challenges in coating durability and scalability, highlighting areas for further optimization in commercial applications. Catalyst formulation played a key role in enhancing process efficiency. The core catalyst used was 5%Fe/Al 2 O 3 (wt%), optimized through wet impregnation, which improved CNT morphology, yielding longer and more uniform CNTs. The catalyst's performance was further enhanced by adding promoters: 2.5 wt% Ni increased methane conversion close to the thermodynamic limit, while 2.5 wt% Mn improved CNT alignment and crystallinity, and 1.5 wt% NaCl boosted CNT morphology but slightly lowered methane conversion. These adjustments allowed the reactor to maintain high methane conversion while producing high-quality CNTs, enabling stable performance over multiple cycles. To address carbon buildup and ensure uninterrupted operation, a pneumatic conveying tube was implemented for effective carbon dislodging in the fluidized bed configuration. CO 2 and H 2 O co-feeds were also introduced to enhance carbon removal, with CO 2 boosting CNT yield by approximately 15%. This setup enabled stable reactor operation across multiple cycles, preventing clogging and minimizing catalyst wear, making the process suitable for industrial scaling. Techno-economic analysis (TEA) projected hydrogen production costs between $\$$1.00 and $\$$1.64 per kilogram, with CNT values assumed at $\$$375/ton and $\$$100/ton. The life cycle assessment showed that CO 2 emissions could be as low as 0.64 kg CO 2 e/kg H 2 at 95% methane conversion assuming an electricity input of 50 kg CO 2 e/MWh. Even at 50% methane conversion, emissions remained below 1 kg CO 2 e/kg H 2 , demonstrating the process's low-emission potential and making it a viable alternative to traditional steam methane reforming. Overall, the results from this project demonstrate the feasibility of a pyrolysis process where carbon is continuously removed from the catalyst surface and hydrogen is continuously produced until a catalyst regeneration step is required to fully clean the catalyst surface and renew catalyst performance. Major open challenges are related to avoiding the loss of catalyst material in the dislodged carbon during fluidized bed conditions, since our best result demonstrated a carbon purity of ~70 wt. % (rest being iron and alumina). A monolith reactor was used to favor dislodgement of carbon compared to fluidized bed conditions but our results do not demonstrate an advantage of the monolith configuration. Catalyst performance was similar to fluidized bed conditions with slower deactivation rates overall, but we could not observe carbon dislodging in any of the tens of experiments that were run at Stanford. Our results show that the most relevant areas of improvement are related to the fundamental understanding of the iron-carbon interface for dislodging, and the development of catalyst that can produce CNTs via a base-growth mechanism such that catalyst is not lost in the dislodgement steps. The final report documents all findings and methodologies in detail, providing a valuable resource for the scientific community. By building on these results, researchers can further advance methane pyrolysis technology, moving toward a more sustainable, scalable pathway for hydrogen production. This work lays the foundation for future research and commercial efforts to reduce emissions in hydrogen production while generating valuable carbon products.

08 HYDROGEN↗

DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗