SEARCH · Engineering Papers
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.
CALAND is a system of algorithms (.r files) and data developed to quantify the impacts of various suites of California State-supported land use and land management strategies on landscape carbon and greenhouse gas emissions (CO2, CH4, and optional black carbon) relative to a baseline scenario for the Draft California 2030 Natural and Working Lands Climate Change Implementation Plan (2019).
Resiliency is one of the key challenges in today's power system. Natural disasters and cyber-attacks both can limit communications between microgrids and the central management system. Thus, having a Distributed Microgrid Coordination (DMC) algorithm can improve the system resiliency, which enables the microgrids to operate without communication with the central management system. Peer-to-Peer Energy Management System for Distributed Microgrid Coordination adopts a primal-dual approach, where each microgrid controller keeps a local estimate of the dual variables. The estimate is updated with local measurements and peer-to-peer communication, leading to a fully distributed algorithm. While the DMC is developed for microgrid coordination, it can be used for general distributed control purpose.
BiG-SLiCE was originally an open source Python-based command line bioinformatics software that offers a highly scalable clustering analysis on biosynthetic gene clusters (BGC) data. It allows a simultaneous analysis of millions of BGCs, exceeding the capability of other existing tools (around one hundred thousands). As a tradeoff, the clustering accuracy is relatively lower and sometimes fall short in corner cases and specific BGC classes such as the RiPPs (Ribosomally-translated, Post-translationally modified Peptides). In BiG-SLiCE V2 (developed in LBNL), the clustering algorithm has been significantly improved to deliver a much accurate result even for RiPPs and other previous corner case classes. Moreover, the speed of the overall pipeline has been improved by 50-100%. Finally, additional features were implemented to support downstream analyses of BiG-SLiCE results, such as customized tabular (TSV/CSV) and columnar (Parquet) outputs.
This HADREC GUI is designed to visualize power system emergency control time-series data and facilitate the comparison of the newly developed meta strategy optimization algorithm-based emergency control with existing out-of-step-based emergency control approach, and a no-action case for reference.
Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.
The multivariate regression interpretation of the Gaussian chain graph model simultaneously parametrizes (i) the direct effects of p predictors on q outcomes and (ii) the residual partial covariances between pairs of outcomes. We introduce a new method for fitting sparse versions of these models with spike-and-slab LASSO (SSL) priors. We develop an Expectation Conditional Maximization algorithm to obtain sparse estimates of the p × q matrix of direct effects and the q × q residual precision matrix. Our algorithm iteratively solves a sequence of penalized maximum likelihood problems with self-adaptive penalties that gradually filter out negligible regression coefficients and partial covariances. Because it adaptively penalizes individual model parameters, our method is seen to outperform fixed-penalty competitors on simulated data. We establish the posterior contraction rate for our model, buttressing our method’s excellent empirical performance with strong theoretical guarantees. Using our method, we estimated the direct effects of diet and residence type on the composition of the gut microbiome of elderly adults.
As part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC), the HELiX uncrewed aircraft system (UAS) was deployed over the sea ice in the central Arctic Ocean during summer 2020. Albedo measurements were obtained with stabilized pyranometers, and melt pond fraction was calculated from orthomosaic imagery from a surface-imaging multispectral camera. This study analyzed HELiX flight data to provide insights on the temporal and spatial evolution of albedo and melt pond fraction of the MOSAiC floe during the melt season as it drifted south through Fram Strait. The surface albedo distributions showed peak values changing from high albedo (0.55–0.6) to lower values (0.3) as the season advanced. Inspired by methods developed for satellite data, an algorithm was established to retrieve melt pond fraction from the orthomosaic images. We demonstrate that the near-surface observations of melt pond fraction were highly dependent on sample area, offering insight into the influence of subgrid scale features and spatial heterogeneity in satellite observations. Vertical observations conducted with the HELiX were used to quantify the influence of melt pond scales on observed surface albedo as a function of sensor footprint. These scaling results were used to link surface-based measurements collected during MOSAiC to broader-scale satellite data to investigate the influence of surface features on observed albedo. Albedo values blend underlying features within the sensor footprint, as determined by the melt pond size and concentration. This study framed the downscaling (upscaling) problem related to the airborne (surface) observations of surface albedo across a variety of spatial scales.
We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.
Heating hot water distribution systems are typically used in commercial buildings to condition spaces to provide occupant thermal comfort. However, recent research shows significant distribution losses within these systems that drive down the overall hot water plant efficiency. This research focuses on detecting passing valves in reheat coils found in variable air volume (VAV) terminal units to reduce distribution losses. A passing valve allows hot water flow when the actuator on the valve is commanded to be closed. The fluid causes unintentional heating or cooling to occur, causing comfort and control issues, and wasting energy. We developed the passing valve detection algorithm using a framework based on the Brick schema and Mortar platform to ensure that the application is portable and can scale to many buildings. We applied the same application to analyze 1,335 VAV reheat terminal units in 20 buildings. The diversity found in these large datasets increases confidence that any building with VAV reheat terminal units with the required sensors and Brick data model can run our open-source algorithm with little or no modification. In aggregate, 5% of VAV units analyzed were categorized as having a sensor fault, 14% with potential passing valve fault, and 81% with no faults detected. However, there is a significant variation in the proportion of VAV units with a passing valve detected (1% to 83%) of each building’s analyzed units.
The performance of the open-source thermochemical software library Thermochimica has been enhanced by developing a re-initialization algorithm to make use of data from previous calculations to reduce the number of convergence steps in the Gibbs energy minimization procedure. This algorithm has been tested in the context of stand-alone Thermochimica calculations, and speedups in the range of 2x-3x achieved for cases with chemistries resembling those of irradiated nuclear fuels. Routines to make use of this re-initialization procedure have been implemented in Bison, in which each node uses the results of the previous calculation at that node as the initial conditions for the following Thermochimica call. Examples based on calculating the diffusion of oxygen in UO$_2$ LWR fuel in 1D and 3D have been tested, with speedups up to 7.38x demonstrated.
Across multi-tenant commercial office and multifamily buildings, centrally metered fuel use represents a substantial fraction of whole-building energy use. Energy audit practitioners understand that improving heating distribution efficiency is typically more of an opportunity than combustion efficiency and that differing thermal comfort preferences between tenants are the bane of operators across these building typologies. There is an unmet market need for retrofit technologies that allow for the delivery of the right amount of heat to the right spaces, at the right time. The Energy Management and Information System (EMIS) package fills this gap through enhanced controls and metering, incorporating low-cost sensors and wireless communication infrastructure to provide a platform for ongoing commissioning and tenant feedback, including heat cost allocation. With support from the US DOE Building Technologies Office, Steven Winter Associates, Inc. (SWA) partnered with Sentient Buildings, E Source, building owners, and utility and policy stakeholders, to demonstrate a market viable EMIS that achieves a reduction in space heating energy use by reducing heating load, improving control, and positively impacting behavior while providing an acceptable financial return. In this study, EMIS packages were implemented in two New York City multifamily rental buildings. Both buildings conducted basic mechanical work (e.g., repairing steam traps) to ensure the heating system was operating well before any tenant feedback was layered in. Heating Energy Use Reports (HEUR) were created to provide tenants with social comparisons and energy savings tips to influence their behavior; these were provided monthly to all tenants in both buildings. Additionally, one building allocated heating costs to a portion of the tenants. Heat cost allocation (HCA) has a long history in the European Union (EU), although it is not common in the US or in steam-heated buildings. SWA leveraged existing EU best practices and stakeholder feedback to develop a Heat Cost Allocation algorithm that was considered equitable and intuitive. Energy use and tenant behavior impacts were tracked throughout the study. The basic mechanical repair work saved between 11-20% of heating energy. Those savings rose to 17-24% with the addition of tenant feedback. While it may not be possible to precisely determine the impact of COVID-19 on research studies like this, there may have been additional savings realized had the study taken place in a period of normal occupancy patterns. These types of central heating systems have been a blind spot for utilities, who have traditionally had little visibility into detailed behind-the-meter gas usage. Heating energy savings stayed consistent during the coldest months, indicating the potential for utilities to utilize EMIS packages for peak gas demand reductions or demand response programs. Tenant comfort was also improved. Post installation, room temperatures more closely matched thermostat set points. Perhaps due to this greater level of control, the vast majority of tenants being billed for heating were accepting of the allocation costs. And tenants receiving heat cost allocations were more likely to reduce their thermostat setpoints than tenants receiving behavioral feedback without financial impacts were. Variation in building specifics makes it difficult to provide precise energy and financial savings estimates. But within the range of expected conditions, the study identified a few key variables that can have the greatest impact on financial returns: the cost of fuel, the ability and willingness to allocate heating costs to tenants, and a well-functioning heating system as a starting point. This study focused on two multifamily buildings, but additional use cases, such as commercial buildings and affordable housing, should be explored to better understand the full market potential. While this type of upgrade has the potential for deep energy reductions and cost savings, future projects should take into account the balance of costs and benefits between owners and tenants, especially in the affordable, regulated, or other low-to-moderate income (LMI) segments of the market. Rent credits, utility allowances, or a shared savings program are possible options to accelerate adoption of this strategy in these market segments.
This document describes Offline Software and Computing for the Deep Underground Neutrino Experiment (DUNE) experiment, in particular, the conceptual design of the offline computing needed to accomplish its physics goals. Our emphasis in this document is the development of the computing infrastructure needed to acquire, catalog, reconstruct, simulate and analyze the data from the DUNE experiment and its prototypes. In this effort, we concentrate on developing the tools and systems that facilitate the development and deployment of advanced algorithms. Rather than prescribing particular algorithms, our goal is to provide resources that are flexible and accessible enough to support creative software solutions as HEP computing evolves and to provide computing that achieves the physics goals of the DUNE experiment.
The space environment is fraught with complex plasmas spanning a wide range of densities and temperatures. Much of space plasma research has focused on large-scale changes in the ambient plasma and fields from the Sun to the Earth’s atmosphere, driven primarily by the solar wind. Yet one particular type of space plasma, which forms from a hypervelocity impact (HVI), remains poorly understood. Hypervelocity impactors include both meteoroids and space debris. Meteoroids are naturally occurring objects in space that travel between 11 and 72.8 km/s and originate primarily from comets and asteroids. In contrast, space debris are human-made objects with speeds typically < 11 km/s. Hypervelocity impactors routinely hit spacecraft, yet the physics behind the formation of the plasma and the dynamics of its expansion remain largely unknown. The complexity of this phenomenon necessitates a research approach that includes both experimental studies and numerical simulation in order to understand the underlying physical processes that occur upon formation and expansion of the impact plasma. Our research has focused on providing a comprehensive understanding of plasma generated by hypervelocity impacts by meteoroids and space debris on spacecraft in order to characterize the behavior of the expanding plasma and its interactions with the ambient environment. Previously. we conducted experimental campaigns at a dust accelerator facility that can accelerate particles up to 60 km/s, which is representative of meteoroid speeds, and at a light gas gun facility that can accelerate larger projectiles up to 7 km/s, which is representative of orbital debris. The experiments included plasma, optical and radio frequency (RF) sensors in order to understand the dynamics and associated RF emission resulting from hypervelocity plasmas. For this research, we developed and applied machine learning algorithms to identify which type of impactor would produce RF and physics-based models to determine the source of the RF. This was a one-year research program that resulted in 2 refereed journal publications (uploaded separately).
Molecular- and nano- junctions, surfaces, interfaces, and interfacial processes have been the subjects of this DOE supported program for about 20 years. We have investigated electronic and magnetic structure, spin-dependent charge transport, and emerging phenomena at surfaces and interfaces in a number of physical systems ranging from nano-clusters, to nanowires, to two-dimensional (2D) crystals and their interactions, to bulk matter. Eighty-four refereed papers were published [1-84], with a majority effort from the group led by the PI, and over a hundred presentations were made at international conferences. We use first-principles methods in the framework of density functional theory in conjunction with nano-equilibrium Green function techniques, semi-empirical Boltzmann transport theory, tight-binding models, and the Hubbard model and many-body wavefunction approaches. In pursuing our scientific goals, we have also developed a number of computational algorithms and computational packages.
This project is focused on performing the best science possible on current generation quantum computers. As such, we have focused on developing techniques that involve robust algorithms that are either insensitive to noise, or are extremely low depth. We also have worked on determining topological properties and employing the robustness of topology to noise to allow for robust quantum computation. The ultimate goal is to be able to do new science on quantum computers that cannot be achieved on classical computers. The main effort of the work we have completed is on time-evolution of quantum systems and robust calculation of properties of systems enabled by being able to perform time evolution.
This project explored the role of tropical oceans in regulating cloud populations within the Madden-Julian oscillation (MJO)–a large-scale tropical weather disturbance that influences extreme weather events across the globe by affecting the position of the jet stream and patterns of moisture transport. Our project sought to 1) document ocean feedback strength to the MJO in climate models, 2) understand how air-sea interactions affect the intensity, organization, and propagation of MJO convection, 3) utilize existing observations to identify sources of model bias for tropical air-sea interactions, and 4) better understand the effects of slow and fast ocean variability on MJO cloudiness. We found that transfers of heat and moisture across the air-sea interface (i.e., surface fluxes) in most climate models are over-supportive in maintaining MJO convection and under-supportive in propagating MJO convection. This can lead to misrepresentations in the frequency of extreme weather events across the globe through erroneous MJO feedbacks to jet stream position and moisture transport. Furthermore, using observations collected by ocean moorings across the tropical oceans, we showed that most models greatly overestimate surface evaporation in the subtropics, and that this bias is likely one contributor to a pervasive bias among climate models in their ability to simulate the mean state of rainfall over the tropical oceans (i.e., the “double ITCZ bias”). Replacing the default surface flux algorithm in two climate models with a state-of-the-art surface flux algorithm modestly reduced the double ITCZ bias in those models. Our work produced two model diagnostics packages that are now publicly available through DOE and NOAA model diagnostics archives. The first, the Analysis of Scales of Precipitation (ASoP) diagnoses scales of precipitation variability across space and time, with temporal resolution as fine as the model integration timestep. The second, the air-sea flux diagnostic diagnoses model surface fluxes and allows direct comparison of model fluxes to observations without having to recompute fluxes using model output and a state-of-the-art flux algorithm. These diagnostics, which were developed using climate models, can be readily applied to regional and global forecast model output to improve forecasts for the benefit of society.
The purpose of the Field Validation of Cloud Properties Sensor – SAIL campaign was to deploy Aerodyne’s extended-wavelength TWST cloud optical properties sensor (TWST-CPS) in an operationally relevant environment with co-located, externally validated sensors. Colocation with The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s deployment of its second Mobile Facility (AMF2) on the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign at Crested Butte, Colorado was ideal because of the variety of operational sensors that can measure some of the same cloud properties using different modalities. While cloud property sensors have long existed, they tend to be costly to produce and maintain. Our sensor measures absolute spectral radiance over 440-1700 nm and retrieves cloud optical depth (COD), droplet effective radius (Reff), and thermodynamic phase (Phase). Our prototype is built predominantly from off-the-shelf components and uses uncooled spectrometers. A low-cost, easy-to-use sensor such as this could allow deployment at many more sites for greater spatial coverage. Instrument deployment to test robust operations and a new stray-light baffle design, and the continued development of analysis and retrieval algorithms using the data from this deployment, are central technical objectives of Aerodyne Research’s U.S. Department of Energy Small Business Innovative Research Phase 2 contract (DE-SC0020473, Low-cost shortwave spectroradiometer for retrieval of cloud properties).
The purpose of the Field Validation of Cloud Properties Sensor – SAIL campaign was to deploy Aerodyne’s extended-wavelength TWST cloud optical properties sensor (TWST-CPS) in an operationally relevant environment with co-located, externally validated sensors. Colocation with The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s deployment of its second Mobile Facility (AMF2) on the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign at Crested Butte, Colorado was ideal because of the variety of operational sensors that can measure some of the same cloud properties using different modalities. While cloud property sensors have long existed, they tend to be costly to produce and maintain. Our sensor measures absolute spectral radiance over 440-1700 nm and retrieves cloud optical depth (COD), droplet effective radius (Reff), and thermodynamic phase (Phase). Our prototype is built predominantly from off-the-shelf components and uses uncooled spectrometers. A low-cost, easy-to-use sensor such as this could allow deployment at many more sites for greater spatial coverage. Instrument deployment to test robust operations and a new stray-light baffle design, and the continued development of analysis and retrieval algorithms using the data from this deployment, are central technical objectives of Aerodyne Research’s U.S. Department of Energy Small Business Innovative Research Phase 2 contract (DE-SC0020473, Low-cost shortwave spectroradiometer for retrieval of cloud properties).