Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable importance”

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 415 records · Page 23

In-situ imaging of flexure-induced fracture in tape-laminate composites using high-resolution X-ray computed tomography

Here, this study presents development of an in-situ flexure test for imaging progressive inter- and intralaminar fracture in tape-laminate composites using X-ray computed tomography (CT). The intent of this test is to provide detailed experimental observations of ply-level damage that can be used to validate existing, and develop new, progressive damage analysis (PDA) tools. The test consists of a vertically mounted specimen which is flexed using two eccentric compressive loads using an in-situ uniaxial load stage. The flexure specimen contains a starter notch mid-span which promotes initiation of composite failure within the X-ray field of view. Specimens with two different laminate stacking sequences were tested, imaged, analyzed. For a quasi-isotropic laminate with large angle changes between adjacent plies, there was near simultaneous growth of transverse cracks and delaminations below the midplane of the laminate. For a laminate with small angles between adjacent plies, there was extensive formation of transverse crack networks which penetrated the laminate thickness without delamination growth. In addition to imaging fracture, the X-ray CT data from both specimen types were used to quantify the variability in thickness and analyze the local orientation of individual plies. Overall, the proposed test and the image-data analysis methodology provided an important insight into the fracture processes in tape laminates and highlights the inherent ply-level geometrical variabilities that should be accounted for in PDA simulations.

42 ENGINEERING↗

Dominant Controls on Preferential Flow and Their Implications for Future Soil Water Fluxes

Abstract Soil water flow, particularly preferential flow (PF), is a critical control on hydrological and biogeochemical processes, including groundwater recharge, contaminant transport, and carbon cycling. However, it remains challenging to predict PF occurrence across large environmental gradients. Here, we developed a deep learning (DL) model to estimate event‐scale soil water flow velocity and the probability of PF occurrence using high‐frequency soil moisture and precipitation data from 33 sites across the National Ecological Observatory Network. The model demonstrated high skill in predicting the binary occurrence of PF (91% F1‐score; 85% accuracy) but the performance was limited in predicting soil water velocity ( R 2 = 0.31). We found that precipitation characteristics (duration, volume, and intensity) were the most important predictors for soil water velocity. Among the non‐precipitation event variables, sand content showed relatively high predictive skill, though differences among non‐event climate variables were generally modest. Lower sand content was associated with increased predicted soil water velocity, a finding that highlights the role of soil structure in producing more non‐uniform flow, which contrasts with traditional uniform flow models. Projecting a reduced DL model under both moderate and high‐emissions future climate scenarios (2060–2099 Representative Concentration Pathways 4.5 and 8.5), we found ∼7.3% increase under RCP4.5 and ∼15% under RCP8.5 of soil water velocities compared to the historical simulation, while modeled likelihood of PF changed little. These findings suggest climate change is not making PF more frequent, but it is making existing PF pathways more efficient with important consequences for associated nutrient and contaminant transport under climate change. Plain Language Summary Water movement in soil is critical for water quality. While often modeled as a uniform flow process, in reality water moves rapidly through cracks and burrows in what is called “preferential flow” (PF), which limits natural filtration and can transport pollutants. We developed a deep learning model, trained on data from 33 U.S. sites, to predict when and how fast this PF occurs based on precipitation, soil, and climate data. The model showed that precipitation characteristics (duration, intensity, volume) were the most important predictors of PF. Lower soil sand content/higher clay content was associated with faster water flow, likely due to clay soils forming aggregates and cracks that water moves through rather than infiltrating uniformly. Further analyses based on climate projections suggest that the speed at which PF occurs will become more rapid under future climate scenarios compared to historical simulation. This highlights the need to represent PF in soil water models when assessing future water quality. Key Points The effect of precipitation peak intensity on soil water velocities declined with increasing precipitation intensity Antecedent soil moisture failed to predict preferential flow (PF), contrasting the high predictive power of sand content Climate predictions suggest that soil water velocities through PF paths will increase ∼15% by 2099

Li, Bonan↗

Interfacial Exchange Phenomena Driven by Ferromagnetic Domains

Abstract Interfacial proximity effects in antiferromagnetic/ferromagnetic (AFM/FM) bilayers control the exchange‐bias (EB) phenomena exploited in most spintronic devices, although still is lack of full understanding. Discordant results, including different exchange‐bias field ( H E ), coercivity ( H C ), or blocking temperature ( T B ) found even in similar systems, are usually ascribed to uncontrolled parameters, namely dissimilar interfacial defects, structure, and thicknesses. Here, it is shown in the very same sample that the magnetic domain structure during the magnetization reversal of the FM layer controls those mentioned effects. Simultaneous transport and vectorial‐resolved magnetic measurements performed in a V 2 O 3 /Co system during warming after different field cooling (FC) procedures exhibit a strong dependence on the FC angle and the domain structure of the FM layer. Remarkably, magnetization reversal analysis reveals 35 K of variation in T B and up to a factor of two in H E . These observations can be explained within the random‐field model for the interfacial exchange coupling with a fixed AFM domain structure in contact with a variable (angle‐dependent) FM domain structure. The results highlight the importance of the domain structure and magnetization reversal of the FM layer (not previously considered) in the EB phenomena, with potential to tailor interfacial effects in future spintronic devices.

36 MATERIALS SCIENCE↗

Seasonal variation in methane ebullition in a temperate hydropower reservoir impacts annualized estimates

Estimates of ebullitive methane (CH4) emissions from reservoirs are highly uncertain, often based on spatial and temporal upscaling from a limited number of sampling sites and periods. While some of these estimates may adjust for seasonal variation in temperature, they rarely account for drawdown and the resulting changes in site depth and the areal extent of the reservoir. In Douglas Reservoir (Tennessee, USA), we conducted monthly sampling of 20 sites over 1 year to capture temporal and spatial variation in CH4 ebullition. We then compared various approaches for obtaining annualized, whole-reservoir ebullition estimates. CH4 ebullition fluxes were zero-inflated and highly skewed, with fluxes rising appreciably above zero in July and peaking in autumn. Estimates that did not account for drawdown overestimated annualized rates by 20%, though with high uncertainty. Seasonal variability and hotspots strongly influence annualized ebullition estimates and are important to account for in inventories and global upscaling.

Cardall, Anna [Brigham Young University]↗

Elucidating texture and grain morphology contributions to the micromechanical response of additively manufactured Inconel 625

Microstructural variation of additively manufactured (AM) metal components in comparison to wrought counterparts makes certification for critical applications a challenge. Microscale simulations leveraging modern computational tools may be used to supplement testing of AM microstructures, thus accelerating certification by reducing the number of experiments needed. However, as micromechanical response is closely tied to critical properties like fatigue-life and fracture, utilization of these simulations with macroscale experimental data alone is insufficient. One means to attain microscale experimental data is in situ diffraction data collected from synchrotron X-ray sources. In this work, such data were collected during in situ compression of AM Inconel 625 superalloy. Interpretation of experimental results was assisted by massive (8M element) complementary micromechanical simulations performed on sets of virtual microstructures generated using cellular automata. Together, micromechanical data from diffraction experiments and simulations were used to probe the effects of textured “track” microstructures generated during laser powder bed fusion and directional strength-to-stiffness on micromechanical response. Though fiber-averaged directional strength-to-stiffness ratios were expected to dominate given the high elastic anisotropy of the material, the combination of small variations in texture and specific grain configurations unique to AM microstructures lead to significant variability in micromechanical response after yield. The findings emphasize the importance of high-fidelity microstructural representation that captures key texture components and AM-specific morphology for property prediction of AM metals.

36 MATERIALS SCIENCE↗

Hydropower operation in future power grid with various renewable power integration

Hydropower generation may play an increasingly important role in the power grid under increasing contribution of variable renewable sources such as wind and solar. An improved understanding of the changes to hydropower dispatch under future higher VRE grid conditions reveals research gap that should be informed power grid planning and reservoir water releases policies considering multiple other water uses and varying hydrologic condition. Here this study aims to understand the role of hydropower in a changing power grid by employing a production cost model, PLEXOS, across future power system scenarios, planning horizons, and regions. We explore optimized hydropower dispatch to understand its potential role in minimizing the system cost and renewable curtailment. We also examine the sensitivity of hydropower revenue under various grid scenarios of the Eastern U.S. and hydrology conditions. Results indicate hydropower generation follows net load and compensates for the variability of solar and wind generation. Although energy prices are lower during some periods in the future grid scenarios, there is a potential for higher revenue for hydropower by providing both energy and ancillary services during times of stress. Additionally, hydropower revenue is sensitive to hydrology in the SERC region, which we considered as an example. The feasibility of hydropower dispatching with higher ramps between low and high hourly-capacity factors, as indicated in the optimization model, requires further study to consider other water use and ecology constraints.

13 HYDRO ENERGY↗

Demons in the North Atlantic: Variability of Deep Ocean Ventilation

Abstract Translation of atmospheric forcing variability into the ocean interior via ocean ventilation is an important aspect of transient climate change. On a seasonal timescale in the subtropics, this translation is mediated by a so‐called “Demon” that prevents access to all except late‐winter mixed‐layer water. Here, we use an eddy‐permitting numerical circulation model to investigate a similar process operating on longer (interannual) timescales in the subpolar North Atlantic. We find that variations in atmospheric forcing are mediated in their translation to the ocean interior, with year‐to‐year changes in the late‐winter mixed layer depth being the critical factor. The signature of persistent strong atmospheric forcing driving deep mixed layers is preferentially ventilated to the interior when the forcing is ceased. Susceptibility to this effect depends on the location and density of subduction—with the rate at which newly ventilated water escapes its region of subduction being the crucial factor.

MacGilchrist, G. A.↗

Robust production of monovalent bispecific IgG antibodies through novel electrostatic steering mutations at the C H 1-C λ interface

Bispecific antibodies represent an increasingly large fraction of biologics in therapeutic development due to their expanded scope in functional capabilities. Asymmetric monovalent bispecific IgGs (bsIgGs) have the additional advantage of maintaining a native antibody-like structure, which can provide favorable pharmacology and pharmacokinetic profiles. The production of correctly assembled asymmetric monovalent bsIgGs, however, is a complex engineering endeavor due to the propensity for non-cognate heavy and light chains to mis-pair. Previously, we introduced the DuetMab platform as a general solution for the production of bsIgGs, which utilizes an engineered interchain disulfide bond in one of the C H 1-C L domains to promote orthogonal chain pairing between heavy and light chains. While highly effective in promoting cognate heavy and light chain pairing, residual chain mispairing could be detected for specific combinations of Fv pairs. Here, we present enhancements to the DuetMab design that improve chain pairing and production through the introduction of novel electrostatic steering mutations at the C H 1-C L interface with lambda light chains (C H 1-C λ ). These mutations work together with previously established charge-pair mutations at the C H 1-C L interface with kappa light chains (C H 1-C κ ) and Fab disulfide engineering to promote cognate heavy and light chain pairing and enable the reliable production of bsIgGs. Importantly, these enhanced DuetMabs do not require engineering of the variable domains and are robust when applied to a panel of bsIgGs with diverse Fv sequences. We present a comprehensive biochemical, biophysical, and functional characterization of the resulting DuetMabs to demonstrate compatibility with industrial production benchmarks. Overall, this enhanced DuetMab platform substantially streamlines process development of these disruptive biotherapeutics.

60 APPLIED LIFE SCIENCES↗

Calibrating the Classical Hardness of the Quantum Approximate Optimization Algorithm

The trading of fidelity for scale enables approximate classical simulators such as matrix product states (MPSs) to run quantum circuits beyond exact methods. A control parameter, the so-called bond dimension $\mathcal{χ}$ for MPSs, governs the allocated computational resources and the output fidelity. Here, we characterize the fidelity for the quantum approximate optimization algorithm by the expectation value of the cost function that it seeks to minimize and find that it follows a scaling law $\mathscr{F}$(ln $\mathcal{χ}$/N), where N is the number of qubits. With ln $\mathcal{χ}$ amounting to the entanglement that a MPS can encode, we show that the relevant variable for investigating the fidelity is the entanglement per qubit. Importantly, our results calibrate the classical computational power required to achieve the desired fidelity and benchmark the performance of quantum hardware in a realistic setup. For instance, we quantify the hardness of performing better classically than a noisy superconducting quantum processor by readily matching its output to the scaling function. Moreover, we relate the global fidelity to that of individual operations and establish its relationship with $\mathcal{χ}$ and N. We sharpen the requirements for noisy quantum computers to outperform classical techniques at running a quantum optimization algorithm in speed, size, and fidelity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Building-Integrated Microscale Sensors for CO 2 -Level Monitoring

This research effort enabled the development of new classes of small-scale sensing systems that use mass and/or electrochemical sensors to detect the presence of carbon dioxide (CO 2 ). The concentration of this gas within the indoor environment is an important metric, as it can help enable the use of variable speed ventilation fans in commercial and governmental buildings, thus saving a significant amount of energy. Furthermore, there is a pressing need for enhanced CO 2 sensing to improve the comfort and productivity of people in commercial and government buildings, including academic spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conceptual Design of Pulverized Coal Electricity Generating Units for Flexible Operation

This presentation was given by Marc Turner, NETL Support Contractor, at the International Center for Sustainable Carbon's "The Energy Transition - the Role for Sustainable Carbon" workshop held in Sardinia, Italy on November 16-18, 2022. The presentation addresses the NETL report titled "Conceptual Design of Pulverized Coal Electricity Generating Units for Flexible Operation" published on September 6, 2022. This study presents an independent conceptual assessment of the features, cost, and performance of pulverized coal power plants (without carbon dioxide capture) intended for flexible, rather than baseload, operation. The objective was to leverage the already substantial knowledge base about flexible operation and develop state-of-the-art pulverized coal power plant cases that are designed to accommodate flexible operation within an electric grid having a high level of renewable generation. This is of significant importance in coal-dependent regions of the world as variable renewable energy becomes more prevalent.

Turner, Marc↗

maslov-group/FUN-PROSE

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE—a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

Liu, Simon↗

FUN-PROSE: A deep learning approach to predict condition-specific gene expression in fungi

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE—a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

59 BASIC BIOLOGICAL SCIENCES↗

Earth as an Exoplanet. I. Time Variable Thermal Emission Using Spatially Resolved Moderate Imaging Spectroradiometer Data

Among the more than 4000 exoplanets known today, some terrestrial planets have been detected in the so-called habitable zone of their host stars and their number is expected to increase in the near future, energizing a drive to understand and interpret the eagerly awaited wealth of data to identify signs of life beyond our solar system. So far, Earth remains the best and only example of a habitable (and inhabited) world. Although, it seems extremely unlikely that any other exoplanets will be true Earth twins, it is important to explore and understand the full range of spectral signatures and variability of Earth in order to inform the design of future instruments and missions, and understand their diagnostic power as well as potential limitations. In this work we use Earth observation data collected by the MODIS instrument aboard the Aqua satellite. The complete data set comprises 15 years of thermal emission observations in the 3.66–14.40 μm range for five different locations on Earth (Amazon Rainforest, Antarctica, Arctic, Indian Ocean, and the Sahara Desert). We then determine flux levels and variations as a function of wavelength and surface type (i.e., climate zone and surface thermal properties) and investigate whether periodic signals indicating Earth’s tilted rotation axis can be detected. Our findings suggest that (1) viewing geometry plays an important role when thermal emission data is analyzed as Earth’s spectrum varies by a factor of three and more depending on the dominant surface type underneath; (2) typically strong absorption bands from CO 2 (15 μm) and O 3 (9.65 μm) are significantly less pronounced and partially absent in data from the polar regions implying that estimating correct abundance levels for these molecules might be challenging in these cases; and (3) the time-resolved thermal emission spectrum encodes information about seasons/planetary obliquity, but the significance depends on the viewing geometry and spectral band considered.

79 ASTRONOMY AND ASTROPHYSICS↗

Universal Compiling and (No-)Free-Lunch Theorems for Continuous-Variable Quantum Learning

Quantum compiling, where a parameterized quantum circuit is trained to learn a target unitary, is an important primitive for quantum computing that can be used as a subroutine to obtain optimal circuits or as a tomographic tool to study the dynamics of an experimental system. While much attention has been paid to quantum compiling on discrete-variable hardware, less has been paid to compiling in the continuous-variable paradigm. Here we motivate several, closely related, short-depth continuous-variable algorithms for quantum compilation. We analyze the trainability of our proposed cost functions and numerically demonstrate our algorithms by learning arbitrary Gaussian operations and Kerr nonlinearities. We further make connections between this framework and quantum learning theory in the continuous-variable setting by deriving no-free-lunch theorems. These generalization bounds demonstrate a linear resource reduction for learning Gaussian unitaries using entangled coherent-Fock states and an exponential resource reduction for learning arbitrary unitaries using two-mode-squeezed states.

97 MATHEMATICS AND COMPUTING↗

A surrogate-model-based approach for estimating the first and second-order moments of offshore wind power

Power curve, the functional relationship that governs the process of converting a set of weather variables experienced by a wind turbine into electric power, is widely used in the wind industry to estimate power output for planning and operational purposes. Existing methods for power curve estimation have three main limitations: (i) they mostly rely on wind speed as the sole input, thus ignoring the secondary, yet possibly significant effects of other environmental factors, (ii) they largely overlook the complex marine environment in which offshore turbines operate, potentially compromising their value in offshore wind energy applications, and (ii) they solely focus on the first-order properties of wind power, with little (or null) information about the variation around the mean behavior, which is important for ensuring reliable grid integration, asset health monitoring, and energy storage, among others. In light of that, this study investigates the impact of several wind-and wave-related factors on offshore wind power variability, with the ultimate goal of accurately predicting its first two moments. Further, our approach couples OpenFAST—a multi-physics wind turbine simulator—with Gaussian Process (GP) regression to reveal the underlying relationships governing offshore weather-to-power conversion. We first find that a multi-input power curve which captures the combined impact of wind speed, direction, and air density, can provide double-digit improvements, in terms of prediction accuracy, relative to univariate methods which rely on wind speed as the sole explanatory variable (e.g. the standard method of bins). Wave-related variables are found not important for predicting the average power output, but interestingly, appear to be extremely relevant in describing the fluctuation of the offshore power around its mean. Tested on real-world data collected at the New York/New Jersey bight, our proposed multi-input models demonstrate a high explanatory power in predicting the first two moments of offshore wind generation, testifying their potential value to the offshore wind industry.

17 WIND ENERGY↗

Improved Simulations of Atmospheric River Climatology and Variability in High-Resolution CESM

Atmospheric rivers (ARs) can trigger extreme rainfall and snowfall events, causing high risk of flooding. Meanwhile, they are sometimes responsible for the occurrence and break of drought events. So, it is of great importance for weather and climate forecast models to accurately represent ARs. Studies have shown that AR variability is linked to various climate modes, such as El Niño/Southern Oscillation, Pacific Decadal Oscillation, Pacific-North America pattern, providing the potential for predicting ARs at subseasonal-to-decadal timescales. In this study, we examine the impact of climate model resolution on AR simulations and their relationship with major climate modes by comparing a pair of high-resolution and a low-resolution (LR) Community Earth System Model (CESM) simulations against the observations. The results reveal major deficiencies in LR that include underestimated climatological AR strength and associated precipitation, and inaccurate representations of ARs relationships with climate modes. Significant improvements are observed by increasing CESM horizontal resolution, suggesting a path for more accurate prediction of ARs.

54 ENVIRONMENTAL SCIENCES↗

A Novel Process for Carbon Dioxide Conversion to Fuel

In this project, TDA developed a new mixed metal oxide-based sorbent that converts CO2 (captured from coal fired power plants) to CO, which can then be combined with renewable H2 generated by water electrolysis or H2 from steam methane reforming to produce different liquid fuels. TDA’s absorbent-based CO2 conversion process uses a redox process, which splits the catalytic reforming of methane with CO2 reaction into two stages: CO2 reduction to CO and CH4 reforming into H2 and CO which eliminates the equilibrium limitations. The CO produced in the two-stage reactor system can then be further reacted with renewable H2 to produce methanol, naphtha, diesel, or gasoline. We worked with the Advanced Power & Energy Program (APEP) of University of California, Irvine (UCI) to design and develop the liquid fuel synthesis process that is built around this new material. We demonstrated the techno-economic viability of the new sorbent based redox process to convert CO2 into synthesis gas by demonstrating continuous carbon dioxide reduction in a prototype test system for over 585 hours while converting up to 10 kg CO2/day. With the successful completion of the R&D effort, the technology is now ready for a larger pilot-scale demonstration and the technology readiness has been raised from TRL 3 to TRL 5. In collaboration with UCI, we completed a high-fidelity process design and economic analysis. The required sale price (RSP) for gasoline (Case 1 NG-MTG) is $4.91/gal and naphtha and diesel (Case 2 NG-FT) are $4.23/gal and $6.07/gal, respectively, on a 2011 dollar basis. To put these costs in perspective, the California prices in current dollars (with its strict specifications) for regular grade gasoline from last year to current year have varied from a low of $3.10/gal in January 2021 to a high of $5.76/gal in March 2022, while prices for diesel from last year to current year have varied from a low of $3.40/gal in January 2021 to a high of $6.41/gal in May 2022 according to the U.S. Energy Information Agency data. It should be noted that the gasoline and diesel produced by these designs of Case 1 and Case 2 would be of very high quality and both nitrogen and sulfur free. These RSPs are based on a cost of imported electricity of $64/MWh based on the low-end current wind generated electricity cost (Genevieve 2011). This cost is by far the largest component of the variable costs used in computing the RSPs. A sensitivity analysis of these RSPs to the cost of imported electricity shows that the cost of the imported electricity has a significant effect on the RSPs. The life cycle analysis (LCA) shows that the total cradle-to-gate CO2 emissions for both liquid fuels (diesel and gasoline) were negative, indicating that overall more CO2 is consumed than released during production of the fuel from CO2 feed stack for both cases (-296 kgCO2 per MT gasoline for Case 1 and -705 kgCO2 per MT diesel for Case 2). On a cradle-to-grave comparison, the use of diesel produced using TDA’s process (2,457 kgCO2 per MT diesel) would release 37.6% less CO2 compared to petroleum based diesel (3,937 kgCO2 per MT diesel) while the use of gasoline produced using TDA’s process (2,792 kgCO2 per MT gasoline) would release 29.3% less CO2 compared to petroleum based gasoline (3,946 kgCO2 per MT gasoline). With the successful completion of the R&D effort, the technology is now ready for a bench-scale demonstration and the technology readiness has been raised from TRL 3 (Analytical and experimental critical function and/or characteristic proof of concept) to TRL 5 (Laboratory scale similar system validation in relevant environment).

20 FOSSIL-FUELED POWER PLANTS↗