Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservative data transfer”

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 217 records · Page 12

Ride Pingo to Transit

In this project, we developed an on-demand microtransit first- and last-mile service. To integrate the service with fixed-route transit, we developed a feature called Transit Connect that prioritized riders’ on-time arrival at the transit station over other service requirements. We first prototyped service and related algorithms in a simulated environment, and then piloted the service in the city of Kent, Washington. Our algorithm incorporates request-specific hard drop-off deadlines to ensure timely arrivals for transit transfers. In the pilot, these constraints were obtained from GTFS Realtime data to accurately determine the schedule of the transit and the location of the stations. This approach introduced the ability to accept or decline new requests based on the timing of transit connections for these new requests and connection status of onboarding customers. The pilot (called “Ride Pingo to Transit”) deployed a fleet of three 14-person vans, ran from September 2021 to March 2023, and served a total of 21,329 trips. Transit Connect was offered for drop-offs at both Kent Station and the Kent Valley hub. In total, 2,844 such trips were completed. This dataset was collected from our pilot, which includes the following: - Requests: List of all trip requests, including those that were actually served and those not materialized. - Fleet: Daily vehicle service logs. - Service details: Daily vehicle stop logs (boarding and alighting). - Trip types: First mile, last mile, or point-to-point. ![pingo to transit](pingo-to-transit.jpg)

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Deploying Traffic Smoothing Cruise Controllers Learned from Trajectory Data

Autonomous vehicle-based traffic smoothing con- trollers are often not transferred to real-world use due to challenges in calibrating many-agent traffic simulators. We show a pipeline to sidestep such calibration issues by collecting trajectory data and learning controllers directly from trajectory data that are then deployed zero-shot onto the highway. We construct a dataset of 772.3 kilometers of recorded drives on the I-24. We then construct a simple simulator using the recorded drives as the lead vehicle in front of a simulated platoon consisting of one autonomous vehicle and five human followers. Using policy-gradient methods with an asymmetric critic to learn the controller, we show that we are able to improve average MPG by 11% in simulation on congested trajectories. We deploy this controller to a mixed platoon of 4 autonomous Toyota RAV-4’s and 7 human drivers in a validation experiment and demonstrate that the expected time-gap of the controller is maintained in the real world test. Finally, we release the driving dataset [1], the simulator, and the trained controller at https://github.com/nathanlct/trajectory-training-icra.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

PemNet: A Transfer Learning-Based Modeling Approach of High-Temperature Polymer Electrolyte Membrane Electrochemical Systems

Widespread adoption of high-temperature electrochemical systems such as polymer electrolyte membrane fuel cells (HT-PEMFCs) requires models and computational tools for accurate optimization and guiding new materials for enhancing fuel cell performance and durability. Furthermore, while robust and better suited for extrapolation, knowledge-based modeling has limitations as it is time-consuming and requires information about the system that is not always available (e.g., material properties and interfacial behavior between different materials). Data-driven modeling, on the other hand, is easier to implement but often necessitates large datasets that could be difficult to obtain. In this contribution, knowledge-based modeling and data-driven modeling are combined by implementing a few-shot learning (FSL) approach. A knowledge-based model originally developed for a HT-PEMFCs was used to generate simulated data (887,735 points) and used to pretrain a neural network source model tuned via a genetic algorithm-based AutoML. Then, experimental datasets from HT-PEMFCs with different materials and operating conditions (~50 points each) were used to train six target models via FSL. Models for the unseen data reached high accuracies in all cases (rRMSE < 10%).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

Characterizing Current THD’s Dependency on Solar Irradiance and Supraharmonics Profiling for a Grid-Tied Photovoltaic Power Plant

The rapidly increasing distributed energy resources (DERs) in power systems are now getting interconnected to set community grid structures, where power quality will be a major concern. The grid-to-grid (G2G) bidirectional power transfer among the distribution microgrid will not be considered commercially feasible unless the upstream harmonics are under the limits. The aggregation of such harmonics, measured as total harmonic distortion (THD), is feared to be beyond tolerable limits with the progression of rooftop grid-tied PV-like installations. Hence, this THD needs to be characterized with DER generation end variables. In this work, the photovoltaic (PV) DERs’ dependency on environment variables such as irradiance was profiled in the context of generating and injecting harmonics into the grid. A mathematical model of a grid-tied three-phase PV DER was developed as part of this correlation characterization, matching the fundamental unit structure of a 1.4 MW solar canopy located on the Florida International University (FIU) Miami campus. To determine the qualitative association with produced THD patterns, the model was evaluated with various irradiance settings. A real-time digital simulation (RTDS) platform was used to verify it. Following this confirmation, sets of data from power quality meters at the point of common coupling and FIU field sensors were utilized to validate further the correlation model. The results showed that the grid current’s THD exhibited a high correlation with the irradiance profile and its variation over time. The early morning and late afternoon periods of the day, associated with a low irradiance, constantly had higher harmonics generated from the PV DER. The midday THD was rather rational with partial shadings, hence a geolocation-dependent factor. These findings were verified by an RTDS and validated by real field data. In quantifying the THD injected by a single DER at a high-frequency (2–150 kHz) supraharmonics (SH) level, a 3% peak increment in magnitude was observed from the high-fixed to the low-fixed irradiance profile. The correlation characteristics depicted that the hybrid microgrid suffered from a daytime-dependent harmonic insertion from the grid-tied DER. This is a global problem unless specific measures are taken to mitigate the harmonics. The electrically notorious higher-frequency SH was found to increase proportionally. The G2G power transfer can be limited because of the higher THD in the early morning and late afternoon, which will also worsen because the numbers of grid-tied PV DERs (i.e., rooftop solar and industrial solar) are likely to increase rapidly soon. The community grid structure can thus have a controlled harmonics filtration setup purposefully designed to address the findings of this work, which also fall within the scope of our future research.

14 SOLAR ENERGY↗

Daylight simulation workflows incorporating measured bidirectional scattering distribution functions

Daylight predictions of architectural spaces depend on good estimates of light transfer through skylights, windows and other fenestration systems. For clear glazing and painted surfaces, parametric transmission and reflection models have proven adequate, but there are many cases where light-scattering, semi-specular shading and daylighting materials defy simple characterization. Something as commonplace as fabric roller shades and venetian blinds may turn daylight prediction into guesswork, and numerous advanced systems on the market tuned specifically to enhance daylight are not sufficiently characterized to distinguish their performance. In this paper, we describe new tools available to handle novel and specialized fabrics, materials, and devices using data-driven modelling of bi-directional scattering distribution functions (BSDFs). These representations are usually tabulated at constant or adjustable angular resolution for efficient point-in-time and annual daylight simulations. We describe a variety of BSDF simulation workflows, including some of the tools and methods that make advanced analysis possible, and highlight some of the current challenges. We conclude with a discussion of future work and how such data might be created and shared worldwide.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Impacts of Sub‐Grid Topographic Representations on Surface Energy Balance and Boundary Conditions in the E3SM Land Model: A Case Study in Sierra Nevada

Abstract Sub‐grid topographic heterogeneity has large impacts on surface energy balance and land‐atmosphere interactions. However, the impacts of representing sub‐grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub‐grid topographic representations on surface energy balance, turbulent heat flux, and scalar (co‐)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub‐grid topographic representations in ELM were compared: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in mean values and standard deviations become small when aggregated to the grid scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and boundary conditions. These results underline the importance of representing sub‐grid topographic heterogeneities in LSMs and motivate future research to understand the sub‐grid topographic effects on land‐atmosphere interactions over mountainous areas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Building envelope anomaly characterization and simulation using drone time-lapse thermography

Defects in building envelopes deteriorate over time without being visible to the human eye, while significantly impacting energy performance due to unaccounted heat transfer. Defects can be characterized in the infrared (IR) spectrum. However, IR readings are typically recorded at singular points in time, when in several cases anomalies can only be revealed at specific times of the day, possibly in different seasons of the year. This paper presents a novel workflow for 3D envelope defect characterization and modeling using aerial time-lapse IR data collection using drones. A comprehensive envelope thermal profile is developed for a case study building employing the photogrammetry software Agisoft Photoscan, which generates temporal IR inspections of building skins using multiple thermography orthomosaics. Point-cloud data is then translated into a CAD model and thermal zones for whole Building Energy Modeling (BEM) using Honeybee as a frontend to EnergyPlus to showcase the potential of inclusion of detailed 4D data. Envelope contributions in this case study’s anomalies showed heat losses of 6447.6 kWh, and Energy Use Intensity (EUI) differences of ~2 kWh/m 2 /year from the baseline. Finally, why there is currently little translation of this work in BEM software is discussed, while identifying limitations and future research in the employment of time-lapse thermography using drones for more accurate building envelope inspection and modeling.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Experimental Characterization of High-Surface Area Thermal Energy Storage

There is growing interest in energy storage technologies due to the expansion of renewable energy sources that are inherently intermittent and the increasing frequency of extreme weather events that disturb the power grid. Power consumption in buildings makes up approximately 76% of all electricity usage on the grid and is primarily used for thermal applications such as space conditioning, hot water, and cooking. This makes thermal energy storage (TES) an ideal solution for many of these applications. Many TES technologies rely on latent energy storage, which utilizes the melting/solidification of phase change materials (PCM) to store energy. Typically, TES designs suffer from low power density due to their low inherent thermal conductivity. This limitation makes the deployment of TES in active applications difficult as the ease of access to energy is essential for effective use. Common routes for improving power density include high thermal conductivity additives or extended features such as fins that increase cost. This study presents an alternative approach to improving performance through increasing the overall surface area to volume ratio of the device, to increase the available area for convection to occur between the working fluid and PCM. In the study a commercial PCM was selected with a transition temperature ideal for space heating applications. The heat exchanger design utilizes a unique application of triply periodic minimal surfaces (TPMS) for macro-encapsulation of the PCM. The use of TPMS for heat exchangers has been growing in interest due to their high-surface area to volume ratios, which were previously unmanufacturable until the development of additive manufacturing. A modular system was designed and manufactured with a stereolithography resin printing system that is then backfilled with PCM. An experiment test set-up is designed to test the charge and discharge performance of the thermal storage using a conditioned air stream. The pressure drop of the design is tested across a variety of flow rates. When compared to existing experimental data within literature, there is excellent agreement based on the Reynolds number at similar hydraulic diameters. Several inlet temperatures are tested at consistent temperature differences from the phase change temperature for both charging and discharging. Additionally, the volumetric flow rate is varied for each temperature set point. It was found that increasing flow rate had diminishing returns in reducing the overall charge time of the TES. The temperature delta from the melting temperature was the primary contributor to the change in average heat flux with limited variation in average heat transfer rate between charging and discharging at similar inlet temperatures and flow rates. The TPMS heat exchanger design has a high air-side pressure drop but it provides high heat transfer rates. This helps maintain a high outlet temperature during discharge, which is important to thermal comfort applications. The design, manufacturing, and experimental characterization of the TES device will be presented as part of this study.

25 ENERGY STORAGE↗

Enhancement of phase change material hysteresis model: A case study of modeling building envelope in EnergyPlus

Nowadays, buildings are expected to offer demand side services to the power grid to enhance the electrical load flexibility, which leads to the concepts of grid-interactive efficient buildings (GEBs). Phase change material (PCM)-based thermal energy storage has seen increasing attention in recent years for peak load shifting of grid-interactive efficient buildings (GEBs). Numerical models are critical tools for design and evaluation of PCM-integrated systems. Most industrial-grade PCMs are reported to melt/freeze over a temperature range instead of at a unique temperature. Such thermal hysteresis effect significantly affects the reliability of simulation results because not only the heat transfer process depends on melting and freezing temperatures, the PCM thermal properties change significantly during the phase change process as well. This study is aimed to develop a model for the PCMs used in the building envelope with the capability to accurately simulate hysteretic behaviors. Further, this model is based on a two-phase assumption and is implemented in a whole building energy performance simulation program (i.e., EnergyPlus). A comparison between numerical results and experimental data shows that during a complete phase transition, the two-phase model could achieve a good agreement with the experimental data. During a partial phase transition, the two-phase model could lead to significant improvements compared to other alternative PCM models, including the existing PCM model in EnergyPlus. Last, whole building simulations were performed to study this model's performance regarding heating/cooling loads and zone mean air temperature of a given building. The results show that the difference in hourly heating/cooling loads introduced by the models was less than 1% in design conditions, while significant changes were observed in both hourly heating/cooling loads and zone mean air temperature when the PCM envelope underwent partial phase transition processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Energy Analysis of Combi Heat Pump System Configurations for Space Conditioning and Domestic Hot Water Heating in Residential Buildings

Combi heat pump systems, also referred to multifunctional variable refrigerant flow heat recovery (MF-VRFHR) systems, are specifically designed for residential applications to manage both space conditioning and domestic hot water (DHW). They have attracted attention due to their potential for energy conservation through heat recovery. The incorporation of a hot water tank introduces various system configurations, each characterized by distinct pros and cons related to energy efficiency, system stability, and maintenance. Despite this, a critical gap exists as the specific energy performance remains unquantified under diverse operational modes (e.g., heating mode and heat recovery mode). This paper aims to bridge this gap by conducting a comprehensive comparative analysis of two prevalent system configurations while considering feasible proposed control logics. Configuration 1 integrates a separate hot water tank and a refrigerant-to-water heat exchanger (HEX), also known as a Hydro Kit while Configuration 2 incorporates a refrigerant-wrapped hot water tank. To facilitate this analysis, we developed high-fidelity system models for both configurations in Modelica, capturing system dynamics and detailed control sequences effectively. These system models were built upon the TIL library for HVAC equipment components and the Buildings library for residential building thermal load calculations. The validation of the simulation testbed utilized data from experiments conducted in the PNNL lab home for Configuration 1. To establish the simulation testbed for Configuration 2, we extended the modeling setup derived from Configuration 1. This extension specifically involved substituting the separate hot water tank and Hydro Kit with a refrigerant-wrapped hot water tank of similar sizing sourced from an actual product. The simulation analysis of heating-only and heat recovery modes reveals that Configuration 2 not only saves energy and maintains warmer tank temperatures but also demonstrates faster water heating capabilities. This is attributed to decreased energy loss and improved heat transfer. The study encompasses a wide range of scenarios, considering diverse thermal loads and water usage patterns across heating and heat recovery modes. Overall, the comprehensive results indicate that Configuration 2 achieves energy savings ranging from 3.5% to 12.2% compared to Configuration 1, depending on factors such as water usage patterns, thermal loads, and operational modes.

Configuration, Comparison, Multi-functional, Resid↗

The N-terminus of varicella-zoster virus glycoprotein B has a functional role in fusion

Varicella-zoster virus (VZV) is a medically important alphaherpesvirus that induces fusion of the virion envelope and the cell membrane during entry, and between cells to form polykaryocytes within infected tissues during pathogenesis. All members of the Herpesviridae , including VZV, have a conserved core fusion complex composed of glycoproteins, gB, gH and gL. The ectodomain of the primary fusogen, gB, has five domains, DI-V, of which DI contains the fusion loops needed for fusion function. We recently demonstrated that DIV is critical for fusion initiation, which was revealed by a 2.8Å structure of a VZV neutralizing mAb, 93k, bound to gB and mutagenesis of the gB-93k interface. To further assess the mechanism of mAb 93k neutralization, the binding site of a non-neutralizing mAb to gB, SG2, was compared to mAb 93k using single particle cryogenic electron microscopy (cryo-EM). The gB-SG2 interface partially overlapped with that of gB-93k but, unlike mAb 93k, mAb SG2 did not interact with the gB N-terminus, suggesting a potential role for the gB N-terminus in membrane fusion. The gB ectodomain structure in the absence of antibody was defined at near atomic resolution by single particle cryo-EM (3.9Å) of native, full-length gB purified from infected cells and by X-ray crystallography (2.4Å) of the transiently expressed ectodomain. Both structures revealed that the VZV gB N-terminus (aa72-114) was flexible based on the absence of visible structures in the cryo-EM or X-ray crystallography data but the presence of gB N-terminal peptides were confirmed by mass spectrometry. Notably, N-terminal residues 109 KSQD 112 were predicted to form a small α-helix and alanine substitution of these residues abolished cell-cell fusion in a virus-free assay. Importantly, transferring the 109 AAAA 112 mutation into the VZV genome significantly impaired viral propagation. These data establish a functional role for the gB N-terminus in membrane fusion broadly relevant to the Herpesviridae .

59 BASIC BIOLOGICAL SCIENCES↗

Comparative Analysis of HEATNETS for Geothermal Network Performance

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification and validation process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study focuses on a validation of HEATNETS using results from TRNSYS. HEATNETS is not intended to replace other modeling tools, but this work demonstrates, via a comparison with an industry standard code, that HEATNETS can be a unique, high-level and rapid modeling tool for estimating the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗

Energy Exascale Earth System Model v2.1.0

First release of version 2.1 of the Energy Exascale Earth System Model. [ATM] The atmosphere component remains EAM. There are no major changes in the default configuration since 2.0. New features include: A semi-lagrangian tracer transport for theta-l dycore, a new algorithm for finding the tropopause, new RRM mesh configurations. Add and update SSP370 and SSP585 cases. Restore the FIDEAL case. [LAND] The land component is ELM. There are no major changes in the default configuration since 2.0. Several option features have been added including: implementation of topography-based subgrid structure (topounits) and accompanying parameterizations and atmospheric forcing downscaling methods; a new plant hydraulics scheme; two-way land-river hydrological coupling through the infiltration of floodplain water; an implementation of perennial crops; updates to the SNICAR-AD snow radiative transfer model; and implementation of soil erosion and sediment yield in ELM-Erosion. Each of these new changes is modular in design and can be turned on or off as the user specifies; they are currently being tested in different “BGC” configurations. [OCEAN] The ocean component remains MPAS-Ocean. Major change since version 2.0 include the addition of the Fox-Kemper et al. 2011 parameterization for submesocale eddies, a correction for barotropic thickness consistency that reduces divergence noise, and the addition of an ocean carbon conservation analysis member. [SEAICE] The sea-ice component remains MPAS-Seaice. Major changes since version 2.0 include: A correction to how shortwave parameters are interpolated in the snicar-ad 5-band radiation scheme, the addition of a sea ice carbon conservation analysis member, updates to the default sea ice biogeochemistry namelist parameters to be consistent with version 2.0 improvements to nitrogen cycling and a correction in the ice-ocean dissolved organic nitrogen coupling. [LAND ICE] The land-ice component remains MPAS-Albany-landIce (MALI). Major changes since 2.0 include an update to the MALI version and the Greenland mesh.[RIVER] The river model is MOSART. There are no major changes in the default configuration since 2.0. A major new optional feature is two-way river-ocean hydrological coupling between MOSART and MPAS-O. This change can be turned on or off as the user specifies, and is being tested in different configurations. [COUPLER] The coupler remains cpl7/MCT. Major changes since version 2.0 include: Carbon budget calculated when heat/water budgets active. Fix a bug in land-atm fluxes for tri-grid configurations. [OTHER] a small bug in the zenith angle calculation was fixed in the data models.

ECP↗

Can a variable alpha induce limit cycle behavior and exponential luminosity decay in transient soft x ray sources?

There has been, recently, a revival of the stability problem of accretion disks. Much of this renewed interest is due to recent observational data on transient soft X-ray novae, which are low-mass X-ray binaries. It is widely believed that nonsteady mass transfer from the secondary onto the compact primary, through an accretion disk, is the reason for the observed spectacular events in the form of often repetitive outbursts, with recurrence times ranging from 1 to 60 yr and duration time on the scale of months. Though not having reached yet a consensus about the nature of the mechanism that regulates the mass transfer, the disk thermal instability model seems to be favored by the fact that the rise in the hard X-ray luminosity is prior to the rise in the soft X-ray luminosity, while the mass transfer instability model seems to be hindered by the fact that the luminosity during quiescence is unable to trigger the thermal instability. However, it should be stressed that, remarkably, the X-ray light curves of these X-ray novae all show overall exponential decays, a feature quite difficult to reproduce in the framework of the viscous disk model, which yields powerlike luminosity decay. Taking into account this observational constraint, we have studied the temporal evolution of perturbations in the accretion rate, under the assumption that alpha is radial and parameter dependent. The chosen dependence is such that the model can reproduce limit cycle behavior (the system is locally unstable but globally stable). However, the kind of dependence we are looking for in alpha does not allow us to use the usual Shakura and Sunyaev procedure in the sense that we no longer can obtain a linearized continuity equation without explicit dependence on the accretion rate. This is so because now we cannot eliminate the accretion rate by using the angular momentum conservation equation.

Meirellesfilho, C.↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗

Design Requirements and Software Specification for the Autonomous Energy Management Software System for Small Commercial Buildings

Commercial buildings are responsible for approximately 20 percent of the total United States energy consumption and greenhouse gas emissions. Over 85 percent of these buildings lack building automation systems. Many of these buildings are small (<50,000 square feet), underserved, and use rooftop units for heating, ventilation, and air-conditioning needs. Because these buildings lack proper energy management systems, they have several operational deficiencies that lead to excess energy consumption. Studies have shown that managing the rooftop units heating and cooling set points, schedules, setbacks, and optimal start can result in 20 to 25 percent reduction in electricity consumption in small commercial buildings. In addition, improving demand flexibility of these buildings will result additional cost savings for the building owner. Therefore, the Department of Energy’s Building Technologies Office approved a project to address the needs for small commercial buildings. The project is led by Pacific Northwest National Laboratory (PNNL) with Intellimation LLC as the cooperative research and development agreement partner. The primary goal of the project is to develop and validate an autonomous energy management software (AEMS) system that will continuously optimize small commercial building operations by minimizing energy consumption and cost, while providing a solution for maximizing decarbonization benefits from electrification of buildings. The work will leverage the vast experience of PNNL research and development staff who have over two decades of experience in developing and successfully transferring software technologies to the private sector. This solution will be jointly developed with Intellimation, a company that plans to use it to scale their building energy efficiency (EE) and grid services offering. Widespread deployment of the AEMS system will improve the EE and demand flexibility of the building commercial building stock. It should also support cities and states in meeting their climate change mitigation goals. This document describes the various EE and grid service features of the AEMS system, infrastructure and data required to implement those features, and how the features should be automated. It also details how the various features will be tested and validated, including field validation. The document also details what flexibility the users have and how they will be able to leverage those capabilities exercise those. The intent is to create an AEMS system that would support scalable deployment, requires minimal configuration, and is easy to maintain over its expected lifespan. The initial alpha release of AEMS system is planned for March 2023, and the beta release is planned for the summer of 2023. The final release is planned for March 2024. Section 2 of the report documents the relevant building types that AEMS is suitable for. Section 3 documents EE features that will be supported. It will also include the data requirements, hardware requirements, implementation details, and how EE features will be tested and validated. Grid service features will be documented in section 4, including data requirements, hardware requirements, implementation details, and how the services will be tested and validated. Planned next steps are described in section 5.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Air Source Heat Pumps in Very Cold Climates

Air source heat pumps (ASHPs) in heating mode use a refrigeration cycle to remove heat from the outside air and transfer it into a building. Therefore, ASHPs have the potential to save energy compared to conventional heat sources that create the heat as opposed to transferring it from outside. Many cold climate ASHPs available on the market today can operate at outside temperatures of -25 degrees C (-13 degrees F) or even lower. As a result of technological advances and other factors, there is growing interest in ASHPs in cold climates, including very cold regions such as Alaska. However, guidance on the selection, installation, and operation of ASHPs in these very cold climates is limited, as significant data gaps exist regarding the performance of ASHPs in very cold environments. To address these data gaps and guide future innovations, our research team has studied the field performance of several ASHP installations in Alaska and done lab evaluations of several ASHP models using a cold chamber. While a cold chamber cannot fully reflect field conditions, it allowed for changing one variable at a time and gaining additional understanding of the behavior of ASHPs that would be difficult to gain from field studies only. This study focused on two main variables: the temperature in the chamber (representing the outdoor temperature) and the level of thermal loading of the ASHP. The results from the field as well as the lab show that ASHPs can operate with relatively high efficiency even in very cold climates if used in appropriate situations and in an appropriate way. It was found that not only the outside temperature, but also the level of thermal loading is a significant factor affecting the ASHP efficiency and needs to be carefully considered when sizing and operating ASHPs in very cold climates.

air source heat pump↗