Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete element modeling”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗
DOE OSTI2023

Finding simplicity: unsupervised discovery of features, patterns, and order parameters via shift-invariant variational autoencoders *

Abstract Recent advances in scanning tunneling and transmission electron microscopies (STM and STEM) have allowed routine generation of large volumes of imaging data containing information on the structure and functionality of materials. The experimental data sets contain signatures of long-range phenomena such as physical order parameter fields, polarization, and strain gradients in STEM, or standing electronic waves and carrier-mediated exchange interactions in STM, all superimposed onto scanning system distortions and gradual changes of contrast due to drift and/or mis-tilt effects. Correspondingly, while the human eye can readily identify certain patterns in the images such as lattice periodicities, repeating structural elements, or microstructures, their automatic extraction and classification are highly non-trivial and universal pathways to accomplish such analyses are absent. We pose that the most distinctive elements of the patterns observed in STM and (S)TEM images are similarity and (almost-) periodicity, behaviors stemming directly from the parsimony of elementary atomic structures, superimposed on the gradual changes reflective of order parameter distributions. However, the discovery of these elements via global Fourier methods is non-trivial due to variability and lack of ideal discrete translation symmetry. To address this problem, we explore the shift-invariant variational autoencoders (shift-VAEs) that allow disentangling characteristic repeating features in the images, their variations, and shifts that inevitably occur when randomly sampling the image space. Shift-VAEs balance the uncertainty in the position of the object of interest with the uncertainty in shape reconstruction. This approach is illustrated for model 1D data, and further extended to synthetic and experimental STM and STEM 2D data. We further introduce an approach for training shift-VAEs that allows finding the latent variables that comport to known physical behavior. In this specific case, the condition is that the latent variable maps should be smooth on the length scale of the atomic lattice (as expected for physical order parameters), but other conditions can be imposed. The opportunities and limitations of the shift VAE analysis for pattern discovery are elucidated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Machine learning pipeline to predict defect behavior in metallic alloy systems

The interaction between defect and solute atoms is critical to the thermodynamic and kinetic behavior of metallic alloys under exposure to high-energy radiation, causing irradiation damage in materials. Radiation can generate non-equilibrium concentrations of point defects such as vacancies and interstitials. The excess point defects not only accelerate diffusional processes such as precipitation that cause radiation embrittlement, but also change the pathway of phase transformations, including nucleation processes. Understanding these defect behaviors is complicated by the challenge and complexity of addressing each possible local and discrete distribution of environments and chemical interactions around targeted defects-solute or solute-solute complexes. To resolve the challenge, machine learning regression techniques have emerged as powerful tools that can train and construct an energy model to accurately describe the chemical interactions of solutes and defects. In Fiscal Year 2022, the work focused on the workflow development and demonstration using machine learning regression, density functional theory, cluster expansion, and Monte Carlo simulation to predict the effects of ternary solute elements (e.g., aluminum and molybdenum) and point defects on the Cr-rich $\alpha^{\prime}$ precipitation in multicomponent FeCr model alloys. The computational outcomes include the prediction of the ternary phase diagram, vacancy formation energy for different compositions, and the effect of vacancies on the nucleation of Cr-rich clusters. The simulations predict a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from $\alpha^{\prime}$ precipitates. Additionally, the simulations show the formation of Cr-vacancy clusters as the initial nuclei for stable nucleation and growth of $\alpha^{\prime}$ particles. The results demonstrate important outcomes and applications of using machine learning pipeline to study model or commercial alloys with multicomponent solute species and point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

The deconstruction of flavor in the privately democratic Higgs sector

The Standard Model (SM) of particle physics fails to explain the observed hierarchy in fermion masses or the origin of fermion-flavor structure. We construct a model to explain these observations in the quark sector. We introduce a spectrum of new particles consisting of six of each — massive singlet vector-like quarks (VLQs), singlet scalars, and SU(2)-doublet scalars. SM quark masses are generated when the neutral components of the SU(2)-doublet scalars acquire non-zero vacuum expectation values (VEVs). We impose global symmetries to ensure that Yukawa couplings stay roughly flavor diagonal and democratic (of the same order), as well as to suppress tree-level flavor-changing neutral currents. Quark-mass hierarchy then follows from a hierarchy in scalar VEVs. The singlet scalars also acquire weak-scale VEVs. Together with the VLQs, they act as messengers between different generations of quarks in the SM. These messenger particles are responsible for generating the elements of the Cabibbo-Kobayashi-Masakawa (CKM) matrix which depend on the ratios of the singlet VEVs and VLQ masses. Constructed this way, the CKM matrix is found to be independent of the SM fermion masses. Using the measured values of the CKM matrix elements and assuming order-one couplings, we derive constraints on the masses of the VLQs and discuss prospects for probing our model in the near future.

CKM Parameters↗
DOE OSTI2026

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-fidelity Pebble Bed Reactor Depletion Based on Pebble Tracking Transport in Griffin

The pebble tracking transport (PTT) method is a high-fidelity, heterogeneous deterministic transport technique for pebble bed reactor analysis. It discretizes the broad-group neutron transport equation in space and angle with the discontinuous finite element and the discrete ordinates method, and utilizes various solving techniques, including mesh sweeping and diffusion acceleration, to provide pebble- wise reaction rates. This work presents the extension of the PTT method to enable fuel depletion capability in the Griffin code. We discuss the implementation details of the PTT-based high-fidelity depletion where isotope inventory of all individual pebbles is tracked through pre-determined pebble flow paths in the core. The implementation is verified with a generic pebble bed reactor model. Some preliminary equilibrium core results are included. Future works are also discussed.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗
DOE OSTI2021

Gluon distributions and their applications to Ioffe-time distributions

We investigate unpolarized and polarized gluon distributions and their applications to the Ioffe-time distributions, which are related to lattice QCD calculations of parton distribution functions. Guided by the counting rules based on the perturbative QCD at large momentum fraction x and the color coherence of gluon couplings at small x, we parametrize gluon distributions in the helicity basis. By fitting the unpolarized gluon distribution, the inferred polarized gluon distribution from our parametrization agrees with the one from global analysis. A simultaneous fit to both unpolarized and polarized gluon distributions is also performed to explore the model uncertainty. The agreement with the global analysis supports the (1 - x) power suppression of the helicity-antialigned distribution relative to the helicity-aligned distribution. The corresponding Ioffe-time distributions and their asymptotic expansions are calculated from the gluon distributions. Our results of the Ioffe-time distributions can provide guidance to the extrapolation of lattice QCD data to the region lacking precise gluonic matrix elements. Therefore, they can help regulate the ill-posed inverse problem associated with extracting the gluon distributions from discrete data from first-principle calculations, which are available in a limited range of the nucleon momentum and the spatial separation between the gluonic currents. Given various limitations in obtaining lattice QCD data at large Ioffe time, phenomenological approaches can provide important complementary information to extract the gluon distributions in the entire momentum fraction region, especially at small x. The possibility of investigating higher-twist effects and other systematic uncertainties in the contemporary first-principle calculations of parton distributions from phenomenologically well-determined Ioffe-time distributions in the large Ioffe-time region is also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The arbitrary‐order virtual element method for linear elastodynamics models: convergence, stability and dispersion‐dissipation analysis

Abstract We design the conforming virtual element method for the numerical approximation of the two‐dimensional elastodynamics problem. We prove stability and convergence of the semidiscrete approximation and derive optimal error estimates under h ‐ and p ‐refinement in both the energy and the L 2 norms. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including nonconvex cells up to order four in the h ‐refinement setting. Exponential convergence is also experimentally observed under p ‐refinement. Finally, we present a dispersion‐dissipation analysis for both the semidiscrete and fully discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion‐dissipation properties.

Antonietti, Paola F.↗
DOE OSTI2023

DFT (and DMFT) for DFT users

The interesting properties of f -electron systems comes from a simultaneous existence of free and discrete level electrons. The challenge we have is to describe free electrons and discrete level electrons equally well in a unified picture. The ELF index can be used to find regions in real space where discrete level physics needs to be taken into account. We have identified the HO gas as a model system that can be used to gain insight about this kind of physics. We will use the HO gas model system for creating a functional suitable for these systems via the subsystem functional scheme. With a predictive capability for the rare earths and lanthindes, we can identify new materials with new/improved properties. This is crucial for, e.g., developing material and technology substitutes for addressing the challenges associated with the use of rare earth elements and other critical materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Verification Suite of Test Cases for the Barotropic Solver of Ocean Models

Abstract The development of any atmosphere or ocean model warrants a suite of test cases (TCs) to verify its spatial and temporal discretizations, order of accuracy, stability, reproducibility, portability, scalability, etc. In this paper, we present a suite of shallow water TCs designed to verify the barotropic solver of atmosphere and ocean models. These include the non‐dispersive coastal Kelvin wave; the dispersive inertia‐gravity wave; the dispersive planetary and topographic Rossby waves; the barotropic tide; and a non‐linear manufactured solution. These TCs check the implementation of the linear pressure gradient term; the linear constant or variable‐coefficient Coriolis and bathymetry terms; and the non‐linear advection terms. Simulation results are presented for a variety of time‐stepping methods as well as two spatial discretizations: a mimetic finite volume method based on the TRiSK scheme, and a high‐order discontinuous Galerkin spectral element method. The experimental procedure for conducting these numerical experiments is detailed. It underscores several key considerations that vary depending on the chosen spatial discretization method. Finally, convergence studies of every TC are conducted with refinement in both space and time, only in space, and only in time. The convergence slopes match the expected theoretical predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Optimal Design of Intensified Towers for CO2 Capture with Internal, Printed Heat Exchangers

Solvent-based carbon capture processes typically suffer from the temperature rise of the solvent due to the heat of absorption of CO2. This increased temperature is not thermodynamically favor-able and results in a significant reduction in performance in the absorber column. As opposed to interstage coolers, which only remove, cool, and return the solvent at discrete locations in the column, internal coolers that are integrated with the packing can cool the process inline, which can result in improved efficiency. This work presents the modeling of these internal coolers within an existing generic, equation-oriented absorber column model that can cool the process while allowing for simultaneous mass transfer. Optimization of this model is also performed, which is capable of optimally choosing the best locations to place these devices, such that heat removal and mass transfer area are balanced. Results of the optimization have shown that optimally placed cooling elements result in a significant increase in the capture efficiency of the process, compared to a similar column with no internal cooling, with a common trend being the cooling of the column in the temperature bulge region. It is observed that by optimally placing an internal cooler, the solvent flow rate can be decreased, and the CO2 lean loading can be increased while still maintaining the same efficiency. These process changes can lead to a substantial reduction in costs due to lower reboiler duty.

Summits, Stephen↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2022

Learning in continuous action space for developing high dimensional potential energy models

Reinforcement learning (RL) approaches that combine a tree search with deep learning have found remarkable success in searching exorbitantly large, albeit discrete action spaces, as in chess, Shogi and Go. Many real-world materials discovery and design applications, however, involve multi-dimensional search problems and learning domains that have continuous action spaces. Exploring high-dimensional potential energy models of materials is an example. Traditionally, these searches are time consuming (often several years for a single bulk system) and driven by human intuition and/or expertise and more recently by global/local optimization searches that have issues with convergence and/or do not scale well with the search dimensionality. Here, in a departure from discrete action and other gradient-based approaches, we introduce a RL strategy based on decision trees that incorporates modified rewards for improved exploration, efficient sampling during playouts and a “window scaling scheme" for enhanced exploitation, to enable efficient and scalable search for continuous action space problems. Using high-dimensional artificial landscapes and control RL problems, we successfully benchmark our approach against popular global optimization schemes and state of the art policy gradient methods, respectively. We demonstrate its efficacy to parameterize potential models (physics based and high-dimensional neural networks) for 54 different elemental systems across the periodic table as well as alloys. We analyze error trends across different elements in the latent space and trace their origin to elemental structural diversity and the smoothness of the element energy surface. Broadly, our RL strategy will be applicable to many other physical science problems involving search over continuous action spaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.

97 MATHEMATICS AND COMPUTING↗