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.
An interactive map showing Populus trichocarpa GWAS population structure estimated by ADMIXTURE (k=3, selected as optimal from k=2-11). Sampling locations are colored by their predominant ancestry proportion among the three inferred populations and geographic origins are searchable by genotype or river system using the search bar.
The increasing frequency of extreme weather events poses significant risks to power distribution systems, leading to widespread outages and severe economic and social consequences. This paper presents a novel simulation framework for assessing and enhancing the resilience of power distribution networks under such conditions. Resilience is estimated through Monte Carlo simulations, which simulate extreme weather scenarios and evaluate the impact on infrastructure fragility. Due to the proprietary nature of power network topology, a distribution network is synthesized using publicly available data. To generate the weather scenarios, an extreme weather generation method is developed. To enhance resilience, renewable resources such as solar panels and energy storage systems (batteries in this study) are incorporated. A customized Genetic Algorithm is proposed to determine the optimal locations and capacities for solar panels and battery installations, maximizing resilience while balancing cost constraints. Experiment results demonstrate that on a large-scale synthetic distribution network with more than 300,000 nodes and 300,000 edges, the proposed framework can efficiently evaluate the resilience, and enhance the resilience through the installations of distributed energy resources (DERs), providing utilities with valuable insights for community-level power system resilience estimation and enhancement.
The progression of electrification in the building and transportation sectors brings new opportunities for energy decarbonization. With higher dependence on the grid power supply, the variation of the grid carbon emission intensity can be utilized to reduce the carbon emissions from the two sectors. Existing coordinated control methods for buildings with distributed energy resources (DERs) either consider electricity price or renewable energy generation as the input signal, or adopt optimization in the decision-making, which is difficult to implement in the real-world environment. This paper aims to propose and validate an easy-to-deploy rule-based carbon responsive control framework that facilitates coordination between all-electric buildings and electric vehicles (EVs). The signals of the grid carbon emission intensity and the local photovoltaics (PV) generation are used for shifting the controllable loads. Extensive simulations were conducted using a model of an all-electric mixed-use community in a cold climate to validate the control performance with metrics such as emissions, energy consumption, peak demand, and EV end-of-day state-of-charge (SOC). Our study identifies that 4.5% to 27.1% of annual emission reduction can be achieved with limited impact on energy costs, peak demand, and thermal comfort. Additionally, up to 32.7% of EV emission reduction can be obtained if the EV owners reduce the target SOC by less than 21.2%.
Three-dimensional characterization of materials at the nano- and meso-scale has become possible with transmission and scanning transmission electron microscopes (S/TEM). Its importance has extended to a wide class of nanomaterials such as hydrogen fuel cells, solar cells, industrial catalysts, new battery materials and semiconductor devices, as well as spanning high-tech industry, universities, and national labs. While capable instrumentation is abundant, this rapidly expanding demand for high-resolution tomography is bottlenecked by software that is instead tailored for lower-dose, biological applications and not optimized for higher-resolution materials applications. Existing tomography fails to utilize the chemical information provided by S/TEM spectrometers. To address this problem, this project delivered a scalable, fully functional, freely-distributable, open Source materials tomography package with a modern user interface that enables automated acquisition, alignment, and real-time reconstruction of raw tomography data, and provides advanced segmentation, three-dimensional chemical visualization and analysis optimized for materials applications. It has established an extendable framework capable of automation for high-throughput for the tomography of materials from data acquisition to visualization. Phase I and II delivered a full-featured, cross-platform, clean and integrated application for S/TEM materials tomography. It can read projection data from the microscope, with graphical tools for alignment of data, tomographic reconstruction, segmentation, and visualization of the reconstructed 3D volume. The entire pipeline can be saved to an XML state file, enabling fully reproducible data processing and analysis with a full Python environment for the development of custom algorithms, processing, and analysis requirements. Phase IIB extended the application to provide real-time tomographic reconstruction, "live updates" as data is processed, analysis of big data, and multi-channel tomography. Here, quantitative assessment of nanomaterials can occur as data is being recorded on a 3D visualization platform that accommodates additional information from multiple channels. We provide unmatched high-throughput tomography, where the complete tomographic pipeline from measurement to 3D visualization can occur rapidly. The project was extended to offer capabilities for micro-CT, X-ray, neutron, focused ion beam, atomic electron tomography and atom probe tomography. With over 600 transmission electron microscopes worldwide and approximately 50 coming online each year, the demand and impact of an open-source tomography tool is large. Significant opportunities exist in high-tech industry, universities, and national labs to enable or enhance three-dimensional imaging at the nanoscale and bring automated high-throughput approaches that will accelerate progress in materials characterization and metrology. The project supports a service based business model by enabling lab-specific acquisition and processing customization and integration-support and development that will be provided into Phase III and beyond.
The use of Gaussian processes (GPs) as models for astronomical time series data sets has recently become almost ubiquitous, given their ease of use and flexibility. In particular, GPs excel at marginalization over the stellar signal when the variability due to starspots is treated as a nuisance, as in exoplanet transit modeling. However, these effective models are less useful in cases where the starspot signal is of primary interest, since it is not obvious how the parameters of the GP relate to physical parameters like the spot size, contrast, and latitudinal distribution. Instead, it is common practice to explicitly model the effect of individual starspots on the light curve and attempt to infer their properties via optimization or posterior inference. Unfortunately, this process is ill-posed and often computationally intractable when applied to stars with more than a few spots and/or to ensembles of many stars. Here we derive a closed-form expression for a GP that describes the light curve of a rotating, evolving stellar surface conditioned on a given distribution of starspot sizes, contrasts, and latitudes. We demonstrate that this model is correctly calibrated, allowing one to robustly infer physical parameters of interest from one or more light curves, including the typical spot radii and latitudes. Our GP has far-ranging implications for understanding the variability and magnetic activity of stars from light curves and radial velocity measurements, as well as for modeling correlated noise in exoplanet searches. Our implementation is efficient, user-friendly, and open-source, available in the package starry-process.
This is a parallel variable population multi-objective optimizer with an adaptive unified differential evolution algorithm or a genetic algorithm. It can also be used for single objective optimization. Some features of this code include: 1) The population size varies from generation to generation to save the total # of objective function evaluations. 2) The population is uniformly distributed to a number of parallel processors for simultaneous objective function evaluation. 3) The objective function evaluation can be attained from an external simulation program with control variables in its input file and objectives calculated from its output files. 4) The optimizer includes an adaptive unified differential evolution algorithm and a real value genetic algorithm. The parameters in the unified differential evolution algorithm can be chosen to attain any mutation schemes in the published literature.
This paper proposes a novel Soft-Actor-Critic (SAC) based Deep Reinforcement Learning (DRL) method for optimizing the cost of microgrid operation by leveraging load flexibility. The proposed SAC-DRL method is designed to coordinate the control of distributed energy resources (DERs) and flexible load, addressing practical energy billing formation by power distribution utilities. Key contributions include an innovative reward function to mitigate sparse reward challenges and a mixed control strategy for discrete and continuous variables, ensuring radial network topology and minimizing power loss. We evaluate the proposed method on the model of a real microgrid located in Southern California, U.S.. The SAC-DRL model is tested to demonstrate its efficacy in reducing grid dependence, optimizing resource use, and minimizing costs. The results highlight the potential of DRL in modern energy systems, offering a sustainable and economically efficient solution for energy management in microgrids.
Utility service providers are often challenged with the synchronization of thermostatically controlled loads. Load synchronization, as a result of naturally occurring and demand-response events, has the potential to damage power distribution equipment. Because thermostatically controlled loads constitute most of the power consumed by the grid at any given time, the proper control of such devices can lead to significant energy savings and improved grid stability. The contribution of this paper is the development of an optimal control algorithm for commonly used variable speed heat pumps. By means of selective peer-to-peer communication, our control architecture allows for the regulation of home temperatures while simultaneously minimizing aggregate power consumption, and aggregate load volatility. An optimal centralized controller is also explored and compared against its decentralized counterpart.
Residual stresses are often overlooked in friction stir processing (FSP), but their significant impact on fatigue performance necessitates their consideration in optimizing processing parameters. The first step in this effort is understanding how process conditions influence residual stress distributions, especially across different alloys. This study focuses on determining and explaining the through-thickness residual stress variations and the effect of process temperature on the residual stress magnitude in wrought AA7075 and cast AA380.0 alloys. Additionally, for AA380.0, the impact of a second FSP pass was investigated. To achieve this, hole-drilling electronic speckle pattern interferometry (ESPI) and the thermal pseudo-mechanical (TPM) model within finite element analysis were employed to study the 3D distributions of in-plane residual stresses in processed samples under various conditions. A key finding was the varying impact of process temperatures on residual stress magnitudes. Higher process temperatures reduced stresses in AA380.0 but increased them in AA7075. Additionally, the through-thickness stress distributions differed between the two alloys. Further analysis revealed that yield stresses are crucial in explaining these phenomena and the effects of additional FSP passes. Further, this fundamental understanding will be vital in guiding the efforts to mitigate residual stresses and assess their impact on the performance of FSP aluminum alloys.
Abstract The method for direct injection of fuel in the cylinder of an IC engines is important to high-efficiency and low-emission performance. Optical spray diagnostics plays an important role in understanding plume movement and interaction for multi-hole injectors, and providing baseline understanding used for computational optimization of fuel delivery. Traditional planar or line-of-sight diagnostics fail to capture the liquid distribution because of optical thickness concerns. This work proposes a high-speed (67 kHz) extinction imaging technique at various injector rotations coupled to computed tomography (CT) for time-resolved reconstruction of liquid volume fraction in three dimensions. The number of views selected and processing were based on synthetic (modeled) liquid volume fraction data where extinction and CT adequately reconstructed each plume. The exercise showed that for an 8-hole, symmetric-design injector (ECN Spray G), only three different views are enough to reproduce the direction of each plume, and particularly the mean plume direction. Therefore, the number of views was minimized for experiments to save expense. Measurements applying this limited-view technique confirm plume–plume variations also detected with mechanical patternation, while providing better spatial and temporal resolution than achieved previously. Uncertainties due to the limited view within pressurized spray chambers, the droplet size, and optically thick regions are discussed. Graphic abstract
The microbe, Zymomonas mobilis, can efficiently convert sugars to 2,3 butanediol (BDO), an important intermediate for downstream chemical products, only in a well-controlled microaerated environment. However, controlling oxygen distribution in industrial-scale bioreactors is challenging, and at-scale BDO production is hard to optimize using traditional engineering methods. This study takes a step towards addressing this problem through a computational model for reacting multiphase flows in large-scale bioreactors. A phenomenological metabolic model was first developed and validated against experiments, and then it was coupled to a multiphase computational fluid dynamics (CFD) solver using a subcycling algorithm for simulating long conversion times (~30 hours). Large-scale (~500 m 3 ) bubble column simulations using this coupled approach demonstrated a 25% improvement in BDO yield with low sparging rates (0.002 m/s) and low height-to-diameter ratio (0.875) compared to a baseline bubble column (35 m height, 5 m diameter) with a superficial gas velocity of 0.06 m/s.
We introduce NeuroSA, a neuromorphic architecture specifically designed to ensure asymptotic convergence to the ground state of an Ising problem using a Fowler-Nordheim quantum mechanical tunneling based threshold-annealing process. The core component of NeuroSA consists of a pair of asynchronous ON-OFF neurons, which effectively map classical simulated annealing dynamics onto a network of integrate-and-fire neurons. The threshold of each ON-OFF neuron pair is adaptively adjusted by an FN annealer and the resulting spiking dynamics replicates the optimal escape mechanism and convergence of SA, particularly at low-temperatures. To validate the effectiveness of our neuromorphic Ising machine, we systematically solved benchmark combinatorial optimization problems such as MAX-CUT and Max Independent Set. Across multiple runs, NeuroSA consistently generates distribution of solutions that are concentrated around the state-of-the-art results (within 99%) or surpass the current state-of-the-art solutions for Max Independent Set benchmarks. Furthermore, NeuroSA is able to achieve these superior distributions without any graph-specific hyperparameter tuning. For practical illustration, we present results from an implementation of NeuroSA on the SpiNNaker2 platform, highlighting the feasibility of mapping our proposed architecture onto a standard neuromorphic accelerator platform.
Coin cells are used extensively as test devices in battery research for evaluation of new materials and optimization of cycling protocols. Here, in situ X-ray diffraction profilometry is used to characterize spatial distribution of the active materials, lithiation, and phase distribution in electrodes of NCM523/graphite coin cells. The X-ray data indicate uneven areal compression of the electrode assembly in such cells, which we trace to a specific design feature that leads to elastic deformation of a metal spacer. Steep lithiation gradients observed in the electrodes imply radially-dependent resistivity, for which uneven compression of the separator is a likely cause. Electrochemical model calculations suggest that variable porosity of the polymer separator would account for the salient features of spatial profiles observed in these coin cells.
Significant developments made in quantum hardware and error correction recently have been driving quantum computing toward practical utility. However, gaps remain between abstract quantum algorithmic development and practical applications in computational sciences. In this perspective article, we propose several properties that scalable quantum computational science methods should possess. We further discuss how block-encodings and polynomial transformations can potentially serve as a unified framework with the desired properties. Recent advancements on these topics are presented, including the construction and assembly of block-encodings, and various generalizations of quantum signal processing (QSP) algorithms to perform polynomial transformations. The scalability of QSP methods on parallel and distributed quantum architectures is also highlighted. Promising applications in simulation and observable estimation in chemistry, physics, and optimization problems are presented. We hope this perspective serves as a gentle introduction to state-of-the-art quantum algorithms for the computational science community and inspires future development of scalable quantum computational science methodologies that bridge theory and practice.
In hybrid quantum systems, the interface between quantum and classical domains is essential for the generation, control, and measurement of quantum states. Quantum-classical interfaces (QCIs) are ubiquitous in devices such as optical modulators, quantum sensors, and signal processors, where classical signals influence quantum dynamics. In this paper, we employ the Lindblad master equation to simulate the evolution of a quantum system interacting with a classical control system. Our model captures both linear and nonlinear interactions by incorporating first- and second-order susceptibilities, and it quantifies the influence of externally applied control parameters on decoherence and state evolution. As an illustrative example, we analyze an optical modulator and demonstrate how variations in material response and drive conditions affect photon statistics, coherence, and phase-space distributions. In conclusion, the findings offer a path to an all-encompassing model for understanding and optimizing QCIs, with wide-ranging implications for the performance, design, and robustness of next-generation quantum devices.
Compact heat exchangers using supercritical fluids such as CO2 are preferred due to their high heat transfer capacity and smaller footprint. Three-dimensional conjugate forced convection heat transfer analysis was performed on several shell-and-tube counter-flow microchannel heat exchangers. Numerical simulations were conducted to test effect of change in mass flow rate, hydraulic diameter and various cross sections on the heat transfer. Increasing mass flow rate improved heat transfer up to a maximum value and then decreased downstream with increasing turbulence. Maximum heat transfer was obtained for the micro channel with the smallest hydraulic diameter. Amongst the cross sections analyzed (circular, square, circular with radial ribs, and square with radial ribs), the most uniform distribution of temperature and maximum heat transfer were obtained for circular cross section with radial ribs. An optimally efficient operation of such a heat exchanger can be attained by considering these factors during multi-objective constrained optimization of geometric parameters and requirements for additive manufacturing of such compact heat exchangers.
Photon Doppler velocimetry (PDV) is a ubiquitous diagnostic method in dynamic compression experiments. It typically requires the roughening of the test surface to ensure scattering of the probe light back into the PDV probe at non-normal angles of incidence. Although surface scattering plays an essential role in the quality of PDV data, surface scattering properties for common materials and preparations have not previously been studied with specific relevance to PDV. Here, we present detailed measurements of the bidirectional reflectance distribution function for four different surface preparations (Scotch-Brite hand and drill roughened, ball-rolled, and grit blasted) and four different materials (copper, aluminum, stainless steel, and tantalum). These measurements employ a conventional PDV probe and obtain a diffraction limited angular resolution with 10 pW accuracy and 2.6% repeatability. In addition to scattering data, we employ scattering theory and simulations to accurately emulate the measured data. We also present a straightforward method to derive the average scattering distribution from surface profilometry and observe several qualitative aspects of the scattering data that may help to optimize PDV signals.
Modern data scientists are tasked to analyze ever-growing data sets with increasingly complex relationships. Tensor decompositions have come to play a central role in identifying underlying latent structures in higher-order data. The problem of fitting tensor models to different distributions is complicated by the combinations of size, dimensionality, and sparsity present in real world data. The situation demands efficient algorithms designed for shared-memory and distributed systems. This work will present new research that tackles these challenges on several different fronts, leveraging optimizations in numerical algorithms and sparse tensor representations in heterogeneous high performance computing environments.