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 1,261 records · Page 70

Exploring the Landscape of Distributed Graph Clustering on Leadership Supercomputers

The rapid growth of large-scale datasets in fields like biology and social networks has driven the need for advanced graph analytics techniques. Community detection, a fundamental task in graph analytics, identifies closely connected groups of nodes within a network, providing valuable insights across various disciplines. This study focuses on two classic community detection methods, the Louvain algorithm and Markov Clustering (MCL), and evaluates the performance of two prominent distributed community detection algorithms: HiPDPL-GPU, our prior implementation, and HipMCL. We conduct experiments on GPU-accelerated heterogeneous HPC systems, Summit and Frontier, to assess their performance under varying conditions. Our objective is to identify the strengths and weaknesses of these algorithms in terms of scalability, and quality of solutions. We evaluate these algorithms on a diverse set of 70+ networks spanning 13 domains, with sizes ranging up to 4.2 billion edges. Our results demonstrate that HiPDPL-GPU consistently outperforms HipMCL, especially for large-scale networks. HiPDPL-GPU achieves significantly faster runtimes (47x to 1439x), higher modularity scores, and improved scalability. These findings highlight HiPDPL-GPU as a promising solution for efficient and effective large-scale graph analytics in diverse application domains, and provide insights into the feasibility of using MCL-based approaches for certain application domains.

Community detection, graph algorithms↗

Guest Editorial: Planning and operation of resilient distribution system for integrated multi-energy

Resilience is the ability of power systems to prepare for and adapt to low-probability, high-impact incidents and withstand and recover rapidly from disruptions. With the ageing of electricity distribution infrastructure and increasing threats of weather-related incidents and natural disasters, the need to effectively enhance the resilience of the electricity distribution system has become urgent and has attracted worldwide attention. Although there are an increasing number of publications related to enhancing resilience strategies, resilience is an emerging concept in power systems. Existing practices are mostly focused on deploying distributed energy resources (DERs) and microgrids, hardening the existing infrastructures and building redundant capacities. However, from a broader perspective, resilience enhancement of power distribution system is a systematic engineering, involving long-term system planning and upgrading (e.g., deployment of smart grid technologies and intelligent switches), short-term proactive scheduling, real-time robust and resilient control of DERs, and post-event restoration and recovery strategies. Based on this point, this Special Issue in IET Energy System Integration focuses on soliciting the most recent and original technologies, scheduling and control strategies for improving the resilience of power distribution system. Eight papers are presented in this Special Issue, covering various aspects related to resilience enhancement of power distribution system, including fault-tolerant frequency measurement, robust scheduling of integrated electricity and district heating systems, robust control of DERs, efficient methods for safety verification as well as novel graph theory-based approach to restore the distribution systems after multiple simultaneous faults. A brief introduction of these 8 papers is provided below.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient Measurement of Length Distribution of 1D Nanoparticles in Solution via Optical Polarimetry

The efficient measurement of the length distribution of nanotubes, nanowires, and other one-dimensional (1D) nanoparticles in solution is important to enable their incorporation into materials and devices and to optimize their processing for properties of interest, such as thermal/electrical conductivity or mechanical strength, in suspensions and composites. We report an electric-field (E-field)-assisted optical-polarimetry technique to measure the length distribution of ensembles of high-aspect-ratio particles in dilute suspension. The degree of alignment of polarizable 1D particles suspended in a fluid under Brownian motion explicitly depends on the E-field strength and the particle length. We show that it is possible to extract the length distribution of 1D nanoparticles suspended in an insulating fluid by applying a range of E-fields and using optical polarimetry to measure the corresponding alignment order parameter. Notably, the method is relatively insensitive to the diameter of the 1D particles, which can be poorly known or vary within a sample. The technique is validated with silver nanowires and carbon nanotubes of known lengths, as well as polymer-depletion-length-separated single-wall carbon nanotube samples with length distributions independently measured with analytical ultracentrifugation. Finally, we demonstrate the ability of the optical-polarimetry technique to quantify changes in the length distribution of ultranarrow, sub-nanometer-diameter single-wall carbon nanotubes under different types and durations of ultrasonication. Within its range of applicability (polarizable 1D nanoparticles in the 0.5 to 15 μm length range, constrained by the voltage stability of the media and the suspended particles), the E-field-assisted optical-polarimetry method is a particularly efficient and accurate method to measure the length distribution of nanowires and nanotubes in suspension.

1D nanoparticles↗

Utilizing the maximum likelihood estimator for flow analysis

We explore the possibility of evaluating flow harmonics by employing the maximum likelihood estimator (MLE). For a given finite multiplicity, the MLE simultaneously furnishes estimations for all the parameters of the underlying distribution function while efficiently suppressing the variance of measures. Also, the method provides a means to assess a specific class of mixed harmonics, which is not straightforwardly feasible by the approaches primarily based on particle correlations. The results are analyzed using the Wald, likelihood ratio, and score tests of hypotheses. Besides, the resultant flow harmonics obtained using MLE are compared with those derived using particle correlations and event plane methods. Here, the dependencies of extracted flow harmonics on the multiplicity of individual events and the total number of events are analyzed. It is shown that the proposed approach works efficiently to deal with the deficiency in detector acceptability. Moreover, we elaborate on a fictitious scenario where the event plane is not a well-defined quantity in the distribution function. For the latter case, the MLE is shown to largely perform better than the two-particle correlation estimator. In this regard, one concludes that the MLE furnishes a meaningful alternative to the existing approaches for flow analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spatial heterogeneity analyses of pore structure and mineral composition of Barnett Shale using X-ray scattering techniques

The spatial heterogeneity of pore structure and mineral composition at mu m to mm scales in mm to cm-sized shale samples complicates the laboratory measurements of porosity & pore size distribution and consequently reservoir quality evaluation. Traditional fluid invasion characterization methods typically take hours to days to get data and can only test relatively small samples (usually 2.54 cm diameter core plugs and mm-sized rock chips). Furthermore, the Xray scattering techniques at the Advanced Photon Source of the Argonne National Laboratory have the capability to non-destructively determine the pore structure and mineral composition of mm x mm area in a short time of seconds and to hold several to hundreds of mm 2 sized sample to conduct large area heterogeneity analysis in the short duration of hours. With the use of ultra-small angle X-ray scattering (USAXS), small-angle X-ray scattering (SAXS), and wide-angle X-ray scattering (WAXS), the porosity, surface area and pore size distribution of two Barnett Shale samples have been mapped out and the mineral compositions are determined qualitatively. Overall, the porosity and pore surface area of the Barnett Shale Sample A over a 36 cm2 area ranges from 1.70 % to 8.13 % and 8.87 m 2 /g to 40.2 m 2 /g respectively. The porosity and surface area of Sample B over a 24 cm 2 area varies from 3.15 % to 11.7 % and 18.3 m 2 /g to 60.2 m 2 /g. In Sample A, the mineral composition shows a shift from carbonate-rich to siliciclastic-rich over a distance of several millimeters. For Sample B, the mineral compositions are mainly siliciclastic and do not show an obvious lithological change. The direct observations of high spatial heterogeneity of pore structure and mineral composition in shale illustrate the utility of integrated X-ray scattering techniques to provide valuable insights into shale studies and reservoir evaluation.

58 GEOSCIENCES↗

Transformer Leakage Inductance Design Methodology

The leakage inductance exhibited by a transformer depends on its winding geometry, which generally involves the selection of several key design parameters in addition to the winding structure and the interleaving configuration. With few resources explaining the effects of these design choices on the observed leakage inductance, numerous trial-and-error iterations become necessary to realize the desired leakage inductance. This paper explores more than a hundred winding geometries feasible in a 2-winding transformer comprising the same magnetic core, number of turns, and wire gauge, and finds the leakage inductance for each unique design using 2-D finite element method (FEM) simulations in association with the semi-analytical double- 2-D model. These leakage inductances are plotted and further analyzed to understand the effects of different design parameters on the effective leakage inductance. The results presented herein and the conclusions drawn from this research can serve as a valuable resource for future design practitioners from both industry and academia.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Accelerating astronomical and cosmological inference with preconditioned Monte Carlo

ABSTRACT We introduce preconditioned Monte Carlo (PMC), a novel Monte Carlo method for Bayesian inference that facilitates efficient sampling of probability distributions with non-trivial geometry. PMC utilizes a Normalizing Flow (NF) in order to decorrelate the parameters of the distribution and then proceeds by sampling from the preconditioned target distribution using an adaptive Sequential Monte Carlo (SMC) scheme. The results produced by PMC include samples from the posterior distribution and an estimate of the model evidence that can be used for parameter inference and model comparison, respectively. The aforementioned framework has been thoroughly tested in a variety of challenging target distributions achieving state-of-the-art sampling performance. In the cases of primordial feature analysis and gravitational wave inference, PMC is approximately 50 and 25 times faster, respectively, than nested sampling (NS). We found that in higher dimensional applications, the acceleration is even greater. Finally, PMC is directly parallelisable, manifesting linear scaling up to thousands of CPUs.

79 ASTRONOMY AND ASTROPHYSICS↗

Roadmap for Advancement of Low-Voltage Secondary Distribution Network Protection

Downtown low-voltage (LV) distribution networks are generally protected with network protectors that detect faults by restricting reverse power flow out of the network. This creates protection challenges for protecting the system as new smart grid technologies and distributed generation are installed. This report summarizes well-established methods for the control and protection of LV secondary network systems and spot networks, including operating features of network relays. Some current challenges and findings are presented from interviews with three utilities, PHI PEPCO, Oncor Energy Delivery, and Consolidated Edison Company of New York. Opportunities for technical exploration are presented with an assessment of the importance or value and the difficulty or cost. Finally, this leads to some recommendations for research to improve protection in secondary networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Sparge Sampling of Molten Salts for Online Monitoring via Laser-Induced Breakdown Spectroscopy

A method was developed to sample molten salts by sparging to generate and transport aerosols to an isolated instrument for compositional analysis by laser-induced breakdown spectroscopy (LIBS). Real-time monitoring of molten salt composition is critical to developing molten salt nuclear reactors, which offer enhanced safety and efficiency. In this article, the sparge sampling method is described and compared with sampling using a Collison nebulizer. The size distribution and transport of aerosols produced from molten eutectic NaNO 3 –KNO 3 salt were compared for multiple gas flow rates (75–1200 mL min –1 ) and transport distances (0.68–2.61 m). Both methods produced aerosols ranging from 0.5 to 5.0 μm determined using a cascade impactor. Aerosols were effectively transported without pre- or trace-heating of gas lines, but transport efficiency was reduced by the formation of agglomerates. Sparge sampling was found to use less sample and less gas than a Collison nebulizer while producing a more concentrated aerosol stream (up to 5 μg L –1 ). The effects of laser energy and delay time on the signal quality of LIBS measurements of these aerosols were also studied. High energy and short delay times were found to enhance signal and repeatability, whereas signal-to-background and signal-to-noise ratios were highest at low energy and longer delay times. The capabilities of this system for online monitoring of molten salts were demonstrated with calibrations for Sr and Li with relative standard deviations of 2.6% and 1.5% and limits of detection of 380 and 180 μg g –1 , respectively.

Aerosols↗

Pair distribution function analysis of nano-object assemblies

The atomic pair distribution function (aPDF) analysis technique, also known as the total scattering method, which considers both Bragg and diffuse scattering, has been used extensively to probe local atomic arrangements in crystalline and disordered materials. In contrast, there have been limited applications of the PDF in self-assembled nanomaterials, which represent a class of materials built from nanoscale objects, such as nano-colloids, micelles and proteins. As distinguished from atoms, nano-objects have polydispersity in size and shape, and such form-factor effects complicate the application of PDF analysis to nano-systems. Here, in this paper, the application of the PDF is extended to spherical nano-object assemblies and the formulae for the nano-PDF (nPDF) are derived, showing some differences from the aPDF. By numerical simulations, the properties of the nPDF (peak broadening and pattern profile) are studied systematically as a function of structural features, such as nano-object parameters (size and size polydispersity) and assembly structural features (size, shape, structure type and lattice disorder), and of data processing parameters ( q cut-off and `missing' data in ultra-small-angle regions). The nPDF analysis method is found to provide an effective route to revealing not only nanoscale but also mesoscale structural properties, for example the morphology of a nano-assembly.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Characterization of Nanoscale Pores in Tight Gas Sandstones Using Complex Techniques: A Case Study of a Linxing Tight Gas Sandstone Reservoir

Pore structures with rich nanopores and permeability in tight gas reservoirs are poorly understood up to date. Advanced techniques are needed to be employed to accurately characterize pore structures, especially tiny pores which include micron and nanopores. In this study, various experimental techniques such as scanning electron microscopy (SEM), nuclear magnetic resonance (NMR) T 2 , nitrogen adsorption method, and NMR cryoporometry (NMRC) are combined to interrogate the complex pore systems of the tight gas reservoir in the Linxing formation, Ordos Basin, China. Results show that tight gas sandstones are primarily comprised of residual interparticle and clay-dominated pores. Clay and quartz are two dominate minerals while pyrite occupies a nontrivial amount as well. The permeability of tight gas sandstones is very low, exhibiting an extremely poor positive correlation with porosity. While pore types and relative pore contents are more influential factors on the permeability, accurate characterization of pore size distribution is critical for the permeability of tight gas sandstones. Therefore, complementary characterization methods are carried out, indicating that neither small pores with radii < 100 nm (around peak 1 in NMR T 2 distribution) nor large pores with radii > 5 μ m (around peak 3 in NMR T 2 distribution) control the permeability by analyzing the connectivity of the pores in various size ranges, but rather pores averaging approximately 350 ± X nm (around peak 2 in NMR T 2 distribution) have sufficient connectivity to host and transmit hydrocarbons. The pore size of tight gas sandstones is dominated by the clay-rich mineral assemblage. The study shows that the NMRC technique can be a very promising method, especially when referred to as a promising “roadmap” on how to interrogate tight formations such as the tight gas sands or even shale especially for the nanopore characterization.

15 GEOTHERMAL ENERGY↗

Harmonic cancellation within AC low voltage distribution for a realistic office environment

An increase of non-linear loads, primarily from power electronics, has substantially increased current harmonics in commercial buildings, which contributes to decreased transformer efficiency / lifespan and poor power quality. This study uses recorded power consumption data from common miscellaneous electric loads (MELs) seen in offices, combined with detailed characterizations of example MELs, to simulate harmonic cancellation within building circuits. Typically, harmonic cancellation studies assume that AC converters operate across their rated power range. However, this study finds that common MELs operate below 40% of rated power the majority of the time when not quiescent; 89% of sampled devices never operated above 60% of rated power. Simulations using these more realistic power levels indicate current-harmonic cancellation (3rd to 13th harmonic) is significantly lower than that predicted when using full-range power assumptions, resulting in minor errors for low-order harmonics and larger errors for higher order harmonics. Furthermore, increased MELs load diversity increases harmonic cancellation, but insufficiently to eliminate errors. In contrast, blending lighting loads with MELs on the secondaries of distribution transformers improves harmonic cancellation to near those predicted by traditional methods. These results indicate that realistic power levels, as well as better characterization of harmonics from typical MELs, should be used to estimate harmonic cancellation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Asynchronous distributed-memory task-parallel algorithm for compressible flows on unstructured 3D Eulerian grids

Here, we discuss the implementation of a finite element method, used to numerically solve the Euler equations of compressible flows, using an asynchronous runtime system (RTS). The algorithm is implemented for distributed-memory machines, using stationary unstructured 3D meshes, combining data-, and task-parallelism on top of the Charm++ RTS. Charm++’s execution model is asynchronous by default, allowing arbitrary overlap of computation and communication. Task-parallelism allows scheduling parts of an algorithm independently of, or dependent on, each other. Built-in automatic load balancing enables continuous redistribution of computational load by migration of work units based on real-time CPU load measurement. The RTS also features automatic checkpointing, fault tolerance, resilience against hardware failure, and supports power-, and energy-aware computation. We demonstrate scalability up to 25 x 10 9 cells at $\mathscr{O}$10 4 compute cores and the benefits of automatic load balancing for irregular workloads. The full source code with documentation is available at https://quinoacomputing.org.

42 ENGINEERING↗

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A method for predicting failure statistics for steady state elevated temperature structural components

This paper presents the initial development of a high temperature life prediction method that accounts for the variability in the material properties of Grade 91 steel. The method accounts for material variability by fitting a variable 3-parameter Weibull distribution to experimental rupture data and accounts for the variability of creep deformation on the steady-state stresses via a Monte Carlo approach. To ensure reasonable computational times, the model represents the material as an extremely viscous Stokes fluid with a non-Newtonian viscosity, therefore solving the stress relaxation problem with a steady, static, instead of transient, analysis. Furthermore, the complete statistical analysis combines this model for creep deformation with a probabilistic model for creep rupture to evaluate the probability of premature failure for a set of sample problems, comparing the predicted failure statistics to the design life predicted by the ASME Boiler and Pressure Vessel Code rules.

42 ENGINEERING↗

Microfracture-pore structure characterization and water-rock interaction in three lithofacies of the Lower Eagle Ford Formation

The microfracture-pore structures and water-rock interactions of eight samples of three typical lithofacies (three limestone, three wackestone, and two mudstone) from the Lower Eagle Ford Formation have been investigated using an integrated methodology. The methods used were X-ray diffraction (XRD), total organic carbon (TOC) content, pyrolysis, thin-section petrography, scanning electron microscopy (SEM), mercury intrusion porosimetry (MIP), spontaneous water imbibition, and water vapor & nitrogen physisorption. Petrophysical properties such as organic richness, maturity, mineral composition, types and patterns of microfractures and pores, porosity, pore-throat size distribution, pore size distribution, pore connectivity, water adsorption and desorption behavior were determined from these methods, and the results compared among the three lithofacies. For example, the three limestone EAS samples and three wackestone DM samples show a low porosity range from 2.02–3.35%, whereas the two mudstone VRU samples exhibit a much large average porosity of 11.7% by MIP results. Nitrogen physisorption results show that pore size distributions have a good similarity in N2 adsorption volume for the samples with similar lithology. Microfractures are morphologically divided into seven types, with their definitions & pictorial examples presented and possible formation mechanisms proposed. In addition, the water vapor physisorption behavior and their controlling factors are discussed. Here, the water adsorption capacity shows a positive relationship (R 2 = 0.90) with clay content and a negative relationship (R 2 = 0.99) with carbonate content. In general, this integrated study of microfracture-pore structure characterization and water-rock interaction provides some insights of water vapor adsorption behavior in different lithofacies and whether microfractures could be an important pathway for fluid flow.

58 GEOSCIENCES↗

Hierarchical, Grid-Aware, and Economically Optimal Coordination of Distributed Energy Resources in Realistic Distribution Systems

Renewable portfolio standards are targeting high levels of variable solar photovoltaics (PV) in electric distribution systems, which makes reliability more challenging to maintain for distribution system operators (DSOs). Distributed energy resources (DERs), including smart, connected appliances and PV inverters, represent responsive grid resources that can provide flexibility to support the DSO in actively managing their networks to facilitate reliability under extreme levels of solar PV. This flexibility can also be used to optimize system operations with respect to economic signals from wholesale energy and ancillary service markets. Here, we present a novel hierarchical scheme that actively controls behind-the-meter DERs to reliably manage each unbalanced distribution feeder and exploits the available flexibility to ensure reliable operation and economically optimizes the entire distribution network. Each layer of the scheme employs advanced optimization methods at different timescales to ensure that the system operates within both grid and device limits. The hierarchy is validated in a large-scale realistic simulation based on data from the industry. Simulation results show that coordination of flexibility improves both system reliability and economics, and enables greater penetration of solar PV. Discussion is also provided on the practical viability of the required communications and controls to implement the presented scheme within a large DSO.

distribution network↗