Search NASA⌕ Search

SEARCH · Search NASA

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 73 records · Page 4

Real-Time Coupling of Geographically Distributed Research Infrastructures: Taxonomy, Overview, and Real-World Smart Grid Applications

Novel concepts enabling a resilient future power system and their subsequent experimental evaluation are experiencing a steadily growing challenge: large scale complexity and questionable scalability. The requirements on a research infrastructure (RI) to cope with the trends of such a dynamic system therefore grow in size, diversity and costs, making the feasibility of rigorous advancements questionable by a single RI. Analysis of large scale system complexity has been made possible by the real-time coupling of geographically separated RIs undertaking geographically distributed simulations (GDS), the concept of which brings the equipment, models and expertise of independent RIs, in combination, to optimally address the challenge. This article presents the outputs of IEEE PES Task Force on Interfacing Techniques for Simulation Tools towards standardization of GDS as a concept. First, the taxonomy for setups utilized for GDS is established followed by a comprehensive overview of the advancements in real-time couplings reported in literature. The overview encompasses fundamental technological design considerations for GDS. The article further presents four application oriented case studies (real-world implementations) where GDS setups have been utilized, demonstrating their practicality and potential in enabling the analysis of future complex power systems.

distributed laboratories↗

Verification of Dose Distribution in Carbon Ion Radiation Therapy for Stage I Lung Cancer

Purpose: To evaluate robustness of dose distribution of carbon-ion radiation therapy (C-ion RT) in non-small cell lung cancer (NSCLC) and to identify factors affecting the dose distribution by simulated dose distribution. Methods and Materials: Eighty irradiation fields for delivery of C-ion RT were analyzed in 20 patients with stage I NSCLC. Computed tomography images were obtained twice before treatment initiation. Simulated dose distribution was reconstructed on computed tomography for confirmation under the same settings as actual treatment with respiratory gating and bony structure matching. Dose-volume histogram parameters, such as %D95 (percentage of D95 relative to the prescribed dose), were calculated. Patients with any field for which the %D95 of gross tumor volume (GTV) was below 90% were classified as unacceptable for treatment, and the optimal target margin for such cases was examined. Results: Five patients with a total of 8 fields (10% of total number of fields analyzed) were classified as unacceptable according to %D95 of GTV, although most patients showed no remarkable change in the dose-volume histogram parameters. Receiver operating characteristic curve analysis showed that tumor displacement and change in water-equivalent pathlength were significant predictive factors of unacceptable cases (P<.001 and P=.002, respectively). The main cause of degradation of the dose distribution was tumor displacement in 7 of the 8 unacceptable fields. A 6-mm planning target volume margin ensured a GTV %D95 of >90%, except in 1 extremely unacceptable field. Conclusions: According to this simulation analysis of C-ion RT for stage I NSCLC, a few fields were reported as unacceptable and required resetting of body position and reconfirmation. In addition, tumor displacement and change in water-equivalent pathlength (bone shift and/or chest wall thickness) were identified as factors influencing the robustness of dose distribution. Such uncertainties should be regarded in planning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Understanding processes that control dust spatial distributions with global climate models and satellite observations

Dust aerosol is important in modulating the climate system at local and global scales, yet its spatiotemporal distributions simulated by global climate models (GCMs) are highly uncertain. In this study, we evaluate the spatiotemporal variations of dust extinction profiles and dust optical depth (DOD) simulated by the Community Earth System Model version 1 (CESM1) and version 2 (CESM2), the Energy Exascale Earth System Model version 1 (E3SMv1), and the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) against satellite retrievals from Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), Moderate Resolution Imaging Spectroradiometer (MODIS), and Multi-angle Imaging SpectroRadiometer (MISR). We find that CESM1, CESM2, and E3SMv1 underestimate dust transport to remote regions. E3SMv1 performs better than CESM1 and CESM2 in simulating dust transport and the northern hemispheric DOD due to its higher mass fraction of fine dust. CESM2 performs the worst in the Northern Hemisphere due to its lower dust emission than in the other two models but has a better dust simulation over the Southern Ocean due to the overestimation of dust emission in the Southern Hemisphere. DOD from MERRA-2 agrees well with CALIOP DOD in remote regions due to its higher mass fraction of fine dust and the assimilation of aerosol optical depth. The large disagreements in the dust extinction profiles and DOD among CALIOP, MODIS, and MISR retrievals make the model evaluation of dust spatial distributions challenging. Our study indicates the importance of representing dust emission, dry/wet deposition, and size distribution in GCMs in correctly simulating dust spatiotemporal distributions.

54 ENVIRONMENTAL SCIENCES↗

A dynamic parameterization of sulfuric acid–dimethylamine nucleation and its application in three-dimensional modeling

Sulfuric acid (SA) is a governing gaseous precursor for atmospheric new particle formation (NPF), a major source of global ultrafine particles, in environments studied around the world. In polluted urban atmospheres with high condensation sinks (CSs), the formation of stable SA–amine clusters, such as SA–dimethylamine (DMA) clusters, usually initializes intense NPF events. Coagulation scavenging and cluster evaporation are dominant sink processes of SA–amine clusters in urban atmospheres, yet these loss processes are not quantitatively included in the present parameterizations of SA–amine nucleation. We herein report a parameterization of SA–DMA nucleation, based on cluster dynamic simulations and quantum chemistry calculations, with certain simplifications to greatly reduce the computational costs. Compared with previous SA–DMA nucleation parameterizations, this new parameterization was able to reproduce the dependences of particle formation rates on temperature and CSs. We then incorporated it in a three-dimensional (3-D) chemical transport model to simulate the evolution of the particle number size distributions. Simulation results showed good consistency with the observations in the occurrence of NPF events and particle number size distributions in wintertime Beijing and represented a significant improvement compared to that using a parameterization without coagulation scavenging. Quantitative analysis shows that SA–DMA nucleation contributes significantly to nucleation rates and aerosol population during the 3-D simulations in Beijing (>99 % and >60 %, respectively). These results broaden the understanding of NPF in urban atmospheres and stress the necessity of including the effects of coagulation scavenging and cluster stability in simulating SA–DMA nucleation in 3-D simulations. Representing these processes is thus likely to improve model performance in particle source apportionment and quantification of aerosol effects on air quality, human health, and climate.

Li, Yuyang↗

Ring Asymmetry and Spin in M87*

Event Horizon Telescope (EHT) images of the supermassive black hole M87* depict an asymmetric ring of emission. General relativistic magnetohydrodynamic (GRMHD) models of M87* and its accretion disk predict that the amplitude and location of the ring’s peak brightness asymmetry should fluctuate due to turbulence in the source plasma. We compare the observed distribution of brightness asymmetry amplitudes to the simulated distribution in GRMHD models, across varying black hole spin a * . We show that, for strongly magnetized (MAD) models, three epochs of EHT data marginally disfavor ∣a * ∣ ≲ 0.2. This is consistent with the Blandford–Znajek model for M87’s jet, which predicts that M87* should have nonzero spin. We show quantitatively how future observations could improve spin constraints and discuss how improved spin constraints could distinguish between differing jet-launching mechanisms and black hole growth scenarios.

Bernshteyn, Vadim [Univ. of Illinois at Urbana-Cha↗

Dynamics modeling of molten salt reactor with reduced and expanded representations of delayed neutron precursors

Molten salt reactors (MSRs) present unique challenges in dynamic behavior due to the mobility of their fuel. In these reactors, delayed neutron precursors (DNPs) drift with the fuel circulation through the primary loop. As a result, a fraction of DNPs decays outside the core, effectively reducing the available delayed neutron population for reactivity control. Consequently, precise modeling of the distribution and behavior of DNPs is critical for accurate reactor dynamics simulations. In this study, the System Dynamics Analysis Tool (SDAT) was used to simulate a thermal-spectrum MSR under steady-state conditions and following transients. The effects of using reduced and expanded representations of DNPs with fewer or more groups than the conventional 6-group model were investigated. Their impact on the simulated distribution of precursors in the primary loop, reactivity loss value, and reactor response to transients was analyzed. Simulation results showed that reduced models lead to the loss of the actual DNPs distribution data, resulting in less accurate estimates of reactivity loss. Reactor power predictions using these reduced models showed significant deviations compared to those using the conventional 6-group model in transient simulations. Expanded models offered a more accurate representation of the distribution of DNPs and reactivity loss estimates. Reactor power predictions using expanded models showed minimal deviation from the conventional 6-group model during the simulated transients.

analysis↗

Quantitative Power System Resilience Metrics and Evaluation Approach

Power system resilience is an emerging topic and plays an essential role in helping the 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↗

A Scalable Multi-Modal Framework for High-Fidelity Distributed Human Mobility Simulations

The development of data-driven models for human mobility in urban settings requires access to substantial and diverse real-world data. However, existing historical data often presents challenges such as limited volume, variety, and veracity, as well as missing data and privacy preservation concerns. Also, urban mobility modeling is inherently time-variant, complex, and multi-modal, encompassing everything from individual walking and running to private road travel and large-scale public transportation. These challenges call for innovative solutions to overcome data limitations and compute needs to model mobility behaviors accurately. To address these challenges, we propose a distributed, co-simulation-based architecture DURMOSim that integrates real-world data with scalable, high-fidelity simulations, demonstrating distributed co-simulation feasibility with existing mobility models. DURMOSim underpins a modular integration that would enable using any available mobility simulators for greater extensibility and scalability in performing various urban scenarios. In this paper, we present the design, implementation, and performance evaluation of DURMOSim, highlighting its capability to model population-scale mobility patterns. Our initial results show its ability to dynamically synchronize multiple simulation models at runtime with negligible computational overhead. We believe DURMOSim could be a robust tool for advancing urban mobility research and intelligent transportation systems.

Yoginath, Srikanth [ORNL] (ORCID:0000000184236050)↗

`SkyPy`: A package for modelling the Universe

SkyPy is an open-source Python package for simulating the astrophysical sky. It comprises a library of physical and empirical models across a range of observables and a command-line script to run end-to-end simulations. The library provides functions that sample realisations of sources and their associated properties from probability distributions. Simulation pipelines are constructed from these models using a YAML-based configuration syntax, while task scheduling and data dependencies are handled internally and the modular design allows users to interface with external software. SkyPy is developed and maintained by a diverse community of domain experts with a focus on software sustainability and interoperability. By fostering development, it provides a framework for correlated simulations of a range of cosmological probes including galaxy populations, large scale structure, the cosmic microwave background, supernovae and gravitational waves. Version 0.4 implements functions that model various properties of galaxies including luminosity functions, redshift distributions and optical photometry from spectral energy distribution templates. Future releases will provide additional modules, for example, to simulate populations of dark matter halos and model the galaxy-halo connection, making use of existing software packages from the astrophysics community where appropriate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Bimodality in simulated precipitation frequency distributions and its relationship with convective parameterizations

Abstract Bimodality in precipitation frequency distributions is often evident in atmospheric models, but rarely in observations. This study i) proposes a metric to objectively quantify the bimodality in precipitation distributions, ii) evaluates model simulations contributed to the Coupled Model Intercomparison Project (CMIP) phase 5 (CMIP5), phase 6 (CMIP6), and the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project by comparing them to satellite-based and reanalysis precipitation products, and iii) investigates possible origins of bimodal precipitation distributions. Our results reveal that about 83% (20 out of 24) of CMIP5 and 70% (21 out of 30) of CMIP6 models used in this study exhibit bimodal distributions. The few DYAMOND models that use a deep convective parameterization also show bimodal distributions, while most DYAMOND models do not. Predictably, the bimodality originates from the separation of precipitation process between resolved grid-scale and parameterized subgrid-scale. However, in a larger number of models bimodality arises from the parameterized subgrid-scale convective precipitation alone.

54 ENVIRONMENTAL SCIENCES↗

Cyber-Physical Events Emulation Based Transmission and Distribution Co-Simulation for Situation Awareness and Grid Anomaly (SAGA) Detection: Preprint

Energy management of transmission and distribution networks is becoming more challenged with the accelerated increasing of distributed energy resources (DERs) such as distributed photovoltaic (PV) generation and distributed energy storage. To better analyze the impacts of DERs on both transmission and distribution systems, a comprehensive transmission and distribution co-simulation platform should be developed. Furthermore, with DERs more actively participated in system operation such as providing real time grid services, their cyber vulnerability should be better understood to maintain system reliability. This paper discussed a cyber-physical events emulation based transmission & distribution co-simulation platform to perform different cyber events emulation and analyze the impacts of cyber physical events happened in distribution system on the T&D system operation. The case studies with both a transmission network and a synthetic distribution network data validate that the proposed T&D co-simulation platform can perform comprehensive cyber physical events emulation. Therefore, with extensive simulation using the proposed model, the system operator can accumulate adequate training data for the system situation awareness and grid anomaly detection purpose.

31 CESER - Office of Cybersecurity, Energy Securit↗

Anthropogenic extremely low volatility organics (ELVOCs) Govern the Growth of Molecular Clusters over the Southern Great Plains during the Springtime

New particle formation (NPF) and growth govern cloud condensation nuclei (CCN) concentrations in many regions. The mechanisms governing the nucleation of molecular clusters vary substantially in different regions of the atmosphere. Additionally, the growth of these clusters from ~2 to 20 nm sizes is often governed by the availability of extremely low volatility organic vapours (ELVOCs). While the pathways to ELVOC formation from the oxidation of biogenic monoterpenes with ozone is better understood, the chemical and mechanistic pathways for ELVOC formation from oxidation of anthropogenic organics are not well understood. We integrate measurements and three-dimensional regional model simulations with the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to understand the processes governing new particle formation and growth and secondary organic aerosol (SOA) formation during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign at the Southern Great Plains (SGP) observatory in Oklahoma, and contrast it with a site within the Bankhead National Forest (BNF), Alabama in Southeast USA, where 5-year long measurements will begin in 2024. Simulations show that nucleation rates are at least an order of magnitude higher at SGP compared to BNF during the springtime days (April 28 and May 14, 2016), largely due to lower H2SO4 concentrations at BNF, which are needed for nucleation. In addition, the larger CS at BNF (compared to SGP) increase the loss of molecular clusters by coagulation to pre-existing particles. Among the 8 different nucleation mechanisms in WRF-Chem, we find that the amine+H2SO4 nucleation mechanism dominates at the SGP site, while the pure organic ion induced nucleation mechanism dominates over BNF. Through various WRF-Chem sensitivity simulations, we find that anthropogenic ELVOCs are critical for explaining the growth of newly formed particles and the resulting number size distribution observed near the surface at the SGP site during the daytime. In addition, we show that treating organic particles as semisolid, with strong diffusion-limited uptake of organic vapours, brings model predictions into closer agreement with the observed evolution of particle size distribution. Simulations also predict that anthropogenic SOA, formed by the oxidation of aromatic volatile organic compounds (VOCs), is the dominant organic aerosol component at SGP, while biogenic SOA dominates particle composition at the BNF site in Southeast USA on these days.

Shrivastava, ManishKumar B.↗

Multivariate quadrature for representing cloud condensation nuclei activity of aerosol populations

Here, sparse representations of atmospheric aerosols are needed for efficient regional- and global-scale chemical transport models. Here we introduce a new framework for representing aerosol distributions, based on the quadrature method of moments. Given a set of moment constraints, we show how linear programming, combined with an entropy-inspired cost function, can be used to construct optimized quadrature representations of aerosol distributions. The sparse representations derived from this approach accurately reproduce cloud condensation nuclei (CCN) activity for realistically complex distributions simulated by a particleresolved model. Additionally, the linear programming techniques described in this study can be used to bound key aerosol properties, such as the number concentration of CCN. Unlike the commonly used sparse representations, such as modal and sectional schemes, the maximum-entropy approach described here is not constrained to pre-determined size bins or assumed distribution shapes. This study is a first step toward a particle-based aerosol scheme that will track multivariate aerosol distributions with sufficient computational efficiency for large-scale simulations.

54 ENVIRONMENTAL SCIENCES↗

Exploring a Variable-Resolution Approach for Simulating Regional Climate in the Rocky Mountain Region Using the VR-CESM: VR-CESM Simulation in the Rocky Mountains

The reliability of climate simulations and projections, particularly in the regions with complex terrains, is greatly limited by the model resolution. Here, we evaluate the variable-resolution Community Earth System Model (VR-CESM) with a high-resolution (0.125°) refinement over the Rocky Mountain region. The VR-CESM results are compared with observations, as well as CESM simulation at a quasi-uniform 1° resolution (UNIF) and Canadian Regional Climate Model version 5 (CRCM5) simulation at a 0.11° resolution. We find that VR-CESM is effective at capturing the observed spatial patterns of temperature, precipitation, and snowpack in the Rocky Mountains with the performance comparable to CRCM5, while UNIF is unable to do so. VR-CESM and CRCM5 simulate better the seasonal variations of precipitation than UNIF, although VR-CESM still overestimates winter precipitation whereas CRCM5 and UNIF underestimate it. All simulations distribute more winter precipitation along the windward (west) flanks of mountain ridges with the greatest overestimation in VR-CESM. VR-CESM simulates much greater snow water equivalent peaks than CRCM5 and UNIF, although the peaks are still 10–40% less than observations. Moreover, the frequency of heavy precipitation events (daily precipitation ≥ 25 mm) in VR-CESM and CRCM5 is comparable to observations, whereas the same events in UNIF are an order of magnitude less frequent. In addition, VR-CESM captures the observed occurrence frequency and seasonal variation of rain-on-snow days and performs better than UNIF and CRCM5. Lastly, these results demonstrate the VR-CESM's capability in regional climate modeling over the mountainous regions and its promising applications for climate change studies.

54 ENVIRONMENTAL SCIENCES↗

An aerosol activation metamodel of v1.2.0 of the pyrcel cloud parcel model: development and offline assessment for use in an aerosol–climate model

We describe an emulator of a detailed cloud parcel model which has been trained to assess droplet nucleation from a complex, multimodal aerosol size distribution simulated by a global aerosol–climate model. The emulator is constructed using a sensitivity analysis approach (polynomial chaos expansion) which reproduces the behavior of the targeted parcel model across the full range of aerosol properties and meteorology simulated by the parent climate model. An iterative technique using aerosol fields sampled from a global model is used to identify the critical aerosol size distribution parameters necessary for accurately predicting activation. Across the large parameter space used to train them, the emulators estimate cloud droplet number concentration (CDNC) with a mean relative error of 9.2% for aerosol populations without giant cloud condensation nuclei (CCN) and 6.9% when including them. Versus a parcel model driven by those same aerosol fields, the best-performing emulator has a mean relative error of 4.6%, which is comparable with two commonly used activation schemes also evaluated here (which have mean relative errors of 2.9 and 6.7%, respectively). We identify the potential for regional biases in modeled CDNC, particularly in oceanic regimes, where our best-performing emulator tends to overpredict by 7%, whereas the reference activation schemes range in mean relative error from -3 to 7%. The emulators which include the effects of giant CCN are more accurate in continental regimes (mean relative error of 0.3%) but strongly overestimate CDNC in oceanic regimes by up to 22%, particularly in the Southern Ocean. Finally, the biases in CDNC resulting from the subjective choice of activation scheme could potentially influence the magnitude of the indirect effect diagnosed from the model incorporating it.

54 ENVIRONMENTAL SCIENCES↗

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗