skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Optimistic Simulations of Physical Systems using Reverse Computation

Abstract

Efficient computer simulation of complex physical phenomena has long been challenging due to their multi-physics and multi-scale nature. In contrast to traditional time-stepped execution methods, we describe an approach using optimistic parallel discrete event simulation (PDES) and reverse computation techniques to execute plasma physics codes. We show that reverse computation-based optimistic parallel execution can significantly reduce the execution time of an example plasma simulation without requiring a significant amount of additional memory compared to conservative execution techniques. We describe an application-level reverse computation technique that is efficient and suitable for complex scientific simulations.

Authors:
 [1];  [2];  [2];  [3];  [3];  [3]
  1. Georgia Institute of Technology
  2. ORNL
  3. SciberQuest Inc.
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE Laboratory Directed Research and Development (LDRD) Program
OSTI Identifier:
1003413
DOE Contract Number:
DE-AC05-00OR22725
Resource Type:
Journal Article
Resource Relation:
Journal Name: SIMULATION: Transactions of The Society for Modeling and Simulation International; Journal Volume: 82; Journal Issue: 1
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICAL METHODS AND COMPUTING; COMPUTERIZED SIMULATION; PLASMA SIMULATION; CALCULATION METHODS

Citation Formats

Tang, Yarong, Perumalla, Kalyan S, Fujimoto, Richard, Karimabadi, Dr. Homa, Driscoll, Jonathan, and Omelchenko, Yuri. Optimistic Simulations of Physical Systems using Reverse Computation. United States: N. p., 2006. Web. doi:10.1177/0037549706065481.
Tang, Yarong, Perumalla, Kalyan S, Fujimoto, Richard, Karimabadi, Dr. Homa, Driscoll, Jonathan, & Omelchenko, Yuri. Optimistic Simulations of Physical Systems using Reverse Computation. United States. doi:10.1177/0037549706065481.
Tang, Yarong, Perumalla, Kalyan S, Fujimoto, Richard, Karimabadi, Dr. Homa, Driscoll, Jonathan, and Omelchenko, Yuri. Sun . "Optimistic Simulations of Physical Systems using Reverse Computation". United States. doi:10.1177/0037549706065481.
@article{osti_1003413,
title = {Optimistic Simulations of Physical Systems using Reverse Computation},
author = {Tang, Yarong and Perumalla, Kalyan S and Fujimoto, Richard and Karimabadi, Dr. Homa and Driscoll, Jonathan and Omelchenko, Yuri},
abstractNote = {Efficient computer simulation of complex physical phenomena has long been challenging due to their multi-physics and multi-scale nature. In contrast to traditional time-stepped execution methods, we describe an approach using optimistic parallel discrete event simulation (PDES) and reverse computation techniques to execute plasma physics codes. We show that reverse computation-based optimistic parallel execution can significantly reduce the execution time of an example plasma simulation without requiring a significant amount of additional memory compared to conservative execution techniques. We describe an application-level reverse computation technique that is efficient and suitable for complex scientific simulations.},
doi = {10.1177/0037549706065481},
journal = {SIMULATION: Transactions of The Society for Modeling and Simulation International},
number = 1,
volume = 82,
place = {United States},
year = {Sun Jan 01 00:00:00 EST 2006},
month = {Sun Jan 01 00:00:00 EST 2006}
}
  • Vehicular traffic simulations are useful in applications such as emergency management and homeland security planning tools. High speed of traffic simulations translates directly to speed of response and level of resilience in those applications. Here, a parallel traffic simulation approach is presented that is aimed at reducing the time for simulating emergency vehicular traffic scenarios. Three unique aspects of this effort are: (1) exploration of optimistic simulation applied to vehicular traffic simulation (2) addressing reverse computation challenges specific to optimistic vehicular traffic simulation (3) achieving absolute (as opposed to self-relative) speedup with a sequential speed equal to that of amore » fast, de facto standard sequential simulator for emergency traffic. The design and development of the parallel simulation system is presented, along with a performance study that demonstrates excellent sequential performance as well as parallel performance.« less
  • Reverse computation is presented here as an important future direction in addressing the challenge of fault tolerant execution on very large cluster platforms for parallel computing. As the scale of parallel jobs increases, traditional checkpointing approaches suffer scalability problems ranging from computational slowdowns to high congestion at the persistent stores for checkpoints. Reverse computation can overcome such problems and is also better suited for parallel computing on newer architectures with smaller, cheaper or energy-efficient memories and file systems. Initial evidence for the feasibility of reverse computation in large systems is presented with detailed performance data from a particle simulation scalingmore » to 65,536 processor cores and 950 accelerators (GPUs). Reverse computation is observed to deliver very large gains relative to checkpointing schemes when nodes rely on their host processors/memory to tolerate faults at their accelerators. A comparison between reverse computation and checkpointing with measurements such as cache miss ratios, TLB misses and memory usage indicates that reverse computation is hard to ignore as a future alternative to be pursued in emerging architectures.« less
  • In this article, we investigated structural and thermal properties of an amorphous alloy of the Ga–Se system. The amorphous GaSe{sub 9} alloy was produced by mechanical alloying and it was studied using EXAFS spectroscopy and cumulant expansion method. We also made reverse Monte Carlo simulations using the total structure factor S(K) obtained from x-ray diffraction and the EXAFS χ(k) oscillations on Se and Ga K edges as input data. Several parameters, such as average coordination numbers and interatomic distances, structural and thermal disorders, asymmetry of the partial distribution functions g{sub ij}(r), and Einstein and Debye temperatures, were determined. The g{submore » ij}{sup E}(r) functions were reconstructed from the cumulants C{sub 1}, C{sub 2}, and C{sub 3} obtained from the Einstein model, and they were compared to the g{sub ij}{sup RMC}(r) functions obtained from the simulations. The simulations also furnished the partial bond angle distribution functions Θ{sub ijℓ}(cosθ), which describe the angular distribution of bonds between first neighbors, and give information about the kind of structural units present in the alloy.« less
  • System descriptors are given for key hydrogen economy components. Hypothetical transition scenarios are given for three generalized social climates with regard to hydrogen: strong government intervention, optimistic and realistic. All transition scenarios are characterized by five eras related to the density of the hydrogen distribution network. (auth)
  • In this article we compare the halo mass function predicted by the excursion set theory with a drifting diffusive barrier against the results of N-body simulations for several cosmological models. This includes the standard ΛCDM case for a large range of halo masses, models with different types of primordial non-Gaussianity, and the Ratra-Peebles quintessence model of Dark Energy. We show that in all those cosmological scenarios, the abundance of dark matter halos can be described by a drifting diffusive barrier, where the two parameters describing the barrier have physical content. In the case of the Gaussian ΛCDM, the statistics aremore » precise enough to actually predict those parameters at different redshifts from the initial conditions. Furthermore, we found that the stochasticity in the barrier is non-negligible making the simple deterministic spherical collapse model a bad approximation even at very high halo masses. We also show that using the standard excursion set approach with a barrier inspired by peak patches leads to inconsistent predictions of the halo mass function.« less