Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Twinning”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Assessing Strain Rate Sensitivity of Nanotwinned Al–Zr Alloys through Nanoindentation

Nanotwinned metals have exhibited many enhanced physical and mechanical properties. Twin boundaries have recently been introduced into sputtered Al alloys in spite of their high stacking fault energy. These twinned Al alloys possess unique microstructures composed of vertically aligned Σ3(112) incoherent twin boundaries (ITBs) and have demonstrated remarkable mechanical strengths and thermal stability. However, their strain rate sensitivity has not been fully assessed. A modified nanoindentation method has been employed here to accurately determine the strain rate sensitivity of nanotwinned Al–Zr alloys. The hardness of these alloys reaches 4.2 GPa while simultaneously exhibiting an improved strain rate sensitivity. The nanotwinned Al–Zr alloys have shown grain size-dependent strain rate sensitivity, consistent with previous findings in the literature. This work provides insight into a previously unstudied aspect of nanotwinned Al alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Fast Growth of Ga2O3 on Highly Offcut (100)-Oriented Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. More importantly, the impact ionization coefficients along the [100] direction are low, leading to the highest possible critical fields. This is advantageous compared to the current state of the art, (001), due to reduced surface defects and increased possible breakdown voltage. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption, and on-axis (100) growth easily forms twin domains. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction. These films show step-flow growth from (20-1) step-edges and high electron mobility due to suppressed twins. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (>5 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3 and critical breakdown field in Schottky barrier diodes comparable with state-of-the-art (001) Ga2O3 without edge termination.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimizing Batch Crystallization with Model-based Design of Experiments

Adaptive and self-optimizing intelligent systems such as digital twins are increasingly important in science and engineering. Digital twins utilize mathematical models to provide added precision to decision-making. However, physics-informed models are challenging to build, calibrate, and validate with existing data science methods. Model-based design of experiments (MBDoE) is a popular framework for optimizing data collection to maximize parameter precision in mathematical models and digital twins. In this work, we apply MBDoE, facilitated by the open-source package Pyomo.DoE, to train and validate mathematical models for batch crystallization. We quantitatively examined the estimability of the model parameters for experiments with different cooling rates. This analysis provides a quantitative explanation for the heuristic of using multiple experiments at different cooling rates.

Lynch, Hailey↗
DOE OSTI2021

Geothermal Operational Optimization with Machine Learning

The Geothermal Operational Optimization with Machine Learning (GOOML) project has developed a generic and extensible component-based system modeling framework to study complex geothermal fields using a data-driven approach. Through building a digital twin of a geothermal steam field with the GOOML modeling framework, operators can analyze historical and forecasted power production, explore possible steam field configurations, and optimize real world operations, all in a cost-effective digital environment. The GOOML modeling software is based on a historical data-assimilation framework that uses first-principal thermodynamics to model steam field components using historical data, and a forecast framework that uses machine-learning-driven models of steam field components to predict future operations. This modeling framework creates countless new opportunities for digital exploration of steam field design and operations. To date, digital twins have been developed for several steam fields in New Zealand and the United States. These digital twins have been validated by comparing hindcast predictions against historical production data. Field design and operations have been explored using genetic optimization and reinforcement learning. Initial results show compelling and often surprising opportunities for improved design and operation of fields with 2 to 5 percent improvements in annual energy production. GOOML is driving a step-change in geothermal operations by applying state-of-the-art machine learning algorithms, comprehensive data analytics, and a first-of-its-kind intelligent geothermal systems model.

40 EE - Geothermal Technologies Office (EE-4G)↗
DOE OSTI2025

Validation of Pronghorn for Natural-Circulation Molten Salt Loops

This paper presents the development and validation of a high-fidelity thermal-hydraulic model of a molten salt natural circulation flow loop, designed for integration within a digital twin framework. The study evaluates the performance of Idaho National Laboratory’s Pronghorn against experimental data from Texas A&M University Molten Salt Flow Loop (MSFL) four Hitec-salt test benchmark data. Natural circulation of high-Prandtl-number fluids exhibits complex, counter-intuitive flow patterns that make pointwise thermocouple readings unreliable. Experimental work at TAMU’s MSFL provides benchmark data, including flow visualization at a test-section and centerline steady-state temperature measurements along the loop. Validation includes four single-phase natural circulation test cases with Hitec salt. Key metrics include flow profile agreement and steady-state temperature accuracy. Pronghorn results for two-dimensional single-phase agree qualitatively with the experimental flow profile. This paper illustrates the importance of Computational Fluid Dynamics (CFD) in elucidating the behavior of high-Prandtl-number thermal-hydraulics, along with how misleading centerline temperature measurements can be. Pronghorn reproduces the axial and radial stratification that makes single thermocouple readings unreliable. Future research will focus on reduced-order modeling techniques to enable rapid simulation suitable for real-time digital twin applications. The validated cases provide a basis for developing reduced-order surrogates aimed at real-time digital-twin applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Plastic deformation of single crystals of the equiatomic Cr-Mn-Fe-Co-Ni high-entropy alloy in tension and compression from 10 K to 1273 K

In this report the plastic deformation behavior of single crystals of the quinary, equiatomic Cr-Mn-Fe-Co-Ni high-entropy alloy (HEA) with the face-centered cubic structure has been investigated in tension and compression as a function of crystal orientation and temperature from 10 K to 1373 K. The critical resolved shear stress (CRSS) for {111}<110> slip at room temperature is 42-45 MPa. It does not depend much on crystal orientation (i.e., the Schmid law holds true) and the sense (tension vs. compression) of the applied load. The CRSS for {111}<110> slip increases with the decrease in temperature, without showing any significant inertia effects at cryogenic temperatures below 77 K. Extrapolation from the measured yield stresses down to 10 K yields a CRSS value at 0 K of 168 MPa. At cryogenic temperatures, the measured strain-rate sensitivity of flow stress is consistent with a very small activation volume. The concept of stress equivalence holds true both for the temperature dependence of CRSS and the stress dependence of activation volume, indicating that solid-solution hardening is the major strengthening mechanism. Deformation twinning occurs at 77 K but not at room temperature, resulting in higher tensile elongation to failure at 77 K than at room temperature. Deformation twinning at 77 K occurs at a shear stress of 378 MPa on conjugate ($\bar1\bar11$) planes in the form of Lüders deformation after large plastic strain (about 85%) achieved by the stage I (easy glide) and stage II (linear work-hardening) deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A nudged hybrid analysis and modeling approach for realtime wake-vortex transport and decay prediction

Here, we put forth a long short-term memory (LSTM) nudging framework for the enhancement of reduced order models (ROMs) of fluid flows utilizing noisy measurements for air traffic improvements. Toward emerging applications of digital twins in aviation, the proposed approach allows for constructing a realtime predictive tool for wake-vortex transport and decay systems. We build on the fact that in realistic application, there are uncertainties in initial and boundary conditions, model parameters, as well as measurements. Moreover, conventional nonlinear ROMs based on Galerkin projection (GROMs) suffer from imperfection and solution instabilities, especially for advection-dominated flows with slow decay in the Kolmogorov n-width. In the presented LSTM nudging (LSTM-N) approach, we fuse forecasts from a combination of imperfect GROM and uncertain state estimates, with sparse Eulerian sensor measurements to provide more reliable predictions in a dynamical data assimilation framework. We illustrate our concept by solving the two-dimensional vorticity transport equation. We investigate the effects of measurements noise and state estimate uncertainty on the performance of the LSTM-N behavior. We also demonstrate that it can sufficiently handle different levels of temporal and spatial measurement sparsity, and offer a huge potential in developing next-generation digital twin technologies for aerospace applications.

42 ENGINEERING↗
DOE OSTI2020

Exploiting ultra-large linear elasticity over a wide temperature range in nanocrystalline NiTi alloy

Many shape memory alloys can support large recoverable strains of a few percent by reversible stress-induced martensite transformation, yet they behave non-linear within a narrow operating temperature range. Developing the bulk metallic materials with ultra-large linear elasticity over a wide temperature range has proven to be difficult. In this work, a material design concept was proposed, that is true elastic deformation and reversible twinning-detwinning deformation run in parallel to overcome this challenge. By engineering the residual internal stress to realize the concurrency of true elastic deformation and twinning-detwinning deformation, a bulk nanocrystalline NiTi that possesses an ultra-large linear elastic strain up to 5.1 % and a high yield stress of 2.16 GPa over a wide temperature range of 270 °C was developed. This study offers a new avenue for developing the metallic materials with ultra-large linear elasticity over a wide temperature range of 270 °C (from 70 °C to -197 °C).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Comparing the deformation mechanism in 316 L stainless steel fabricated by hybrid and additive manufacturing

This study rationalizes the relatively higher elongation with minimal strength drop of 316L stainless steel (316LSS) manufactured using the hybrid approach in comparison to a fully additive approach followed by machining. It was found that twinning was the dominant deformation mechanism in 316LSS produced by both additive and hybrid approaches. Additionally, the relatively lower fraction of twins in the as-fabricated state of hybrid approach compared to a fully additive approach resulted in delaying the onset of steady-state strain hardening rate to higher strain levels, thereby resulting in an overall increase in elongation to failure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Cryogenic mechanical behaviors of CrMnFeCoNi high-entropy alloy

The CrMnFeCoNi high-entropy alloy (HEA) exhibits higher yield strength, ultimate strength and ductility at lower temperature. To further clarify the effect of the testing temperature on microstructure evolution, in this study in-situ synchrotron-based high-energy X-ray diffraction tensile tests were carried out from 298 K down to 123 K. The enhanced yield strength of the alloy at cryogenic temperatures can be attributed to the greater lattice distortion prior to plastic deformation. Higher strain hardening rate leads to the simultaneously enhanced strength and ductility of the studied HEA below room temperature. Both dynamic Hall-Petch hardening (twinning) and dislocation hardening provide high work hardening capacity for this alloy during the plastic deformation at cryogenic temperatures. The increased dislocation density and nano-twins at cryogenic temperatures can be attributed to the decrease in the stacking fault energy as the deformation temperature decreases. These studies could provide an in-depth understanding for the strengthening mechanisms of the HEA in different temperature conditions and guide the exploration of HEAs with superb mechanical properties at cryogenic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In situ investigation of the deformation behaviors of Fe 20 Co 30 Cr 25 Ni 25 and Fe 20 Co 30 Cr 30 Ni 20 high entropy alloys by high-energy X-ray diffraction

In situ synchrotron-based high-energy X-ray diffraction (HE-XRD) technique was employed to investigate the mechanical behaviors and microstructural evolution of face-centered cubic (FCC) Fe 20 Co 30 Cr 25 Ni 25 and Fe 20 Co 30 Cr 30 Ni 20 high entropy alloys (HEAs) during tensile deformation. Fe 20 Co 30 Cr 30 Ni 20 HEA has a good combination of strength (ultimate tensile strength of 864 ± 35 MPa) and ductility (elongation of 0.627 ± 0.021). The HE-XRD investigation reveals that Fe 20 Co 30 Cr 30 Ni 20 HEA has the transformation-induced plasticity (TRIP) effect, which starts at a critical stress of ~555 MPa. Furthermore, transmission electron microscopy confirmed this deformation-induced new phase is hexagonal-close-packed structured ε-martensite, which follows an orientation relationship of {111} γ //(0001) ε , and (110) γ //[11 2 ¯ 0] ε with the FCC γ matrix. The observation of deformation twins in the deformed samples of the two studied HEAs proves that twinning-induced plasticity (TWIP) effect occurs in both HEAs. The combination of TRIP and TWIP effects lead to the high strength, large ductility and improved strain hardening behavior of Fe 20 Co 30 Cr 30 Ni 20 HEA.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Orientation dependence of the effect of short-range ordering on the plastic deformation of a medium entropy alloy

Multi-principal-element alloys (also known as medium and high entropy alloys) offer a much larger and richer design space than conventional alloys, providing opportunities for discovering new functionalities and their governing physics. Some of these alloys exhibit an outstanding combination of high strength and ductility, linked to the activation of various deformation modes triggered by low-energy stacking faults. However, a pressing question remains: Is the plasticity of medium- and high-entropy alloys governed only by stacking fault energy, or does atomic short-range order (SRO) play a role? Here, to answer this, we investigated how SRO affects the deformation in single-crystalline NiCoCr, with previous contradictory findings. First, we established unique experimental evidence for SRO formation in bulk single crystals using high-energy synchrotron transmission X-Ray Diffraction. By tuning the degree of SRO by aging at high temperatures, twinning density and strain-induced martensitic phase transformation can be significantly increased in the [110] and [111] orientations under tension, increasing the tensile ductility; yet, no increase was observed along the [001] orientation due to lack of TWinning-Induced Plasticity (TWIP) or TRansformation-Induced Plasticity (TRIP), indicating a strong crystallographic orientation dependence. Our first-principles thermodynamic calculations unequivocally show SRO exists and governs the observed microstructural evolution and deformation hardening behavior. Here we find direct proof that SRO triggers a simultaneous TWIP and TRIP in NiCoCr, a rare microstructural evolution path. Our findings establish that the interplay of SRO and plasticity could be exploited to alter deformation modes and yield unprecedented mechanical response in medium- and high-entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rapid Solidification Effects in Additively Manufactured Si and SiGe Compositions

SiGe alloys have a proven track record as robust high-temperature thermoelectric materials, powering NASA missions like SNAP-10A, LES-9, and Voyager 1 and 2. Enhancing thermoelectric efficiency hinges on minimizing thermal conductivity while preserving electrical conductivity. Rapid solidification via LPBF can generate microstructural features such as subgrain cellular boundaries and twinning, which may help reduce thermal conductivity while preserving semiconducting behavior. Here, this study investigates the potential of laser powder bed fusion (LPBF) additive manufacturing to fabricate nanostructured Si and SiGe thermoelectric materials. High cooling rates (10 5 to 10 7 K/s) rates inherent to the LPBF process are conducive to forming such nanostructures. Moreover, this fabrication technique could also be suitable for fabricating complex geometries needed to achieve improved device level performance. Process mapping of commercial Si powder with irregular morphology was first performed to understand the LPBF processing behavior of this semiconductor material. Subsequent studies included in-house synthesized B-doped (p-type) Si7 8 Ge 22 spherical powder that was produced via ultrasonic atomization. Scan strategies involved multiple laser exposures to mitigate solidification cracking: a high density of >98% was achieved, but solidification cracking could not be fully eliminated. Subsequently, a high electrical resistivity (i.e., low conductivity) was observed, but the measured Seebeck coefficient, ∼230 μV/K @ 500 °C, proved that good semiconductor material was being fabricated. A subgrain cellular structure (5–10 μm) was observed as defined by Ge segregation to the intercellular boundaries. The remelting strategies helped lower the cooling rates in processing SiGe, but this still resulted in high residual stresses, which induced a remarkably high density of twins (78–95%) to accommodate the deformation. This unique grain structure offers an avenue for phonon scattering and potential improvements in thermoelectric performance.

figure of merit↗
DOE OSTI2022

Solutes that reduce yield strength anisotropies in magnesium from first principles

Using Labusch-type solid solution strengthening models parameterized with DFT-computed solute-dislocation interaction energies, we perform a computational search for 63 solutes across the periodic table to find those that lower anisotropy ratios (non-basal to basal CRSS) of magnesium potentially increasing its ductility per the von Mises criterion. For this purpose, we compute changes in strength for solutes as a function of composition and temperature, and compute anisotropy ratios for solutes that include both rare earth and non-rare earth elements. Here we specifically focus on solute-dislocation interaction energies in the following DFT-optimized dislocations as representative of three non-basal plastic deformation modes: $\langle c + a \rangle$ edge, (10$\bar1$2) tension twinning edge, and the (10$\bar{1}$1) compression twinning edge. We find that solute-induced changes in non-basal deformation modes can be approximated using a second-order polynomial in the size misfit of the solutes, which permits rapid screening of solutes. Our approach to identify solutes known to improve strengthening incorporates solute solubility, and suggests other solutes that not have been previously explored for strengthening. The 8 rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios at their optimal concentrations, are: Gd, Tb, Dy, Nd, Ho, Er, Tm, and Yb. The 12 non-rare-earth solutes that our method suggests as the best, ordered by increasing anisotropy ratios, are: Y, Mn, Sc, Pb, Ca, Ag, Bi, Tl, Zn, Li, Ga, and Al. Of these, Gd, Nd, Er, Yb, Y, Mn, Ca, Zn, Li, and Al are used in commercial Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structure and glide of Lomer and Lomer-Cottrell dislocations: Atomistic simulations for model concentrated alloy solid solutions

Lomer (L) and Lomer-Cottrell (LC) dislocations have long been considered to be central to work hardening in face-centered cubic (FCC) metals and alloys. These dislocations act as barriers of motion for other dislocations, and can serve as sites for twin nucleation. Recent focus on multicomponent concentrated FCC solid solution alloys has resulted in many reported observations of LC dislocations. While these and L dislocations are expected to have a role in the mechanical behavior of these alloys, little is understood about how variations in composition and associated fault energies change the response of these dislocations under stress. Here we present atomistic simulations of L and LC dislocations in a model Cu-Ni system and find that changes in composition and applied stress conditions result in a wide variety of responses, including changes in core configuration and (100) glide. The results are compared to and extend previous literature related to the nature of L/LC core structures and how they vary with respect to intrinsic materials properties and stress states. This study also provides insights into mechanisms such as twin nucleation that could have important implications for work hardening in FCC solid-solution alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

DEVELOPMENT AND DEMONSTRATION TESTBED FOR THE REMOTE OPERATIONS AND MONITORING OF MICROREACTORS

The nuclear industry is rapidly developing many advanced-reactor concepts for near-term deployment in both traditional and non-traditional nuclear-powered applications. One such category of advanced reactor is the microreactor, a class of reactor with less than 20MWth power output, intended for applications where the economics or logistics of traditional power sources are difficult. This includes applications such as remote communities, mining sites, defense installations, or humanitarian and disaster-relief missions. One key enabling feature for the successful deployment of microreactors is a remote operations capability. Remote operations provide monitoring and control capabilities which can significantly reduce staffing costs by eliminating the need for licensed operators at each reactor facility and improve the economic viability for microreactor deployment. A remote concept of operations is not currently an established capability in the nuclear industry. In addition, no demonstration microreactor is expected to complete construction or go critical until at least 2026. This leaves two major capability gaps: the successful demonstration of a remote concept of operations for microreactors and a test bed suitable for said demonstration. Both gaps must be addressed in order to advance the remote concepts of nuclear operation and, more broadly, microreactors themselves from paper to reality. This paper aims to fill these gaps and describes a test bed that would support development and deployment of a remote concept of nuclear operations, initial experimental results from that test bed, and the application of the test bed and experimental results for a digital-twin-based remote concept of operations underdevelopment at Idaho National Laboratory (INL). The platform chosen as a remote concept of nuclear operations test bed is the Single Primary Heat Extraction and Removal Emulator, known as SPHERE, located at INL. SPHERE is a small-scale non-nuclear test bed that emulates thermal behavior of a microreactor. The small-scale and non-nuclear nature of SHPERE limit safety concerns associated with remote operations while still providing the physical response representative of a microreactor. A network connection was added to SPHERE that enables remote-monitoring capability. This allows for real-time data streaming to networked workstations, data historians, and human-machine interfaces (HMIs). These are all critical components in a remote concept of operations, thus providing a robust development and demonstration platform. An initial experiment was performed using the SPHERE remote operations testbed. This included running a comprehensively instrumented SPHERE through a series of steady-state and transient operating scenarios in both normal and abnormal operating conditions, all while streaming live test data to a remote HMI and data warehouse. This initial experiment served three purposes: (1) characterizing the response of SPHERE, (2) demonstrating the remote connection to SPHERE, and (3) providing a baseline data set for development of a digital-twin-based remote concept of operations that is under development at INL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Transient Efficiency Flexibility and Reliability Optimization of Coal-Fired Power Plants (Model-Based Estimation Library Report)

This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a model-based estimation (MBE) library for generating digital twins of coal-fired power plants (CFPPs). The models used by the MBE are transient and physics-based. Model classes in the library that correspond to different components in CFPPs are presented with physical explanations and representative results. MBE library provides Extended Kalman Filter (EKF), and Unscented Kalman Filter (UKF) algorithms. The MBE library is generic and can be used to represent a variety of CFPPs. A digital twin of a CFPP will be created as the reduced-order model of that CFPP is embedded in the developed MBE architecture.

20 FOSSIL-FUELED POWER PLANTS↗