Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed method”

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 577 records · Page 32

Optimal PV Inverter Control in Distribution Systems via Data-Driven Distributionally Robust Optimization

Distribution systems with high penetration of uncertain solar generation call for advanced control strategies of photovoltaics (PVs) inverters. This paper proposes a data-driven distributionally robust optimization (DDDRO) approach to optimally controlling the PV inverters to improve the system operation performance under solar power uncertainties. In the proposed DDDRO approach, a Wasserstein ball-based method is proposed to construct the distributional ambiguity set to model the uncertainties of PV generation through partial observations of historical data without knowing exact probability distributions. We further reformulate the computationally intractable DDDRO model to a mixed integer second order cone programming (MISOCP) problem. The effectiveness and out-of-sample performance of the proposed approach have been demonstrated on a modified IEEE 33-node system. We conduct a comparative study to compare the proposed method with traditional chance constrained programming (CCP). It shows that the proposed DDDRO approach can provide a less conservative yet robust solution to minimize the worse-case expectation of the total network loss while maintaining nodal voltages in a secure range.

Xue, Yaosuo↗

Practical considerations for measuring global spin density matrix elements of vector mesons in heavy-ion collisions

The STAR Collaboration has reported a significant 𝜙-meson global spin alignment (𝜌 00 ) signal in Au+Au collisions at $\sqrt{s_{NN}}$ ≤ 62 GeV by measuring the polar angle distribution of 𝜙-meson daughters with respect to the orbital angular momentum direction of the collision system. Here, in this paper, a new method is explored for studying vector-meson global spin alignment in heavy-ion collisions by examining the two-dimensional (2D) polar and azimuthal angle distribution. This method allows simultaneous extraction of 𝜌 00 and off-diagonal spin density matrix elements (SDMEs), providing unique access to local quark-antiquark spin correlations and spin hydrodynamics in quark-gluon plasma. The new 2D method also removes potential biases from nonzero off-diagonal SDMEs on 𝜌 00 with the 1D method. A detailed procedure to correct for detector acceptance and resolution effects is also presented and validated by simulation studies.

Wilks, Gavin [University of Illinois, Chicago, IL ↗

EGS Stimulation Design with Uncertainty Quantification at the EGS Collab Site

Engineering a robust hydraulic connection between wells is one of the most difficult aspects of enhanced geothermal systems (EGS). Designing and constructing such hydraulic connections requires and understanding of the in situ state of stress and the heterogeneities and discontinuities that naturally exist and may control the stimulation. Even with comprehensive stress and formation characterization programs substantial uncertainty remains in these key parameters. This is especially the case in high-temperature EGS environments where drilling conditions are often difficult and a far fewer logging and testing options are available. This paper presents a new approach for explicitly quantifying the uncertainties in the state of stress using a Bayesian Markov Chain Monte Carlo method. This approach produces a probability distribution for the stress tensor, including a general 3D orientation, that reflects the uncertainties in all the observations or indicators used to constrain the stress state. This method is demonstrated on the characterization data for the EGS Collab Experiment 2 site. The output of the analysis is used to guide the design of the planned stimulations. In the case of research projects like EGS Collab, explicitly quantifying the uncertainties in the stress state allow for more rigorous hypothesis testing by allowing conclusions drawn from the experiments to be interpreted in the context of the uncertain knowledge about conditions in the test bed.

Burghardt, Jeffrey A.↗

Using Median Point in Keeling Plot to Reduce the Uncertainty of the Isotopic Composition of Evapotranspiration

Abstract The isotopic composition of evapotranspirationδ ET is a crucial parameter in isotope-based evapotranspiration (ET) partitioning and moisture recycling studies. The Keeling plot method is the most prevalent method to calculateδ ET , though it contains large extrapolated uncertainties from the least squares regression. Traditional Keeling regression uses the mean point of individual measurements. Here, a modified Keeling plot framework was proposed using the median point of individual measurements. We tested theδ ET uncertainty using the mean point [σ ET (mean)] and median point [σ ET (median)]. Multiple resolutions of input and output data from six independent sites were used to test the performance of the two methods. Theσ ET (mean) would be greater thanσ ET (median) when the mean value of inverse vapor concentration () is greater than the median value of inverse vapor concentration []. When applying the filter ofr 2 > 0.8, around 70% ofσ ET (mean) was greater thanσ ET (median). This phenomenon might be due to the normality of the vapor concentrationC υ producing the asymmetric distribution of 1/C υ . The median method could perform significantly better than the mean method when inputting high-resolution measurements (e.g., 1 Hz) and when the water vapor concentrationC υ is relatively low. Compared to the mean method, applying the median method could on average reduce 6.88% of ET partitioning uncertainties and could on average reduce 9.00% of moisture recycling uncertainties. This study provided a new insight of the Keeling plot method and emphasized handling model output uncertainty from multiple perspectives instead of only from input parameters.

Meteorology & Atmospheric Sciences↗

How “hot” are hotspots: Statistically localizing the high-activity areas on soil and rhizosphere images

The topic of microbial hotspots in soil requires not only visualizing their spatial distribution and biochemical analyses, but also statistical approaches to identify these hotspots and separate them from the surrounding activities (background). We hypothesized that each hotspot type (e.g. enzyme activities in the rhizosphere, root exudation, localization of herbicide accumulation) is a result of local process driven by biotic and/or abiotic factors, and the process rates in the hotspots are much faster than those in the soil background. We further hypothesized that the background and hotspot activities in soil belong to different statistical distributions. Consequently, hotspot determination should be based on statistical separation of activities significantly higher than the background. We analyzed for the statistical distributions of grey values on three groups of published images: 1) 14 C images of carbon input by roots into the rhizosphere, 2) 14 C glyphosate accumulation in the plant, and 3) zymogram of leucine aminopeptidase activity in rooted soil. The two Gaussian distributions were fit (the first representing the background, the second the hotspots) to the distribution of grey values in the images, the parameters (means and standard deviations, SD) of the fitted distributions were calculated, and the background was removed. Thus, we identified hotspots as areas outside of the Mean+2SD image intensity (corresponding to the upper ~ 2.5% of activity, being over 97.5% of background values) and finally, visualized images of solely hotspot locations. Finally, these results were compared with previously used decisions on hotspot intensity thresholding (i.e. Top-25% and 17 standard thresholding approaches in ImageJ) and discussed the advantages of the Mean+2SD as well as Mean+3SD approaches. These advantages include: i) simple unification of the thresholding approach for several imaging methods with various principles of activity distribution, ii) identification of hotspots with various activity levels, iii) analysis of “time-specific” hotspots in temporal sequences of images. Compared with 17 standard thresholding methods, we concluded that objectively elucidating and separating the hotspots should be based on statistical distribution analysis, e.g. using the Mean+2SD or Mean+3SD approaches. Furthermore, this simple Mean+2SD approach delivered suitable results for three groups of images and so, helps to understand the processes responsible for the highest activities and elucidate hotspots.

59 BASIC BIOLOGICAL SCIENCES↗

Sample preparation with sucrose cryoprotection dramatically alters Zn distribution in the rodent hippocampus, as revealed by elemental mapping

Transition metal ions (Fe, Mn, Cu, Zn) are essential for healthy brain function, but altered concentration, distribution, or chemical form of the metal ions has been implicated in numerous brain pathologies. Currently, it is not possible to image the cellular or sub-cellular distribution of metal ions in vivo and therefore, studying brain-metal homeostasis largely relies on ex vivo in situ elemental mapping. Sample preparation methods that accurately preserve the in vivo elemental distribution are essential if one wishes to translate the knowledge of elemental distributions measured ex vivo toward increased understanding of chemical and physiological pathways of brain disease. Furthermore, the choice of sample preparation is particularly important for metal ions that exist in a labile or mobile form, for which the in vivo distribution could be easily distorted by inappropriate sample preparation. One of the most widely studied brain structures, the hippocampus, contains a large pool of labile and mobile Zn. Herein, we describe how sucrose cryoprotection, the gold standard method of preparing tissues for immuno-histochemistry or immuno-fluorescence, which is also often used as a sample preparation method for elemental mapping studies, drastically alters hippocampal Zn distribution. Based on the results of this study, in combination with a comparison against the strong body of published literature that has used either rapid plunge freezing of brain tissue, or sucrose cryo-protection, we strongly urge investigators in the future to cease using sucrose cryoprotection as a method of sample preparation for elemental mapping, especially if Zn is an analyte of interest.

59 BASIC BIOLOGICAL SCIENCES↗

Developing ML/AI Methods for High-Throughput Characterization of Multiple-Sensor Streams of Tokamak Dynamics for High-Speed Control (Final Report)

This project evaluated and developed new mathematical and algorithmic techniques capable of handling (in real-time) the growing amounts of data generated by modern fusion research. While existing numerical linear algebra (NLA) methods provide the backbone to classical data analysis and algorithms, these methods fundamentally do not port to distributed architectures nor do they allow low-latency data reduction for control. Motivated by the needs for modern fusion reactors, this project explored and implemented new numerical methods to characterize plasma dynamics, respond in real-time to discharge evolution, and to process massive-scale data accurately and rapidly more fully. This project links expertise in multiple-sensor diagnostics of tokamak plasma dynamics from Columbia University’s Plasma Physics Laboratory with expertise in massive-scale data reduction and extreme data control algorithms at Columbia University’s Data Science Institute. This interdisciplinary project (i) applied machine learning methods, (ii) implemented a properly-trained neural-network for very fast processing of high-speed plasma videography, and (ii) developed the applied mathematical methods, based on randomized-NLA (rNLA) routines, for data analysis, reduction, and real-time control. The Columbia University High Beta Tokamak-Extended Pulse (HBT-EP) facility provided data to test new algorithms and partnership with Columbia University's Data Sciences Institute evaluated the broader use of new algorithms for many challenging control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neural Image Compression: Generalization, Robustness, and Spectral Biases

Recent advances in neural image compression (NIC) have resulted in models which are starting to outperform traditional codecs. While this has led to growing excitement about using these methods in real-world applications, the successful adoption of any machine learning system (including NIC) in the wild requires it to generalize (and be robust) to unseen distribution shifts at deployment time. Unfortunately, current research lacks comprehensive datasets and informative tools to evaluate and understand compression performance in real-world settings. To bridge this crucial gap, first, this paper presents a comprehensive benchmark suite to evaluate the out-of-distribution (OOD) performance of image compression methods. Specifically, we design CLIC-C and Kodak-C by introducing 15 common corruptions to popular CLIC and Kodak benchmarks. Next, we propose spectrally inspired introspection tools to gain a deeper understanding of errors introduced by image compression methods as well as their OOD performance. To this end, we carry out a detailed performance comparison of the classical codec with various variants of NIC (e.g., original, variable rate, pruned), revealing intriguing findings that challenge our current understanding of the strengths and limitations of NIC. Finally, we corroborate our empirical findings with theoretical analysis, providing an in-depth view of the OOD performance of NIC. Our benchmarks, spectral introspection tools, and findings provide a crucial bridge to the real-world adoption of NIC. We hope that our work will propel future efforts in designing more robust and generalizable NIC methods.

neural networks, Variational Autoencoder, robustne↗

Establishing metrics to quantify spatial similarity in spherical and red blood cell distributions

As computational power increases and systems with millions of red blood cells can be simulated, it is important to note that varying spatial distributions of cells may affect simulation outcomes. Since a single simulation may not represent the ensemble behavior, many different configurations may need to be sampled to adequately assess the entire collection of potential cell arrangements. In order to determine both the number of distributions needed and which ones to run, we must first establish methods to identify well-generated, randomly placed cell distributions and to quantify distinct cell configurations. We utilize metrics to assess (1) the presence of any underlying structure to the initial cell distribution and (2) similarity between cell configurations. We propose the use of the radial distribution function to identify long-range structure in a cell configuration and apply it to a randomly distributed and structured set of red blood cells. To quantify spatial similarity between two configurations, we make use of the Jaccard index, and characterize sets of red blood cell and sphere initializations. As an extension to our work submitted to the International Conference on Computational Science, we significantly increase our data set size from 72 to 1048 cells, include a similar set of studies using spheres, compare the effects of varying sphere size, and utilize the Jaccard index distribution to probe sets of extremely similar configurations. Our results show that the radial distribution function can be used as a metric to determine long-range structure in both distributions of spheres and RBCs. We determine that the ideal case of spheres within a cube versus bi-concave shaped cells within a cylinder affects the shape of the Jaccard index distributions, as well as the range of Jaccard values, showing that both the shape of particle and the domain may play a role. Furthermore, we also find that the distribution is able to capture very similar configurations through Jaccard index values greater than 95% when appending several nearly identical configurations into the data set.

59 BASIC BIOLOGICAL SCIENCES↗

Mass agnostic jet taggers

Searching for new physics in large data sets needs a balance between two competing effects—signal identification vs background distortion. In this work, we perform a systematic study of both single variable and multivariate jet tagging methods that aim for this balance. The methods preserve the shape of the background distribution by either augmenting the training procedure or the data itself. Multiple quantitative metrics to compare the methods are considered, for tagging 2-, 3-, or 4-prong jets from the QCD background. This is the first study to show that the data augmentation techniques of Planing and PCA based scaling deliver similar performance as the augmented training techniques of Adversarial NN and uBoost, but are both easier to implement and computationally cheaper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rapid mapping of electrochemical processes in energy-conversion devices

Electrochemical impedance spectroscopy (EIS) is ubiquitously applied to identify physicochemical processes governing the performance of energy-conversion devices. However, deconvolution and interpretation of impedance phenomena are limited by measurement throughput and a dearth of scalable analysis methods. Here, we demonstrate an approach to quickly collect and coherently analyze large volumes of electrochemical data. In this study, we accelerate impedance characterization by combining rapid measurements in time and frequency domains, which are interpretably transformed using the distribution of relaxation times (DRT) and a new distribution of phasances (DOP) model. This method provides excellent agreement with EIS and decreases measurement time by an order of magnitude. High-throughput spectra are then distilled into detailed electrochemical maps. This approach is applied to a Li-ion battery and a protonic ceramic electrochemical cell as practical case studies, demonstrating how mapping can richly characterize physicochemical relationships that are difficult to decipher with conventional measurement and analysis methods.

25 ENERGY STORAGE↗

Mechanical and Electrochemical Implications of Drying Temperature on Lithium‐Ion Battery Electrodes

Abstract Lithium‐ion battery (LIB) electrodes are typically produced with n‐methyl‐2‐pyrrolidone, a toxic solvent that is a known carcinogen and reproductive hazard. Accordingly, aqueous processing has been an expanding area of research interest in the field of LIB manufacturing. Although aqueous processing has been widely successful in anode processing, serious challenges remain in processing the cathode. In this work, the drying mechanics of cathode processed with both solvents is investigated though implementation of a chemical‐engineering‐based model to better understand the utilization of heat provided by experimentally determining the heat and mass transfer coefficients. Electrochemical performance is also evaluated to determine the impact of drying temperature on cycling performance. Binder distribution is determined via various methods to confirm differences in binder homogeneity as a function of both solvent and drying temperature. Identified is the large difference in the efficiency in which the heat is used as well as an ideal drying temperature for both aqueous and non‐aqueous processed cathodes. Also identified is the increased sensitivity to processing temperature for aqueous processed electrodes compared to non‐aqueous processed counterparts, pointing to the possibility of tuned drying regimes which would capitalize on the potential cost savings of aqueous processing for cathodes.

25 ENERGY STORAGE↗

Understanding the Value of Energy Storage for Power System Reliability and Resilience Applications

Abstract Purpose of Review The need for energy storage in the electrical grid has grown in recent years in response to a reduced reliance on fossil fuel baseload power, added intermittent renewable investment, and expanded adoption of distributed energy resources. While the methods and models for valuing storage use cases have advanced significantly in recent years, the value of enhanced resilience remains an open research question. Recent Findings The findings of the recent research indicate that energy storage provides significant value to the grid, with median benefit values for specific use cases ranging from under $10/kW-year for voltage support to roughly $100/kW-year for capacity and frequency regulation services. While the value of lost load is used widely to estimate the benefits of mitigating short-duration outages, reaching as high as $719/kilowatt-year, there is no consensus when it comes to monetizing the value of improving grid resilience. Summary This paper presents a use case taxonomy for energy storage and uses the taxonomy to conduct a meta-analysis of an extensive set of energy storage valuation studies. It reviews several approaches for monetizing reliability and resiliency services and presents a proposed approach for valuing resiliency for energy storage investments.

25 ENERGY STORAGE↗

RFB Side Effects - Shunt Currents

The various models related to fluid flow, electrochemistry and shunt current, and their interactions are reviewed, and gaps identified in the development of flow battery holistic models. In a flow battery stack, the cells are electrically connected through the manifolds via the electrolyte. This results in shunt current through the electrolyte, which is a source of loss in the battery stack. This issue is examined in depth with a mathematical method of calculating shunt current distribution, the associated losses, and strategies to mitigate shunt current losses, along with the trade-offs.

Crawford, Aladsair J.↗

A FETI approach to domain decomposition for meshfree discretizations of nonlocal problems

We propose a domain decomposition method for the efficient simulation of nonlocal problems. Our approach is based on a multi-domain formulation of a nonlocal diffusion problem where the subdomains share “nonlocal” interfaces of the size of the nonlocal horizon. This system of nonlocal equations is first rewritten in terms of minimization of a nonlocal energy, then discretized with a meshfree approximation and finally solved via a Lagrange multiplier approach in a way that resembles the finite element tearing and interconnect method. Specifically, we propose a distributed projected gradient algorithm for the solution of the Lagrange multiplier system, whose unknowns determine the nonlocal interface conditions between subdomains. Several two-dimensional numerical tests on problems as large as 191 million unknowns illustrate the strong and the weak scalability of our algorithm, which outperforms the standard approach to the distributed numerical solution of the problem. Finally, this work is the first rigorous numerical study in a two-dimensional multi-domain setting for nonlocal operators with finite horizon and, as such, it is a fundamental step towards increasing the use of nonlocal models in large scale simulations.

42 ENGINEERING↗

Advocating for view and daylight in buildings: Next steps

With the exponential growth in population and commensurate increased density in urban cities, access to daylight and views to nature is being severely curtailed in buildings. In parallel, increasingly urgent demands to sharply reduce building energy use and associated greenhouse gas emissions are being made to mitigate climate change. There are many challenges and performance tradeoffs associated with the building facade (i.e., daylight and view versus solar and glare control); increased prioritization of health and well-being as a fundamental human requirement could adversely affect building energy-efficiency. Given the current state of knowledge on the effects of daylight and view on health and well-being in buildings, we identify critical needs in research, tools and technologies that if satisfied may enable more effective use of daylight and view in buildings within the constraints of climate change. Lack of knowledge regarding the complex causal mechanisms of window views on human factors is a severely limiting factor in forward progress. Current models and methods to derive bidirectional scattering distribution functions (BSDFs) will need to be modified. Developers of energy-efficient window technologies will need more guidance to shape product development. Advanced window technologies and integrated design can enable attainment of both health and well-being and net zero energy goals, but considerable work will be needed to make such options turnkey and broadly available.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗

Seismic monitoring of CO 2 geosequestration using multi-well 4D DAS VSP: Stage 3 of the CO2CRC Otway project

An important part of any CO 2 geosequestration project is to ensure CO 2 containment and conformance in the subsurface. This is generally done by implementing a comprehensive, risk-based Measurement, Monitoring and Verification plan, a key element of which is active time-lapse seismic monitoring. However, high cost and environmental impact of the standard surface seismic monitoring dictate the need for a cost-effective and environmentally friendly alternative. An opportunity to develop such method emerges with advances in distributed acoustic sensing (DAS) technology, which turns an optical fibre into a seismic sensor with dense spatial sampling. DAS can be permanently deployed in multiple wells across the geosequestration site providing a robust and non-intrusive network of seismic receivers. This approach was developed and tested in the CO2CRC Otway project, where injection of 15 kt of CO 2 at 1.5 km depth was monitored with a 4D vertical seismic profiling (VSP) using five borehole DAS arrays and mobile vibroseis sources. The 4D DAS VSP in each of the five wells provides broadly consistent images of the CO 2 plume with some differences due to different illumination of the target horizon, lateral variation of velocities, and seismic anisotropy. When the newly injected CO 2 reaches a CO 2 plume created as a result of an earlier injection into the same formation ~600 m updip, 4D DAS VSP shows a change in reflectivity in that area and beyond. Furthermore, this shows a potential of 4D DAS VSP for monitoring gas injection into gas-saturated reservoirs.

58 GEOSCIENCES↗