Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and Modeling Tools”

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 991 records · Page 55

Active and Transfer Learning of High-Dimensional Neural Network Potentials for Transition Metals

Classical molecular dynamics (MD) simulations represent a very popular and powerful tool for materials modeling and design. The predictive power of MD hinges on the ability of the interatomic potential to capture the underlying physics and chemistry. There have been decades of seminal work on developing interatomic potentials, albeit with a focus predominantly on capturing the properties of bulk materials. Such physics-based models, while extensively deployed for predicting the dynamics and properties of nanoscale systems over the past two decades, tend to perform poorly in predicting nanoscale potential energy surfaces (PESs) when compared to high-fidelity first-principles calculations. These limitations stem from the lack of flexibility in such models, which rely on a predefined functional form. Machine learning (ML) models and approaches have emerged as a viable alternative to capture the diverse size-dependent cluster geometries, nanoscale dynamics, and the complex nanoscale PESs, without sacrificing the bulk properties. Here, in this study, we introduce an ML workflow that combines transfer and active learning strategies to develop high-dimensional neural networks (NNs) for capturing the cluster and bulk properties for several different transition metals with applications in catalysis, microelectronics, and energy storage, to name a few. Our NN first learns the bulk PES from the high-quality physics-based models in literature and subsequently augments this learning via retraining with a higher-fidelity first-principles training data set to concurrently capture both the nanoscale and bulk PES. Our workflow departs from status-quo in its ability to learn from a sparsely sampled data set that nonetheless covers a diverse range of cluster configurations from near-equilibrium to highly nonequilibrium as well as learning strategies that iteratively improve the fingerprinting depending on model fidelity. All the developed models are rigorously tested against an extensive first-principles data set of energies and forces of cluster configurations as well as several properties of bulk configurations for 10 different transition metals. Our approach is material agnostic and provides a methodology to transfer and build upon the learnings from decades of seminal work in molecular simulations on to a new generation of ML-trained potentials to accelerate materials discovery and design.

36 MATERIALS SCIENCE↗

PV Hosting Capacity Estimation: Experiences with Scalable Framework; Preprint

Hosting capacity is an indication of the amount of photovoltaics (PV) can be hosted in a distribution system. This paper presents a framework for estimating distributed PV hosting capacity at scale. We first analyze the key challenges of performing relevant large scale simulation including computational and modeling challenges. Then, we develop two python-based software tools in order to conduct repeatable distribution analyses: Distribution Integration Solution Cost Options (DISCO) for configuring and analyzing simulations, and JADE for parallelizing jobs on HPC clusters. A case study of hosting capacity estimation for SMART-DS SFO 2000+ synthetic feeders is used to demonstrate the capability of the developed framework and tools. The framework and tools can help utilities assess the overall hosting capacity of their service territory, through which the overall upgrade cost can be better planned in order to integrate more PV in the future.

distributed energy resources↗

Future ATM Concepts Evaluation Tool (FACET) Interface Control Document

This Interface Control Document (ICD) documents the airspace adaptation and air traffic inputs of NASA's Future ATM Concepts and Evaluation Tool (FACET). Its intended audience is the project manager, project team, development team, and stakeholders interested in interfacing with the system. FACET equips Air Traffic Management (ATM) researchers and service providers with a way to explore, develop and evaluate advanced air transportation concepts before they are field-tested and eventually deployed. FACET is a flexible software tool that is capable of quickly generating and analyzing thousands of aircraft trajectories. It provides researchers with a simulation environment for preliminary testing of advanced ATM concepts. Using aircraft performance profiles, airspace models, weather data, and flight schedules, the tool models trajectories for the climb, cruise, and descent phases of flight for each type of aircraft. An advanced graphical interface displays traffic patterns in two and three dimensions, under various current and projected conditions for specific airspace regions or over the entire continental United States. The system is able to simulate a full day's dynamic national airspace system (NAS) operations, model system uncertainty, measure the impact of different decision-makers in the NAS, and provide analysis of the results in graphical form, including sector, airport, fix, and airway usage statistics. NASA researchers test and analyze the system-wide impact of new traffic flow management algorithms under anticipated air traffic growth projections on the nation's air traffic system. In addition to modeling the airspace system for NASA research, FACET has also successfully transitioned into a valuable tool for operational use. Federal Aviation Administration (FAA) traffic flow managers and commercial airline dispatchers have used FACET technology for real-time operations planning. FACET integrates live air traffic data from FAA radar systems and weather data from the National Weather Service to summarize NAS performance. This information allows system operators to reroute flights around congested airspace and severe weather to maintain safety and minimize delay. FACET also supports the planning and post-operational evaluation of reroute strategies at the national level to maximize system efficiency. For the commercial airline passenger, strategic planning with FACET can result in fewer flight delays and cancellations. The performance capabilities of FACET are largely due to its architecture, which strikes a balance between flexibility and fidelity. FACET is capable of modeling the airspace operations for the continental United States, processing thousands of aircraft on a single computer. FACET was written in Java and C, enabling the portability of its software to a variety of operating systems. In addition, FACET was designed with a modular software architecture to facilitate rapid prototyping of diverse ATM concepts. Several advanced ATM concepts have already been implemented in FACET, including aircraft self-separation, prediction of aircraft demand and sector congestion, system-wide impact assessment of traffic flow management constraints, and wind-optimal routing.

FACET↗

Low-order nonlinear dynamic model of IC engine-variable pitch propeller system for general aviation aircraft

This paper presents a dynamic model of an internal combustion engine coupled to a variable pitch propeller. The low-order, nonlinear time-dependent model is useful for simulating the propulsion system of general aviation single-engine light aircraft. This model is suitable for investigating engine diagnostics and monitoring and for control design and development. Furthermore, the model may be extended to provide a tool for the study of engine emissions, fuel economy, component effects, alternative fuels, alternative engine cycles, flight simulators, sensors, and actuators. Results show that the model provides a reasonable representation of the propulsion system dynamics from zero to 10 Hertz.

Richard, Jacques C.↗

A Web-Based Negotiation Tool for Conflict Resolution in Upper Class E Traffic Management

In Upper Class E airspace 60,000 feet (or Flight level / FL600), vehicles such as High Altitude Long Endurance (HALE) balloons and slow fixed-wing gliders have diverse vehicle characteristics and limited maneuverability. Due to these unique characteristics of high-altitude operations, a new type of strategic negotiation-based conflict resolution method has been proposed to avoid potential conflict between vehicles in Upper Class E airspace, well in advance of the conflict point and with ample time for bilateral negotiation. This paper introduces the first real-time web-based negotiation tool for high-altitude operations. To facilitate negotiation, this tool incorporates a bilateral negotiation model with an algorithm to assess conflict risks and generate new flight trajectories with a given flight path deviation. This tool allows users to make decisions during negotiation manually and automatically while enforcing the needed constraints for negotiation to be successful. Operators can make decisions at every step while they interact with each other on separate devices. To conduct sensitivity analysis, a mechanism that automates human inputs to the user interface is also developed for this web-based negotiation tool, such that fast-time simulations can be constructed and used to explore various scenarios to gain insights into this web-based negotiation model. Using this tool, a study was conducted to evaluate the impact of different negotiation strategies, vehicle types, vehicle crossing angles, and operator response times during the negotiations on metrics such as total negotiation completion time, number of negotiation rounds, and extra flight distance to avoid the conflict due to negotiation, compared to ones without negotiation. The results suggest that operators benefit from using negotiated flight paths with quick response times, across various crossing angles and vehicle types. Various Negotiation Strategies are investigated to simulate the behaviors associated with different types of negotiations. The findings demonstrate that, in comparison to the conventional method where a single operator assumes full responsibility, negotiation-based strategic deconfliction reduces the total extra flight distance by an average of 24% - 27%. On an individual basis, each operator may be able to save an average of 65% of their extra flight distance.

Upper Class E Traffic Management (ETM)↗

Ablation and Heating During Atmospheric Entry and Its Effect on Airburst Risk

Large meteoroids and asteroids entering the atmosphere endure tremendous heating from the shock heated air, and thereby lose a significant fraction of their mass during atmospheric entry a process known as ablation. The predicted evolution of the asteroids mass as it passes through the atmosphere can affect both the predicted energy deposition profile relevant to an airburst event, or the residual mass that strikes the ground in the case of an impact event. This presentation is divided roughly into two parts. In the first part, an overview of traditional models for heat transfer and ablation that are historically used in the meteor physics community is presented, and the validity in the asteroid entry regime discussed. Sensitivity analyses performed using the recently developed Fragment-Cloud Model (FCM) will be presented which show illustrate the range of sizes and entry parameters for which the predicted asteroid threat is most sensitive to the models for ablation and heat transfer. The second part of the presentation shall focus on recent work done under NASAs Asteroid Threat Assessment Project (ATAP) to develop new models for heat transfer and ablation using high-fidelity numerical simulation in concert with state-of-the-art experiments. Coupled computational fluid dynamics (CFD)radiation transport simulations preformed using the state-of-the-art entry modeling tools at NASA show that, for large meteoroids and asteroids, there can significant attenuation of the heat transfer to the surface (95 in some cases) by the products of ablation. In addition to the heat transfer, new models for the material response and ablation of asteroidal material have been developed [cite]. In the current work, we present finding from recent novel experiments performed in the arc jet facility at NASA Ames, which allows us to, in part, simulate the extreme environment experienced by the asteroid during entry. Briefly, the experimental set-up was comprised of a 1.5 conical article of machined H5 chondrite, which was exposed to a high-enthalpy flow resulting in approximately 4 kWcm2 of heating to the surface. A still frame capture from high-speed video taken during this experiment can be seen in Figure 1. In this figure, we can observe some of the major mechanisms for meteoroid ablation, such as melt flow, spallation (mechanical removal of material), and vaporization. Major findings from this, and other experiments will be discussed, as well progress on utilizing the data from the experiments to inform and develop improved models for ablation.

Asteroids↗

What Machine Learning Can and Cannot Do for Inertial Confinement Fusion

Machine learning methodologies have played remarkable roles in solving complex systems with large data, well-defined input–output pairs, and clearly definable goals and metrics. The methodologies are effective in image analysis, classification, and systems without long chains of logic. Recently, machine-learning methodologies have been widely applied to inertial confinement fusion (ICF) capsules and the design optimization of OMEGA (Omega Laser Facility) capsule implosion and NIF (National Ignition Facility) ignition capsules, leading to significant progress. As machine learning is being increasingly applied, concerns arise regarding its capabilities and limitations in the context of ICF. ICF is a complicated physical system that relies on physics knowledge and human judgment to guide machine learning. Additionally, the experimental database for ICF ignition is not large enough to provide credible training data. Most researchers in the field of ICF use simulations, or a mix of simulations and experimental results, instead of real data to train machine learning models and related tools. They then use the trained learning model to predict future events. This methodology can be successful, subject to a careful choice of data and simulations. However, because of the extreme sensitivity of the neutron yield to the input implosion parameters, physics-guided machine learning for ICF is extremely important and necessary, especially when the database is small, the uncertain-domain knowledge is large, and the physical capabilities of the learning models are still being developed. In this work, we identify problems in ICF that are suitable for machine learning and circumstances where machine learning is less likely to be successful. This study investigates the applications of machine learning and highlights fundamental research challenges and directions associated with machine learning in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Single Primary Heat Extraction and Removal Emulator (SPHERE) Long Duration Testing

For the development of heat-pipe cooled microreactors, it is crucial to thoroughly understand the characteristics and functioning of heat pipes across a wide spectrum of operating conditions. Passive heat removal and its long-term performance stability are critical factors in this context. Enhanced experimental data is vital for evaluating the operational lifespan of alkali metal heat pipes. Idaho National Laboratory (INL) has successfully conducted an extended duration test on a high-performance sodium-filled heat pipe, closely monitoring the axial temperature profile, power supplied by the heaters, and heat removed by a gas-gap calorimeter. The results from this testing provide valuable data that are instrumental in supporting heat pipe validation efforts. Specifically, this data aids in the development and validation of Sockeye, the Multiphysics Object-Oriented Simulation Environment (MOOSE) tool under the US-DOE NEAMS program designed for heat pipe modeling. By comparing experimental results with Sockeye’s predictions, the tool's accuracy and reliability can be assessed and improved, thereby enhancing its capability to simulate heat pipe operations under various conditions.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation of the ERO2.0 code using W7-X and JET experiments and predictions for ITER operation

Abstract The paper provides an overview of recent modelling of global material erosion and deposition in the fusion devices Wendelstein 7-X (W7-X), JET and ITER using the Monte-Carlo code ERO2.0. For validating the modelling tool in a three-dimensional environment, W7-X simulations are performed to describe carbon erosion from the graphite test divertor units, which were equipped in operational phase OP 1.2 and analysed post-mortem. Synthetic spectroscopy of carbon line emission is compared with experimental results from the divertor spectrometer measurement system, showing a good agreement in the e-folding lengths in the radial intensity profiles of carbon. In the case of metallic wall materials, earlier modelling of the Be/W environment in JET and ITER is revisited and extended with an updated set of sputtering and reflection data, as well as including the mixing model for describing the Be/W dynamics in the divertor. Motivated by recent H/D/T isotope experiments in JET, limited and diverted configuration pulses are modelled, showing the expected trend of both Be and W erosion increasing with isotope mass. For the JET diverted configuration pulses, it is shown that Be migrates predominantly to the upper part of the inner divertor where it initially leads to strong W erosion. With longer exposure time, the growth of a Be deposited layer leads to a reduction of W erosion in that region. A similar trend is observed in simulations of the ITER baseline Q = 10 scenario, however with a more symmetric Be migration pattern leading to deposition also on the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Battery State of Health Estimator: Cooperative Research and Development Final Report

NREL has developed a software tool to enable Renewance to estimate the degradation of batteries from basic information such as the type of battery and the application of that battery during its first life, so that used batteries may be evaluated for potential repurposing at low cost. This software tool utilizes NREL's BLAST-Lite battery degradation modeling code, which was updated with additional models for commercially produced lithium-ion batteries as a part of this CRADA. The software tool enables users to input details such as battery type and application so that lifetime estimates can be made without any programming or expert battery knowledge. The application input loads in saved values for parameters such as cycles per year, depth-of-discharge, and other battery operating parameters from a file defined by Renewance. These parameters may be modified to refine simulations for specific batteries. The software tool also incorporates a degradation model optimization tool, whereby existing battery degradation models may be tuned according to measured battery health. This ensures that new models still predict degradation behaviors expected from a certain battery chemistry, but with the overall degradation rate tuned to a specific battery make and model. The new model can then be saved for estimating the degradation of other similar batteries. An additional task was planned to utilize machine-learning to enable battery health diagnosis from rapid EIS measurements to accelerate the screening of used batteries. This task was not completed due to lack of available data for training a machine-learning model. CRADA benefit to DOE, Participant, and US Taxpayer: Further development of open-source software tool BLAST-Lite for predicting the lifetime of commercially produced Lithium-ion batteries (NREL SWR-22-69).

25 ENERGY STORAGE↗

Assimilation of Multiscale Data into Multifidelity Biogeochemical Models (Final Report)

Quantitative predictions of subsurface processes rely on computational models that capture, with different degrees of fidelity, complex interactions between hydrologic and biogeochemical processes. Molecular- and pore-scale models provide a high-fidelity representation of these processes but are impractical at the field scale. Reduced complexity (e.g., field-scale or data-driven) models sacrifice some degree of fidelity in favor of computational efficiency. Uncertainty and assimilation of data into model predictions, pose a question of model selection: Given a significant difference in computational cost, can a lower-fidelity model be preferable to its higher-fidelity counterpart? Since a predictive model must be computable in reasonable time it has to operate at the field scale, with molecular- and pore-scale data (obtained either from simulations or measurements) determining both the model's structure and parameters. Availability of such multi-resolution data raises a question hitherto undressed in hydrogeology: Can a coarse model dynamically ``learn'' its own structure (e.g., adjust the reaction pathways in its transport module) as more fine-scale data become available during simulations? Our results in multifidelity simulations, data assimilation, and machine learning led us to hypothesize that the answer to these questions is ``yes''. The overarching goal of this project was to confirm this hypothesis by developing scalable, computationally efficient tools for assimilation of data into multiscale models, in which fine-scale data and simulations dynamically inform and autonomously modify coarse-scale models. Our research led to nine manuscripts, three of which have been published and the other six are currently under review.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Control of machining-induced residual stress via tool geometry and process parameter modification

Distortion generated in machined, monolithic, thin-walled aerospace components due to residual stresses leads to significant material and economic waste in the manufacturing industry. Inherent residual stress (IRS) present in stock materials combines with machining-induced residual stress (MIRS) to influence the final machined part distortion. It is hypothesized that MIRS can be controlled, based on the part geometry, through deliberate cutting tool geometry and process parameter modifications to negate the effect of IRS on distortion, consequently resulting in distortion-free parts. A finite element (FE) orthogonal cutting model is developed to study how tool geometry and process parameters influence machining-induced residual stress (MIRS). Orthogonal cutting experiments are performed on Al 7075-T651 samples to measure cutting forces and MIRS. A cutting force dynamometer is used to measure forces during cutting and a novel digital image correlation (DIC) based hole drilling technique is employed to measure the near-surface residual stress (RS) in the cut samples. These data are subsequently utilized to validate the FE prediction model. Various cases of cutting simulations involving different depths of cut, tool tip radii, and rake angles are performed to study their effects on RS. Similar to prior literature, increasing the depth of cut, tool tip radius, or rake angle is found to promote the formation of near-surface tensile stresses. The competing effects of material plowing and temperature are shown to determine the type of RS at the end of the cut. Moreover, a window of variation of RS (up to ± 400 MPa) is estimated within the given range of conditions, allowing for the control of MIRS through tool and process modification.

Mathews, Ritin [ORNL] (ORCID:0000000301440828)↗

Multibody Based Digital Astronaut Dynamics Simulation

BACKGROUND: This study provides the Software, Robotics, & Simulation Division at the NASA Johnson Space Center with a verification tool for multibody dynamics simulation requiring human motion. The motivation stems from current studies of several Vibration Isolation & Stabilization(VIS)system designs that attenuate the moments and forces which would be transmitted to a spacecraft during an exercise. A multibody dynamics model for a proposed VIS was available previously[1], therefore modeling of the VIS was not needed for this work. The interest here is in creating the multibody dynamics model of an astronaut in motion which may be utilized independently or while attached to a mechanism. An existing simulation [2] that utilizes OpenSim [3,4] and an in-house multibody dynamics package (MBDyn) [5] is used in order to verify the astronaut model. The main advantage this model will have over the existing simulation is that everything will be processed in one tool. METHODS AND RESULTS: Creating the simulation required; estimation of Body Segment Inertial Parameters (BSIP),a multibody model of the human-VIS system, joint acceleration profiles, and input files for MBDyn, which is used for this analysis. The scaling factors provided by Dumas et al. [6] are utilized in estimating the BSIP. Anthropometric data are used for estimating these parameters, the Anthropometric Survey of US Army Personnel (ANSUR II) [7] was the source. The astronaut model consists of 15 bodies, 14 joints and 32 degrees of freedom, with the dynamics topology generated using the center of mass locations and anthropometric data. MBDyn has an option for prescribed joint motion (PJM), which requires joint acceleration data as input. The joint angle data is first obtained from a motion capture system and then processed through code that has been created to generate approximate joint acceleration profiles. The topology tree for the astronaut model begins at the right foot up to the pelvis where there is one branch for going down the left leg and another for the torso. The torso branch leads to branches for the arms and a leaf body for the head/neck segment. For attachment to the VIS, the heel of the right foot is connected to the VIS platform through a fixed joint, resembling a foot restraint. The left foot does not attach to the platform in order to prevent a system with a closed loop. Topology and symmetry of the astronaut model were verified through kinematic analysis. Further verification of the forces and moments transmitted to the VIS were verified against the existing simulation. There was a satisfactory level of agreement when testing a simple motion, for example, rocking back and forth. Full exercise motions are to be tested soon. The main outcome has been a novel application of MBDyn for biomechanics modeling that is now available for dynamic simulations involving human motion. The estimation of BSIP was another useful result of this study, requiring only 15 inputs for generating mass properties of a theoretical astronaut model. Expansion on this work is possible by going through an alternative route in obtaining the joint motion data. Instead of high-tech and often expensive motion capture systems, an individual may watch videos with high focus and at a slow motion for each individual segment in order to determine the initial and final time and angle for that specific degree of freedom. Synthetic trajectories may also be created if there is no video reference available.

F N Matari↗

Biorefinery upgrading of herbaceous biomass to renewable hydrocarbon fuels, Part 1: Process modeling and mass balance analysis

Process simulation has long been a well-established tool to track key operational, design, and mass and energy balance metrics for pre-commercial technologies such as advanced lignocellulosic biofuels. While tools such as this are well-documented in the public literature around 2nd-generation cellulosic ethanol technologies (which have been scaled up to commercial deployment to some degree over the past decade), such models and analysis information remain more sparse for more complex biorefinery pathways focused on producing drop-in hydrocarbon fuels and blend-stocks, particularly regarding information required to support air emissions or other environmental analysis. In this work, we summarize key details for an established "design case" modeling the conversion of herbaceous lignocellulosic biomass into a renewable diesel hydrocarbon blend-stock based on a representative lipid pathway from oleaginous yeast. The process is based on a biochemical deconstruction and upgrading approach utilizing deacetylation and dilute acid pretreatment, followed by enzymatic hydrolysis, fermentation, and catalytic upgrading of hydrolysate sugars to fuels. We provide key mass and energy balance outputs from the process models, with accompanying stream tables and component-level flowrates. A total of 12 model scenarios are presented spanning two feedstocks, three biorefinery scales, and two processing approaches for the lignin/residual solids waste streams. This "Part 1" manuscript presents the resulting impacts across the 12 cases on fuel yields and key output streams, focused here on direct biorefinery air emissions for selected components including CO 2 as well as sulfur (SOx) and nitrogen oxides (NOx). In the context of cleaner production, the latter focus on selected biorefinery air emission outputs establishes an initial baseline estimate and accompanying framework of the model cases, upon which an accompanying "Part 2" study will build to refine the values for these and other air pollutants across these scenarios, also considering mitigation options to comply with applicable regulatory standards. We also highlight further optimization opportunities based on potential tradeoffs identified here between air emissions versus life-cycle greenhouse gas profiles attributed to the disposition of lignin/residual solids.

09 BIOMASS FUELS↗

BiRD (BioReactorDesign) [SWR-24-35]

Numerical simulations can accelerate scale up and optimization of bioreactors. This project provides a toolbox to facilitate such numerical simulations in OpenFOAM. The tool box contains 1) mesh generations tools for several typical bio reactor types in OpenFOAM 2) Inverse modeling tools to calibrate models using available experimental data 3) post processing tools that accelerate and streamline the analysis of the numerical simulations.

Hassanaly, Malik↗

Simulation of the Impact of Point Defects and Edge Dislocations on X-Ray Diffraction in Hexagonal (Ni,Co) 1+2 x Ti 1– x O 3 Thin Films

In this work, a computer code for simulating high-resolution X-ray diffraction (XRD) data from disordered crystals with arbitrary spatial composition and local lattice parameters is developed. Simulated patterns are compared with the experimental data collected on a single phase, highly crystalline (Ni 0.42 Co 0.58 ) 2.22 Ti 0.39 O 3 solid-solution thin film exhibiting a large number of subsidiary minima. As a case study, the edge dislocations in hexagonal (Ni,Co) 3 O 3 thin films with Burgers vector parallel and perpendicular to the a-axis and c-axis, respectively, are modeled. No peak profiles are assumed and thus issues related to profile fitting are avoided. Both macroscopic features, such as film thickness, and atomic-scale structure, such as dislocations, are simultaneously modeled. Simulations are run on a desktop machine. Commonly applied database-based phase identification routines can result in wrong phase identification, unless intensities are properly modeled. Modeling tools of diffractometer manufacturers are compared with the present approach. An example of how simulated reciprocal lattice patterns can be used to choose relevant measurement geometries in defected thin films is given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Methodology to evaluate design modifications intended to eliminate frosting and high discharge temperatures in air-source heat pumps (ASHPs) in cold climates

Air-source heat pumps (ASHPs) operating in cold climates experience problems with frosting and high refrigerant temperatures. These problems increase energy consumption, and their severity depends on the climatic conditions. In the present paper, a methodology for identifying the prevailing problem between frosting and high discharge temperatures is presented. Three performance indices, the frosting index (FI), the discharge index (DI), and the total loss index (TLI), are proposed to quantify the impacts of frosting and high discharge temperatures on the annual performance of ASHPs in different climatic conditions. The FI and DI show which problem (frosting or high discharge temperature) dominates, and the TLI indicates the combined effect of frosting and high discharge temperatures on the performance of an ASHP. A thermodynamic model of an ASHP coupled with the TRNSYS building simulation tool is used to estimate the performance of an ASHP and the proposed loss indices to estimate the impact of both frosting and high discharge temperatures for 45 cities in Canada. The results can be extended to other parts of the world that experience similar climatic conditions The results reveal that in cities in ASHRAE climatic zones 5 and 6 (classified as cold regions) where the ambient air temperatures are predominantly between -15 °C to 6 °C, ASHPs are heavily impacted by frosting. The problem of high discharge temperatures in ASHPs is predominant in cities in climate zones 7 and 8 (classified as very cold and subarctic regions) where the temperatures are frequently below -20 °C in winter. Among the cities considered, St. John, NL has the highest fraction of heating hours experiencing frosting (90 %), where the annual increase in energy consumption due to frosting is 13.5 % of the annual heating energy consumption. The highest annual increase in energy consumption due to high discharge temperatures is in Isachsen, NU (zone 8), where the increase is 30 % of the annual heating energy consumption. Based on the proposed indices, another index called the performance gain index (PGI) is created, which can be used as a first step to assess the energy-saving potential of design modifications applied to ASHPs to solve the problems of frosting and high discharge temperatures. The PGI will aid in developing climate specific ASHPs. One possible design modification is the use of a two-stage ASHP with an economizer. It is observed that the two-stage ASHP with economizer can mitigate high discharge temperatures and improve performance in very cold and subarctic regions (zones 7 and 8). However, it is not as beneficial in zones 5 and 6, where the impact of high discharge temperatures on performance is minimal and frosting dominates. Finally, a case study, using the PGI to evaluate the economic and environmental effectiveness of a two-stage ASHP with economizer is presented for the city of Saskatoon.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗