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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantitative Power System Resilience Metrics and Evaluation Approach: Preprint

Power system resilience is an emerging topic and plays an essential role in helping power industry understand and respond to the increasing threats of extreme weather events. The first step of power system resilience analysis is to introduce metrics to quantify the resilience reasonably. Existing resilience metrics are typically restrained by the limited data for extreme event modeling and fall short in terms of physical interpretation and comparability. This paper develops novel quantitative metrics to evaluate power system resilience in pre- and post-event contexts. The developed metrics illustrate clear physical meanings and can be effectively used to compare resilience across different systems under different extreme events. Moreover, the developed metrics can be applied to both transmission and distribution systems. Simulation on a distribution system is employed to validate the effectiveness of the proposed resilience metrics and resilience evaluation approach.

power system resilience↗
DOE OSTI2025

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗
DOE OSTI2015

Comparison of film measurements and Monte Carlo simulations of dose delivered with very high-energy electron beams in a polystyrene phantom

Purpose: To measure radiation dose in a water-equivalent medium from very high-energy electron (VHEE) beams and make comparisons to Monte Carlo (MC) simulation results. Methods: Dose in a polystyrene phantom delivered by an experimental VHEE beam line was measured with Gafchromic films for three 50 MeV and two 70 MeV Gaussian beams of 4.0–6.9 mm FWHM and compared to corresponding MC-simulated dose distributions. MC dose in the polystyrene phantom was calculated with the EGSnrc/BEAMnrc and DOSXYZnrc codes based on the experimental setup. Additionally, the effect of 2% beam energy measurement uncertainty and possible non-zero beam angular spread on MC dose distributions was evaluated. Results: MC simulated percentage depth dose (PDD) curves agreed with measurements within 4% for all beam sizes at both 50 and 70 MeV VHEE beams. Central axis PDD at 8 cm depth ranged from 14% to 19% for the 5.4–6.9 mm 50 MeV beams and it ranged from 14% to 18% for the 4.0–4.5 mm 70 MeV beams. MC simulated relative beam profiles of regularly shaped Gaussian beams evaluated at depths of 0.64 to 7.46 cm agreed with measurements to within 5%. A 2% beam energy uncertainty and 0.286° beam angular spread corresponded to a maximum 3.0% and 3.8% difference in depth dose curves of the 50 and 70 MeV electron beams, respectively. Absolute dose differences between MC simulations and film measurements of regularly shaped Gaussian beams were between 10% and 42%. Conclusions: The authors demonstrate that relative dose distributions for VHEE beams of 50–70 MeV can be measured with Gafchromic films and modeled with Monte Carlo simulations to an accuracy of 5%. The reported absolute dose differences likely caused by imperfect beam steering and subsequent charge loss revealed the importance of accurate VHEE beam control and diagnostics.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Investigating Characteristic Droplet Size Distributions in Large Eddy Simulations of Stratocumulus Clouds

Cloud processes relevant to radiative and precipitation properties depend on the shape of the cloud droplet size distribution. Recent holographic observations revealed that cloud droplet populations do not have the same size distribution shapes throughout but form regions of characteristic distributions with similar microphysical properties. We investigate the existence and properties of these characteristic distributions within Large‐Eddy Simulations of stratocumulus clouds using Lagrangian and bin microphysics schemes. Distribution types are identified, revealing localized characteristic distributions that vary on the scale of the largest convective cell for simulations with bin microphysics. The results from the Lagrangian microphysics scheme hint at similar behavior. Compared to observations, the simulated clouds are much more uniform. Analysis of the LES results suggests a connection to the local entrainment rate, so the poorly resolved entrainment interface in LES may be a cause of the uniformity. The uniformity of the large‐scale forcing could also be a factor.

cloud droplet size distributions↗
DOE OSTI2016

ASSESSING ASTROPHYSICAL UNCERTAINTIES IN DIRECT DETECTION WITH GALAXY SIMULATIONS

We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions that are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Furthermore, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulations. We show that the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way’s estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models that are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Penumbra: A spatially distributed, mechanistic model for simulating ground-level incident solar energy across heterogeneous landscapes

Landscape solar energy is a significant environmental driver, yet it remains complicated to model well. Several solar radiation models simplify the complexity of light by estimating it at discrete point locations or by averaging values over larger areas. These modeling approaches may be useful in certain cases, but they are unable to provide spatially distributed and temporally dynamic representations of solar energy across entire landscapes. We created a landscape-scale ground-level shade and solar energy model called Penumbra to address this deficiency. Penumbra simulates spatially distributed ground-level shade and incident solar energy at user-defined timescales by modeling local and distant topographic shading and vegetative shading. Spatially resolved inputs of a digital elevation model, a normalized digital surface model, and landscape object transmittance are used to estimate spatial variations in solar energy at user-defined temporal timesteps. The research goals for Penumbra included: 1) simulations of spatiotemporal variations of shade and solar energy caused by both objects and topographic features, 2) minimal user burden and parameterization, 3) flexible user defined temporal parameters, and 4) flexible external model coupling. We test Penumbra’s predictive skill by comparing the model’s predictions with monitored open and forested sites, and achieve calibrated mean errors ranging from -17.3 to 148.1 μmoles/m 2 /s. Penumbra is a dynamic model that can produce spatial and temporal representations of shade percentage and ground-level solar energy. Outputs from Penumbra can be used with other ecological models to better understand the health and resilience of aquatic, near stream terrestrial, and upland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2011

Margin Evaluation in the Presence of Deformation, Rotation, and Translation in Prostate and Entire Seminal Vesicle Irradiation With Daily Marker-Based Setup Corrections

Purpose: To develop a method for margin evaluation accounting for all measured displacements during treatment of prostate cancer. Methods and Materials: For 21 patients treated with stereographic targeting marker-based online translation corrections, dose distributions with varying margins and gradients were created. Sets of possible cumulative delivered dose distributions were simulated by moving voxels and accumulating dose per voxel. Voxel motion was simulated consistent with measured distributions of systematic and random displacements due to stereographic targeting inaccuracies, deformation, rotation, and intrafraction motion. The method of simulation maintained measured correlation of voxel motions due to organ deformation. Results: For the clinical target volume including prostate and seminal vesicles (SV), the probability that some part receives <95% of the prescribed dose, the changes in minimum dose, and volume receiving 95% of prescription dose compared with planning were 80.5% {+-} 19.2%, 9.0 {+-} 6.8 Gy, and 3.0% {+-} 3.7%, respectively, for the smallest studied margins (3 mm prostate, 5 mm SV) and steepest dose gradients. Corresponding values for largest margins (5 mm prostate, 8 mm SV) with a clinical intensity-modulated radiotherapy dose distribution were 46.5% {+-} 34.7%, 6.7 {+-} 5.8 Gy, and 1.6% {+-} 2.3%. For prostate-only clinical target volume, the values were 51.8% {+-} 17.7%, 3.3 {+-} 1.6 Gy, and 0.6% {+-} 0.5% with the smallest margins and 5.2% {+-} 7.4%, 1.8 {+-} 0.9 Gy, and 0.1% {+-} 0.1% for the largest margins. Addition of three-dimensional rotation corrections only improved these values slightly. All rectal planning constraints were met in the actual reconstructed doses for all studied margins. Conclusion: We developed a system for margin validation in the presence of deformations. In our population, a 5-mm margin provided sufficient dosimetric coverage for the prostate. In contrast, an 8-mm SV margin was still insufficient owing to deformations. Addition of three-dimensional rotation corrections was of minor influence.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Study of In-Trap Ion Clouds by Ion Trajectory Simulations

Gaussian distribution has been utilized to describe the global number density distribution of ion cloud in the Paul trap, which is known as the thermal equilibrium theory and widely used in theoretical modeling of ion clouds in the ion traps. Using ion trajectory simulations, however, the ion clouds can now also be treated as a dynamic ion flow field and the location-dependent features could now be characterized. This study was carried out to better understand the in-trap ion cloud properties, such as the local particle velocity and temperature. The local ion number densities were found to be heterogeneously distributed in terms of mean and distribution width; the velocity and temperature of the ion flow varied with pressure depending on the flow type of the neutral molecules; and the “quasi-static” equilibrium status can only be achieved after a certain number of collisions, for which the time period is pressure-dependent. This work provides new insights of the ion clouds that are globally stable but subjected to local rf heating and collisional cooling. .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Impact of DER Communication Delay in AGC: Cyber-Physical Dynamic Co-Simulation: Preprint

Distributed energy resource (DER) frequency regulations are promising technologies for future grid operation. Unlike conventional generators, DERs might require open communication networks to exchange signals with control centers, possibly through DER aggregators; therefore, the impacts of the communication variations including the latency on the system frequency stability need to be investigated. This paper develops a cyber-physical dynamic simulation (CPDS) model based on the Hierarchical Engine for Large-Scale Co-Simulation (HELICS) to evaluate the impact of the communication variations, such as delays in DER frequency regulations. The feasible delay range can be obtained under different AGC parameter settings. The results show that the risk of instability generally increases with the communication delay and that the difference of AGC provided from DER and conventional generation units should be recognized when designing frequency regulation controllers.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 2: Stratospheric and free-tropospheric response

The paper constitutes Part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric aerosol injection (SAI) at various single latitudes in the tropics, as simulated by three state-of-the-art Earth system models – CESM2-WACCM6, UKESM1.0, and GISS-E2.1-G. Building on Part 1 (Visioni et al., 2023) we demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics in driving the inter-model differences in the simulated sulfate distributions. We then characterize the simulated changes in stratospheric and free-tropospheric temperatures, ozone, water vapor, and large-scale circulation, elucidating the role of the above aspects in the surface SAI responses discussed in Part 1. We show that the differences in the aerosol spatial distribution can be explained by the significantly faster shallow branches of the Brewer–Dobson circulation in CESM2, a relatively isolated tropical pipe and older tropical age of air in UKESM, and smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age of air) in the two GISS versions used. We also find a large spread in the magnitudes of the tropical lower-stratospheric warming amongst the models, driven by microphysical, chemical, and dynamical differences. These lead to large differences in stratospheric water vapor responses, with significant increases in stratospheric water vapor under SAI in CESM2 and GISS that were largely not reproduced in UKESM. For ozone, good agreement was found in the tropical stratosphere amongst the models with more complex microphysics, with lower stratospheric ozone changes consistent with the SAI-induced modulation of the large-scale circulation and the resulting changes in transport. In contrast, we find a large inter-model spread in the Antarctic ozone responses that can largely be explained by the differences in the simulated latitudinal distributions of aerosols as well as the degree of implementation of heterogeneous halogen chemistry on sulfate in the models. The use of GISS runs with bulk microphysics demonstrates the importance of more detailed treatment of aerosol processes, with contrastingly different stratospheric SAI responses to the models using the two-moment aerosol treatment; however, some problems in halogen chemistry in GISS are also identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as identifying and narrowing the uncertainty in model projections of climate impacts from SAI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Simulation Model of Microsphere Distribution for Selective Internal Radiation Therapy Agrees With Observations

Purpose: To perform a detailed analysis of microsphere distribution in biopsy material from a patient treated with {sup 90}Y-labeled resin spheres and characterize microsphere distribution in the hepatic artery tree, and to construct a novel dichotomous bifurcation model for microsphere deposits and evaluate its accuracy in simulating the observed microsphere deposits. Methods and Materials: Our virtual model consisted of arteries that successively branched into 2 new generations of arteries at 20 nodes. The artery diameter exponentially decreased from the lowest generation to the highest generation. Three variable parameters were optimized to obtain concordance between simulations and measure microsphere distributions: an artery coefficient of variation (ACV) for the diameter of all artery generations and the microsphere flow distribution at the nodes; a hepatic tree distribution volume (HDV) for the artery tree; and an artery diameter reduction (ADR) parameter. The model was tested against previously measured activity concentrations in 84 biopsies from the liver of 1 patient. In 16 of 84 biopsies, the microsphere distribution regarding cluster size and localization in the artery tree was determined via light microscopy of 30-μm sections (mean concentration, 14 microspheres/mg; distributions divided into 3 groups with mean microsphere concentrations of 4.6, 14, and 28 microspheres/mg). Results: Single spheres and small clusters were observed in terminal arterioles, whereas large clusters, up to 450 microspheres, were observed in larger arterioles. For 14 microspheres/mg, the optimized parameter values were ACV=0.35, HDV = 50 cm{sup 3}, and ADR=6 μm. For 4.6 microspheres/mg, ACV and ADR decreased to 0.26 and 0 μm, respectively, whereas HDV increased to 130 cm{sup 3}. The opposite trend was observed for 28 microspheres/mg: ACV = 0.49, HDV = 20 cm{sup 3}, and ADR = 8 μm. Conclusion: Simulations and measurements reveal that microsphere clusters are larger and more common in volumes with high microsphere concentrations and indicate that the spatial distribution of the artery tree must be considered in estimates of microsphere distributions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Parallel Simulation of Quantum Networks with Distributed Quantum State Management

Quantum network simulators offer the opportunity to cost-efficiently investigate potential avenues for building networks that scale with the number of users, communication distance, and application demands by simulating alternative hardware designs and control protocols. Several quantum network simulators have been recently developed with these goals in mind. As the size of the simulated networks increases, however, sequential execution becomes time-consuming. Parallel execution presents a suitable method for scalable simulations of large-scale quantum networks, but the unique attributes of quantum information create unexpected challenges. In this work, we identify requirements for parallel simulation of quantum networks and develop the first parallel discrete-event quantum network simulator by modifying the existing serial simulator SeQUeNCe. Our contributions include the design and development of a quantum state manager (QSM) that maintains shared quantum information distributed across multiple processes. We also optimize our parallel code by minimizing the overhead of the QSM and decreasing the amount of synchronization needed among processes. Using these techniques, we observe a speedup of 2 to 25 times when simulating a 1,024-node linear network topology using 2 to 128 processes. We also observe an efficiency greater than 0.5 for up to 32 processes in a linear network topology of the same size and with the same workload. We repeat this evaluation with a randomized workload on a caveman network. We also introduce several methods for partitioning networks by mapping them to different parallel simulation processes. We have released the parallel SeQUeNCe simulator as an open source tool alongside the existing sequential version.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

On electron bunching and stratification of glow discharges

Plasma stratification and excitation of ionization waves is one of the fundamental problems in gas discharge physics. Significant progress in this field is associated with the name of Lev Tsendin. He advocated the need for the kinetic approach to this problem contrary to the traditional hydrodynamic approach, introduced the idea of electron bunching in spatially periodic electric fields, and developed a theory of kinetic resonances for analysis of moving striations in rare gases. The present paper shows how Tsendin's ideas have been further developed and applied for understanding the nature of the well-known S-, P-, and R-striations observed in glow discharges of inert gases at low pressures and currents. We review numerical solutions of a Fokker-Planck kinetic equation in spatially periodic electric fields under the effects of elastic and inelastic collisions of electrons with atoms. We illustrate the formation of kinetic resonances at specific field periods for different shapes of injected Electron Distribution Functions (EDF). Computer simulations illustrate how self-organization of the EDFs occurs under nonlocal conditions and how Gaussian-like peaks moving along resonance trajectories are formed in a certain range of discharge conditions. The calculated EDFs agree well with the experimentally measured EDFs for the S, P, and R striations in noble gases. We discuss how kinetic resonances affect dispersion characteristics of moving striations and mention some non-linear effects associated with glow discharge stratification. We propose further studies of stratification phenomena combining physical kinetics and non-linear physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Diffraction pattern simulation of cellulose fibrils using distributed and quantized pair distances

Intensity simulation of X-ray scattering from large twisted cellulose molecular fibrils is important in understanding the impact of chemical or physical treatments on structural properties such as twisting or coiling. This paper describes a highly efficient method for the simulation of X-ray diffraction patterns from complex fibrils using atom-type-specific pair-distance quantization. Pair distances are sorted into arrays which are labelled by atom type. Histograms of pair distances in each array are computed and binned and the resulting population distributions are used to represent the whole pair-distance data set. These quantized pair-distance arrays are used with a modified and vectorized Debye formula to simulate diffraction patterns. This approach utilizes fewer pair distances in each iteration, and atomic scattering factors are moved outside the iteration since the arrays are labelled by atom type. As a result, this algorithm significantly reduces the computation time while maintaining the accuracy of diffraction pattern simulation, making possible the simulation of diffraction patterns from large twisted fibrils in a relatively short period of time, as is required for model testing and refinement.

09 BIOMASS FUELS↗
DOE OSTI2018

An Assessment of GCM Performance at a Regional Scale Using a Score-Based Method

A multicriteria score-based method was developed to assess the performances of 18 general circulation models (GCMs) in the study region from 1970 to 2005. The results indicate the following. (1) GCMs simulate temperature better than rainfall. The temporal and spatial distributions of simulated temperature performed well compared with those from the observations. In comparison to temperature, the spatial distribution of simulated precipitation performed poorly. Most of the GCMs underestimated temperature and overestimated precipitation. (2) The Grubbs test was used to detect anomalous moving changes in the rank score (RS) results; the inm-cm4 and ipsl-cm5b-lr models were rejected when simulating temperature, while the bnu-esm and canesm2 models performed poorly when simulating precipitation. (3) Adding or removing any criterion does not significantly influence the RS result, which indicates that the multicriteria score-based method is robust. The advantages of using multicriteria score-based method to assess GCMs performance were demonstrated, and this method also provides a more comprehensive assessment when compared with the single-criterion method. The multicriteria method could replace other criteria as the research requirements and could be easily extended to different study regions; the results could be used for better informed regional climate change impact analyses.

Shi, Fangxin↗
DOE OSTI2014

An improved Monte Carlo (MC) dose simulation for charged particle cancer therapy

Heavy-particle therapy such as carbon ion therapy are more popular nowadays because of the nature characteristics of charged particle and almost no side effect to patients. An effective treatment is achieved with high precision of dose calculation, in this research work, Geant4 based Monte Carlo simulation method has been used to calculate the radiation transport and dose distribution. The simulation have the same setting with the treatment room in Heavy Ion Medical Accelerator, HIMAC. The carbon ion beam at the isocentric gantry nozzle for the therapeutic energy of 290 MeV/u was simulated, experimental work was carried out in National Institute of Radiological Sciences, NIRS, Chiba, Japan by using the HIMAC to confirm the accuracy and qualities dose distribution by MC methods. The Geant4 based simulated dose distribution were verified with measurements for Bragg peak and spread out Bragg peak (SOBP) respectively. The verification of results shows that the Bragg peak depth-dose and SOBP distributions in simulation has good agreement with measurements. In overall, the study showed that Geant4 based can be fully applied in the heavy-ion therapy field for simulation, further works need to be carry on to refine and improve the Geant4 MC simulations.

61 RADIATION PROTECTION AND DOSIMETRY↗