Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Validation of MCNP Critical Benchmark Models of PU-MET-FAST-016

A new centralized repository of high-quality MCNP models of critical benchmark experiments is currently under development at Los Alamos National Laboratory (LANL). The repository is the combined effort of the Nuclear Criticality Safety, Nuclear Data, and Monte Carlo code development/application organizations at LANL. The initial set of benchmark models in the repository are derived from the Whisper Suite provided with MCNP6.2 and are evaluated against the benchmark experiments described in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The goal is to build a single LANL benchmark collection that is up to date with the latest ICSBEP revision, has a formal review and revision process, is contained in an open-source repository, and utilizes new Python tools for improved input and output file review. This paper describes the validation of the models associated with PU-MET-FAST-016, “Flooded 3X3X3 Arrays of 3-kg Plutonium (Pu) Metal Cylinders – Phase I”. In this critical experiment, twenty-seven 3-kg Pu metal cylinders with a height of 46.33 ± 0.15 mm and diameter of 65.25 ± 0.05 mm were arranged in a cubic array. The Pu cylinders were sealed in aluminum cans, and three cans were placed in each of nine perforated aluminum sleeves, as shown in Fig. 1. The nine sleeves were supported by an aluminum frame that allowed for control of the lateral spacing of the cylinders, both in the X and Y directions. The benchmark experiment measured thirteen configurations. Of these configurations, six were modeled using MCNP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗
DOE OSTI2022

Materials Dynamics

The Materials Dynamics area of leadership focuses on understanding process-structure-properties-performance (PSPP) relationships for the extreme conditions of dynamic loading. This research encompasses controlled synthesis of materials to meet dynamic performance requirements and entails computational coupling across length and time scales for three-dimensional microstructure modeling. For this leadership area, we define dynamic loading as strain rates ≥ 10 3 /s and often, high pressures. A key grand challenge of this area is to predict and measure the evolution of microstructural phases, defect structures, and electronic structure under dynamic conditions while also measuring local temperature to understand transition states. Solving this challenge will require agile, multi-dimensional data analysis and interpretation capability.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulating Real Atmospheric Boundary Layers at Gray-Zone Resolutions: How Do Currently Available Turbulence Parameterizations Perform?

Recent computational and modeling advances have led a diverse modeling community to experiment with atmospheric boundary layer (ABL) simulations at subkilometer horizontal scales. Accurately parameterizing turbulence at these scales is a complex problem. The modeling solutions proposed to date are still in the development phase and remain largely unvalidated. This work assesses the performance of methods currently available in the Weather Research and Forecasting (WRF) model to represent ABL turbulence at a gray-zone grid spacing of 333 m. We consider three one-dimensional boundary layer parameterizations (MYNN, YSU and Shin-Hong) and coarse large-eddy simulations (LES). The reference dataset consists of five real-case simulations performed with WRF-LES nested down to 25 m. Results reveal that users should refrain from coarse LES and favor the scale-aware, Shin-Hong parameterization over traditional one-dimensional schemes. Overall, the spread in model performance is large for the cellular convection regime corresponding to the majority of our cases, with coarse LES overestimating turbulent energy across scales and YSU underestimating it and failing to reproduce its horizontal structure. Despite yielding the best results, the Shin-Hong scheme overestimates the effect of grid dependence on turbulent transport, highlighting the outstanding need for improved solutions to seamlessly parameterize turbulence across scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Frustrated Heisenberg J 1 – J 2 model within the stretched diamond lattice of LiYbO 2

We investigate the magnetic properties of LiYbO 2 , containing a three-dimensionally frustrated, diamondlike lattice via neutron scattering, magnetization, and heat capacity measurements. The stretched diamond network of Yb 3+ ions in LiYbO 2 enters a long-range incommensurate, helical state with an ordering wave vector k = (0.384,±0.384,0) that “locks-in” to a commensurate k = (1/3,±1/3,0) phase under the application of a magnetic field. The spiral magnetic ground state of LiYbO 2 can be understood in the framework of a Heisenberg J 1 –J 2 Hamiltonian on a stretched diamond lattice, where the propagation vector of the spiral is uniquely determined by the ratio of J 2 /|J 1 |. Here, the pure Heisenberg model, however, fails to account for the relative phasing between the Yb moments on the two sites of the bipartite lattice, and this detail as well as the presence of an intermediate, partially disordered, magnetic state below 1 K suggests interactions beyond the classical Heisenberg description of this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Materials Characterization, Prediction, and Control Project: Summary Report on Material Characterization, Part 1

The Pacific Northwest National Laboratory (PNNL) undertook the Materials Characterization, Prediction, and Control (MCPC) Laboratory Directed Research and Development Project to advance understanding of nuclear material processing and enable multifold acceleration in the development and qualification of new material systems in national security and advanced energy applications (Smith 2021). The MCPC Project executed research across three scientific vertices—material characterization, predictive modeling, and data analytics—with extensive support by a data curation and management team. The central technical objective in the MCPC Project was to improve the prediction and characterization of the process-structure-property relationships within the microstructurally refined region of stainless-steel samples prepared utilizing friction stir processing (FSP). Application of the FSP technique is well established at PNNL within the Solid Phase Processing capability through many years of investment across a range of materials and applications (PNNL 2024). Three distinct rounds of FSP experiments were performed by the experimental team, producing replicate samples utilizing across different nominal processing conditions (Condition IDs) listed in Table 1. The starting material on which FSP was applied was commercially available unprocessed stainless-steel type 316L material. Chosen processing conditions were very diverse, and some were intentionally chosen to produce defects. Several samples experienced tool breakage during experimentation, so a full set of three replicates was not produced for every nominal processing condition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Measurements and models of enhanced recombination following inner-shell vacancies in liquid xenon

Electron-capture decays of 125 Xe and 127 Xe , and double-electron-capture decays of 124 Xe , are backgrounds in searches for weakly interacting massive particles (WIMPs) conducted by dual-phase xenon time projection chambers such as LUX-ZEPLIN (LZ). These decays produce signals with more light and less charge than equivalent-energy 𝛽 decays and correspondingly overlap more with WIMP signals. We measure three electron-capture charge yields in LZ: the 1.1 keV M-shell, 5.2 keV L-shell, and 33.2 keV K-shell at drift fields of 193 and 96.5 V/cm. The LL double-electron-capture decay of 124 Xe exhibits even more pronounced shifts in charge and light. We provide a first model of double-electron-capture charge yields using the link between ionization density and electron-ion recombination, and identify a need for more accurate calculations. Finally, we discuss the implications of the reduced charge yield of these decays and other interactions creating inner-shell vacancies for future dark matter searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Surface heat flux and its association with the MJO in the tropical western Pacific using ARM observations

The Madden-Julian oscillation (MJO) is a major atmospheric phenomenon in the tropics that moves eastward every 20 to 100 days. It brings heavy rain and strong winds, influencing extreme weather events far beyond the tropics – including flooding, hurricanes, tornadoes and heavy snow in the United States. Therefore, an improved understanding of the MJO is critical for enhancing weather forecasts and supporting better decision-making for communities, emergency managers, and the private sector. Past studies, based on short-term observations or long-term model data, have emphasized the dominant role of atmospheric humidity in driving the MJO. Our research, using long-term observations (2000-2014) from three U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) sites in the tropical western Pacific (Manus, Nauru, and Darwin) confirms those past findings. However, we also show that atmospheric temperature, especially in the mid and upper troposphere, also plays a key role in the MJO’s evolution from its quiet (suppressed) to convective (active) phase. These findings provide valuable insights into how the MJO evolves and may help evaluate and improve weather and climate models. Our study also examined how the exchange of heat between the Earth's surface and the atmosphere, called surface heat flux, interacts with the MJO. This flux includes components such as latent heat (related to evaporation), sensible heat (related to temperature difference between surface and atmosphere), and radiation. We found that modulation in MJO convection is well connected to surface heat flux over the tropical western Pacific. One surface heat flux component often overlooked is the sensible heat flux due to precipitation (Q P ). Since falling raindrops are typically cooler than the land or ocean surface, they can cool the surface and affect local weather. To test the impact of Q P on convection during MJO, we incorporated it into a weather model, and ran simulations for two periods: April 2009 (when the MJO was active) and June 2006 (when it was not). Including Q P reduced errors in simulating the daily peak of rainfall—by 83% during the strong MJO and 23% during the inactive phase. It also improved the timing of peak daily rainfall and reduced the overall precipitation error by about 10%. These improvements were especially noticeable during periods of heavy rainfall. Our results suggest that even small heat exchanges from rainfall can play a significant role in shaping local weather. Accounting for these effects, particularly in regions like the tropical islands, can lead to more accurate simulations of the MJO-associated precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multi-GPU porting of a phase-change cascaded lattice Boltzmann method for three-dimensional pool boiling simulations

The Lattice Boltzmann method (LBM) has proven effective in simulating phase-change phenomena, such as melting, solidification, evaporation, and boiling. In this work, we develop a highly parallelized multi-GPU implementation of LBM for three-dimensional pool boiling simulations. The code is based on the OpenACC programming model, which enables the code to be deployed efficiently on multi-core CPUs, GPUs, and potentially other accelerators, without the need for architecture-specific rewrites. To support large-scale simulations, the domain is decomposed and distributed across multiple compute nodes using MPI. We demonstrate that the code exhibits excellent scaling properties, with ideal strong-scaling running with up to 256 GPUs on the MareNostrum5 cluster.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An underestimated negative cloud feedback from cloud lifetime changes

As the atmosphere warms, part of the cloud population shifts from ice and mixed phase (“cold”) to liquid (“warm”) clouds. Because warm clouds are more reflective and longerlived, this phase change reduces the solar flux absorbed by Earth and constitutes a negative radiative feedback. This cooling feedback is weaker in CMIP6 than in CMIP5, contributing to greater greenhouse warming. While this change is often attributed to improvements in simulated cloud phase, another model bias persists: warm clouds precipitate too readily, potentially leading to underestimated negative lifetime feedbacks. In this work, we modify a climate model to better simulate warm rain probability and find it exhibits a cloud lifetime feedback nearly three times larger than the default model. This suggests model errors in cloud-precipitation processes may bias cloud feedbacks as much as CMIP5–CMIP6 climate sensitivity differences. Reliable climate model projections therefore require improved cloud process realism guided by process-oriented observations and observational constraints.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Development of Conceptual Lead Cartridge Design to Perform Irradiation Experiments in VTR

This paper seeks to introduce the latest design of the Extended Length Test Assembly–Cartridge Lead (ELTA-CL) with associated thermal-hydraulic (TH) assessment and related experiment activities to support the critical component development performed by the ELTA-CL team (Los Alamos National Laboratory, Westinghouse Electric Company, and the University of New Mexico). The goal of the ELTA-CL program is to develop and validate an experimental capability to perform irradiation experiments in the Versatile Test Reactor (VTR) addressing Lead Fast Reactor (LFR) technology gaps, in support of the commercial development of advanced lead-cooled fast reactor concepts. Through a design maturation process and parametric study, a conceptual design is proposed to meet the requirements for material and corrosion testing. Thermal-hydraulic characteristics for the conceptual design at desired operating conditions are assessed with systems-level (one-dimensional) and computational fluid dynamics (three-dimensional) simulations. Along with the conceptual design work, experimental activities for the development of critical components such as the pump and flowmeter are undertaken. From both the modeling study and the experimental results, the design requirements of the Phase 1 ELTA-CL (e.g., 500°C and 2 m/s) are achievable with the current conceptual design. Additional design improvements and safety assessments at both steady-state and transient conditions for the final ELTA-CL design will be pursued.

42 ENGINEERING↗
DOE OSTI2021

Frustrated Kondo impurity triangle: A simple model of deconfinement

We report that the concepts of deconfinement and topological order are of great current interest for quantum information science and for our understanding of quantum materials. Here, we introduce a simple model of three antiferromagnetically coupled Kondo impurities, a Kondo triangle, which can be used to further extend the application of these concepts to electronic systems. We show that, by tuning the magnetic frustration, the Kondo triangle undergoes a quantum phase transition between two phases of unbroken symmetry, signaling a phase transition beyond the Landau paradigm. We demonstrate that the frustrated spin liquid phase is described by a three-channel Kondo (3CK) fixed point and thus displays an irrational ground state degeneracy. Using an Abrikosov pseudofermion representation, this quantum state is categorized by an emergent U(1) gauge field and its projective symmetry group. The gauge theory is deconfining in the sense that a miniature Wilson loop orders and topological defects (instantons in the gauge field) are expelled. This phase persists in the presence of moderate Kondo screening until proliferation of topological defects leads to a quantum phase transition to an unfrustrated Fermi liquid phase. Based on this evidence, we propose that the 3CK phase displays topological order in a similar sense as gapless spin liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Numerical Study on the Energy Performance of a Novel Furnace With Acidic Gas Trap Absorbers

Natural gas furnaces are widely used in US residential and commercial building markets. An important issue for natural gas furnaces is serious corrosion and fouling problems caused by acidic gas, such as SOx. An advanced adsorption technology based on acidic gas trap (AGT) absorbers offers the possibility to remove SOx acidic gas from natural gas furnaces with high efficiency and low cost, thereby enabling the development of condensing furnaces without the use of expensive corrosion resistant materials in the heat exchanger. A three-dimensional (3D) computational fluid dynamics (CFD) model has been developed to evaluate the heat transfer performance of a furnace with AGT absorbers and to compare it with a baseline conventional furnace without the AGT. Moreover, an axisymmetric model has been built focusing on the absorbing process in the AGT. The baseline conventional furnace used for the study is a commercial condensing furnace (Rheem 92% AFUE 84,000 BTU Multi-Position Gas Furnace). This furnace was completely disassembled, and the dimensions of each part were carefully measured and used to build a detailed CFD model. A model representing the new furnace, incorporating the AGT absorbers, was developed by adding the AGT system to the conventional furnace model. For the CFD analysis, a mixture model was employed to characterize the heat and mass transfer during the condensing process in the furnace while considering three components—air, water vapor and liquid water. Condensation takes place in the condensing heat exchanger, where water vapor changes phase to liquid water, and the latent heat is thus used in the furnace for useful heating. The simulation results characterize the energy performance of both the conventional furnace and the novel furnace with AGT absorbers, as well as the reactive processing in the AGT. These results provide insightful guidance for the development of the AGT absorber-based furnace from the perspective of its energy performance and will be used to further optimize this novel furnace design.

Laclair, Tim↗
DOE OSTI2021

Improving Data Center Energy Efficiency Through End-to-End Cooling Modeling and Optimization. Final Report

This report summarizes the development and demonstration of the three features developed for the DOE project “Improving Data Center Energy Efficiency through End-to-End Cooling Modeling and Optimization” (Award No. DE-EE0007688). The development of the three software features, Feature 1: Optimization of Cooling System Operation, Feature 2: Optimization of Airflow Management System, and Feature 3: Simultaneous Optimization of Airflow Management and Cooling Systems, are described first. Feature 1 is implemented in Phase 1 of the project and demonstrates 15% energy savings from cooling system optimization by applying the proposed retrofit solutions. Feature 2 is implemented in Phase 2 and the assessments show cooling energy-saving potential between 19.3% and 54.0%, depending on the energy-saving measures employed. In addition, investigating the energy saving potential for a hypothetical simple data center shows about 6.8% energy savings by optimizing the server load placement. Feature 3 is implemented in Phase 3 and the results show energy savings of 53.4% and 74.6% for the two partner data centers, respectively, compared to the baselines defined in Phase 1. Lastly, the project outcomes including software releases and publications are discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Data-Driven Recommendation of Optimal Tuning Scheme for Range-Separated Hybrid Functionals in Solution-Phase UV/Vis Absorption Energy Prediction

Time-dependent density functional theory (TDDFT) combined with range-separated hybrid (RSH) functionals and a tuned range-separation parameter γ offers a computationally economical approach for high-throughput excited- state property predictions. The γ-tuning procedure in the gas phase is well established. However, no agreement on the best γ- tuning procedure has been made when considering the solvent effect with implicit solvent models like the polarizable continuum model (PCM). To answer that question, this study created a diverse dataset with 937 molecules with experimental solutionphase UV/vis absorption spectra. Three γ-tuning methods, the gasphase γ-tuning (GPγT), the partial vertical γ-tuning (PVγT), and the strict vertical γ-tuning (SVγT), were evaluated for the ωPBEh functional over the entire dataset. Additional benchmarks are done for the optimally tuned screened range-separated hybrid combined with the PCM approach (SRSH-PCM) and the solvation-mediated tuning procedure (sol-med-OT). Our findings revealed that the optimal γ-values obtained by the PVγT and the SVγT are significantly smaller than the GPγT. This trend holds consistently across all molecules in our dataset, and we explained the origin of this phenomenon. TDDFT calculations with PVγTand SVγT-tuned γ-values and default global Fock exchange fraction achieve superior performance compared to those using GPγTtuned or default γ and slightly outperform SRSH-PCM and sol-med-OT with similar or lesser computational cost. Furthermore, we found that the smaller γ-values from SVγT captured the expected 1/(εR) asymptotic behavior in the solution phase, resulting in accurate prediction of solution-phase CT excitations, consistent with the screened asymptote behavior encoded in SRSH-PCM. These results show that SVγT is the best scheme for high-throughput UV/vis absorption spectrum calculations using the ωPBEh functional from a data-driven perspective.

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

Faster-than-real-time Simulation with Demonstration for Resilient DER Integration

The US electric grid is facing operational, stability, and security challenges. Transmission system operators need some measure of visibility into distribution system renewable generation. Distribution system generation needs to support transmission system voltage. The grid is experiencing an expansion in measurement systems. How to take full advantage of this expansion and defend against attacks, both cyber and physical, poses additional challenges. The Faster-than-real-time Simulation with demonstration for Resilient DER Integration project set out to do the following: a. Flatten the voltage profile through the feeders and system for cost saving and voltage stabilization needs. b. Increase the amount of intermittent distributed energy resources (IDERs) that could be deployed on a utility feeder and provide 100% or more energy needed for the demands on that feeder, and c. based on an accurate model (Digital Twin) of the utilities system, be able to detect any abnormalities on the utilities distribution system. To manage the voltage and increase IDER penetration (a,b), Graph Trace Analysis is employed in a time-series, optimal power flow to coordinate the time-varying feedback control setpoints of a distribution feeder’s utility control devices. Under the coordinated control are a Load Tap Changing Transformer, a voltage regulator, and five switched capacitor banks. The feeder serves over 2000 customers, the feeder secondaries are modeled, and the feeder has 2.3 MW of PV generation, corresponding to a 17.4% penetration of PV generation. The feeder model has over 12,000 components, where every customer load bus and PV generator are modeled. The accuracy of the power flow solution is compared against historical meter voltage measurements, the improvement in conservation voltage reduction energy savings as a function of the coordinated control desired voltage profile is investigated, and the increase in PV penetration of the coordinated control over the existing control is presented. To achieve improved control performance while observing system operation constraints, bellwether Advanced Metering Infrastructure (AMI) voltage measurements are used to adjust the desired voltage profile used by the optimal power flow analysis. To detect and alleviate or negate attacks or failures on the distribution and transmission utility grids (c) the grid needs to be resilient and self-healing. In this project software was designed to do just that. At the center of the software is an Integrated System Model (ISM) that spans from transmission to secondary distribution. The ISM is employed in real-time abnormality detection, voltage stability forecasting, and multi-mode control. Testing results are presented for: 1—attacks on utility infrastructure; 2—energy savings from optimal control; 3—distribution system control response during a low voltage transmission system event; 4—cyber-attacks on PV inverters, where physical inverters are used in hard-ware-in-the-simulation-loop studies. Contributions of this work include real-time analysis that spans from three-phase transmission through secondary distribution; an approach for detecting abnormalities that employs measurements from three independent measurement systems; and a multi-mode distribution system control that responds to cyber-attacks, physical attacks, equipment failures, and transmission system needs.

Integrated System Model, Graph Trace Analysis, Adv↗