Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed algorithms”

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,621 records · Page 90

Snow Distribution Patterns Revisited: A Physics-Based and Machine Learning Hybrid Approach to Snow Distribution Mapping in the Sub-Arctic

Snowpack distribution in Arctic and alpine landscapes often occurs in repeating, year-to-year patterns due to local topographic, weather, and vegetation characteristics. Previous studies have suggested that with years of observational data, these snow distribution patterns can be statistically integrated into a snow process modeling workflow. Recent advances in snow hydrology and machine learning (ML) have increased our ability to predict snowpack distribution using in-situ observations, remote sensing data sets, and simple landscape characteristics that can be easily obtained for most environments. Here, we propose a hybrid approach to couple a ML snow distribution pattern (MLSDP) map with a physics-based, snow process model. We trained a random forest ML algorithm on tens of thousands of snow survey observations from a subarctic study area on the Seward Peninsula, Alaska, collected during peak snow water equivalent (SWE). We validated hybrid model outputs using in-situ snow depth and SWE observations, as well as a light detection and ranging data set and a distributed temperature profiling sensor data set. When the hybrid results were compared with the physics-based method, the hybrid method more accurately depicted the spatial patterns of the snowpack, areas of drifting snow, and years when no in-situ observations were used in the random forest ML training data set. The hybrid method also showed improvements in root mean squared error at 61% of locations where time-series estimations of snow depth were observed. These results can be applied to any physics-based model to improve the snow distribution patterning to reflect observed conditions in high latitude and high elevation cold region environments.

54 ENVIRONMENTAL SCIENCES↗

Enhancing Gaussian Process Surrogates for Optimization and Posterior Approximation via Random Exploration

This paper proposes novel noise-free Bayesian optimization strategies that rely on a random exploration step to enhance the accuracy of Gaussian process surrogate models. The new algorithms retain the ease of implementation of the classical GP-UCB algorithm, but the additional random exploration step accelerates their convergence, nearly achieving the optimal convergence rate. Furthermore, to facilitate Bayesian inference with intractable likelihoods, we propose to utilize optimization iterates for maximum a posteriori estimation to build a Gaussian process surrogate model for the unnormalized log-posterior density. We provide bounds for the Hellinger distance between the true and the approximate posterior distributions in terms of the number of design points. We demonstrate the effectiveness of our Bayesian optimization algorithms in nonconvex benchmark objective functions, in a machine learning hyperparameter tuning problem, and in a black-box engineering design problem. The effectiveness of our posterior approximation approach is demonstrated in two Bayesian inference problems for parameters of dynamical systems.

Bayesian inference↗

SU-E-J-20: Adaptive Aperture Morphing for Online Correction for Prostate Cancer Radiotherapy

Purpose: Online adaptive aperture morphing is desirable over translational couch shifts to accommodate not only the target position variation but also anatomic changes (rotation, deformation, and relation of target to organ-atrisks). We proposed quick and reliable method for adapting segment aperture leaves for IMRT treatment of prostate. Methods: The proposed method consists of following steps: (1) delineate the contours of prostate, SV, bladder and rectum on kV-CBCT; (2) determine prostate displacement from the rigid body registration of the contoured prostate manifested on the reference CT and the CBCT; (3) adapt the MLC segment apertures obtained from the pre-treatment IMRT planning to accommodate the shifts as well as anatomic changes. The MLC aperture adaptive algorithm involves two steps; first move the whole aperture according to prostate translational/rotational shifts, and secondly fine-tune the aperture shape to maintain the spatial relationship between the planning target contour and the MLC aperture to the daily target contour. Feasibility of this method was evaluated retrospectively on a seven-field IMRT treatment of prostate cancer patient by comparing dose volume histograms of the original plan and the aperture-adjusted plan, with/without additional segments weight optimization (SWO), on two daily treatment CBCTs selected with relative large motion and rotation. Results: For first daily treatment, the prostate rotation was significant (12degree around lateral-axis). With apertureadjusted plan, the D95 to the target was improved 25% and rectum dose (D30, D40) was reduced 20% relative to original plan on daily volumes. For second treatment-fraction, (lateral shift = 6.7mm), after adjustment target D95 improved by 3% and bladder dose (D30, maximum dose) was reduced by 1%. For both cases, extra SWO did not provide significant improvement. Conclusion: The proposed method of adapting segment apertures is promising in treatment position correction, including target translational displacement, rotation and deformation. Additional SWO could improve ROIs dose distribution.

60 APPLIED LIFE SCIENCES↗

Determination of Collision Cross Sections Using a Fourier Transform Electrostatic Linear Ion Trap Mass Spectrometer

Collision cross sections (CCSs) were determined from the frequency-domain linewidths in a Fourier transform electrostatic linear ion trap. With use of an ultrahigh-vacuum precision leak valve and nitrogen gas, transients were recorded as the background pressure in the mass analyzer chamber was varied between 4× 10{sup -8} and 7 × 10{sup -7} Torr. The energetic hard-sphere ion–neutral collision model, described by Xu and coworkers, was used to relate the recorded image charge to the CCS of the molecule. In lieu of our monoisotopically isolating the mass of interest, the known relative isotopic abundances were programmed into the Lorentzian fitting algorithm such that the linewidth was extracted from a sum of Lorentzians. Although this works only if the isotopic distribution is known a priori, it prevents ion loss, preserves the high signal-to-noise ratio, and minimizes the experimental error on our homebuilt instrument. Six tetraalkylammonium cations were used to correlate the CCS measured in the electrostatic linear ion trap with that measured by drift-tube ion mobility spectrometry, for which there was an excellent correlation (R{sup 2} ≈ 0.9999). Although the absolute CCSs derived with our method differ from those reported, the extracted linear correlation can be used to correct the raw CCS. With use of [angiotensin II]{sup 2+} and reserpine, the corrected CCSs (334.9 ± 2.1 and 250.1 ± 0.5, respectively) were in good agreement with the reported ion mobility spectrometry CCSs (335 and 254.3, respectively). With sufficient signal-to-noise ratio, the CCSs determined are reproducible to within a fraction of a percent, comparable to the uncertainties reported on dedicated ion mobility instruments. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Big Data Framework with Machine Learning for D and D Applications - 19108

The nuclear industry is experiencing a steady increase in maintenance costs even though plants are maintained under high levels of safety, capability and reliability. Nuclear power plants are expected to run every unit at maximum capacity at all times, efficiently utilizing assets with minimal downtime. Surveillance and maintenance of nuclear decommissioning infrastructure provides lot challenges with respect to maintenance or decommissioning of the buildings. There is a need for a framework to analyze the huge amount of data generated by the sensors on the nuclear reactor components as well as structures, to monitor the conditions of these building over a period of time. Emerging technologies such as big data analytics have become a requirement in the nuclear industry to improve structural health monitoring and diagnostics. FIU will make use of existing mature technologies in the areas of imaging, robotics, big data, and machine learning/deep learning to assess the structural integrity of aging facilities at DOE sites. As these facilities await decommissioning, there is a need to understand the structural health of these structures. Many of these facilities were built over 50 years ago and in some cases these facilities have gone beyond the operational life expectancy. In other cases, the facilities have been placed in a state of 'cold and dark' and they are sitting unused, awaiting decommissioning. Finally, some older facilities are one-of-a-kind operational/production facilities supporting critical DOE missions and cannot be shut down. In any of these scenarios, the structural integrity of these facilities may be compromised, so it is imperative that adequate inspections and data collection/analysis be performed on a continuous and ongoing basis. The primary goals of the research include collecting various formats of data such as structured, unstructured and streaming data from the various sensors deployed in buildings, collect videos/pictures from various imaging sources, ingest them using a Hadoop distributed file system and process the data using Spark to perform batch processing and real time analytics. Research and development on various machine learning/deep learning algorithms will also be performed to analyze the heterogeneous data collected from nuclear decommission infrastructure. FIU will design the big data framework to ingest, store and process huge amounts of heterogeneous data collected from many sources and optimize the algorithms to provide insights into the data and predict anomalies observed when compared against baseline conditions. Various modules of this framework will include heterogeneous data sources, message broker, Hadoop distributed file system, Spark for stream and batch processing, machine learning/ deep learning, Cassandra - persistent data store and visualization. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Beneficial Integration of Energy Storage and Load Management with Photovoltaic (PV)

In recent years, a number of industry activities have aimed at addressing the integration challenges posed by the variability and uncertainty of higher penetration of renewable generation sources, like solar photovoltaic (PV) – one of the key objectives of the Sustainable and Holistic Integration of Energy Storage and Solar PV (SHINES) program launched by the U.S. Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE). This EPRI led Beneficial Integration of Energy Storge, and Load Management with PV project aimed to design, develop, and demonstrate end-to-end distributed energy resource (DER) integration solution to build on these activities. EPRI led project team designed and implemented a local controller that uses model predictive control (MPC) algorithm to optimally manage DERs on site by planning for a receding horizon while executing the control settings for the first step of its plan. The team has also developed a system controller to interface with the local controller to demonstrate the hierarchical control and how it can leverage DER capabilities to address challenges like over voltage and thermal limit violations which typically limits the DER hosting capacity of distribution feeders. Team has demonstrated how the local controller with optimization algorithm can effectively manage controllable loads like HVAC, water heater, and pool pumps to allow for greater integration of PV with relatively smaller energy storage system requirements. Optimal utilization of the load control can also reduce the depth of discharge of batteries to meet grid export/import limit from behind-the-meter (BTM) DERs. Proper utilization of DER capabilities via local control intelligence, like the one developed and demonstrated in this project can help the industry to address integration challenges of higher penetration of solar PV in economically efficient manner. This can help to accelerate deployment of clean renewable energy systems at lower societal cost.

14 SOLAR ENERGY↗

Disentangling instrumental features of the 130 GeV Fermi line

We study the instrumental features of photons from the peak observed at E{sub γ} = 130 GeV in the spectrum of Fermi-LAT data. We use the sPlots algorithm to reconstruct — seperately for the photons in the peak and for background photons — the distributions of quantities such as incident angles, the recorded time, features of the spacecraft position, the zenith angles, the conversion type and details of the energy and direction reconstruction. The presence of a striking feature or cluster in such a variable would suggest an instrumental cause for the peak. In the publically available data, we find several suggestive features which may inform further studies by instrumental experts, though the size of the signal sample is too small to draw statistically significant conclusions.

79 ASTRONOMY AND ASTROPHYSICS↗

Domi

Domi (Greek for "structure") is a Trilinos package that provides support for multi-dimensional data arrays distributed over multiple processors of a high performance computing system. The software supports data arrays of arbitrary number of dimensions mapped to a multi-dimensional grid of processors. Domi provides local access to the data, as well as the ability to compute global quantities such as norms, maxima and minima. Data arrays can be constructed with boundary and communication padding, necessary for finite difference algorithms, and provides high-level methods for communicating data between processors to update these buffers. The distributed maps and vectors in Domi can be converted to the equivalent data structures in Epetra and Tpetra (both of which assume unstructured data), using view semantics when appropriate, for use with other Trilinos solvers.

Spotz, WilliamF.↗

An advanced joint inversion system for CO 2 storage modeling with large date sets for characterization and real-time monitoring-enhancing storage performance and reducing failure risks under uncertainties

As large-scale, commercial storage projects become operational, the problem of utilizing information from diverse sources becomes more critically important. In this project, we developed, tested, and applied an advanced joint data inversion system for CO 2 storage modeling with large data sets for use in site characterization and real-time monitoring. Emphasis was on the development of advanced and efficient computational algorithms for joint inversion of hydro-geophysical data, coupled with state-of-the-art forward process simulations. The developed system consists of (1) inversion tools using characterization data, such as 3D seismic survey (amplitude images), borehole log and core data, as well as hydraulic, tracer and thermal tests before CO 2 injection, (2) joint inversion tools for updating the geologic model with the distribution of rock properties, thus reducing uncertainty, using hydro-geophysical monitoring data, and (3) highly efficient algorithms for directly solving the dense or sparse linear algebra systems derived from the joint inversion. The system combines methods from stochastic analysis, fast linear algebra, and high performance computing. The developed joint inversion tools have been tested through synthetic CO 2 storage examples.

01 COAL, LIGNITE, AND PEAT↗

Domi

Domi (Greek for "structure") is a Trilinos package that provides support for multi-dimensional data arrays distributed over multiple processors of a high performance computing system. The software supports data arrays of arbitrary number of dimensions mapped to a multi-dimensional grid of processors. Domi provides local access to the data, as well as the ability to compute global quantities such as norms, maxima and minima. Data arrays can be constructed with boundary and communication padding, necessary for finite difference algorithms, and provides high-level methods for communicating data between processors to update these buffers. The distributed maps and vectors in Domi can be converted to the equivalent data structures in Epetra and Tpetra (both of which assume unstructured data), using view semantics when appropriate, for use with other Trilinos solvers.

Spotz, William F.↗

Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics

Abstract The weight function W ( X ) is a diagnostic sensitivity to phase-space variables X that relates the measured signal C to the distribution function F ( X ) through the equation C = ∫ W ( X ) F ( X ) d X . In the present work, an algorithm to calculate W for a diagnostic that measures 3 MeV protons produced in d(d,p)t fusion reactions between a fast ion and a thermal deuteron is developed. The emitted protons escape the tokamak on curved orbits and are detected. These curved orbits constitute effective diagnostic ‘sightlines’. The presented algorithm accounts for the complications associated with these curved sightlines. An initial calculation of time-reversed proton orbits computes effective solid angles and sightlines for the relevant range of incident proton velocity vectors. These precomputed orbits are then used within the framework of FIDASIM (2020 Plasma Phys. Cont. Fusion 62 105008) to calculate the reactivity averaged over the thermal distribution of the ‘target’ deuterons and the probability that a fast ion of specified energy and pitch has a gyroangle that is consistent with the kinematic equations along each of the sightlines. Comparisons with analytic formulas and with independent calculations for the Mega Amp Spherical Tokamak 3 MeV proton diagnostic verify the algorithm.

Physics↗

A BAYESIAN ESTIMATE OF THE CMB–LARGE-SCALE STRUCTURE CROSS-CORRELATION

Evidences for late-time acceleration of the universe are provided by multiple probes, such as Type Ia supernovae, the cosmic microwave background (CMB), and large-scale structure (LSS). In this work, we focus on the integrated Sachs–Wolfe (ISW) effect, i.e., secondary CMB fluctuations generated by evolving gravitational potentials due to the transition between, e.g., the matter and dark energy (DE) dominated phases. Therefore, assuming a flat universe, DE properties can be inferred from ISW detections. We present a Bayesian approach to compute the CMB–LSS cross-correlation signal. The method is based on the estimate of the likelihood for measuring a combined set consisting of a CMB temperature and galaxy contrast maps, provided that we have some information on the statistical properties of the fluctuations affecting these maps. The likelihood is estimated by a sampling algorithm, therefore avoiding the computationally demanding techniques of direct evaluation in either pixel or harmonic space. As local tracers of the matter distribution at large scales, we used the Two Micron All Sky Survey galaxy catalog and, for the CMB temperature fluctuations, the ninth-year data release of the Wilkinson Microwave Anisotropy Probe ( WMAP 9). The results show a dominance of cosmic variance over the weak recovered signal, due mainly to the shallowness of the catalog used, with systematics associated with the sampling algorithm playing a secondary role as sources of uncertainty. When combined with other complementary probes, the method presented in this paper is expected to be a useful tool to late-time acceleration studies in cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-F-T-216: Evaluating Dosimetry Accuracy of a Treatment Planning System On Small Proton Fields

Purpose: This study is aiming to identify the smallest field size for which a treatment planning system (TPS) can accurately calculate the relative dose distribution. The finding would be used as a guideline to choose the smallest proton field for clinical treatment. Methods: Mevion S250™ double scattering proton delivery system and Eclipse™ TPS (Varian) with pencil beam convolution (PBC) dose algorithm were used in this study. Square sized fields of 1 cm, 2 cm, 3 cm, 4 cm, 5 cm, and 10 cm were planned on a cubical water phantom with iso-center placed at 10 cm depth. All beams used the same proton beam option: range 15 cm and modulation 10 cm. Dose in water was calculated without any compensator. Gafchromic™ EBT3 film and diode detectors were used to measure the central axis dose distribution and lateral dose profiles at 5 cm, 10 cm, and 14 cm depth. Results: The preliminary film measurement shows good agreement between Eclipse calculated lateral dose profiles for all tested field sizes. The differences on full width half maximum were ≤ 1 mm while the differences on the penumbras were between 1 mm and 2 mm between Eclipse and film. For the depth dose, Eclipse results matched well with film measurements for field sizes down to 2 cm{sup 2}. With smaller field size of 1 cm{sup 2}, Eclipse was able to predict the decreasing of SOBP due to the lack of lateral charged particle equilibrium in depth. However, it did not match the film measurement. Diode measurement results will be available at the time of presentation. Conclusion: The PBC dose algorithm in Eclipse can accurately calculate relative dose distribution in double scattered proton system for field size down to 2 cm{sup 2}.

60 APPLIED LIFE SCIENCES↗

Field Validated Utility-Scale Battery Storage Control Models for Quasi-Static Analyses

This paper presents a co-simulation based approach to implement and validate advanced battery storage control modes in OpenDSS that replicate behavior of a utility scale vanadium red-ox flow battery (VRFB). Inverter-interfaced utiliy-scale battery storage systems (BSS) have modes of operation to facilitate multiple grid-support functions. In order to enable accurate understanding of the impacts from grid-support functions modes of operation it is required to develop validated models that replicates behavior of the actual control system. This paper presents the pseudo code for the implemented control schemes including peak shaving and capacity firming control modes. Additionally, this paper also provides insight on convergence sensitivity of the presented algorithms to related parameters. The results obtained from simulation studies will be compared against field operated utility-scale BSS measurements to verify the implemented algorithms.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Stochastic robust mathematical programming model for power system optimization

We present a stochastic robust framework for two-stage power system optimization problems with uncertainty. The model optimizes the probabilistic expectation of different worst-case scenarios with different uncertainty sets. In addition, a case study of unit commitment shows the effectiveness of the proposed model and algorithms.

97 MATHEMATICS AND COMPUTING↗

SDSS-IV MaStar: Stellar parameter determination with continuum-supplemented full-spectrum fitting

We present a stellar parameter catalog built to accompany the MaStar Stellar Library, which is a comprehensive collection of empirical, medium-resolution stellar spectra. We constructed this parameter catalog by using a multicomponent χ 2 fitting approach to match MaStar spectra to models generated by interpolating the ATLAS9-based BOSZ model spectra. The total χ 2 for a given model is defined as the sum of components constructed to characterize narrow-band features of observed spectra (e.g., absorption lines) and the broadband continuum shape separately. Extinction and systematics due to flux calibration were taken into account in the fitting. The χ 2 distribution for a given region of model space was sampled using a Markov chain Monte Carlo (MCMC) algorithm, the data from which were then used to extract atmospheric parameter estimates (T eff , log g, [Fe/H], and [α/Fe]), their corresponding uncertainties, and direct extinction measurements. Two methods were used to extract parameters and uncertainties: one that accepts the MCMC’s prescribed minimum-χ 2 result, and one that uses Bayesian inference to compute a likelihood-weighted mean from the χ 2 distribution sampled by the MCMC. Results were evaluated for internal consistency using repeat observations where available and by comparing them with external data sets (e.g., APOGEE-2 and Gaia DR2). Our spectral-fitting exercise reveals possible deficiencies in current theoretical model spectra, illustrating the potential power of MaStar spectra for helping to improve the models. This paper represents an update to the parameters that were originally presented with SDSS-IV DR17. The MaStar parameter catalog containing our BestFit results is available on the SDSS-IV DR17 website as part of version 2 of the MaStar stellar parameter value-added catalog.

79 ASTRONOMY AND ASTROPHYSICS↗

Complex network approach to fractional time series

In order to extract correlation information inherited in stochastic time series, the visibility graph algorithm has been recently proposed, by which a time series can be mapped onto a complex network. We demonstrate that the visibility algorithm is not an appropriate one to study the correlation aspects of a time series. We then employ the horizontal visibility algorithm, as a much simpler one, to map fractional processes onto complex networks. The degree distributions are shown to have parabolic exponential forms with Hurst dependent fitting parameter. Further, we take into account other topological properties such as maximum eigenvalue of the adjacency matrix and the degree assortativity, and show that such topological quantities can also be used to predict the Hurst exponent, with an exception for anti-persistent fractional Gaussian noises. To solve this problem, we take into account the Spearman correlation coefficient between nodes' degrees and their corresponding data values in the original time series.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurement of the splashback feature around SZ-selected Galaxy clusters with DES, SPT, and ACT

We present a detection of the splashback feature around galaxy clusters selected using the Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, rsp, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that rsp inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these cluster samples with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases, potentially ameliorating causes of systematic error for optically selected clusters. We find that the measured r sp for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters with similar mass and redshift distributions, r sp is ~2σ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogues and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy colour, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster.

Shin, T.↗