Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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 937 records · Page 52

Distributed Merge Trees

Improved simulations and sensors are producing datasets whose increasing complexity exhausts our ability to visualize and comprehend them directly. To cope with this problem, we can detect and extract significant features in the data and use them as the basis for subsequent analysis. Topological methods are valuable in this context because they provide robust and general feature definitions. As the growth of serial computational power has stalled, data analysis is becoming increasingly dependent on massively parallel machines. To satisfy the computational demand created by complex datasets, algorithms need to effectively utilize these computer architectures. The main strength of topological methods, their emphasis on global information, turns into an obstacle during parallelization. We present two approaches to alleviate this problem. We develop a distributed representation of the merge tree that avoids computing the global tree on a single processor and lets us parallelize subsequent queries. To account for the increasing number of cores per processor, we develop a new data structure that lets us take advantage of multiple shared-memory cores to parallelize the work on a single node. Finally, we present experiments that illustrate the strengths of our approach as well as help identify future challenges.

97 MATHEMATICS AND COMPUTING↗

Extracting Non-propagating Oscillatory Fields in Concrete to Detect Distributed Cracking

We introduce work to find and locate distributed subsurface cracks in concrete by extracting non-propagating oscillatory fields. The medium of interest is concrete, but the approach also applies to other types of inhomogeneous media. The theoretical basis of the work is first presented through a one-dimensional point-scatterer model that considers the wavefield set up by multiple distinct scatterers. More complex scattering scenarios are then investigated using numerical simulation. The numerical models consider two types of scatterers: elliptic large-scale particles distributed throughout a medium and small-sized cracks localized within a damage zone. The theoretical and numerical analyses show that forward propagating waves undergo distinct scattering behavior within the crack damaged zone: non-propagating resonance-like oscillatory fields are set up within the cracked zone that are distinct from the scatter caused by the large-scale particles. We employ frequency-wavenumber (f-k) domain analysis to extract the energy of non-propagating oscillatory fields and thus to detect and locate zones of distributed cracking. The proposed approach is evaluated using numerical simulation and experimental data collected from concrete specimens that contain simulated distributed cracks. The results reflect the location of distributed crack zones in discrete random media such as concrete can be successfully detected.

36 MATERIALS SCIENCE↗

Community-Based Transactive Coordination Mechanism for Enabling Grid Edge Systems

The changing landscape of the electricity industry, characterized by a surge in distributed energy resources (DERs) and proactive customers, necessitates practical solutions for coordinated operations especially at the distribution-level. This paper introduces a community-based transactive coordination mechanism designed to incentivize customers for providing localized and system-level services reflected through real-time prices. The work presents a bidding approach for communities with DERs, such as solar photovoltaic (PV) and battery energy storage systems (BESS), to formulate their price-responsiveness for retail energy coordination, emphasizing a community-centric model. By sending bidding curves to a third-party, the mechanism enables customers with DER assets to actively participate in localized coordination with the load serving entity (LSE), thereby supplementing each other’s and even the utilities needs through a shared energy economy. The proposed transactive mechanism is implemented leveraging a co-simulation framework that integrates a distribution grid simulator and Python-based agents for performance evaluation. Collaboration with a local utility to access real distribution feeder models and consumption profiles yields simulation results demonstrating the potential to reduce costs by 12\% for communities with DERs like PV and BESS.

Community-based coordination, grid-edge systems, R↗

Development of spatially and temporally resolved electron density measurements for the assessment of radiation hydrodynamics simulations of laboratory X-ray photoionized plasmas

The photoionized plasma gas cell experiment is an established platform we use to make at-parameter (ξ > > 1 ergs cm s –1 ) measurements of plasma properties with application to high-energy astrophysical systems. We model the experiments with 1D radiation hydrodynamics simulations using the HELIOS-CR code to inform our understanding and assist in the interpretation of results. The simulations predict that the bulk of the plasma is in a quasi-uniform and hydrodynamically unperturbed state throughout the duration of the experiment. To evaluate this prediction, we introduced a photonic Doppler velocimetry (PDV) diagnostic to measure spatially and temporally resolved plasma electron density. The initial measurements were successful but had limitations that made model-data comparisons challenging. To address this, we re-designed the gas cell PDV diagnostic and doubled the number of measurement locations to sample across two thirds of the depth of the cell. We also present a comparison of the results from the upgraded PDV diagnostic to the HELIOS-CR simulations for the first time. As a result, the experimental data confirms the prediction of an unperturbed region in the bulk of the plasma but reveals discrepancies in the time evolution and spatial distribution of the simulated electron density.

79 ASTRONOMY AND ASTROPHYSICS↗

Short-Term State Forecasting-Based Optimal Voltage Regulation in Distribution Systems: Preprint

A novel short-term state forecasting-based optimal power flow (OPF) approach for distribution system voltage regulation is proposed in this paper. An extreme learning machine (ELM) based state forecaster is developed to accurately predict system states (voltage magnitudes and angles) in the near future. Based on the forecast system states, a dynamically weighted three-phase AC OPF problem is formulated to minimize the voltage violations with higher penalization on buses which are forecast to have higher voltage violations in the near future. By solving the proposed OPF problem, the controllable resources in the system are optimally coordinated to alleviate the potential severe voltage violations and improve the overall voltage profile. The proposed approach has been tested in a 12-bus distribution system and simulation results are presented to demonstrate the performance of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

I-125 ROPES eye plaque dosimetry: Validation of a commercial 3D ophthalmic brachytherapy treatment planning system and independent dose calculation software with GafChromic{sup ®} EBT3 films

Purpose: The purpose of this study was to measure the dose distributions for different Radiation Oncology Physics and Engineering Services, Australia (ROPES) type eye plaques loaded with I-125 (model 6711) seeds using GafChromic{sup ®} EBT3 films, in order to verify the dose distributions in the Plaque Simulator™ (PS) ophthalmic 3D treatment planning system. The brachytherapy module of RADCALC{sup ®} was used to independently check the dose distributions calculated by PS. Correction factors were derived from the measured data to be used in PS to account for the effect of the stainless steel ROPES plaque backing on the 3D dose distribution.Methods: Using GafChromic{sup ®} EBT3 films inserted in a specially designed Solid Water™ eye ball phantom, dose distributions were measured three-dimensionally both along and perpendicular to I-125 (model 6711) loaded ROPES eye plaque's central axis (CAX) with 2 mm depth increments. Each measurement was performed in full scatter conditions both with and without the stainless steel plaque backing attached to the eye plaque, to assess its effect on the dose distributions. Results were compared to the dose distributions calculated by Plaque Simulator™ and checked independently with RADCALC{sup ®}.Results: The EBT3 film measurements without the stainless steel backing were found to agree with PS and RADCALC{sup ®} to within 2% and 4%, respectively, on the plaque CAX. Also, RADCALC{sup ®} was found to agree with PS to within 2%. The CAX depth doses measured using EBT3 film with the stainless steel backing were observed to result in a 4% decrease relative to when the backing was not present. Within experimental uncertainty, the 4% decrease was found to be constant with depth and independent of plaque size. Using a constant dose correction factor of T= 0.96 in PS, where the calculated dose for the full water scattering medium is reduced by 4% in every voxel in the dose grid, the effect of the plaque backing was accurately modeled in the planning system. Off-axis profiles were also modeled in PS by taking into account the three-dimensional model of the plaque backing.Conclusions: The doses calculated by PS and RADCALC{sup ®} for uniformly loaded ROPES plaques in full and uniform scattering conditions were validated by the EBT3 film measurements. The stainless steel plaque backing was observed to decrease the measured dose by 4%. Through the introduction of a scalar correction factor (0.96) in PS, the dose homogeneity effect of the stainless steel plaque backing was found to agree with the measured EBT3 film measurements.

61 RADIATION PROTECTION AND DOSIMETRY↗

Hydration dynamics in water clusters via quantum molecular dynamics simulations

We have investigated the hydration dynamics in size selected water clusters with n = 66, 104, 200, 500, and 1000 water molecules using molecular dynamics simulations. To study the most fundamental aspects of relaxation phenomena in clusters, we choose one of the simplest, still realistic, quantum mechanically treated test solute, an excess electron. The project focuses on the time evolution of the clusters following two processes, electron attachment to neutral equilibrated water clusters and electron detachment from an equilibrated water cluster anion. The relaxation dynamics is significantly different in the two processes, most notably restoring the equilibrium final state is less effective after electron attachment. Nevertheless, in both scenarios only minor cluster size dependence is observed. Significantly different relaxation patterns characterize electron detachment for interior and surface state clusters, interior state clusters relaxing significantly faster. This observation may indicate a potential way to distinguish surface state and interior state water cluster anion isomers experimentally. A comparison of equilibrium and non-equilibrium trajectories suggests that linear response theory breaks down for electron attachment at 200 K, but the results converge to reasonable agreement at higher temperatures. Relaxation following electron detachment clearly belongs to the linear regime. Cluster relaxation was also investigated using two different computational models, one preferring cavity type interior states for the excess electron in bulk water, while the other simulating non-cavity structure. While the cavity model predicts appearance of several different hydrated electron isomers in agreement with experiment, the non-cavity model locates only cluster anions with interior excess electron distribution. The present simulations show that surface isomers computed with the cavity predicting potential show similar dynamical behavior to the interior clusters of the non-cavity type model. Relaxation associated with cavity collapse presents, however, unique dynamical signatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

COSMIC REIONIZATION ON COMPUTERS. III. THE CLUMPING FACTOR

We use fully self-consistent numerical simulations of cosmic reionization, completed under the Cosmic Reionization On Computers project, to explore how well the recombinations in the ionized intergalactic medium (IGM) can be quantified by the effective “clumping factor.” The density distribution in the simulations (and, presumably, in a real universe) is highly inhomogeneous and more-or-less smoothly varying in space. However, even in highly complex and dynamic environments, the concept of the IGM remains reasonably well-defined; the largest ambiguity comes from the unvirialized regions around galaxies that are over-ionized by the local enhancement in the radiation field (“proximity zones”). That ambiguity precludes computing the IGM clumping factor to better than about 20%. We also discuss a “local clumping factor,” defined over a particular spatial scale, and quantify its scatter on a given scale and its variation as a function of scale.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark Matter Annihilation at the Galactic Center

Observations by the WMAP and PLANCK satellites have provided extraordinarily accurate observations on the densities of baryonic matter, dark matter, and dark energy in the universe. These observations indicate that our universe is composed of approximately five times as much dark matter as baryonic matter. However, efforts to detect a particle responsible for the energy density of dark matter have been unsuccessful. Theoretical models have indicated that a leading candidate for the dark matter is the lightest supersymmetric particle, which may be stable due to a conserved R-parity. This dark matter particle would still be capable of interacting with baryons via weak-force interactions in the early universe, a process which was found to naturally explain the observed relic abundance of dark matter today. These residual annihilations can persist, albeit at a much lower rate, in the present universe, providing a detectable signal from dark matter annihilation events which occur throughout the universe. Simulations calculating the distribution of dark matter in our galaxy almost universally predict the galactic center of the Milky Way Galaxy (GC) to provide the brightest signal from dark matter annihilation due to its relative proximity and large simulated dark matter density. Recent advances in telescope technology have allowed for the first multiwavelength analysis of the GC, with suitable effective exposure, angular resolution, and energy resolution in order to detect dark matter particles with properties similar to those predicted by the WIMP miracle. In this work, I describe ongoing efforts which have successfully detected an excess in gamma-ray emission from the region immediately surrounding the GC, which is difficult to describe in terms of standard diffuse emission predicted in the GC region. While the jury is still out on any dark matter interpretation of this excess, I describe several related observations which may indicate a dark matter origin. Finally, I discuss the role of future telescopes in differentiating a dark matter model from astrophysical emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring thermal profiles in high explosives using neural networks

We present a new method for calculating the temperature profile of high explosive (HE) material using a Convolutional Neural Network (CNN). To train/test the CNN, we have developed a hybrid experiment/simulation method for collecting acoustic and temperature data. We experimentally heat cylindrical containers of HE material until detonation/deflagration, where we continuously measure the acoustic bursts through the HE using multiple acoustic transducers lined around the exterior container circumference. However, measuring the temperature profile in the HE in an experiment would require inserting a large number of thermal probes, which would disrupt the heating process. Thus, we use two thermal probes, one at the HE center and one at the wall. We then use numerical simulation of the heating process to calculate the temperature distribution and correct the simulated temperatures based on the experimental center and wall temperatures. We calculate temperature errors on the order of 15 °C, which is ∼12% of the range of temperatures in the experiment. We also investigate how the algorithm’s accuracy is affected by the number of acoustic receivers used to collect each measurement and the resolution of the temperature prediction. This work provides a means of assessing the safety status of HE material, which cannot be achieved using existing temperature measurement methods. In addition, it has implications for a range of other applications where internal temperature profile measurements would provide critical information. These applications include detecting chemical reactions, observing thermodynamic processes such as combustion, monitoring metal or plastic casting, determining the energy density in thermal storage capsules, and identifying abnormal battery operations.

97 MATHEMATICS AND COMPUTING↗

Simulation of Powder Layer Deposition in Additive Manufacturing Processes Using the Discrete Element Method

This document serves as a final report to a small effort where several improvements were added to a LLNL code GEODYN-­L to develop Discrete Element Method (DEM) algorithms coupled to Lagrangian Finite Element (FE) solvers to investigate powder-­bed formation problems for additive manufacturing. The results from these simulations will be assessed for inclusion as the initial conditions for Direct Metal Laser Sintering (DMLS) simulations performed with ALE3D. The algorithms were written and performed on parallel computing platforms at LLNL. The total funding level was 3-­4 weeks of an FTE split amongst two staff scientists and one post-­doc. The DEM simulations emulated, as much as was feasible, the physical process of depositing a new layer of powder over a bed of existing powder. The DEM simulations utilized truncated size distributions spanning realistic size ranges with a size distribution profile consistent with realistic sample set. A minimum simulation sample size on the order of 40-­particles square by 10-­particles deep was utilized in these scoping studies in order to evaluate the potential effects of size segregation variation with distance displaced in front of a screed blade. A reasonable method for evaluating the problem was developed and validated. Several simulations were performed to show the viability of the approach. Future investigations will focus on running various simulations investigating powder particle sizing and screen geometries.

36 MATERIALS SCIENCE↗

Measuring thermal profiles in high explosives using neural networks

We present a new method for calculating the temperature profile of high explosive (HE) material using a Convolutional Neural Network (CNN). To train/test the CNN, we have developed a hybrid experiment/simulation method for collecting acoustic and temperature data. We experimentally heat cylindrical containers of HE material until detonation/deflagration, where we continuously measure the acoustic bursts through the HE using multiple acoustic transducers lined around the exterior container circumference. However, measuring the temperature profile in the HE in an experiment would require inserting a large number of thermal probes, which would disrupt the heating process. Thus, we use two thermal probes, one at the HE center and one at the wall. We then use numerical simulation of the heating process to calculate the temperature distribution and correct the simulated temperatures based on the experimental center and wall temperatures. We calculate temperature errors on the order of 15 °C, which is ~12% of the range of temperatures in the experiment. We also investigate how the algorithm’s accuracy is affected by the number of acoustic receivers used to collect each measurement and the resolution of the temperature prediction. This work provides a means of assessing the safety status of HE material, which cannot be achieved using existing temperature measurement methods. In addition, it has implications for a range of other applications where internal temperature profile measurements would provide critical information. These applications include detecting chemical reactions, observing thermodynamic processes such as combustion, monitoring metal or plastic casting, determining the energy density in thermal storage capsules, and identifying abnormal battery operations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superresolution Full-polarimetric Imaging for Radio Interferometry with Sparse Modeling

We propose a new technique for radio interferometry to obtain superresolution full-polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the ℓ {sub 1} norm and the total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope. We confirm that ℓ {sub 1}+TV regularization can achieve an optimal resolution of ~25%–30% of the diffraction limit λ/D{sub max}, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e., Stokes I) and linear polarizations (i.e., Stokes Q and U). This optimal resolution is better than that obtained from the widely used Cotton–Schwab CLEAN algorithm or from using ℓ {sub 1} or TV regularizations alone. Furthermore, we find that ℓ {sub 1}+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the superresolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectral Effects in Albedo and Rearside Irradiance Measurment for Bifacial Performance Estimation: Preprint

Albedo data are essential for accurate prediction of bifacial PV module performance. However, spectral response mismatch between PV modules and irradiance sensors used in albedometers can limit the accuracy of performance predictions. In order to provide quantitative assessment of this effect, we investigate via simulation the differences in spectrally responsive albedo measured with thermopile pyranometers and crystallinesilicon PV reference cells in comparison to a representative crystalline-silicon bifacial PV module for nine different representative ground surface materials. Calculations are performed using simulated solar spectra together with catalogued spectral reflectivity data distributed with the SMARTS simulation software. For the specific materials considered, the results show that albedo measurement using thermopile pyranometers could over or under-estimate the ground-reflected radiation usable by a bifacial PV module by up to 10%, versus only approximately 4% total range of variation for a PV reference cell.

Albedo↗

Output power fluctuations due to different weights of macro particles used in particle-in-cell simulations of Cerenkov devices

The output power fluctuations caused by weights of macro particles used in particle-in-cell (PIC) simulations of a backward wave oscillator and a travelling wave tube are statistically analyzed. It is found that the velocities of electrons passed a specific slow-wave structure form a specific electron velocity distribution. The electron velocity distribution obtained in PIC simulation with a relative small weight of macro particles is considered as an initial distribution. By analyzing this initial distribution with a statistical method, the estimations of the output power fluctuations caused by different weights of macro particles are obtained. The statistical method is verified by comparing the estimations with the simulation results. The fluctuations become stronger with increasing weight of macro particles, which can also be determined reversely from estimations of the output power fluctuations. With the weights of macro particles optimized by the statistical method, the output power fluctuations in PIC simulations are relatively small and acceptable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy Storage Sizing, Placement and Valuation for Distribution Networks

SAND2023-05682O The Energy Storage Sizing, Placement, and Valuation for Distribution software allows a user to simulate sizing, placement, and operation of one or more storage assets on a distribution grid. Among its features are measuring differences to determine the best combinations, and options for configuring and comparing results. This new simulation capability leverages the Open Distribution System Simulator and the Hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS). Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Tamrakar, Ujjwol↗

Density relaxation and particle motion characteristics in a non-ionic deep eutectic solvent (acetamide + urea): Time-resolved fluorescence measurements and all-atom molecular dynamics simulations

Temperature dependent relaxation dynamics, particle motion characteristics, and heterogeneity aspects of deep eutectic solvents (DESs) made of acetamide (CH{sub 3}CONH{sub 2}) and urea (NH{sub 2}CONH{sub 2}) have been investigated by employing time-resolved fluorescence measurements and all-atom molecular dynamics simulations. Three different compositions (f) for the mixture [fCH{sub 3}CONH{sub 2} + (1 − f)NH{sub 2}CONH{sub 2}] have been studied in a temperature range of 328-353 K which is ∼120-145 K above the measured glass transition temperatures (∼207 K) of these DESs but much lower than the individual melting temperature of either of the constituents. Steady state fluorescence emission measurements using probe solutes with sharply different lifetimes do not indicate any dependence on excitation wavelength in these metastable molten systems. Time-resolved fluorescence anisotropy measurements reveal near-hydrodynamic coupling between medium viscosity and rotation of a dissolved dipolar solute. Stokes shift dynamics have been found to be too fast to be detected by the time-resolution (∼70 ps) employed, suggesting extremely rapid medium polarization relaxation. All-atom simulations reveal Gaussian distribution for particle displacements and van Hove correlations, and significant overlap between non-Gaussian (α{sub 2}) and new non-Gaussian (γ) heterogeneity parameters. In addition, no stretched exponential relaxations have been detected in the simulated wavenumber dependent acetamide dynamic structure factors. All these results are in sharp contrast to earlier observations for ionic deep eutectics with acetamide [Guchhait et al., J. Chem. Phys. 140, 104514 (2014)] and suggest a fundamental difference in interaction and dynamics between ionic and non-ionic deep eutectic solvent systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing climate-change-induced flood risk in the Conasauga River watershed: an application of ensemble hydrodynamic inundation modeling

Abstract. This study evaluates the impact of potential future climate change on flood regimes, floodplain protection, and electricity infrastructures across the Conasauga River watershed in the southeastern United States through ensemble hydrodynamic inundation modeling. The ensemble streamflow scenarios were simulated by the Distributed Hydrology Soil Vegetation Model (DHSVM) driven by (1) 1981–2012 Daymet meteorological observations and (2) 11 sets of downscaled global climate models (GCMs) during the 1966–2005 historical and 2011–2050 future periods. Surface inundation was simulated using a GPU-accelerated Two-dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) hydrodynamic model. A total of 9 out of the 11 GCMs exhibit an increase in the mean ensemble flood inundation areas. Moreover, at the 1 % annual exceedance probability level, the flood inundation frequency curves indicate a ∼ 16 km2 increase in floodplain area. The assessment also shows that even after flood-proofing, four of the substations could still be affected in the projected future period. The increase in floodplain area and substation vulnerability highlights the need to account for climate change in floodplain management. Overall, this study provides a proof-of-concept demonstration of how the computationally intensive hydrodynamic inundation modeling can be used to enhance flood frequency maps and vulnerability assessment under the changing climatic conditions.

54 ENVIRONMENTAL SCIENCES↗