Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scenario optimization”

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 37 records · Page 2

Long duration battery sizing, siting, and operation under wildfire risk using progressive hedging

Battery sizing and siting problems are computationally challenging due to the need to make long-term planning decisions that are cognizant of short-term operational decisions. This paper considers sizing, siting, and operating batteries in a power grid to maximize their benefits, including price arbitrage and load shed mitigation, during both normal operations and periods with high wildfire ignition risk. Here we formulate a multi-scenario optimization problem for long duration battery storage while considering the possibility of load shedding during Public Safety Power Shutoff (PSPS) events that de-energize lines to mitigate severe wildfire ignition risk. To enable a computationally scalable solution of this problem with many scenarios of wildfire risk and power injection variability, we develop a customized temporal decomposition method based on a progressive hedging framework. Extending traditional progressive hedging techniques, we consider coupling in both placement variables across all scenarios and state-of-charge variables at temporal boundaries. This enforces consistency across scenarios while enabling parallel computations despite both spatial and temporal coupling. The proposed decomposition facilitates efficient and scalable modeling of a full year of hourly operational decisions to inform the sizing and siting of batteries. With this decomposition, we model a year of hourly operational decisions to inform optimal battery placement for a 240-bus WECC model in under 70 min of wall-clock time.

25 ENERGY STORAGE↗

Scalable Pattern Matching in Metadata Graphs via Constraint Checking

Pattern matching is a fundamental tool for answering complex graph queries. Unfortunately, existing solutions have limited capabilities: They do not scale to process large graphs and/or support only a restricted set of search templates or usage scenarios. Moreover, the algorithms at the core of the existing techniques are not suitable for today’s graph processing infrastructures relying on horizontal scalability and shared-nothing clusters, as most of these algorithms are inherently sequential and difficult to parallelize. In this article we present an algorithmic pipeline that bases pattern matching on constraint checking. The key intuition is that each vertex and edge participating in a match has to meet a set of constraints implicitly specified by the search template. These constraints can be verified independently and typically are less expensive to compute than searching the full template. The pipeline we propose generates these constraints and iterates over them to eliminate all the vertices and edges that do not participate in any match, thus reducing the background graph to a subgraph that is the union of all template matches—the complete set of all vertices and edges that participate in at least one match. Additional analysis can be performed on this annotated, reduced graph, such as full match enumeration, match counting, or computing vertex/edge centrality. Furthermore, a vertex-centric formulation for constraint checking algorithms exists, and this makes it possible to harness existing high-performance, vertex-centric graph processing frameworks. This technique (i) enables highly scalable pattern matching in metadata (labeled) graphs; (ii) supports arbitrary patterns with 100% precision; (iii) enables tradeoffs between precision and time-to-solution, while always selects all vertices and edges that participate in matches, thus offering 100% recall; and (iv) supports a set of popular data analytics scenarios. We implement our approach on top of HavoqGT, an open-source asynchronous graph processing framework, and demonstrate its advantages through strong and weak scaling experiments on massive scale real-world (up to 257 billion edges) and synthetic (up to 4.4 trillion edges) labeled graphs, respectively, and at scales (1,024 nodes / 36,864 cores), orders of magnitude larger than used in the past for similar problems. This article serves two purposes: First, it synthesises the knowledge accumulated during a long-term project. Second, it presents new system features, usage scenarios, optimizations, and comparisons with related work that strengthen the confidence that pattern matching based on iterative pruning via constraint checking is an effective and scalable approach in practice. The new contributions include the following: (i) We demonstrate the ability of the constraint checking approach to efficiently support two additional search scenarios that often emerge in practice, interactive incremental search and exploratory search. (ii) We empirically compare our solution with two additional state-of-the-art systems, Arabsque and TriAD. (iii) We show the ability of our solution to accommodate a more diverse range of datasets with varying properties, e.g., scale, skewness, label distribution, and match frequency. (iv) We introduce or extend a number of system features (e.g., work aggregation, load balancing, and the ability to cap the generated traffic) and design optimizations and demonstrate their advantages with respect to improving performance and scalability. (v) We present bottleneck analysis and insights into artifacts that influence performance. (vi) We present a theoretical complexity argument that motivates the performance gains we observe.

97 MATHEMATICS AND COMPUTING↗

Hierarchical Gaussian process-based Bayesian optimization for materials discovery in high entropy alloy spaces

Bayesian optimization (BO) is a powerful and data-efficient method for iterative materials discovery and design, particularly valuable when prior knowledge is limited, underlying functional relationships are complex or unknown, and the cost of querying the materials space is significant. Traditional BO methodologies typically utilize conventional Gaussian Processes (cGPs) to model the relationships between material inputs and properties, as well as correlations within the input space. However, cGP-BO approaches often fall short in multi-objective optimization scenarios, where they are unable to fully exploit correlations between distinct material properties. Leveraging these correlations can significantly enhance the discovery process, as information about one property can inform and improve predictions about others. Here, this study addresses this limitation by employing advanced kernel structures to capture and model multi-dimensional property correlations through multi-task (MTGPs) or deep Gaussian Processes (DGPs), thus accelerating the discovery process. We demonstrate the effectiveness of MTGP-BO and DGP-BO in rapidly and robustly solving complex materials design challenges that occur within the context of complex multi-objective optimization over FCC FeCrNiCoCu high entropy alloy (HEA) spaces, where traditional cGP-BO approaches fail. Furthermore, we highlight how the differential costs associated with querying various material properties can be strategically leveraged to make the materials discovery process more cost-efficient.

36 MATERIALS SCIENCE↗

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSTX-U theory, modeling and analysis results

Here, the mission of the low aspect ratio spherical tokamak NSTX-U is to advance the physics basis and technical solutions required for optimizing the configuration of next-step steady-state tokamak fusion devices. NSTX-U will ultimately operate at up to 2 MA of plasma current and 1 T toroidal field on axis for 5 s, and has available up to 15 MW of neutral beam injection power at different tangency radii and 6 MW of high harmonic fast wave heating. With these capabilities NSTX-U will develop the physics understanding and control tools to ramp-up and sustain high performance fully non-inductive plasmas with large bootstrap fraction and enhanced confinement enabled via the low aspect ratio, high beta configuration. With its unique capabilities, NSTX-U research also supports ITER and other critical fusion development needs. Super-Alfvénic ions in beam-heated NSTX-U plasmas access energetic particle (EP) parameter space that is relevant for both α-heated conventional and low aspect ratio burning plasmas. NSTX-U can also generate very large target heat fluxes to test conventional and innovative plasma exhaust and plasma facing component solutions. This paper summarizes recent analysis, theory and modelling progress to advance the tokamak physics basis in the areas of macrostability and 3D fields, EP stability and fast ion transport, thermal transport and pedestal structure, boundary and plasma material interaction, RF heating, scenario optimization and real-time control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nature-inspired lotus-shaped fins combined with hybrid nanoparticles and metal foam for high-performance latent heat thermal energy storage

Latent heat thermal energy storage (LHTES) systems play a critical role in renewable energy integration by providing high energy density and nearly isothermal operation during phase transitions. However, their performance is often limited by slow melting/charging rates, which motivates the search for enhanced heat transfer designs. This study investigates the melting behavior of RT-82 phase change material (PCM) using novel lotus-shaped fins combined with copper metal foam and conductive graphene nanoparticles and carbon nanotubes. A two-dimensional enthalpy-porosity model in ANSYS Fluent was developed to simulate the charging/melting process, capturing non-thermal equilibrium between the foam and PCM/nano-PCM. In this study, effects of fin geometry, nanoparticle concentration, and foam porosity on melting dynamics and cost-performance trade-offs were investigated. Results showed that natural convection accelerated melting by ~12% compared to conduction-only scenarios. Optimized lotus-shaped fins with higher fin density (T3F4 and T3F10) achieved up to 63% faster melting relative to sparse configurations. Graphene nanoparticles improved thermal conductivity, with a 6% volume fraction, by reducing melting time by ~6.9%, while their combination with 75% porosity foam achieved a maximum reduction in the melting time of ~51% compared to pure PCM. Cost-performance analysis identified T3F4 as the most balanced design, offering rapid thermal response without excessive material costs, while moderate-density designs like T3S6 provided economical alternatives with acceptable performance. These results highlight the performance enhancement that can be achieved by integrating bio-inspired fins, nanoparticles, and foams, into compact and efficient LHTES for solar heating, building thermal management, and industrial waste-heat recovery applications.

25 ENERGY STORAGE↗

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗

Measurement of helicon waves with phase contrast imaging on DIII-D – A theoretical feasibility study

A DIII-D high-beta H-mode discharge, with I p = 850 kA, B t = 2.1 T n e = 4 10 19 m –3 , has been designed to validate full wave modeling of helicon waves by optimizing the expected response of the Phase Contrast Imaging diagnostic. Helicon waves have been predicted to have high current drive efficiency off-axis without facing the accessibility issues of lower-hybrid waves. To test these predictions experimentally, DIII-D has recently commissioned a high-power helicon antenna. To confidently predict the behavior of helicon waves in future devices, measurements of their fundamental properties and validation against models will be essential. Phase contrast imaging (PCI) is an absolutely calibrated internal reference interferometer able to measure density fluctuations with radial wavenumbers k R between 1.5 cm –1 and 20 cm –1 . For helicon waves 2 cm –1 < k R < 10 cm –1 is expected, allowing PCI to measure their envelope and wavenumber spectrum. This makes PCI a powerful tool for the validation of state-of-the-art models, like the AORSA full wave code. AORSA is used to compute the density perturbations measured by the PCI with 2D calculations corresponding to 11 different toroidal mode numbers combined to resolve the trajectory of the helicon wave in 3D. This is necessary because the waves travel 120 degrees toroidally from the antenna to the PCI. Here, a cold plasma finite element model (CPFEM) [4] is used to predict propagation through the scrape-off layer. The result of the CPFEM model is connected to AORSA by creating an artificial Gaussian antenna on the last closed flux surface. PCI shows best results for waves with small vertical wavenumbers k z . Modeling the helicon waves for several past DIII-D experiments shows that k z is minimized if the intersection of the helicon and the PCI laser beams occurs in the midplane. For such an optimized scenario the predicted signal level is two orders of magnitude larger than the background density fluctuations arising from broadband turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental study of the influence of gas puff locations on H-mode boundary plasmas with argon seeding on EAST

To investigate the optimal scenario of impurity seeding to obtain divertor plasma detachment for target protection, experiments with Ar&D 2 seeding from two different poloidal locations, the upper outer (UO) divertor target and lower outer (LO) target, were carried out on EAST. Partial energy detachment (the electron temperature near the strike point T e,spt ≤ 10 eV) were obtained with Ar&D 2 mixture puffing from the UO target and, for the first time, from the LO target into H-mode plasmas in the upper single null (USN) configuration. The peak heat flux q t on the UO target was significantly reduced (by ~80%). The rollover of ion flux density j s did not appear probably due to insufficient momentum loss, which is independent of the puff locations. The poloidal asymmetries of particle and heat fluxes on the targets have also been investigated. The UO-dominant asymmetry of particle flux was reversed, while the UO-dominant asymmetry of heat flux was mitigated but not reversed. The plasma confinement dropped by only 14% and 8.2% in the LO-puff case and UO-puff case, respectively, during detachment. The high level of C may contributed greatly to the higher radiation in the bulk plasma region and the greater decline in W MHD in the LO-puff case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotope dependence of transport in ST40 hot ion mode plasmas

The ST40 compact, high-field spherical tokamak, operating at 2.1 T and with 1.8 MW of neutral beam heating, achieved central carbon impurity ion temperatures in excess of 10 keV, surpassing their business milestone of 100 M °C (8.6 keV). The high temperature discharges were in the hot ion mode, with $T_i \gg T_e$, and they were achieved in both hydrogen and deuterium plasmas with deuterium neutral beam injection. In order to achieve the high temperature scenarios, careful wall conditioning and scenario optimization was carried out in $H^0 \rightarrow H^+$, $D^0 \rightarrow H^+$, and finally in $D^0 \rightarrow D^+$ plasmas. The TRANSP transport code was employed to study the dependence of confinement and transport on isotopic mass, and the conditions that led to the high measured central ion temperatures in the $D^0 \rightarrow H^+$ and $D^0 \rightarrow D^+$ plasmas. The kinetic profiles input into TRANSP were inferred from line-of-sight and limited radial measurements as well as consistency with a number of other experimental constraints. The TRANSP results first showed that the main species central ion temperature was only 0.5–1 keV lower than the measured carbon impurity temperature, and thus in the high performance plasmas also surpassed the 100 M °C level. TRANSP also showed that while the electron thermal conduction loss was dominant, reductions in central ion transport, and an effective decoupling of the ions from the electrons, led to an increase in central ion temperatures with increasing plasma mass. In fact, the associated confinement times exhibited a strong dependence on the isotopic mass of the thermal plasma. These preliminary results are foundations for dedicated experiments with full kinetic profile measurements on ST40 in the next run campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predicting core transport in ITER baseline discharges with neon injections

Achieving self-consistent performance predictions for ITER requires integrated modeling of core transport and divertor power exhaust under realistic impurity conditions. We present results from a systematic power-flow and impurity-content study for the ITER 15 MA baseline scenario constrained directly by existing SOLPS-ITER neon-seeded divertor solutions. Using the OMFIT STEP workflow, stationary temperature and density profiles are predicted with TGYRO for $1.5 \unicode{x2A7D} Z_\textrm{eff} \unicode{x2A7D} 2.5$, and the corresponding power crossing the separatrix $P_\textrm{sep}$ is evaluated. We find that $P_\textrm{sep}$ varies by more than a factor of 1.7 across this scan and matches the ${\sim}100$ MW SOLPS-ITER prediction when $Z_\textrm{eff} \simeq 1.6$ or when auxiliary heating is reduced to ${\sim}75\%$ of nominal. Rotation-sensitivity studies show that plausible variations in toroidal flow magnitude modify $P_\textrm{sep}$ by $\lesssim 20\%$, while AURORA modeling confirms that charge-exchange radiation inside the separatrix is dynamically negligible under predicted ITER neutral densities. These results identify a restricted compatibility window, $Z_\textrm{eff} \approx 1.6$ –1.75 and $0.75 \lesssim f_{P_\textrm{aux}} \unicode{x2A7D} 1.0$, in which core transport predictions remain aligned with neon-seeded divertor protection targets. This self-consistent, model-constrained framework provides actionable guidance for impurity control and auxiliary-heating scheduling in early ITER operation and supports future whole-device scenario optimization.

ITER↗

Predictive modeling of NSTX discharges with the updated multi-mode anomalous transport module

Abstract The objective of this study is twofold: firstly, to demonstrate the consistency between the anomalous transport results produced by updated Multi-Mode Model (MMM) version 9.0.4 and those obtained through gyrokinetic simulations; and secondly, to showcase MMM’s ability to predict electron and ion temperature profiles in low aspect ratio, high beta NSTX discharges. MMM encompasses a range of transport mechanisms driven by electron and ion temperature gradients, trapped electrons, kinetic ballooning, peeling, microtearing, and drift resistive inertial ballooning modes. These modes within MMM are being verified through corresponding gyrokinetic results. The modes that potentially contribute to ion thermal transport are stable in MMM, aligning with both experimental data and findings from linear CGYRO simulations. The isotope effects on these modes are also studied and higher mass is found to be stabilizing, consistent with the experimental trend. The electron thermal power across the flux surface is computed within MMM and compared to experimental measurements and nonlinear CGYRO simulation results. Specifically, the electron temperature gradient modes (ETGM) within MMM account for 2.0 MW of thermal power, consistent with experimental findings. It is noteworthy that the ETGM model requires approximately 5.0 ms of computation time on a standard desktop, while nonlinear CGYRO simulations necessitate 8.0 h on 8 K cores. MMM proves to be highly computationally efficient, a crucial attribute for various applications, including real-time control, tokamak scenario optimization, and uncertainty quantification of experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗

Radiation-front position control with scenario adjustments during the first tungsten divertor operation on KSTAR

This study reports the first experimental demonstration of a radiation front-based detachment control system during the initial tungsten divertor campaign in KSTAR. The system employs real-time infrared video bolometry with tomographic reconstruction to track the inboard radiation-front position, which served as a feedback variable for impurity seeding. In a baseline scenario without pre-D 2 fueling, tungsten accumulation triggered core radiation surges and repetitive H–L back transitions. In contrast, with pre-D 2 fueling, detachment was stabilized by maintaining higher edge-localized mode frequency, suppressing core tungsten buildup, and reducing sputtering, thereby enabling sustained feedback control. The suitability of the inboard front position as a detachment control metric was confirmed by its correlation with the degree of detachment: with D 2 fueling the relation exhibited a broad hysteresis loop characteristic of momentum-loss-dominated detachment, whereas with N 2 fueling, it became narrow and nearly vertical, consistent with radiation-loss-dominated pinning. These results provide a proof-of-principle for radiation-front-based detachment control, demonstrating its viability as a control approach. At the same time, these results clarify that the detachment outcomes were further conditioned by impurity transport characteristics of tungsten divertor operation in KSTAR, providing guidance for future scenario optimization and controller development in tungsten devices.

KSTAR↗

Virtual Tokamak for Test and Development of Plasma Control Applied to NSTX-U

Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deconstructing laws of accessibility and facility distribution in cities

The era of the automobile has seriously degraded the quality of urban life through costly travel and visible environmental effects. A new urban planning paradigm must be at the heart of our road map for the years to come, the one where, within minutes, inhabitants can access their basic living needs by bike or by foot. In this work, we present novel insights of the interplay between the distributions of facilities and population that maximize accessibility over the existing road networks. Results in six cities reveal that travel costs could be reduced in half through redistributing facilities. In the optimal scenario, the average travel distance can be modeled as a functional form of the number of facilities and the population density. As an application of this finding, it is possible to estimate the number of facilities needed for reaching a desired average travel distance given the population distribution in a city.

42 ENGINEERING↗

Hybrid Power Plants for Energy Resilience: A Case Study

As renewable energy technologies are increasingly adopted, they pose an opportunity to improve the sustainability and resilience of distributed grids, especially when their design and operation is coordinated as a hybrid power plant. When included in hybrid power plants, distributed wind turbines in particular have the potential to enhance the resilience of distributed grids in areas with good wind resource, due to their ability to provide more consistent generation and ancillary services as compared to photo-voltaic (PV) solar panels. Despite this benefit, U.S. distributed wind adoption is lower than other comparable renewable energy technologies. In this study, we seek to demonstrate how hybrid power plants that include distributed wind turbines can contribute to distribution grid resilience by meeting loads (especially critical loads) more consistently, increasing reserve capacity, and providing value to customers during outages. To demonstrate these contributions, we integrate three separate frameworks and apply them to a case study in a rural electric cooperative in Iowa. Through this case study, we simulate and compare hybrid power plant design and operation during two hazard events: a tornado that causes a 48-hour distribution outage and a winter weather event that causes a 6-hour generation outage. The inclusion of a hybrid power plant that leverages 1) increased battery duration and 2) advanced forecasting and dispatch strategies that reserve capacity leading up to a hazard event best reduce lost loads as well as diesel consumption that would otherwise be used to meet those loads during short- and long-duration hazard events. Depending on the hybrid power plant capacity and operation, we find that the outage mitigation value of a hybrid power plant (measured in value to customers to avoid an outage and avoided lost revenues for the utility) is significant in both hazard events; adding wind, solar, and battery assets to the existing system adds about $50-$100M in avoided lost load and at least $4-$8k in utility value in the tornado hazard event, and $570k-$2.2M in avoided lost load and at least $220-$650 in utility value in the winter hazard scenario. In both the tornado and winter hazard scenarios, optimizing the operation of the hybrid system for resilience can lend similar value as increasing battery duration by 5 MWh for the lower capacity systems considered.

17 WIND ENERGY↗

Global particle buildup simulations with gas puff scan: application to WEST discharge

This paper deals with the distribution of sources, transport, and exhaust of particles in a tokamak. Knowledge and understanding of all the physical phenomena involved in the global particle buildup are necessary to study and predict density regimes and subsequently to develop optimized scenarios for tokamak operation in order to control heat and particle exhaust. Neutral particles and their interactions with plasma are central in this perspective. This paper discusses the impact of varying the intensity of particle fueling in 2D transport simulations of a WEST discharge. Simulations are performed with an updated version of SOLEDGE-HDG that allows a more realistic transport of neutrals using a self-consistent diffusive model based on charge exchange and ionization processes. New code capabilities allow the entire WEST poloidal cross section to be simulated in a realistic configuration for both geometry and the range of control parameters. A gas puff scan illustrates the main features of the sheath-limited, high-recycling, and detached regimes, such as the buildup of the temperature gradient and the pressure drop in the scrape-off layer (SOL), the target temperature falling to 1 eV, and the ionization source moving away from the targets, as well as the particle flux rollover. A crude estimate of wall erosion is also provided, showing the respective role of each plasma wall component in each of these regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗