Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble simulations”

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

Pulsed studies of intervalley transfer in Al 0.35 In 0.65 As: A paradigm for valley photovoltaics

Al 0.35 In 0.65 As is a direct semiconductor whose multivalley conduction-band structure has been proposed for use in so-called valley photovoltaics. In such hot-carrier solar cells, energetic (hot) photocarriers are stored in satellite valleys away from the Γ point, allowing them to be extracted prior to thermalization and to thereby increase power-conversion efficiency. While prior theoretical work has highlighted the potential of Al 0.35 In 0.65 As—a widely used barrier material in electronic and optoelectronic devices, for use in valley photovoltaics—surprisingly little is known about its electrical properties, especially how these are impacted by the application of high fields. In this work, we therefore undertake a detailed characterization of the electrical properties of Te-doped ($\textit{n}$-type) Al 0.35 In 0.65 As, over wide ranges of temperature (3–400 K) and electric field (<50 kV/cm). Using pulsed measurements to suppress the influence of Joule heating, we reveal the presence of clear negative-differential conductance in the current-voltage characteristics of the films, suggestive of the intervalley transfer of hot electrons. This conclusion is supported by the results of ensemble Monte Carlo simulations of the hot-carrier action, which confirm the connection of the observed negative-differential conductance to hot-electron transfer from the conduction-band (Γ) minimum, to the side valleys at the $\textit{L}$ point. Further, the quantitative features of the experimentally determined velocity-field curves are found to be in good agreement with the results of these calculations, providing further confidence in the role of the implied intervalley transfer mechanism. Overall, these results confirm the excellent potential of Al 0.35 In 0.65 As for use as the absorber material in hot-carrier solar cell technology.

36 MATERIALS SCIENCE↗

Parallel Dislocation Simulator

ParaDiS, or Parallel Dislocation Simulator, is a simulation tool that performs direct numerical simulation of dislocation ensembles, the carriers of plasticity, to predict the strength in crystalline materials from the fundamental physics of defect motion, evolution, and interaction. The code has been successfully deployed on high performance computing architectures and used to study the origins of strength and strain hardening for cubic crystals, the strength of micro-pillars, and irradiated materials at LLNL. The ParaDiS code has been successfully deployed on more than one hundred thousand CPU's with over ten million active degrees of freedom.

Bulatov, Vasily↗

A Hybrid Climate Modeling System Using AI-assisted Process Emulators

This white paper addresses Focus Area II. We advocate developing a hybrid modeling system to improve the understanding of decadal- and longer-scale predictability of high impact water cycle components. This hybrid model combines a partial differential equation (PDE)-based dynamic core with AI/ML based emulators to represent many of the computationally expensive processes in Earth’s climate models. The hybrid modeling system has the potential to exploit emerging graphics processing unit (GPU)-accelerated architectures and allows for the generation of large ensemble (~1000’s) simulations to better characterize the model uncertainty and understand predictability.

58 GEOSCIENCES↗

Data-Driven Specification of Earth's Magnetic Field

As our society is increasingly dependent on technologies in space, predicting space weather becomes ever more important. We have come to rely on technologies such as cellphones, GPS, the Internet, and other commercial and military assets. In order to protect these systems a high priority is to develop better understanding of the harmful conditions of space weather encountered in the near-Earth radiation environment which are mostly dictated by the dynamics of energetic charged particles. The main objective of our proposal was to improve the specification and forecast of the particle distribution in the Earth’s inner magnetosphere using data assimilation together with machine learning (ML) algorithms in our ring current-atmosphere interactions model with self-consistent magnetic and electric fields (RAM-SCBE). Specifically, we explored the application of machine learning algorithms to emulate and improve, through infusion of observational data, the magnetic field (SCBE) of our ring current model. Our work advanced towards the characterization of the RAM-SCBE particle flux field and its potentially different simulation scenarios. We explored how the model behaves in a number of simulation scenarios in order to better understand how the model reacts to perturbations in the solar wind. We were able to capture this behavior through the simulation of an ensemble of RAM-SCB model solutions. The ensemble showed the important correlations between the particle fluxes and the magnetic field. We were able to start the training of a neural network (Autoencoder) to identify the dominant features of the magnetic field, and the main correlations (linear and non-linear) between the particle fluxes and the magnetic field. Our next steps consists of training a neural network to capture the time evolution of the dominant features and obtain a consistent surrogate model for the RAM-SCB.

42 ENGINEERING↗

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗

NOW-23: the 2023 National Offshore Wind Data Set

In this report, we present the latest wind resource data specifically tailored for offshore regions in the United States. The data set, known as the 2023 National Offshore Wind data set (NOW-23), has been developed by the National Renewable Energy Laboratory (NREL) and its partners, and surpasses the previous resource data set, the Wind Integration National Dataset (WIND) Toolkit, which was released approximately ten years ago, in its offshore component. The WIND Toolkit has been widely utilized by stakeholders involved in wind resource assessments across the continental United States. However, with significant advancements in numerical weather prediction modeling over the past decade, the NOW-23 data set incorporates the latest research and development progress to provide stakeholders with an updated and cutting-edge resource for offshore wind analysis. The NOW-23 data set is created using the Weather Research and Forecasting (WRF) model and its output is available, as for its predecessor, at 5-minute time resolution and 2-kilometer horizontal spatial resolution. However, the NOW-23 data set improves upon the WIND Toolkit through: 1. A modeling period of at least 20 years (and as long as 23 years in selected regions), starting in 2000 (compared to the 7-year 2007–2013 modeling period in the WIND Toolkit). 2. For several offshore regions, a region-specific sensitivity analysis, driven by an ensemble of WRF simulations, to assess the most adequate region-specific WRF setup. 3. An updated WRF model, from Version 3.4 used in the WIND Toolkit to Version 4.2.1 used for the NOW-23 data set, which incorporates significant research advancements. 4. The use of the state-of-the-art reanalysis product ERA5 (which supersedes the older ERA-Interim used in the WIND Toolkit) to provide atmospheric forcing at the WRF domain boundaries. Figure 1 shows the NOW-23 mean wind speed at 160 m above sea level (asl), which we use as proxy for hub-height of a commercial offshore wind turbine in this report across all regions.

17 WIND ENERGY↗

Protein Conformational States—A First Principles Bayesian Method

Automated identification of protein conformational states from simulation of an ensemble of structures is a hard problem because it requires teaching a computer to recognize shapes. We adapt the naïve Bayes classifier from the machine learning community for use on atom-to-atom pairwise contacts. The result is an unsupervised learning algorithm that samples a ‘distribution’ over potential classification schemes. We apply the classifier to a series of test structures and one real protein, showing that it identifies the conformational transition with >95% accuracy in most cases. A nontrivial feature of our adaptation is a new connection to information entropy that allows us to vary the level of structural detail without spoiling the categorization. This is confirmed by comparing results as the number of atoms and time-samples are varied over 1.5 orders of magnitude. Further, the method’s derivation from Bayesian analysis on the set of inter-atomic contacts makes it easy to understand and extend to more complex cases.

97 MATHEMATICS AND COMPUTING↗

Challenges for Implementing FAIR Digital Objects with High Performance Workflows

New types of workflows are being used in science that couple traditional distributed and high-performance computing (HPC) with data-intensive approaches, and orchestrate ensembles of numerical simulations and artificial intelligence (AI) models. Such workflows may use AI models to supplement computation where numerical simulations may be too computationally expensive, to automate trivial yet time consuming operations, to perform preliminary selections among intractable numbers of combinations in domains as diverse as protein binding, fine-grid climate simulations, and drug discovery.

97 MATHEMATICS AND COMPUTING↗

The Role of Coupled Feedbacks in the Decadal Variability of the Southern Hemisphere Eddy-Driven Jet

The Southern Hemisphere summertime eddy-driven jet and storm tracks have shifted poleward over the recent few decades. In previous studies, explanations have mainly stressed the influence of external forcing in driving this trend. Here we examine the role of internal tropical SST variability in controlling the austral summer jet’s poleward migration, with a focus on interdecadal time scales. The role of external forcing and internal variability are isolated by using a hierarchy of Community Earth System Model version 1 (CESM1) simulations, including the pre-industrial control, large ensemble, and pacemaker runs. Model simulations suggest that in the early twenty-first century, both external forcing and internal tropical Pacific SST variability are important in driving a positive southern annular mode (SAM) phase and a poleward migration of the eddy-driven jet. Additionally, Tropical Pacific SST variability, associated with the negative phase of the interdecadal Pacific oscillation (IPO), acts to shift the jet poleward over the southern Indian and southwestern Pacific Oceans and intensify the jet in the southeastern Pacific basin, while external forcing drives a significant poleward jet shift in the South Atlantic basin. In response to both external forcing and decadal Pacific SST variability, the transient eddy momentum flux convergence belt in the middle latitudes experiences a poleward migration due to the enhanced meridional temperature gradient, leading to a zonally symmetric southward migration of the eddy-driven jet. This mechanism distinguishes the influence of the IPO on the midlatitude circulation from the dynamical impact of ENSO, with the latter mainly promoting the subtropical wave-breaking critical latitude poleward and pushing the midlatitude jet to higher latitudes.

54 ENVIRONMENTAL SCIENCES↗

Role of Tropical Variability in Driving Decadal Shifts in the Southern Hemisphere Summertime Eddy-Driven Jet

The Southern Hemisphere summertime eddy-driven jet and storm tracks have shifted poleward over the recent few decades. In previous studies, explanations have mainly stressed the influence of external forcing in driving this trend. Here we examine the role of internal tropical SST variability in controlling the austral summer jet’s poleward migration, with a focus on interdecadal time scales. The role of external forcing and internal variability are isolated by using a hierarchy of Community Earth System Model version 1 (CESM1) simulations, including the pre-industrial control, large ensemble, and pacemaker runs. Model simulations suggest that in the early twenty-first century, both external forcing and internal tropical Pacific SST variability are important in driving a positive southern annular mode (SAM) phase and a poleward migration of the eddy-driven jet. Tropical Pacific SST variability, associated with the negative phase of the interdecadal Pacific oscillation (IPO), acts to shift the jet poleward over the southern Indian and southwestern Pacific Oceans and intensify the jet in the southeastern Pacific basin, while external forcing drives a significant poleward jet shift in the South Atlantic basin. In response to both external forcing and decadal Pacific SST variability, the transient eddy momentum flux convergence belt in the middle latitudes experiences a poleward migration due to the enhanced meridional temperature gradient, leading to a zonally symmetric southward migration of the eddy-driven jet. This mechanism distinguishes the influence of the IPO on the midlatitude circulation from the dynamical impact of ENSO, with the latter mainly promoting the subtropical wave-breaking critical latitude poleward and pushing the midlatitude jet to higher latitudes.

54 ENVIRONMENTAL SCIENCES↗

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Meteorology & Atmospheric Sciences↗

Ettringite hysteresis under sorption from molecular simulations

Highlights: • Hybrid Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) simulations are performed. • Anisotropic large strain coefficient of crystal shrinkage of ettringite is provided. • Entropy increases when ettringite converts into metaettringite. • Sorption-induced stress may lead to columns collapse upon drying. • H-bonds involving OH groups in columns play a role in hysteresis. The response of ettringite under sorption is critical for the utilization of this mineral as an energy storage material and to understand the concrete durability problems associated with ettringite formation. We report molecular simulations directly sampling the osmotic ensemble to understand the sorption processes and structural changes in ettringite-metaettringite transition. Desorption branch and sorption-induced volume changes from simulations agree with the experiments. The structural changes at low relative humidity are analyzed, revealing for the first time details of metaettringite structure. The reversibility under sorption observed experimentally is captured by simulation in the ettringite domain. The reversibility in the conversion of metaettringite into ettringite is not captured by the direct simulations in the osmotic ensemble due to the large volume changes associated with ettringite desorption. Finally, we discuss the role of hydrogen bonds on the hysteresis in ettringite.

36 MATERIALS SCIENCE↗

Cyclone Fani: the tug-of-war between regional warming and anthropogenic aerosol effects

Before Cyclone Amphan took place in 2020, Cyclone Fani (May 2019) is the strongest pre-monsoon cyclone in the Bay of Bengal (BOB) since 1991, killing 90 people in eastern India and Bangladesh while causing US$1.81 billion of damages. Fani developed during a period of high concentration of anthropogenic aerosols in the BOB with abnormally high sea surface temperature (SST), thereby presenting an opportunity to understand the compound effects of atmospheric aerosols and regional climate warming on a tropical cyclone. A quantitative attribution analysis was conducted using the Weather Research and Forecasting model with chemistry (WRF-Chem) run at the convection-permitting (4 km) grid spacing, accompanied by an ensemble of coarser-resolution simulations to quantify the uncertainty. The removal of post-1990 trends in the tropospheric variables and SST from WRF-Chem's initial conditions (IC) and boundary conditions (BC, including the lateral and lower boundary conditions) resulted in a reduction of cyclone precipitation by about 51% during the 5 d of April 28-May 2. The removal of tropospheric warming shows approximately twice as strong an effect on Fani (39% reduction in precipitation) as that of SST warming (22% reduction). When aerosol's direct and indirect effects were removed from the simulations, i.e., no aerosol influence on radiation and cloud microphysics, Fani initially strengthened but later weakened, as measured by geopotential height and precipitation amounts. These results suggest that aerosol and its interaction with the atmosphere acted to mitigate the strengthening effect of anthropogenic warming on Fani, but was not strong enough to entirely counteract it. Although the ensemble of coarser simulations appears to overestimate Cyclone Fani in terms of precipitation, the direction of the effects is in agreement with that obtained from the 4 km simulations. Given the increasing anthropogenic aerosols in the BOB, future attribution studies using more sophisticated dynamical aerosol models on BOB tropical cyclones are urged.

54 ENVIRONMENTAL SCIENCES↗

High-Resolution Comonomer Sequencing of Blocky Brominated Syndiotactic Polystyrene Copolymers Using 13C NMR Spectroscopy and Computer Simulations

This work demonstrates the first high-resolution comonomer sequencing of Blocky brominated syndiotactic polystyrene (sPS-co-sPS-Br) copolymers based on pentad assign- ments of the quaternary carbon region of the nuclear magnetic resonance spectrum. Copolymers containing p-bromostyrene (Br-Sty) units were prepared in matched sets using postpolymerization bromination methods carried out in the heterogeneous gel state (Blocky) and homogeneous solution state (Random). Quantitative information from the quaternary carbon spectra, heteronuclear multiple bond correlation spectroscopy, electronic structure calculations, and simulated statistically random copolymers was correlated to confirm the carbon resonance assignments for all 20 possible pentad comonomer sequences. Using the experimental pentad sequence prevalences, a computer code was developed to simulate chains with microstructures typical of each sample as a means to visually represent the copolymer blockiness with quantitative precision. Based on the microstructure and distribution of run lengths in these chains, the simulations revealed that the Blocky copolymers contain a high degree of blockiness. By comparing the run lengths in the simulated chains to the average number of styrene units in a crystalline segment of sPS (found by small-angle X-ray scattering), copolymer crystallizability was predicted. For the simulated Blocky B-21% (21 mol % Br-Sty) chain, the probability of randomly selecting a styrene unit in a crystallizable block was 25.8%, while that in the simulated Random R-18% was zero, in excellent agreement with the experimental crystallization behavior measured by differential scanning calorimetry. Additionally, these predictions confirmed that the simulated chains accurately represent the ensemble of chains in their respective copolymer samples. Furthermore, each simulated Blocky chain contained one or more long sPS blocks that paralleled the measured 38-40 styrene units spanning a crystalline segment within the sPS/CCl4 gel. This finding affirmed that the long sPS segments originated from the precise lamellar structure within the heterogeneous gel morphology (i.e., block length is correlated with lamellar thickness). Overall, the ability to tailor the copolymer microstructure through control of the semicrystalline gel morphology opens the door to synthesizing ordered copolymers by postpolymerization functionalization processes with unprecedented levels of compositional control.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗

CERESMIP: a climate modeling protocol to investigate recent trends in the Earth's Energy Imbalance

The Clouds and the Earth's Radiant Energy System (CERES) project has now produced over two decades of observed data on the Earth's Energy Imbalance (EEI) and has revealed substantive trends in both the reflected shortwave and outgoing longwave top-of-atmosphere radiation components. Available climate model simulations suggest that these trends are incompatible with purely internal variability, but that the full magnitude and breakdown of the trends are outside of the model ranges. Unfortunately, the Coupled Model Intercomparison Project (Phase 6) (CMIP6) protocol only uses observed forcings to 2014 (and Shared Socioeconomic Pathways (SSP) projections thereafter), and furthermore, many of the ‘observed' drivers have been updated substantially since the CMIP6 inputs were defined. Most notably, the sea surface temperature (SST) estimates have been revised and now show up to 50% greater trends since 1979, particularly in the southern hemisphere. Additionally, estimates of short-lived aerosol and gas-phase emissions have been substantially updated. These revisions will likely have material impacts on the model-simulated EEI. We therefore propose a new, relatively low-cost, model intercomparison, CERESMIP, that would target the CERES period (2000-present), with updated forcings to at least the end of 2021. The focus will be on atmosphere-only simulations, using updated SST, forcings and emissions from 1990 to 2021. The key metrics of interest will be the EEI and atmospheric feedbacks, and so the analysis will benefit from output from satellite cloud observation simulators. The Tier 1 request would consist only of an ensemble of AMIP-style simulations, while the Tier 2 request would encompass uncertainties in the applied forcing, atmospheric composition, single and all-but-one forcing responses. We present some preliminary results and invite participation from a wide group of models.

54 ENVIRONMENTAL SCIENCES↗

Visual HPC Workflows for the Analysis of System Dynamics Models

Visual analytics supported by high performance computing (HPC) accelerates and enhances the discovery, exploration, and analysis of causal patterns in complex system dynamics (SD) models. We present a suite of visualization-assisted ensemble-based techniques for hypothesis generation and testing, and for sensitivity analysis. By employing HPC to provide parallel, on-demand simulation of SD models, one can “steer” an ensemble of simulated scenarios in real time as one first formulates and then informally tests those hypotheses: this provides rapid feedback for analysts to refine their understanding of the causal relationships emergent from a model. Such understandings can be followed and augmented by rigorous application of statistical methods, namely global variance-based sensitivity analysis, Monte-Carlo filtering, adaptive regional sensitivity analysis, and self-organized maps: here timely computation relies on HPC, while effective presentation emphasizes high-dimensional multivariate data visualization. Immersive visualization in virtual 3D environments provides an excellent adjunct to the traditional 2D graphics typically used for SD models, as it generates an embodied understanding of model behavior and facilitates an active, collaborative critique of model structure and output. Finally, we summarize prospects for HPC-enabled visual analytics applied to SD modeling.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗