Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation System”

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 1,279 records · Page 71

Feasibility of a Novel Density Functional Method Outside the Kohn-Sham Framework for Modeling Global Potential Energy Surfaces of Molecular Chemical Reactions (Final Technical Report)

The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DWT-Based Aggregated Load Modeling and Evaluation for Quasi-Static Time-Series Simulation on Distribution Feeders: Preprint

This paper presents an extension of a previously reported discrete wavelet transform (DWT) based load modeling methodology which targets at modeling time-series load profiles to enable the more effective and realistic quasi-static time-series simulation of distribution feeders. The time-series load model is composed of two major parts: 1) The low-resolution field-measured load data which usually have a resolution of 30 or 15 min., 2) The high-resolution load variability model data extracted from the established variability database which can have resolution up to 1 sec. A load aggregation methodology is developed to aggregate the load profiles so that load profiles at different transformer ratings can be effectively modeled. Validation, evaluation, and analysis of the developed load modeling approach has been performed on the IEEE-123 feeder and an actual utility feeder from California. The analysis completed on the two feeders have demonstrated the effectiveness and revealed the value of the developed model for distribution feeder quasi-static time-series simulation at high temporal resolution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Uncertainty quantification tools for multiphase gas-solid flow simulations using MFIX

Computational fluid dynamics (CFD) has been widely studied and used in the scientific community and in the industry. Various models were proposed to solve problems in different areas. However, all models deviate from reality. Uncertainty quantification (UQ) process evaluates the overall uncertainties associated with the prediction of quantities of interest. In particular it studies the propagation of input uncertainties to the outputs of the models so that confidence intervals can be provided for the simulation results. In the present work, a non-intrusive quadrature-based uncertainty quantification (QBUQ) approach is proposed. The probability distribution function (PDF) of the system response can be then reconstructed using extended quadrature method of moments (EQMOM) and extended conditional quadrature method of moments (ECQMOM). The report first explains the theory of QBUQ approach, including methods to generate samples for problems with single or multiple uncertain input parameters, low order statistics, and required number of samples. Then methods for univariate PDF reconstruction (EQMOM) and multivariate PDF reconstruction (ECQMOM) are explained. The implementation of QBUQ approach into the open-source CFD code MFIX is discussed next. At last, QBUQ approach is demonstrated in several applications. The method is first applied to two examples: a developing flow in a channel with uncertain viscosity, and an oblique shock problem with uncertain upstream Mach number. The error in the prediction of the moment response is studied as a function of the number of samples, and the accuracy of the moments required to reconstruct the PDF of the system response is discussed. The QBUQ approach is then demonstrated by considering a bubbling fluidized bed as example application. The mean particle size is assumed to be the uncertain input parameter. The system is simulated with a standard two-fluid model with kinetic theory closures for the particulate phase implemented into MFIX. The effect of uncertainty on the disperse-phase volume fraction, on the phase velocities and on the pressure drop inside the fluidized bed are examined, and the reconstructed PDFs are provided for the three quantities studied. Then the approach is applied to a bubbling fluidized bed with two uncertain parameters, particle-particle and particle-wall restitution coefficients. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities and gas pressure are provided. The PDFs of the response are reconstructed using EQMOM with appropriate kernel density functions. The simulation results are compared to experimental data provided by the 2013 NETL small-scale challenge problem. Lastly, the proposed procedure is demonstrated by considering a riser of a circulating fluidized bed as an example application. The mean particle size is considered to be the uncertain input parameter. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities, and granular temperature are provided. Mean values and confidence intervals of the quantities of interest are compared to the experiment results. The univariate and bivariate PDF reconstructions of the system response are performed using EQMOM and ECQMOM.

97 MATHEMATICS AND COMPUTING↗

Numerical discretization of completely integrable nonlinear Hamiltonian systems

Completely integrable nonlinear Hamiltonian systems offer attractive models for novel accelerator systems that are being developed to push the frontiers of high intensity and high power beam applications. However, all known numerical methods used to discretize and simulate these systems spoil the integrability property of the original, time-continuous, theoretically integrable system. Since the hallmark of integrability is the existence of sufficiently many functionally independent first integrals in involution, we investigated several methods’ performance in terms of invariant preservation over longer timescales. We conclude that, while novel structure preserving integrators would be ideal, symplectic methods under mild restrictions offer adequate invariant preservation. We especially highlight the very good performance of the Stormer-Verlet method at low accuracy and a modified Lobatto method at high accuracy for the Integrable Optics Test Accelerator (IOTA). However, we show that symplectic integrators are limited by numerically induced resonances, which make interpretation of some simulation results challenging. These limitations and workarounds are discussed.

43 PARTICLE ACCELERATORS↗

Reducing Long‐Standing Surface Ozone Overestimation in Earth System Modeling by High‐Resolution Simulation and Dry Deposition Improvement

The overestimation of surface ozone concentration in low‐resolution global atmospheric chemistry and climate models has been a long‐standing issue. We first update the ozone dry deposition scheme in both high‐ (0.25°) and low‐resolution (1°) Community Earth System Model (CESM) version 1.3 runs, by adding the effects of leaf area index and correcting the sunlit and shaded fractions of stomatal resistances. With this update, 5‐year‐long summer simulations (2015–2019) using the low‐resolution CESM still exhibit substantial ozone overestimation (by 6.0–16.2 ppbv) over the U.S., Europe, eastern China, and ozone pollution hotspots. The ozone dry deposition scheme is further improved by adjusting the leaf cuticle conductance, reducing the mean ozone bias by 19%, and increasing the model resolution further reduces the ozone overestimation by 43%. We elucidate the mechanism by which model grid spacing influences simulated ozone, revealing distinctive pathways in urban versus rural areas. In rural areas, grid spacing mainly affects daytime ozone levels, where additional NO x emissions from nearby urban areas result in an ozone boost and overestimation in low‐resolution simulations. In contrast, over urban areas, daytime ozone overestimation follows a similar mechanism due to the influence of volatile organic compounds from surrounding rural areas. However, nighttime ozone overestimation is closely linked to weakened NO titration owing to the redistribution of urban NO x to rural areas. Additionally, stratosphere‐troposphere exchange may also contribute to reducing ozone bias in high‐resolution simulations, warranting further investigation. This optimized high‐resolution CESM may enhance understanding of ozone formation mechanisms, sources, and changes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

Exploring Microsimulation Process for Energy Impact Evaluation of Connected and Automated Vehicles

In this paper, the authors present a microsimulation-based methodological approach for evaluating the energy impact of connected and automated vehicles (CAVs). They use an open-source micro-simulator, SUMO, and provide a way to set up a simulation environment that emulates real-world traffic dynamics. They also employ the Intelligent Driver Model to represent human drivers and calibrate its driving behavior using real-world traffic data and driving statistics. The authors conduct extensive simulation studies considering different penetration rates of CAVs, different car-following models, and varying car-following model parameters. Using the state-of-the-art Future Automotive System Technology Simulator (FASTSim), they estimate the fuel economy of each vehicle and analyze the energy impact of the given CAV implementation. Finally, the authors analyze the possible factors affecting the simulation results, and also discuss limitations and future work.

ADVANCED PROPULSION SYSTEMS↗

Power Grid Simulation Testbed for Transactive Energy Management Systems

To effectively engage demand-side and distributed energy resources (DERs) for dynamically maintaining the electric power balance, the challenges of controlling and coordinating building equipment and DERs on a large scale must be overcome. Although several control techniques have been proposed in the literature, a significant obstacle to applying these techniques in practice is having access to an effective testing platform. Performing tests at scale using real equipment is impractical, so simulation offers the only viable route to developmental testing at scales of practical interest. Existing power-grid testbeds are unable to model individual residential end-use devices for developing detailed control formulations for responsive loads and DERs. Furthermore, they cannot simulate the control and communications at subminute timescales. To address these issues, this paper presents a novel power-grid simulation testbed for transactive energy management systems. Detailed models of primary home appliances (e.g., heating and cooling systems, water heaters, photovoltaic panels, energy storage systems) are provided to simulate realistic load behaviors in response to environmental parameters and control commands. The proposed testbed incorporates software as it will be deployed, and enables deployable software to interact with various building equipment models for end-to-end performance evaluation at scale.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Extension of Plant Dynamics Code Capabilities for Simulation of TerraPower Pascal Reactor

Argonne National Laboratory has been developing the Plant Dynamics Code (PDC) for design and transient analysis of supercritical carbon dioxide (sCO 2 ) Brayton cycles. In previous analyses with PDC, only indirect sCO 2 cycles, where the heat is being added through a heat exchanger, such as sodium-to-CO 2 HX, were analyzed. Under the U.S. Department of Energy Technology Commercialization Fund (TCF), Argonne cooperated with TerraPower to bring the Plant Dynamics Code to commercial market. The main focus of the TCF project is to extend the application base and the code usability by developing the capabilities to be able to simulate reactor systems with direct sCO 2 cycles. These new PDC capabilities have been developed in application to the TerraPower’s Pascal reactor concept. This report documents the Pascal reactor modeling with the PDC, including simulation of the Pascal split-expansion cycle, the development of the reactor module in the PDC, modeling of the Pascal’s shutdown heat removal system in PDC, and other updates to the code. The report also describes the results of the steady state and transient demonstration of the newly developed code features.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Recovery act: development of design and simulation tool for hybrid geothermal heat pump system

The ground source heat pump (GSHP) system is one of the most energy efficient HVAC technologies in the current market. However, the heat imbalance may degrade the ability of the ground loop heat exchanger (GLHX) to absorb or reject heat. The hybrid GSHP system, which combines a geothermal well field with a supplemental boiler or cooling tower, can balance the loads imposed on the ground loop heat exchangers to minimize its size while retaining superior energy efficiency. This paper presents a recent simulation-based study with an intention to compare multiple common control strategies used in hybrid GSHP systems, including fixed setpoint, outside air reset, load reset, and wetbulb reset. A small office in Oklahoma City conditioned by a hybrid GSHP system was simulated with the latest version of eQUEST 3.7[1]. The simulation results reveal that the hybrid GSHP system has the excellent capability to meet the cooling and heating setpoints during the occupied hours, balance thermal loads on the ground loop, as well as improve the thermal comfort of the occupants with the undersized well field.

15 GEOTHERMAL ENERGY↗

HIGH-RESOLUTION SIMULATIONS OF THE REIONIZATION OF AN ISOLATED MILKY WAY-M31 GALAXY PAIR

We present the results of a set of numerical simulations aimed at studying reionization at the galactic scale. We use a high-resolution realization of the formation of the Milky Way (MW)-M31 system to simulate the reionization of the Local Group. The reionization calculation was performed with the post-processing radiative transfer code ATON and the underlying cosmological simulation was performed as part of the CLUES project (http://www.clues-project.org). We vary the source models to bracket the range of source properties used in the literature. We investigate the structure and propagation of the galactic ionization fronts by a visual examination of our reionization maps. Within the progenitors, we find that reionization is patchy and proceeds locally inside-out. The process becomes patchier with decreasing source photon output. It is generally dominated by one major H II region and one to four additional isolated smaller bubbles, which eventually overlap. Higher emissivity results in faster and earlier local reionization. In all models, the reionization of the MW and M31 are similar in duration, i.e., between 203 Myr and 22 Myr depending on the source model, placing their z{sub reion} between 8.4 and 13.7. In all models except the most extreme, the MW and M31 progenitors reionize internally, ignoring each other despite being relatively close to each other, even during the epoch of reionization. Only in the case of strong supernova feedback suppressing star formation in halos less massive than 10{sup 9} M{sub ☉}, and using our highest emissivity, do we find that the MW is reionized by M31.

79 ASTRONOMY AND ASTROPHYSICS↗

Real-time x-ray studies of indium island growth kinetics

The authors have investigated the early stages of indium island formation and growth by vapor phase deposition on room temperature sapphire using real-time grazing incidence small angle x-ray scattering (GISAXS), followed by ex-situ atomic force microscopy and scanning electron microscopy. The results are consistent with the formation and coalescence of hemispherical islands, as described by Family and Meakin. Monte Carlo simulations of systems of coalescing islands were used to supplement and quantify the results of GISAXS, and a good agreement is seen between the data and the simulations.

42 ENGINEERING↗

Measurement of Photoneutrons from Depleted Uranium and Comparison Study Using MCNPX-PoliMi

High-energy photon interrogation is a nondestructive technique that is used to detect special nuclear materials and characterize nuclear waste. The development of such systems is complex and requires Monte Carlo simulations to optimize system performance. Monte Carlo simulations rely on various scattering, absorption, and photonuclear cross-section data. While the scattering and absorption cross-section data have been extensively studied and validated with experiments, the results obtained from photonuclear simulations are often found to underpredict measured results, indicating uncertainties in the photonuclear cross sections themselves. Thus, there is a need for new measured results that can be used to quantify underpredictions in simulations using photonuclear cross-section data. In the present work, we interrogated depleted uranium with a 9-MV electron linac and detected photoneutrons with trans-stilbene organic scintillators. The measurement of photoneutrons with organic scintillators is challenging due to the presence of the intense photon flux, which causes issues such as pulse pile-up, detector saturation, and poor signal-to-background ratio. To mitigate these challenges, we used iron and polyethylene shielding of varying thicknesses around the depleted uranium target and a neural network–based digital pulse processing algorithm to recover neutron and photon information from piled-up events. Our goal was to compare the measured photoneutron count rate with the simulated rate obtained using the MCNPX-PoliMi transport code. For a light output window of 0.28 to 2.67 MeVee (1.66- to 6.85-MeV proton recoil energy), we found that the simulated count rate obtained using the ENDF/B-VII photonuclear cross-section library underpredicts the measured rate by 32.8% 3.2%. Additionally, we compared the simulated and measured photoneutron light output distributions. For the least thicknesses of shielding, the simulation was found to underpredict measurements in the 0.70- to 2.67-MeVee light output window. For the greatest thicknesses of shielding, the simulation was found to underpredict the measurement across the entire light output window of 0.28 to 2.67 MeVee.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

STOCHASTIC OPTIMAL POWER FLOW FOR REAL-TIME MANAGEMENT OF DISTRIBUTED RENEWABLE GENERATION AND DEMAND RESPONSE (Final Report)

To meet the grand challenge of a sustainable energy future, there has been a surge of interest in renewable energy. Today, the uncertainty associated with renewable resources is handled by using operating reserves. The high penetration of renewable resources, however, introduces difficult-to-control dynamics and challenges for power system operation. Decision support tools are necessary at the bulk system operational level to recognize and efficiently utilize renewable resources and distributed demand response products in concert with traditional grid resources. It is envisaged that responsive load can potentially have very significant cost advantages over either spinning or non-spinning ramping reserve. Critical decisions are made during hour(s)-ahead and real-time power system operation regarding the commitment and dispatch of generators to ensure power delivery is both reliable and economic. These decisions are typically made by a security constrained optimal flow, which determines future generator commitments, dispatches, and ensures adequate reserves are available in the event of a contingency (unexpected outage) or if future system conditions deviate from forecasts. However, security has been always based on a pre-specified subset of contingency constraints whose enforcement does not guarantee security under all possible future possibilities while also giving little or no weight to the likelihood of each contingent event or the severity of its consequences. Existing tools, which are based exclusively on deterministic optimization models, do not yield optimal operational decisions to address these new challenges, in terms of both reliability and cost-effectiveness. This project has focused on developing a stochastic optimal power flow (SOPF) framework, which integrates renewable resource uncertainty, load uncertainty, distributed storage (DS), demand response (DR) products, in a holistic manner to address the uncertainty associated with ever-increasing renewable resources, along with the inclusion of distributed demand response products in future power systems. A proof-of-concept problem was created using the Pennsylvania-Jersey-Maryland (PJM) power system network. Synthetic wind generation was added to the system to simulate 50% wind penetration. A 1-hour test of SOPF operation indicated more than 6% operational cost savings. The project continued by adding the Midwestern Independent System Operator (MISO) as a partner, with focus shifting from SOPF to Stochastic Look-Ahead Unit Commitment (SLAC). Unlike PJM, MISO is faced with significant renewable energy resources within its footprint and is challenged with substantial uncertainty in its operations. The SLAC distinguishes itself from existing tools that operators use. At best, today’s tools solve two to three cases independently, where one or two system parameters, such as forecasted load level (e.g., a low, base, and high forecast), are varied and the resulting scenarios are analyzed independently. The stochastic-based optimization of SLAC leverages statistical information from an ensemble of potential operational scenarios and their respective likelihood. The SLAC output can be translated into valuable information to the operator such as suggested commitments, optimal scheduling and dispatch of resources, reserve requirements at both locational and zonal resolutions, ramping availability and requirements, availability of demand response including operational guidance concerning the near-term and real-time coordination between distributed energy resources, and utilization of distributed storage resources. The developed SOPF/SLAC tool, a stand-alone tool compatible with existing EMSs, will provide system operators with unprecedented visibility, flexibility and predictability to these resources and operational guidance concerning the real-time coordination between DERs and DR/DS products. The game changing and practical impact of this disruptive technology will be dramatic and will usher in a new era in the electric power industry, wherein green energy concepts are fully embraced, and electric power costs are lowered throughout the nation.

42 ENGINEERING↗

Quantum Information for Fusion Energy Sciences (Final Technical Report)

The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical simulation of NSLS-II fast orbit feedback system

We discuss recent work on the construction of a numerical simulation of the fast orbit feedback system at NSLS-II. The simulation operates in the time domain and also includes the spatial domain, i.e. all beam position monitors (BPMs) and correctors. It can accept inputs of the real beam orbit and the measured orbit response matrix. It can also add errors to every stage of the calculation. We present methods to verify the simulation results by comparing the simulation results with experimental data collected at NSLS-II. The results are in very good agreement. The effects of errors from BPMs and correctors, as well as the amplitude of excitation, on the feedback performance, are also explored. This simulation can be used to predict the behavior and performance of the system for a future upgrade.

36 MATERIALS SCIENCE↗

Shaker-structure interaction modeling and analysis for nonlinear force appropriation testing

Nonlinear force appropriation is an extension of its linear counterpart where sinusoidal excitation is applied to a structure with a modal shaker and phase quadrature is achieved between the excitation and response. While a standard practice in modal testing, modal shaker excitation has the potential to alter the dynamics of the structure under test. Previous studies have been conducted to address several concerns, but this work specifically focuses on a shaker-structure interaction phenomenon which arises during the force appropriation testing of a nonlinear structure. Under pure-tone sinusoidal forcing, a nonlinear structure may respond not only at the fundamental harmonic but also potentially at sub- or superharmonics, or it can even produce aperiodic and chaotic motion in certain cases. Shaker-structure interaction occurs when the response physically pushes back against the shaker attachment, producing non-fundamental harmonic content in the force measured by the load cell, even for pure tone voltage input to the shaker. This work develops a model to replicate these physics and investigates their influence on the response of a nonlinear normal mode of the structure. Experimental evidence is first provided that demonstrates the generation of harmonic content in the measured load cell force during a force appropriation test. This interaction is replicated by developing an electromechanical model of a modal shaker attached to a nonlinear, three-mass dynamical system. Several simulated experiments are conducted both with and without the shaker model in order to identify which effects are specifically due to the presence of the shaker. Finally, the results of these simulations are then compared to the undamped nonlinear normal modes of the structure under test to evaluate the influence of shaker-structure interaction on the identified system’s dynamics.

42 ENGINEERING↗

Core Modeling & Decision Support Capabilities: FASTSim, RouteE, T3CO & OpenPATH

These are the slides on NREL's Core Modeling and Decision Support Capabilities project for presentation at the 2023 Annual Merit Review for the U.S. Department of Energy, Vehicle Technologies Office. The project is funded by the Energy-Efficient Mobility Systems program to support four core tools: the Future Automotive Systems Technology Simulator (FASTSim), the Route Energy Prediction Model (RouteE), the Transportation Technology Total Cost of Ownership (T3CO) tool, and the Open Platform for Agile Trip Heuristics (OpenPATH).

ADVANCED PROPULSION SYSTEMS↗

Electrification of residential and commercial buildings integrated with hybrid renewable energy systems: A techno-economic analysis

Here, this study assesses the techno-economic viability of city-scale building electrification, comparing a base case scenario with conventional HVAC systems to an electrified scenario with heat pump systems. EnergyPlus simulations reveal a 28% increase in annual electricity consumption in the electrified scenario, shifting the peak consumption from summer to winter with a 49% surge. Hybrid renewable energy systems (HRESs) are then evaluated using HOMER Pro, considering renewable penetration from 0% to 100%. The net present cost (NPC) of the base case and electrified scenarios range from $\$10.34$ to $\$24.38$ billion and $\$8.79$ to $\$50.31$ billion, respectively. Sensitivity analysis explores the impact of carbon tax, fuel price, and renewable subsidies on HRESs in the electrified scenario. Results indicate that applying a carbon tax up to $\$120$/tonne of CO 2 increases the optimum renewable fraction from 5% to 26%, while higher fuel prices and renewable subsidies further enhance the feasibility by shifting the optimum renewable fraction to 60% and 54%, respectively, demonstrating the interplay between economic factors and renewable integration in electrified urban contexts.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗