Search NASA⌕ Search

SEARCH · Search NASA

Results for “EVAPORATION MODEL”

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 271 records · Page 15

Influence of electron evaporative cooling on ultracold plasma expansion

The expansion of ultracold neutral plasmas (UCP) is driven primarily by the thermal pressure of the electron component and is therefore sensitive to the electron temperature. For typical UCP spatial extents, evaporative cooling has a significant influence on the UCP expansion rate at lower densities (less than 10{sup 8}/cm{sup 3}). We studied the effect of electron evaporation in this density range. Owing to the low density, the effects of three-body recombination were negligible. We modeled the expansion by taking into account the change in electron temperature owing to evaporation as well as adiabatic expansion and found good agreement with our data. We also developed a simple model for initial evaporation over a range of ultracold plasma densities, sizes, and electron temperatures to determine over what parameter range electron evaporation is expected to have a significant effect. We also report on a signal calibration technique, which relates the signal at our detector to the total number of ions and electrons in the ultracold plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Primordial Black Holes: Observational characteristics of the final evaporation

For many early universe theories predict the creation of Primordial Black Holes (PBHs). PBHs could have masses ranging from the Planck mass to 105 solar masses or higher depending on the size of the universe at formation. A Black Hole (BH) has a Hawking temperature which is inversely proportional to its mass. Hence a sufficiently small BH will quasi-thermally radiate particles at an ever-increasing rate as emission lowers its mass and raises its temperature. Moreover, the final moments of this evaporation phase should be explosive and its description is dependent on the particle physics model. In this work we investigate the final few seconds of BH evaporation, using the Standard Model and incorporating the most recent Large Hadron Collider (LHC) results, and provide a new parameterization for the instantaneous emission spectrum. We calculate for the first time energy-dependent PBH burst light curves in the GeV/TeV energy range. Moreover, we explore PBH burst search methods and potential observational PBH burst signatures. We have found a unique signature in the PBH burst light curves that may be detectable by GeV/TeV gamma-ray observatories such as the High Altitude Water Cerenkov (HAWC) observatory. Finally, the implications of beyond the Standard Model theories on the PBH burst observational characteristics are also discussed, including potential sensitivity of the instantaneous photon detection rate to a squark threshold in the 5–10 TeV range.

79 ASTRONOMY AND ASTROPHYSICS↗

Global variation of transpiration and soil evaporation and the role of their major climate drivers

Although global variation in actual evapotranspiration has been widely investigated, it remains unclear how its two major components, transpiration and soil evaporation, are driven by climate drivers across global land surface. This paper uses a well-validated, process-based model that estimates transpiration and soil evaporation, and for the first time investigates and quantifies how the main global drivers, associated to vegetation process and the water and energy cycle, drive the spatiotemporal variation of the two components. The results show that transpiration and soil evaporation dominate the variance of actual evapotranspiration in wet and dry regions, respectively. Dry southern hemisphere from 13°S to 27°S is highlighted since it contributes to 21% global soil evaporation variance, with only 11% global land area. In wet regions, particularly in the humid tropics, there are strong correlations between transpiration, actual evapotranspiration, and potential evapotranspiration, with precipitation playing a relatively minor role, and available radiative energy is the major contributor to the interannual variability in transpiration and actual evapotranspiration in Amazonia. Conversely in dry regions, there are strong correlations between soil evaporation, actual evapotranspiration, and precipitation. Lastly, our findings highlight that ecohydrological links are highly related to climate regimes, and the small region such as Australia has important contribution to interannual variation in global soil evaporation and evapotranspiration, and anthropogenic activities strongly influence the variances in irrigation regions

54 ENVIRONMENTAL SCIENCES↗

Creation of Multiple Effect Evaporator and Combined Cycle Modelica Modules, and Optimization of Potable Water Generation from Saltwater Sources (Final Technical Report)

Modelica models for a Heat Recovery Steam Generator (HRSG) and Multiple Effect Evaporator have been developed. These models have been validated against results from literature and from an actual cogeneration system on the campus of NC State University. The models are capable of connecting to other components to create multi-input, multi-output energy systems with storage, also known as integrated energy systems. An existing 35 MW gas turbine model is connected to the HRSG model created for this project to produce 9 bar steam. The high pressure steam feeds a backpressure steam turbine from the library to generate additional electric power. 2 bar steam leaves the turbine to supply an 8 effect MEE to generate water for a community. As part of the effort, three research papers were submitted (two published in ANS Proceedings and one submitted to Applied Energy).

20 FOSSIL-FUELED POWER PLANTS↗

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗

Toward a new parameterization of hydraulic conductivity in climate models: Simulation of rapid groundwater fluctuations in N orthern C alifornia

Abstract Preferential flow through weathered bedrock leads to rapid rise of the water table after the first rainstorms and significant water storage (also known as “rock moisture”) in the fractures. We present a new parameterization of hydraulic conductivity that captures the preferential flow and is easy to implement in global climate models. To mimic the naturally varying heterogeneity with depth in the subsurface, the model represents the hydraulic conductivity as a product of the effective saturation and a background hydraulic conductivity K bkg , drawn from a lognormal distribution. The mean of the background K bkg decreases monotonically with depth, while its variance reduces with the effective saturation. Model parameters are derived by assimilating into Richards’ equation 6 years of 30 min observations of precipitation (mm) and water table depths (m), from seven wells along a steep hillslope in the Eel River watershed in Northern California. The results show that the observed rapid penetration of precipitation and the fast rise of the water table from the well locations, after the first winter rains, are well captured with the new stochastic approach in contrast to the standard van Genuchten model of hydraulic conductivity, which requires significantly higher levels of saturated soils to produce the same results. “Rock moisture,” the moisture between the soil mantle and the water table, comprises 30% of the moisture because of the great depth of the weathered bedrock layer and could be a potential source of moisture to sustain trees through extended dry periods. Furthermore, storage of moisture in the soil mantle is smaller, implying less surface runoff and less evaporation, with the proposed new model.

58 GEOSCIENCES↗

Comparative analysis of thermal management systems in electric vehicles at extreme weather conditions: Case study on Nissan Leaf 2019 Plus, Chevrolet Bolt 2020 and Tesla Model 3 2020

With the surge in electric vehicle (EV) adoption and the need for extended driving ranges, optimizing energy efficiency, particularly through thermal management, is critical, especially in extreme weather. Managing the substantial energy needed for cabin climate control and battery temperature regulation can increase energy demands by over 50 %, severely limiting range. This study conducts a comparative analysis of thermal management systems (TMS) in three popular EV vehicles, 2020 Chevrolet Bolt, 2019 Nissan Leaf Plus, and 2020 Tesla Model 3, evaluating their distinct TMS configurations and performance under varied weather conditions. Using both numerical simulations and experimental data collected on a controlled test bench at Argonne National Laboratory, we assess how TMS architecture and operational modes influence energy consumption and range. A comprehensive TMS model was developed, integrating cabin and battery thermal sub-models in the Autonomie software platform, to simulate temperature fluctuations and range impacts. Cabin climate was modeled using a mono-zonal approach, while battery cell temperature distribution was estimated through a 2D nodal structure. Each vehicle's distinct TMS setup was evaluated: the Chevrolet Bolt and Tesla Model 3 use a dual evaporator vapor compression cycle with a PTC heater for the cabin and a coolant loop for battery thermal management; the Nissan Leaf Plus employs a heat pump with a PTC heater for the cabin and air-cooling for the battery. Tests conducted at ambient temperatures of 35°C, 22°C, -7°C, and -18°C reveal significant differences in energy use and range reduction across both configurations and conditions. At 35°C, the Tesla Model 3, Chevrolet Bolt, and Nissan Leaf Plus have a range reduction of 8%, 9%, and 13%, respectively, due to air conditioning. In winter, heating technology is paramount; at -7°C, the Nissan Leaf's heat pump configuration achieves a lower range reduction (19.3%) compared to the Tesla and Chevrolet Bolt PTC heaters, which reduce range by 28.3% and 31%, respectively. Further, this study provides valuable insights for automotive engineers, EV technology researchers, and thermal management system designers aiming to enhance electric vehicle performance by understanding how different weather conditions and TMS architectures impact energy consumption and driving range.

33 ADVANCED PROPULSION SYSTEMS↗

CFD modeling of crystallization fouling with CO 2 desorption incorporated for a falling-film evaporator in thermal desalination

Seawater flowing around horizontal tubes in falling-film evaporators is a common configuration for thermal desalination. Heat is transferred from in-tube steam condensation to the shell-side seawater by conduction through the tube wall and scaling layer, and conduction and convection to the evaporating liquid film. CO 2 is simultaneously released from seawater and mixed with the produced vapor. Fouling coupled with CO 2 desorption has deleterious impacts on both evaporation and condensation. The evaluation of spatiotemporal dependent crystallization fouling and CO 2 desorption is vital to selecting optimal operating conditions. Here, in this work, a CO 2 desorption model is integrated into CFD modeling for predictions of coupled heat and mass transfer, to understand and predict scale formation and CO 2 desorption with local and transient profiles of temperature, carbonate species concentrations, pH, and total alkalinity. The porosity of the scale layer is experimentally determined as 67.7% and invoked in the calculation of effective thermal conductivity. The simulation results reveal that high steam temperatures increase seawater total alkalinity and accelerate scale formation. The scale thickness on the bottom tube reaches 0.30 μm and 5.24 μm for steam temperatures of 60 °C and 80 °C, respectively. High salinity leads to a large CO 2 desorption rates. The CO 2 desorption rate increases 43% when the seawater salinity increases from 35 to 55 g/kg. The effects of operating conditions on carbonate speciation and pH have been compared and analyzed. This model can serve as a comprehensive tool for the design of thermal desalination systems and optimal operation.

42 ENGINEERING↗

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗

MO-A-BRD-10: A Fast and Accurate GPU-Based Proton Transport Monte Carlo Simulation for Validating Proton Therapy Treatment Plans

Purpose: To build a GPU-based Monte Carlo (MC) simulation of proton transport with detailed modeling of elastic and non-elastic (NE) protonnucleus interactions, for use in a very fast and cost-effective proton therapy treatment plan verification system. Methods: Using the CUDA framework, we implemented kernels for the following tasks: (1) Simulation of beam spots from our possible scanning nozzle configurations, (2) Proton propagation through CT geometry, taking into account nuclear elastic and multiple scattering, as well as energy straggling, (3) Bertini-style modeling of the intranuclear cascade stage of NE interactions, and (4) Simulation of nuclear evaporation. To validate our MC, we performed: (1) Secondary particle yield calculations in NE collisions with therapeutically-relevant nuclei, (2) Pencil-beam dose calculations in homogeneous phantoms, (3) A large number of treatment plan dose recalculations, and compared with Geant4.9.6p2/TOPAS. A workflow was devised for calculating plans from a commercially available treatment planning system, with scripts for reading DICOM files and generating inputs for our MC. Results: Yields, energy and angular distributions of secondaries from NE collisions on various nuclei are in good agreement with the Geant4.9.6p2 Bertini and Binary cascade models. The 3D-gamma pass rate at 2%–2mm for 70–230 MeV pencil-beam dose distributions in water, soft tissue, bone and Ti phantoms is 100%. The pass rate at 2%–2mm for treatment plan calculations is typically above 98%. The net computational time on a NVIDIA GTX680 card, including all CPU-GPU data transfers, is around 20s for 1×10{sup 7} proton histories. Conclusion: Our GPU-based proton transport MC is the first of its kind to include a detailed nuclear model to handle NE interactions on any nucleus. Dosimetric calculations demonstrate very good agreement with Geant4.9.6p2/TOPAS. Our MC is being integrated into a framework to perform fast routine clinical QA of pencil-beam treatment plans. Hardware for such a system will cost under $5000. This work was funded in part by a grant from Varian Medical Systems, Inc.

61 RADIATION PROTECTION AND DOSIMETRY↗

Recent Approaches to Modeling Transport of Mercury in Surface Water and Groundwater - Case Study in Upper East Fork Poplar Creek, Oak Ridge, TN - 13349

In this case study, groundwater/surface water modeling was used to determine efficacy of stabilization in place with hydrologic isolation for remediation of mercury contaminated areas in the Upper East Fork Poplar Creek (UEFPC) Watershed in Oak Ridge, TN. The modeling simulates the potential for mercury in soil to contaminate groundwater above industrial use risk standards and to contribute to surface water contamination. The modeling approach is unique in that it couples watershed hydrology with the total mercury transport and provides a tool for analysis of changes in mercury load related to daily precipitation, evaporation, and runoff from storms. The model also allows for simulation of colloidal transport of total mercury in surface water. Previous models for the watershed only simulated average yearly conditions and dissolved concentrations that are not sufficient for predicting mercury flux under variable flow conditions that control colloidal transport of mercury in the watershed. The transport of mercury from groundwater to surface water from mercury sources identified from information in the Oak Ridge Environmental Information System was simulated using a watershed scale model calibrated to match observed daily creek flow, total suspended solids and mercury fluxes. Mercury sources at the former Building 81-10 area, where mercury was previously retorted, were modeled using a telescopic refined mesh with boundary conditions extracted from the watershed model. Modeling on a watershed scale indicated that only source excavation for soils/sediment in the vicinity of UEFPC had any effect on mercury flux in surface water. The simulations showed that colloidal transport contributed 85 percent of the total mercury flux leaving the UEFPC watershed under high flow conditions. Simulation of dissolved mercury transport from liquid elemental mercury and adsorbed sources in soil at former Building 81-10 indicated that dissolved concentrations are orders of magnitude below a target industrial groundwater concentration beneath the source and would not influence concentrations in surface water at Station 17. This analysis addressed only shallow concentrations in soil and the shallow groundwater flow path in soil and unconsolidated sediments to UEFPC. Other mercury sources may occur in bedrock and transport though bedrock to UEFPC may contribute to the mercury flux at Station 17. Generally mercury in the source areas adjacent to the stream and in sediment that is eroding can contribute to the flux of mercury in surface water. Because colloidally adsorbed mercury can be transported in surface water, actions that trap colloids and or hydrologically isolate surface water runoff from source areas would reduce the flux of mercury in surface water. Mercury in soil is highly adsorbed and transport in the groundwater system is very limited under porous media conditions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Growth of sputter-deposited metamagnetic epitaxial Ni-Co-Mn-In films

Metamagnetic thin films represent a promising geometry for more efficient magnetocaloric cooling applications due to a fast heat transfer. Here, we identify suitable growth conditions to obtain epitaxial Ni-Mn-In-Co films with a metamagnetic transition in vicinity of room temperature. We show that both increased substrate temperature and target aging result in loss of indium. This can be attributed to evaporation and preferential sputtering, respectively. We present a model that treats the effect of target aging and temperature dependence of evaporation on the film composition independently and enables predictions of the film composition as a function of initial target composition, target age, and deposition temperature. Furthermore, our analysis reveals that a sufficient degree of chemical B2 order is required for a transformation, in addition to an appropriate film composition.

36 MATERIALS SCIENCE↗

A New Model to Simulate Energy Performance of VRF Systems

This paper presents a new model to simulate energy performance of variable refrigerant flow (VRF) systems in heat pump operation mode (either cooling or heating is provided but not simultaneously). The main improvement of the new model is the introduction of the evaporating and condensing temperature in the indoor and outdoor unit capacity modifier functions. The independent variables in the capacity modifier functions of the existing VRF model in EnergyPlus are mainly room wet-bulb temperature and outdoor dry-bulb temperature in cooling mode and room dry-bulb temperature and outdoor wet-bulb temperature in heating mode. The new approach allows compliance with different specifications of each indoor unit so that the modeling accuracy is improved. The new VRF model was implemented in a custom version of EnergyPlus 7.2. This paper first describes the algorithm for the new VRF model, which is then used to simulate the energy performance of a VRF system in a Prototype House in California that complies with the requirements of Title 24 ? the California Building Energy Efficiency Standards. The VRF system performance is then compared with three other types of HVAC systems: the Title 24-2005 Baseline system, the traditional High Efficiency system, and the EnergyStar Heat Pump system in three typical California climates: Sunnyvale, Pasadena and Fresno. Calculated energy savings from the VRF systems are significant. The HVAC site energy savings range from 51 to 85percent, while the TDV (Time Dependent Valuation) energy savings range from 31 to 66percent compared to the Title 24 Baseline Systems across the three climates. The largest energy savings are in Fresno climate followed by Sunnyvale and Pasadena. The paper discusses various characteristics of the VRF systems contributing to the energy savings. It should be noted that these savings are calculated using the Title 24 prototype House D under standard operating conditions. Actual performance of the VRF systems for real houses under real operating conditions will vary.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

On the Modeling of Wall Heat Flux Partitioning in Subcooled Flow Boiling

Accurate prediction of void distribution in reactor coolant channels at subcooled boiling conditions is an important factor in the analysis of the operating conditions of both pressurized water reactors (PWRs) and boiling water reactors (BWRs). Its importance arises from purely thermal-hydraulic considerations (thermal limits, CRUD formation, etc.) and the impact of void reactivity feedback on core neutronics. A mechanistic multidimensional model of subcooled boiling was originally developed by Kurul and Podowski. Whereas this model has been implemented in several commercial CFD codes and successfully used in computer simulations of flow and phase distributions in the coolant channels of light water reactor, additional new modeling aspects of forced-convection subcooled boiling heat transfer have been investigated in recent studies. The purpose of this paper is to formulate a modified model of wall heat flux partitioning, which allows for capturing and quantifying the effect of simultaneous evaporation and condensation at the heated wall on the net evaporation rate and wall temperature. The new model has been implemented in the NPHASE-CMFD code, numerically tested and validated against experimental data. A mechanistic wall heat flux partitioning model for subcooled flow boiling has been presented. This model has been combined with extensively tested models of interfacial mass, momentum and energy transfer terms. It turns out that the model predictions of the wall temperature and void fraction are in good agreement with the experimental data of reference. Future work will focus on the continuation of model validation against experimental data, and on extending its scope to include the phenomena governing departure from nucleate boiling. (authors)

42 ENGINEERING↗

A mathematical design framework for membrane pre-concentration in energy-efficient recovery of fermentation products

Due to the dilute nature of products manufactured via fermentation and cell-free bioprocessing, dewatering is a common unit operation in downstream processing (DSP) for bioproduct recovery, but it is typically energy intensive. To improve DSP energy efficiency for bio-based small molecules, integrating high-pressure membrane pre-concentration is a promising process option. However, this approach is typically constrained by a tradeoff between concentration factor (CF) and product recovery (PR), namely increasing the CF typically results in greater product loss, and vice versa. Here we developed a model that enables process design guidelines to: (i) identify scenarios in which the additional energy consumption and product loss from membrane pre-concentration are justified for use in DSP, and (ii) determine the optimal CF that minimizes process specific energy consumption. We compared the energy consumption of high-pressure membrane-integrated processes to evaporation-only processes and applied the model to an experimental case study for the separation and purification of butyric acid from Clostridium tyrobutyricum fermentation using an in situ product recovery (ISPR) process. The model estimated that integrating a tangential-flow reverse osmosis (RO) pre-concentration unit could reduce process energy consumption up to 45%. The use of advanced membrane pre-concentration technologies, such as negative rejection membranes and organic solvent reverse osmosis (OSRO), have the potential to further reduce the overall process specific energy consumption up to 96%, projected based on modeling. Overall, membrane pre-concentration, especially when strategically integrated prior to an evaporation step with optimized process conditions, holds significant potential for improving DSP energy efficiency, particularly in applications requiring substantial solvent removal for product recovery from dilute mixtures.

09 BIOMASS FUELS↗

Low global warming potential (GWP) refrigerant supermarket refrigeration system modeling and its application

As an environmentally friendly low global warming potential (GWP) refrigerant, Carbon dioxide (CO 2 ) has continuously gained popularity and research attention as alternative refrigerant for supermarket refrigeration system. In this paper, to fulfill the increasing need of accurate Low-GWP supermarket refrigeration models for development of supervisory level control and optimization strategies, a high fidelity model is developed for CO 2 transcritical supermarket refrigeration system which includes compressor rack of low temperature (LT) compressors and medium temperature (MT) compressors, air-cooled gas cooler, evaporator, expansion valves and other auxiliary equipment. A resistance-capacity model structure is proposed to simulate the display cases. Semi-thermodynamic models are proposed to estimate reciprocating compressors volumetric efficiency and power consumption. Furthermore the zone modeling approach is used for evaporator simulation, and air-cooled gas cooler is modeled with distributed modeling method. The expansion valve simulation is based on orifice flow model. To calibrate these models, both manufacture performance data and experimental data are used. The experiments are conducted with a full instrumental CO 2 supermarket refrigeration system installed in Oak Ridge National Lab Flexible Research Platform (FRP). The simulation model can predict the system performance, including power consumption, cooling capacity, mass flow rate, temperature, and pressure, with high accuracy (within ±4%) compared to experimental data. In addition, this developed model has been used to create the system optimum high side pressure for high side expansion valve control, and to generate wide operating range simulation data for developing the batter-equivalent commercial refrigeration system model which can be used in grid interactive control development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dynamic Simulation of a Sub-Critical Coal Fired Power Plant

In order to address the demanding operating conditions for remaining coal-fired power plants, a dynamic model and a suite of tools have been developed for studying load cycling and to find optimization opportunities. A sub-critical steam cycle power plant was modeled in a flow-sheet modeling tool, APROS™. The model represented the firing system, economizer, evaporator, superheat, and reheat systems. Four loads from 100% TMCR to 25% TMCR were calibrated and tested such that low-to-high cycling could be studied. The model was run through various load cycles; one of which is presented here. This modeling is a prototype for general use in developing cutting edge controls products and for maximizing economic, low-emissions, and efficient operation of the existing coal power fleet.

Braun, Timothy↗