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 235 records · Page 13

Three Dimensional Computer Graphics Federates for the 2012 Smackdown Simulation

The Simulation Interoperability Standards Organization (SISO) Smackdown is a two-year old annual event held at the 2012 Spring Simulation Interoperability Workshop (SIW). A primary objective of the Smackdown event is to provide college students with hands-on experience in developing distributed simulations using High Level Architecture (HLA). Participating for the second time, the University of Alabama in Huntsville (UAHuntsville) deployed four federates, two federates simulated a communications server and a lunar communications satellite with a radio. The other two federates generated 3D computer graphics displays for the communication satellite constellation and for the surface based lunar resupply mission. Using the Light-Weight Java Graphics Library, the satellite display federate presented a lunar-texture mapped sphere of the moon and four Telemetry Data Relay Satellites (TDRS), which received object attributes from the lunar communications satellite federate to drive their motion. The surface mission display federate was an enhanced version of the federate developed by ForwardSim, Inc. for the 2011 Smackdown simulation. Enhancements included a dead-reckoning algorithm and a visual indication of which communication satellite was in line of sight of Hadley Rille. This paper concentrates on these two federates by describing the functions, algorithms, HLA object attributes received from other federates, development experiences and recommendations for future, participating Smackdown teams.

Fordyce, Crystal↗

Simulation of electron distributions within auroral acceleration regions

Plasma data from the Dynamics Explorer 1 satellite are used to investigate the electron distributions found within regions of parallel potential drops. A computer simulation was developed to predict major features of the observed electron distributions found within these regions. The simulation results are discussed and compared with observations. Good agreement is observed only when the potential distribution along the field line is split into separate potential drops, one above and one below the satellite. The simulated electron distribution exhibits a trapped particle region, the formation of the electron hole, and the evacuation of the loss cone.

Gurgiolo, C.↗

DERHCAT (Distributed Energy Resource Hosting Capacity Analysis Tool) [SWR-16-32]

NREL's advanced hosting capacity analysis and System Analysis Reports can help utilities, policymakers, and solar developers better understand the impact of adding new distributed photovoltaic (DPV) systems to the electrical distribution system. Advanced hosting capacity analysis considers the thresholds at which new DPV systems will trigger upgrades or changes to the electrical distribution system and evaluates the cost of different options for expanding the hosting capacity. It is one part of NREL's work to help stakeholders better understand the distribution system costs associated with integrating new DPV systems at different penetration levels. NREL's analysis uses a bottom-up methodology that involves simulating distribution systems. The analysis considers sequential increases in hosting capacity and corresponding upgrade costs from a baseline scenario with no DPV up to penetration levels greater than 100% of peak load.

Ding, Fei↗

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↗

TH-A-19A-04: Latent Uncertainties and Performance of a GPU-Implemented Pre-Calculated Track Monte Carlo Method

Purpose: Assessing the performance and uncertainty of a pre-calculated Monte Carlo (PMC) algorithm for proton and electron transport running on graphics processing units (GPU). While PMC methods have been described in the past, an explicit quantification of the latent uncertainty arising from recycling a limited number of tracks in the pre-generated track bank is missing from the literature. With a proper uncertainty analysis, an optimal pre-generated track bank size can be selected for a desired dose calculation uncertainty. Methods: Particle tracks were pre-generated for electrons and protons using EGSnrc and GEANT4, respectively. The PMC algorithm for track transport was implemented on the CUDA programming framework. GPU-PMC dose distributions were compared to benchmark dose distributions simulated using general-purpose MC codes in the same conditions. A latent uncertainty analysis was performed by comparing GPUPMC dose values to a “ground truth” benchmark while varying the track bank size and primary particle histories. Results: GPU-PMC dose distributions and benchmark doses were within 1% of each other in voxels with dose greater than 50% of Dmax. In proton calculations, a submillimeter distance-to-agreement error was observed at the Bragg Peak. Latent uncertainty followed a Poisson distribution with the number of tracks per energy (TPE) and a track bank of 20,000 TPE produced a latent uncertainty of approximately 1%. Efficiency analysis showed a 937× and 508× gain over a single processor core running DOSXYZnrc for 16 MeV electrons in water and bone, respectively. Conclusion: The GPU-PMC method can calculate dose distributions for electrons and protons to a statistical uncertainty below 1%. The track bank size necessary to achieve an optimal efficiency can be tuned based on the desired uncertainty. Coupled with a model to calculate dose contributions from uncharged particles, GPU-PMC is a candidate for inverse planning of modulated electron radiotherapy and scanned proton beams. This work was supported in part by FRSQ-MSSS (Grant No. 22090), NSERC RG (Grant No. 432290) and CIHR MOP (Grant No. MOP-211360)

07 ISOTOPE AND RADIATION SOURCES↗

Implications of Simulation Conceptual Model Development for Simulation Management and Uncertainty Assessment

A simulation conceptual model is a simulation developers way of translating modeling requirements (i. e., what is to be represented by the simulation or its modification) into a detailed design framework (i. e., how it is to be done), from which the software, hardware, networks (in the case of distributed simulation), and systems/equipment that will make up the simulation can be built or modified. A conceptual model is the collection of information which describes a simulation developers concept about the simulation and its pieces. That information consists of assumptions, algorithms, characteristics, relationships, and data. Taken together, these describe how the simulation developer understands what is to be represented by the simulation (entities, actions, tasks, processes, interactions, etc.) and how that representation will satisfy the requirements to which the simulation responds. Thus the conceptual model is the basis for judgment about simulation fidelity and validity for any condition that is not specifically tested. The more perspicuous and precise the conceptual model, the more likely it is that the simulation development will both fully satisfy requirements and allow demonstration that the requirements are satisfied (i. e., validation). Methods used in simulation conceptual model development have significant implications for simulation management and for assessment of simulation uncertainty. This paper suggests how to develop and document a simulation conceptual model so that the simulation fidelity and validity can be most effectively determined. These ideas for conceptual model development apply to all simulation varieties. The paper relates these ideas to uncertainty assessments as they relate to simulation fidelity and validity. The paper also explores implications for simulation management from conceptual model development methods, especially relative to reuse of simulation components.

Pace, Dale K.↗

Long-Term Simulation of Dust Distribution with the GOCART Model: Correlation with the North Atlantic Oscillation

Global distribution of aeolian dust is simulated from 1981 to 1996 with the Goddard Ozone Chemistry Aerosol Radiation and Transport (GOCART) model. The results are assessed with in-situ measurements and the Total Ozone Mapping Spectrometer (TOMS) aerosol products. The annual budget over the different continents and oceans are analyzed. It is found that there is a maximum of 25% difference of global annual emission from the minimum in 1996 to the maximum in 1988. There is a downward trend of dust emission over Africa and East Asia, of 6 and 2 Tg/yr, respectively. The inter-annual variability of dust distribution is analyzed over the North Atlantic and Africa. It is found that in winter most of the North Atlantic and Africa dust loading is correlated with the North Atlantic Oscillation. The GOCART model indicates that a controlling factor of such correlation can be attributed to dust emission from the Sahel. The Bodele depression is the major dust source in winter and its inter-annual variability is highly correlated with the NAO. However, it is not possible to conclude without further analysis that the North Atlantic Oscillation is forcing the inter-annual variability of dust emission and in-turn dust concentration over the North Atlantic.

Ginoux, P.↗

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↗

VME rollback hardware for time warp multiprocessor systems

The purpose of the research effort is to develop and demonstrate innovative hardware to implement specific rollback and timing functions required for efficient queue management and precision timekeeping in multiprocessor discrete event simulations. The previously completed phase 1 effort demonstrated the technical feasibility of building hardware modules which eliminate the state saving overhead of the Time Warp paradigm used in distributed simulations on multiprocessor systems. The current phase 2 effort will build multiple pre-production rollback hardware modules integrated with a network of Sun workstations, and the integrated system will be tested by executing a Time Warp simulation. The rollback hardware will be designed to interface with the greatest number of multiprocessor systems possible. The authors believe that the rollback hardware will provide for significant speedup of large scale discrete event simulation problems and allow multiprocessors using Time Warp to dramatically increase performance.

Robb, Michael J.↗

SU-E-J-175: Comparison of the Treatment Reproducibility of Tumors Affected by Breathing Motion

Purpose: The aim of the dose distribution simulations was to form a global idea of intensity-modulated radiation therapy (IMRT) realization, by its comparison to three-dimensional conformal radiation therapy (3DCRT) delivery for tumors affected by respiratory motion. Methods: In the group of 10patients both 3DCRT and IMRT plans were prepared.For each field the motion kernel was generated with the largest movement amplitude of 4;6 and 8mm.Additionally,the sets of reference measurements were made in no motion conditions(0 mm).The evaluation of plan delivery,using a diode array placed on moving platform,was based on the Gamma Index analysis with distance to agreement of 3mm and dose difference of 3%. Results: IMRT plans tended to spare doses delivered to lungs compared to 3DCRT.Nonetheless,analyzed volumes showed no significant difference between the static and dynamic techniques,except for the volumes of both lungs receiving 10 and 15Gy.After adding the components associated with the respiratory movement,all IMRT lung parameters evaluated for the ipsilateral,contralateral and both lungs together,revealed considerable differences between the 0vs.6, 0vs.8 and 4vs.8-mm amplitudes.Similar results were obtained for the 3DCRT lung measurements,but without significance between the 0vs.6-mm amplitude.Taking into account the CTV score parameter in 3DCRT and IMRT plans,there was no statistically significant difference between the motion patterns with the smallest amplitudes.The differences were found for the 8-mm amplitude when it was compared both with static conditions and 4-mm amplitude (for 3DCRT) and between 0vs.6, 0vs.8 and 4vs.8-mm amplitudes (for IMRT).All accepted and measured 3DCRT and IMRT doses to spinal cord,esophagus and heart were always below the QUANTEC limits. Conclusion: The application of IMRT technique in lung radiotherapy affords possibilities for reducing the lung doses.For maximal amplitudes of breathing trajectory below 4mm,the disagreement between CTV planned and delivered dose distribution can be neglected.Whereas respiratory motion with an amplitude above 5mm leads to significant changes in IMRT delivery.

60 APPLIED LIFE SCIENCES↗

Technical analysis of prospective photovoltaic systems in Utah

This report explores the technical feasibility of prospective utility-scale photovoltaic system (PV) deployments in Utah. Sandia National Laboratories worked with Rocky Mountain Power (RMP), a division of PacifiCorp operating in Utah, to evaluate prospective 2-megawatt (MW) PV plants in different locations with respect to energy production and possible impact on the RMP system and customers. The study focused on 2-MW{sub AC} nameplate PV systems of different PV technologies and different tracking configurations. Technical feasibility was evaluated at three different potential locations in the RMP distribution system. An advanced distribution simulation tool was used to conduct detailed time-series analysis on each feeder and provide results on the impacts on voltage, demand, voltage regulation equipment operations, and flicker. Annual energy performance was estimated.

42 ENGINEERING↗

High Performance Interactive System Dynamics Visualization

This brochure describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

High Performance Interactive System Dynamics Visualization

This presentation describes a system dynamics simulation (SD) framework that supports an end-to-end analysis workflow that is optimized for deployment on ESIF facilities(Peregrine and the Insight Center). It includes (I) parallel and distributed simulation of SD models, (ii) real-time 3D visualization of running simulations, and (iii) comprehensive database-oriented persistence of simulation metadata, inputs, and outputs.

97 MATHEMATICS AND COMPUTING↗

A Sparse Distributed Gigascale Resolution Material Point Method

In this paper, we present a four-layer distributed simulation system and its adaptation to the Material Point Method (MPM). The system is built upon a performance portable C++ programming model targeting major High-Performance-Computing (HPC) platforms. A key ingredient of our system is a hierarchical block-tile-cell sparse grid data structure that is distributable to an arbitrary number of Message Passing Interface (MPI) ranks. We additionally propose strategies for efficient dynamic load balance optimization to maximize the efficiency of MPI tasks. Our simulation pipeline can easily switch among backend programming models, including OpenMP and CUDA, and can be effortlessly dispatched onto supercomputers and the cloud. Finally, we construct benchmark experiments and ablation studies on supercomputers and consumer workstations in a local network to evaluate the scalability and load balancing criteria. We demonstrate massively parallel, highly scalable, and gigascale resolution MPM simulations of up to 1.01 billion particles for less than 323.25 seconds per frame with 8 OpenSSH-connected workstations.

97 MATHEMATICS AND COMPUTING↗

Range Systems Simulation for the NASA Shuttle: Emphasis on Disaster and Prevention Management During Lift-Off

This article describes a decision-making system composed of a number of safety and environmental models for the launch phase of a NASA Space Shuttle mission. The components of this distributed simulation environment represent the different systems that must collaborate to establish the Expectation of Casualties (E(sub c)) caused by a failed Space Shuttle launch and subsequent explosion (accidental or instructed) of the spacecraft shortly after liftoff. This decision-making tool employs Space Shuttle reliability models, trajectory models, a blast model, weather dissemination systems, population models, amount and type of toxicants, gas dispersion models, human response functions to toxicants, and a geographical information system. Since one of the important features of this proposed simulation environment is to measure blast, toxic, and debris effects, the clear benefits is that it can help safety managers not only estimate the population at risk, but also to help plan evacuations, make sheltering decisions, establish the resources required to provide aid and comfort, and mitigate damages in case of a disaster.

Rabelo, Lisa↗

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗

RETIRED A STARS: THE EFFECT OF STELLAR EVOLUTION ON THE MASS ESTIMATES OF SUBGIANTS

Doppler surveys have shown that the occurrence rate of Jupiter-mass planets appears to increase as a function of stellar mass. However, this result depends on the ability to accurately measure the masses of evolved stars. Recently, Lloyd called into question the masses of subgiant stars targeted by Doppler surveys. Lloyd argues that very few observable subgiants have masses greater than 1.5 M {sub Sun }, and that most of them have masses in the range 1.0-1.2 M {sub Sun }. To investigate this claim, we use Galactic stellar population models to generate an all-sky distribution of stars. We incorporate the effects that make massive subgiants less numerous, such as the initial mass function and differences in stellar evolution timescales. We find that these effects lead to negligibly small systematic errors in stellar mass estimates, in contrast to the Almost-Equal-To 50% errors predicted by Lloyd. Additionally, our simulated target sample does in fact include a significant fraction of stars with masses greater than 1.5 M {sub Sun }, primarily because the inclusion of an apparent magnitude limit results in a Malmquist-like bias toward more massive stars, in contrast to the volume-limited simulations of Lloyd. The magnitude limit shifts the mean of our simulated distribution toward higher masses and results in a relatively smaller number of evolved stars with masses in the range 1.0-1.2 M {sub Sun }. We conclude that, within the context of our present-day understanding of stellar structure and evolution, many of the subgiants observed in Doppler surveys are indeed as massive as main-sequence A stars.

79 ASTRONOMY AND ASTROPHYSICS↗