Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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 109 records · Page 6

Emergency department visits in California associated with wildfire PM 2.5 : differing risk across individuals and communities

The threats to human health from wildfires and wildfire smoke (WFS) in the United States (US) are increasing due to continued climate change. A growing body of literature has documented important adverse health effects of WFS exposure, but there is insufficient evidence regarding how risk related to WFS exposure varies across individual or community level characteristics. To address this evidence gap, we utilized a large nationwide database of healthcare utilization claims for emergency department (ED) visits in California across multiple wildfire seasons (May through November, 2012–2019) and quantified the health impacts of fine particulate matter <2.5 μm (PM 2.5 ) air pollution attributable to WFS, overall and among subgroups of the population. We aggregated daily counts of ED visits to the level of the Zip Code Tabulation Area (ZCTA) and used a time-stratified case-crossover design and distributed lag non-linear models to estimate the association between WFS and relative risk of ED visits. We further assessed how the association with WFS varied across subgroups defined by age, race, social vulnerability, and residential air conditioning (AC) prevalence. Over a 7 day period, PM 2.5 from WFS was associated with elevated risk of ED visits for all causes (1.04% (0.32%, 1.71%)), non-accidental causes (2.93% (2.16%, 3.70%)), and respiratory disease (15.17% (12.86%, 17.52%)), but not with ED visits for cardiovascular diseases (1.06% (–1.88%, 4.08%)). Analysis across subgroups revealed potential differences in susceptibility by age, race, and AC prevalence, but not across subgroups defined by ZCTA-level Social Vulnerability Index scores. These results suggest that PM 2.5 from WFS is associated with higher rates of all cause, non-accidental, and respiratory ED visits with important heterogeneity across certain subgroups. Notably, lower availability of residential AC was associated with higher health risks related to wildfire activity.

54 ENVIRONMENTAL SCIENCES↗

Heat exposure and maternal stress: evidence from the GRAPHS pregnancy cohort in Ghana

Heat exposure has been linked to psychosocial stress, an established antecedent of perinatal depression; however, evidence on heat-related stress during pregnancy in sub-Saharan Africa remains limited. We analyzed psychosocial stress scores and covariate data from the Ghana Randomized Air Pollution and Health Study, linking daily maximum and minimum shaded wet bulb globe temperature (WBGT) metrics to participants’ stress scores derived from the Crisis in Family Systems-Revised Life Events Questionnaire. We evaluated associations using ordinal logistic regression of pregnancy-average and trimester-average exposures and distributed lag non-linear models (DLNMs) to assess time-varying associations across gestation. Higher average maximum WBGT exposure across pregnancy was associated with increased odds of higher psychosocial stress; each 1 °C increase in maximum WBGT was associated with 64% higher odds of belonging to a higher stress category (OR = 1.64; 95% CI = 1.17–2.31; p = 0.0040). In trimester-average models, higher first-trimester maximum WBGT was also associated with higher stress (OR = 1.44; 95% CI = 1.15–1.81; p = 0.0014). DLNMs suggested that relatively cooler daily maximum WBGT values (25th percentile) were associated with decreased odds of stress in early pregnancy, whereas extreme daily maximum WBGT values (99th percentile) showed a pattern consistent with increased odds of stress in mid-to-late gestation. These findings highlight gestational windows in which heat exposure may influence stress, emphasizing the need for further research into underlying mechanisms and effective interventions to protect maternal mental health in heat-vulnerable settings.

White, Lewis [Columbia University] (ORCID:00090005↗

Generalized Analytical Estimation of Sensitivity Matrices in Unbalanced Distribution Networks

Fast and accurate estimation of sensitivity matrices is significant for the enhancement of distribution system modeling and automation. Analytical estimations have mainly focused on voltage magnitude sensitivity to active/reactive power injections for unbalanced networks with Wye-connected loads and neglecting DERs' smart inverter functionality. Hence, this paper enhances the scope of analytical estimation of sensitivity matrices for unbalanced networks with 1- Φ, 2- Φ, and 3- Φ Delta/Wye-connected loads, DERs with smart inverter functionality, and substation/line step-voltage regulators (SVR). A composite bus model comprising of DER, Delta- and Wye-connected load is proposed to represent a generic distribution bus, which can be simplified to load, PV, or voltage-controlled bus as required. Furthermore, the proposed matrix-based analytical method consolidates voltage magnitude and angle sensitivity to active/reactive power injection and tap-position of all SVRs into a single algorithm. Extensive case studies on IEEE and EPRI networks show the accuracy and wide scope of the proposed algorithm compared to the existing benchmark method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

HFTS-2 (Final Report)

This report delves into the research program, identified tasks, various studies, observations, results, and conclusions under the Hydraulic Fracturing Test Site (HFTS-2) project. This test site is located in Loving County close to the Texas, New Mexico border. The program targeted the Wolfcamp formation in Delaware Basin. The aim of the project was to improve our understanding of hydraulic fracturing processes to not only optimize unconventional resource development in US Basins but to minimize environmental footprint as well. Various new diagnostic technologies were demonstrated at this site, including distributed acoustic microseismic, high resolution strain using Rayleigh frequency band (DSS-RFS), and in-fill proppant logging tool. Significant conclusions from this study include: 1) lateral and vertical fracture growth mapping as well as proppant transport behavior, 2) Use of DSS-RFS as a diagnostic tool to understand cluster level behavior, 3) Validation of proppant log, both at core scale as well as in-fill well implementation, and 4) New modeling tools to map far-field strain response using distributed strain sensing data. Specifically, far-field fracture distribution is influenced by linear fracture corridors defined by the stress-state as well as clustering behavior influenced by spatio-temporal effects. Depletion impacts on new stimulation was demonstrated and technically quantified. Significant vertical fracture height growth was observed. However, our understanding is limited by available gauge data availability. Numerical strain modeling techniques have been developed and demonstrated to help understand observed strain response behavior. Finally, proppant log has been further validated as a useful tool to map propped zones within the SRV. It has also shown very strong correlation with observed drawdown behavior at this test site. This report will detail many of these studies as well as technical results and conclusions from said studies. These include general subsurface characterization work, completion designs as implemented, studies looking at fracture geometry, high resolution microseismic study including source mechanisms, advanced DSS-RFS studies including drainage characterization, core characterization results, proppant analysis (core and in-fill child well), geochemistry, etc.

02 PETROLEUM↗

Bias correcting regional scale Earth system model projections: novel approach using empirical mode decomposition

Bias correction is a crucial step in using Earth system model outputs for assessments, as it adjusts systematic errors by comparing the model to observations. However, standard methods – ranging from mean-based linear scaling to distribution-based quantile mapping typically treat bias correction as a single-scale process, overlooking the fact that biases can manifest differently across daily, seasonal, and annual timescales. In this study, we propose a novel, timescale-aware bias-correction approach built on Empirical Mode Decomposition. By decomposing the meteorological signal into multiple oscillatory components and aggregating them to represent distinct timescales, we apply targeted corrections to each component, thereby preserving both short- and long-term structure in the data. Experimental illustrations show that the timescale-aware EMDBC framework matches the performance of conventional quantile-delta mapping (QDM) at the native daily scale and achieves progressively larger bias reductions at bi-weekly, seasonal, and annual scales. As a result, the proposed approach offers a more robust path to accurate and reliable Earth system projections, strengthening their utility for resilience and adaptation planning.

Ganguli, Arkaprabha [Argonne National Laboratory (↗

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. Here, we demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

97 MATHEMATICS AND COMPUTING↗

Data-Driven Generic Turbines for Distributed Wind Modeling, Optimization, and Economic Studies

As distributed energy resources (DER) become less expensive and more popular, utilities, project developers, and customers have an increasing need to model the performance of existing and proposed DER systems. Distributed wind has been shown to have widespread economic potential but is often represented by a simplified model in - or excluded from - DER modeling tools and studies. There is often no economic imperative to extend models and studies to give full consideration to distributed wind. We present a set of data-driven generic turbines derived from 16 years of annual distributed wind market survey data. The proposed methodology can be used to derive generic turbines from separate or updated data sets. Finally, a mixed-integer linear programming approach to optimal distributed wind project sizing is used to demonstrate the generic turbine models. Combined, these models and methods can reduce barriers to considering distributed wind in modeling tools and studies.

Reiman, Andrew P.↗

Gradient Coding With Iterative Block Leverage Score Sampling

Gradient coding is a method for mitigating straggling servers in a centralized computing network that uses erasure-coding techniques to distributively carry out first-order optimization methods. Randomized numerical linear algebra uses randomization to develop improved algorithms for large-scale linear algebra computations. In this study, we propose a method for distributed optimization that combines gradient coding and randomized numerical linear algebra. The proposed method uses a randomized ℓ 2 -subspace embedding and a gradient coding technique to distribute blocks of data to the computational nodes of a centralized network, and at each iteration the central server only requires a small number of computations to obtain the steepest descent update. The novelty of our approach is that the data is replicated according to importance scores, called block leverage scores, in contrast to most gradient coding approaches that uniformly replicate the data blocks. Furthermore, we do not require a decoding step at each iteration, avoiding a bottleneck in previous gradient coding schemes. We show that our approach results in a valid ℓ 2 -subspace embedding, and that our resulting approximation converges to the optimal solution.

97 MATHEMATICS AND COMPUTING↗

Improved Evaluation of Large Network Matrices for Linear Power Flow Within Optimization Problems: Preprint

This work discusses methods for evaluating the Power Transfer Distribution Factor (PTDF) and Line Outage Distribution Factor (LODF) matrices by employing sparse linear algebra for large-scale computing applications. These matrices are critical in many power systems applications, such as the Unit Commitment Problem (UC), pre- and post-contingency power flow analysis, and transmission expansion. These matrices are typically dense, which means they require a significant amount of time and memory to be computed for large networks. However, by analyzing the structure of the matrices and their computation method, it is possible to use reduced memory methods based on sparse matrix operations. This paper shows that sparse linear algebra algorithms are faster and require less memory and time than traditional dense approaches. Additionally, we explore the effect of matrix sparsification by eliminating trailing digits on power flow calculations.

ENERGY PLANNING, POLICY, AND ECONOMY↗

Jet wake from linearized hydrodynamics

We explore how to improve the hybrid model description of the particles originating from the wake that a jet produced in a heavy ion collision leaves in the droplet of quark-gluon plasma (QGP) through which it propagates, using linearized hydrodynamics on a background Bjorken flow. Jet energy and momentum loss described by the hybrid model become currents sourcing linearized hydrodynamics. By solving the linearized hydrodynamic equations numerically, we investigate the development of the wake in the dynamically evolving droplet of QGP, study the effect of viscosity, scrutinize energy-momentum conservation, and check the validity of the linear approximation. We find that linearized hydrodynamics works better in the viscous case because diffusive modes damp the energy-momentum perturbation produced by the jet. We calculate the distribution of particles produced from the jet wake by using the Cooper-Frye prescription and find that both the transverse momentum spectrum and the distribution of particles in azimuthal angle are similar in shape in linearized hydrodynamics and in the hybrid model. Their normalizations are different because the momentum-rapidity distribution in the linearized hydrodynamics analysis is more spread out, due to sound modes. Since the Bjorken flow has no transverse expansion, we explore the effect of transverse flow by using local boosts to add it into the Cooper-Frye formula. After including the effects of transverse flow in this way, the transverse momentum spectrum becomes harder: more particles with transverse momenta bigger than 2 GeV are produced than in the hybrid model. Although we defer implementing this analysis in a jet Monte Carlo, as would be needed to make quantitative comparisons to data, we gain a qualitative sense of how the jet wake may modify jet observables by computing proxies for two example observables: the lost energy recovered in a cone of varying open angle, and the fragmentation function. We find that linearized hydrodynamics with transverse flow effects added improves the description of the jet wake in the hybrid model in just the way that comparison to data indicates is needed. Our study illuminates a path to improving the description of the wake in the hybrid model, highlighting the need to take into account the effects of both transverse flow and the broadening of the energy-momentum perturbation in spacetime rapidity on particle production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parallel String Graph Construction and Transitive Reduction for De Novo Genome Assembly

One of the most computationally intensive tasks in computational biology is de novo genome assembly, the decoding of the sequence of an unknown genome from redundant and erroneous short sequences. A common assembly paradigm identifies overlapping sequences, simplifies their layout, and creates consensus. Despite many algorithms developed in the literature, the efficient assembly of large genomes is still an open problem. In this work, we introduce new distributed-memory parallel algorithms for overlap detection and layout simplification steps of de novo genome assembly, and implement them in the diBELLA 2D pipeline. Our distributed memory algorithms for both overlap detection and layout simplification are based on linear-algebra operations over semirings using 2D distributed sparse matrices. Our layout step consists of performing a transitive reduction from the overlap graph to a string graph. We provide a detailed communication analysis of the main stages of our new algorithms. diBELLA 2D achieves near linear scaling with over 80% parallel efficiency for the human genome, reducing the runtime for overlap detection by 1.2-1.3× for the human genome and 1.5-1.9× for C.elegans compared to the state-of-the-art. Our transitive reduction algorithm outperforms an existing distributed-memory implementation by 10.5-13.3× for the human genome and 18-29× for the C. elegans. Our work paves the way for efficient de novo assembly of large genomes using long reads in distributed memory.

59 BASIC BIOLOGICAL SCIENCES↗

Machine learning based unfolding of x-ray spectra from filter stack spectrometer data

We demonstrate the application of neural networks to perform x-ray spectra unfolding from data collected by filter stack spectrometers. A filter stack spectrometer consists of a series of filter-detector pairs, where the detectors behind each filter measure the energy deposition through each layer as photo-stimulated luminescence (PSL). The network is trained on synthetic data, assuming x-rays of energies < 1 MeV and of two different distribution functions (Maxwellian and Gaussian) and the corresponding measured PSL values obtained from five different filter stack spectrometer designs. Predicted unfolds of single distributions are near identical reproductions of the ground truth spectra, with differences in the values lower than 20% at the higher energy end in some cases. The neural network has also demonstrated robustness to experimental measurement errors of < 5% and some capability of performing unfolds for linear combinations of the two distributions without previous training. The network can perform unfolds at rates > 1 Hz, ideal for application to some high-repetition-rate systems.

47 OTHER INSTRUMENTATION↗

Characterizing degraded visual environments with scanning lenses to improve contrast recovery

Atmospheric aerosols, such as fog and smoke, degrade optical signals via the scattering and absorption of light, making it difficult to recover information from the surrounding environment. Here, we present a technique that measures the angular distribution of light in scattering environments to help recover these degraded optical signals. This is done by passing a collimated beam through the aerosol environment and collecting the scattered and unscattered light with an f-theta scan lens. The lens transforms the angular distribution of the scattered light into the linear domain, mapping the collected light onto a detector. The measured angular scattering distribution is then used to estimate the point spread function of the environment at an arbitrary stand-off distance. This point spread function can be used to deconvolve images blurred by the aerosol environment in which the angular scattering distribution was initially measured. With this approach, we demonstrate improvements in the effective resolution of deblurred images up to 26% over a 120-cm stand-off distance. Experimentation was performed in the Sandia National Laboratories tabletop fog chamber, a platform that generates aerosol environments analogous to real-world fogs. We intend to leverage this technique in future work to develop a low-SWaP system for the in situ characterization of aerosol environments.

Aerosols↗

unWISE tomography of Planck CMB lensing

CMB lensing tomography, or the cross-correlation between CMB lensing maps and large-scale structure tracers over a well-defined redshift range, has the potential to map the amplitude and growth of structure over cosmic time, provide some of the most stringent tests of gravity, and break important degeneracies between cosmological parameters. In this work, we use the unWISE galaxy catalog to provide three samples at median redshifts z ~ 0.6, 1.1 and 1.5, fully spanning the Dark Energy dominated era, together with the most recent Planck CMB lensing maps. We obtain a combined cross-correlation significance S/N = 79.3 over the range of scales 100 < ℓ < 1000. We measure the redshift distribution of unWISE sources by a combination of cross-matching with the COSMOS photometric catalog and cross-correlation with BOSS galaxies and quasars and eBOSS quasars. We also show that magnification bias must be included in our analysis and perform a number of null tests. In a companion paper, we explore the derived cosmological parameters by modeling the non-linearities and propagating the redshift distribution uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-Driven Distributed Algorithms for Estimating Eigenvalues and Eigenvectors of Interconnected Dynamical Systems

Here, the paper presents data-driven algorithms to estimate in a distributed manner the eigenvalues, right and left eigenvectors of an unknown linear (or linearized) interconnected dynamic system. In particular, the proposed algorithms do not require the identification of the system model in advance before performing the estimation. As a first step, we consider interconnected dynamical system with distinct eigenvalues. The proposed strategy first estimates the eigenvalues using the well-known Prony method. The right and left eigenvectors are then estimated by solving distributively a set of linear equations. One important feature of the proposed algorithms is that the topology of communication network used to perform the distributed estimation can be chosen arbitrarily, given that it is connected, and is also independent of the structure or sparsity of the system (state) matrix. The proposed distributed algorithms are demonstrated via a numerical example.

97 MATHEMATICS AND COMPUTING↗

Distributed Data-Driven Power Iteration for Strongly Connected Networks

Here, this paper presents data-driven power iteration to distributively estimate the dominant eigenvalues of an unknown linear time-invariant system. The proposed strategy only requires a single trajectory data or measurements. Furthermore, in order to perform the distributed estimation, the communication network topology can be chosen to be any strongly connected directed graphs. The proposed data-driven power iteration is demonstrated using several numerical examples and is then applied to estimate the generalized algebraic connectivity of cooperative systems and to control the epidemic spreading.

Gusrialdi, Azwirman↗