Search NASA⌕ Search

SEARCH · Search NASA

Results for “Virtual element method”

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 55 records · Page 3

Virtual Design and Analysis of the Balancing Element of an External Gear Machine Considering Cavitation and Mixed Lubrication Effects

This report presents a simulation based approach for sizing the axial balance elements of a pressure compensated external gear machine. The reference hydraulic unit must be able to guarantee operation in two quadrants, as either a pump or a motor. It also needs to operate in a wide operating range of speed and pressure, as it is to be used as primary hydraulic unit for an electro-hydraulic actuator (EHA). The design procedure builds upon past work at the author’s center, but it extends it to the case of multiple quadrant units. Also, a unique method for modeling of both mixed lubrication and cavitation that might occur in the lubricating film is used. After describing the design procedure, the paper discusses the features of the balancing of a reference unit, along with the effects of both cavitation and mixed lubrication.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗

Multi‐fidelity digital twin structural model for a sub‐scale downwind wind turbine rotor blade

Abstract This paper presents the development of a multi‐fidelity digital twin structural model (virtual model) of an as‐built wind turbine blade. The goal is to develop and demonstrate an approach to produce an accurate and detailed model of the as‐built blade for use in verifying the performance of the operating two‐bladed, downwind rotor. The digital twin model development methodology, presented herein, involves a novel calibration process to integrate a wide range of information including design specifications, manufacturing information, and structural testing data (modal and static) to produce a multi‐fidelity digital twin structural model: a detailed high‐fidelity model (i.e., 3D finite element analysis [FEA]) and consistent beam‐type models for aeroelastic simulation. A key element is that the multi‐fidelity structural digital twin method follows the rotor from the stages of design, to manufacturing, then to the ground testing and field operation. The result of this comprehensive approach is an accurate multi‐fidelity digital twin structural model for the geometric, structural, and structural dynamic properties of the as‐built blade within a 1% match in mass properties, 3.2% in blade frequencies, and 6% in deflection. The different stages of processing this information within the methodology are discussed. The rotor examined is the SUMR‐Demonstrator (SUMR‐D), which was installed on the Controls Advanced Research Testbed (CART‐2) wind turbine at the National Wind Technology Center. The digital twin model developed here was utilized to design controllers to safely operate SUMR‐D in field tests, which are providing additional data for further evaluation and development of the multi‐fidelity digital twin structural model.

Chetan, Mayank↗

Comprehensive Study of the Potential of Extracting and Processing Critical Minerals from Coal-Based Resources - Phase I

The Phase I report prepared for the Department of Energy addresses U.S. Executive Order 13817 titled A Federal Strategy to Ensure Secure and Reliable Supplies of Critical Minerals, issued on December 20, 2017, that lists 35 critical minerals that are vulnerable to supply disruption. A comprehensive review of each of the minerals was conducted to determine the criticality based primarily on extractability from coal-based resources. Several other factors were also considered such as gaps in supply and demand, use in current technology, and the existence of viable substitutes. It was determined that the critical minerals that show highest potential for extraction from coal-based resources are lithium, rare earth elements (REEs), cobalt, and manganese. All four of the critical minerals listed serve an important role in the technology industry, have few substitutes, and have a heavy import reliance. Most notably are lithium, which is widely used in the electric vehicle industry, and the REEs which can be found in virtually all electronic devices. The research of critical mineral extraction from coal-based resources was completed using a combination of literature review from public sources as well as cooperation from coal mines and power plants across the United States. Samples collected from six different geographical locations across the U.S. were subjected to sample preparation (i.e., pH measurement, moisture content, particle size analysis) and characterization studies using Inductively Coupled Plasma-Mass Spectroscopy (ICP-MS) and Scanning Electron Microscopy, Energy Dispersive X-Ray Spectroscopy (SEM-EDX) instruments. 27 samples of coal waste materials such as refuse, sludge, and fly ash were tested to characterize the rare earth element concentration by total rare earth elements (TREEs), heavy rare earth elements (HREEs), and light rare earth elements (LREEs). Of the 27 samples tested, 22 contained a TREE concentration higher than the threshold of 300 ppm, which is considered a viable feedstock material. 3 samples contained less than 300 ppm of TREEs; however, they were within 20 ppm of the threshold, and could potentially be considered viable sources in the future pending the advancement of more efficient extraction technologies. 2 of the 27 samples had significantly low TREE concentrations, which does not imply any potential for being a source for REEs. For the minerals identified as most critical in the literature review, a conceptual process flow diagram (PFD) was developed for their extraction from different coal-based feedstocks. The process targets selective recovery of one commodity (i.e., rare earths, lithium, cobalt, and manganese) via several hydrometallurgical separation methods. By identifying potentially extractable coal-based critical mineral resources, a study of the current and future market environments for each critical mineral, and a review of current processing methodologies for critical mineral extraction from coal-based resources, the foundation has been laid to further characterize and explore new resources and extraction techniques. As reliance on technologies in industries such as the production of electronic devices, batteries, and alloys containing critical minerals utilizing critical minerals continues to increase, a sound understanding of our nation’s dependence on and even the global criticality of certain critical minerals, will serve as a catalyst for innovation in virtually all fields of science.

01 COAL, LIGNITE, AND PEAT↗

Gluon unpolarized, polarized, and transversity GPDs from lattice QCD: Lorentz-covariant parametrization

We identify the matrix elements necessary to determine the leading-twist gluon generalized parton distributions (GPDs) H g , E g , H ˜ g , E ˜ g , H g T , E g T , H ˜ g T , E ˜ g T in lattice QCD calculations. We present a method to achieve a Lorentz-covariant parameterization of the matrix elements in terms of a linearly independent basis of tensor structures. This parameterization is crucial for projecting lattice QCD matrix elements onto light cone distributions. For the first time, we determine the corresponding components that project onto the linear combinations of invariant amplitudes, which reduce to the different gluon GPDs in the light cone limit and enable their separation in a lattice QCD calculation for spin-0 and spin- 1 2 hadrons. Hence, this work lays the foundation for the numerical determination of the gluon GPDs from first-principle lattice QCD calculations, directly advancing our understanding of the mass and spin structures and mechanical properties of the nucleon, as well as the physics underlying deeply virtual Compton scattering and deeply virtual meson production in a range of experimental processes. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A VIPA Spectrograph with Ultra-high Resolution and Wavelength Calibration for Astronomical Applications

A prototype spectrograph using a Virtually Imaged Phased Array (VIPA) as the main dispersion element is presented, and its performance is fully examined in our laboratory. The single-mode, fiber-fed spectrograph with simultaneous wavelength calibration possesses a spectral resolution well in excess of R≈1.12×10{sup 6} while the size of the VIPA is several orders of magnitude smaller than that of a conventional échelle with comparable resolution. In laboratory tests, the VIPA-based instrument with a homemade Yb:fiber ring laser frequency comb demonstrates a mode-to-mode tracking stability of 41 cm s{sup −1} over a period of 6 hr. The VIPA spectrograph has promising applications in various astronomical observations in which ultra-high resolution and calibration precision are imperative, such as solar physics research, exoplanet searching with the radial velocity method, and O{sub 2} detection in the atmosphere of Earth-like planets. Ultimately, feasible optimizations for night-sky observations under seeing limited conditions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Rethinking the Price Formation Problem–Part 2: Rewarding Flexibility and Managing Price Risk

In this study, part 1 of this two-part paper describes the impact that uncertainty has on the design and analysis of price formation policies in the non-convex auctions conducted by U.S. wholesale electricity market operators. Using first a toy model and then a large-scale test system, Part 2 demonstrates the difference in prices under the idealized benchmark of ex ante convex hull pricing defined in Part 1 versus existing methods, in particular documenting the potential for suppression of volatility and therefore under-compensation of flexibility by existing methods. The examples highlight that inefficient spot price formation can induce inefficient forward commitments of generators, necessitating out-of-market intervention to restore a reliable and efficient operating plan.Given the potential side effects of existing policies for investment and operation, we suggest two elements in a reoriented approach to the price formation problem: first ensuring that prices exhibit full-strength volatility, and second ensuring that risk-averse market participants have sufficient ability to manage this volatility.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multiphysics Coupling Methods for Molten Salt Reactor Modeling and Simulation in VERA

New multiphysics coupling capabilities for molten salt reactor (MSR) analysis have been developed in the Virtual Environment for Reactor Applications (VERA). This development consisted of two main efforts. First, a generic species transport module was added in the CTF code, which is the thermal-hydraulics (TH) code for VERA. This module uses the velocity fields for which CTF solves during the TH calculation to transport species through the core and around the primary loop. Additionally, a gas sparging model has been added to CTF to model the movement of certain species, namely, fission products such as xenon gas, to transport between the molten salt and gas bubbles present in the salt. The second effort in this development was coupling this capability to VERA’s neutron transport code MPACT. This effort focused on coupling the detailed TH transport models in CTF to MPACT to account for feedback effects in the neutron transport calculations. Finally, the thermochemistry code Thermochimica has also been coupled to VERA. Thermochimica performs pointwise calculations for chemical potential and Gibbs free energy and determines what phases are produced by the temperature, pressure, and elemental concentrations at different locations in the primary loop. These capabilities are demonstrated using a model of the Molten Salt Reactor Experiment (MSRE). Furthermore, this reactor operated at Oak Ridge National Laboratory in the 1960s, providing sources of experimental data that were used to develop the model. Various combinations of species were modeled using VERA’s new multiphysics coupling capabilities. Species distributions and reactivity effects behaved as expected for the MSRE model, demonstrating that the coupling is behaving correctly and causing appropriate feedback. The results of these calculations show the potential for VERA to be used for a wide variety of MSR analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An optimization-based strategy for peridynamic-FEM coupling and for the prescription of nonlocal boundary conditions

We develop and analyze an optimization-based method for the coupling of a static peridynamic (PD) model and a static classical elasticity model. The approach formulates the coupling as a control problem in which the states are the solutions of the PD and classical equations, the objective is to minimize their mismatch on an overlap of the PD and classical domains, and the controls are virtual volume constraints and boundary conditions applied at the local-nonlocal interface. Our numerical tests performed on three-dimensional geometries illustrate the consistency and accuracy of our method, its numerical convergence, and its applicability to realistic engineering geometries. We demonstrate the coupling strategy as a means to reduce computational expense by confining the nonlocal model to a subdomain of interest, and as a means to transmit local (e.g., traction) boundary conditions applied at a surface to a nonlocal model in the bulk of the domain.

97 MATHEMATICS AND COMPUTING↗

Initial Simulations of Empty Room Collapse and Reconsolidation at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Virtual Engineering of Low-Temperature Conversion

In this work, we present the development of an overarching software framework and supporting multiphysics models to simulate the end-to-end process of biomass conversion. This virtual engineering (VE) software is designed with the goal of accelerating research and development and reducing risk for market-relevant biomass conversion processes. We currently support multiple models, computing paradigms, and fidelities representing the steps of feedstock pretreatment, enzymatic hydrolysis, and bioconversion. Although this VE approach was developed to support a biomass workflow, we have designed each component in a way that allows us to easily support new domains, unit models, and feedstocks. We begin by presenting the user-facing aspects of the VE software and highlight how simulated elements are defined and linked by VE functions. We then present an overview of the high-fidelity computational fluid dynamics (CFD) models developed to support our target domain before segueing into our efforts to develop accurate and fast surrogate models, capturing the salient outcomes from the CFD simulations in a significantly less computationally demanding manner. We then present how VE calculations interface with a commercial techno-economic analysis software, Aspen Plus. We conclude by presenting VE case studies that leverage these methods and discuss how our methods can be extended to support a wide variety of accelerated biofuel commercialization pathways in the future.

biofuel↗

Inertia estimation for power grids: A review of methods, challenges, and future prospects

The electric power grid is undergoing a significant transformation, shifting from traditional synchronous generators to inverter-based resources (IBRs) such as solar photovoltaics, wind turbines, and energy storage systems. This evolution leads to a reduction in system inertia, a critical attribute for maintaining frequency stability in response to disturbances. Consequently, the ability to monitor and estimate system inertia has become increasingly essential. This paper provides a comprehensive review of existing inertia estimation methodologies, analyzing them from multiple perspectives, including the types of data utilized, underlying estimation principles, operational modes, and system-wide applicability. A comparative summary table is included to distill commonalities and key characteristics across various studies. In addition, the paper examines practical implementations of inertia estimation across several major power systems worldwide, including the U.S. interconnections, the Nordic power system, and the U.K. grid. Key challenges are identified, particularly in estimating contributions from virtual inertia sources and load-induced inertia in increasingly converter-dominated networks. To address these emerging challenges, the paper proposes an integrated framework for real-time inertia estimation and monitoring. This framework encompasses critical components such as data acquisition, inertia estimation from both synchronous and non-synchronous sources, load-induced effects, optimization techniques, forecasting, and virtual inertia scheduling. Collectively, these elements enable dynamic, system-wide monitoring and adaptive control of grid inertia.

Inertia estimation↗

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Dynamic Networks Experiments: Virtual Experiments to Quantify Gains in Nuclear Explosion Monitoring

We describe an ongoing series of virtual experiments conducted collaboratively by four United States National Laboratories: Sandia National Laboratories, Los Alamos National Laboratory, Lawrence Livermore National Laboratory, and Pacific Northwest National Laboratory. These Dynamic Network Experiments (DNEs) provide an experimental framework to evaluate the potential impact of new research tools on nuclear explosion monitoring. The second DNE (DNE2), completed in 2024, exploited waveform data (seismic, infrasound, and electromagnetic) that was recorded by multi-modal sensors within and near the Nevada National Security Site and synthetic radionuclide signatures over multiple time periods. During the execution of DNE2, we processed and analyzed data through a multi-stage event processing pipeline that ingested raw data, performed quality control, detected signals, built events from these signals, located these events, and characterized the events’ source types and sizes. For each stage and over the entire event processing pipeline, we evaluated performance changes by comparing the performance of new data processing methods, models, and algorithms against a baseline. We also performed an additional execution phase to assess event processing pipeline function, speed, and efficiency against that of an expert analyst, including computational and manual efforts. Finally, we assessed the impact and effort of modern computing infrastructure on the monitoring pipeline. This paper describes key elements of the DNEs, from formulation through execution, as demonstrated in DNE2. The DNEs introduce several novel concepts to quantitatively measure the potential impact of new methods on explosion monitoring, including the collaborative design of multi-modal datasets, performance and logistical metrics, and integrated analyses.

42 ENGINEERING↗

Rancor-HUNTER: A Virtual Plant and Operator Environment for Predicting Human Performance

Advances in simulation capabilities to model physical systems have outpaced the development of simulations for humans using those physical systems. There is an argument that the infinite span of potential human behaviors inherently render human modeling more challenging than physical systems. Despite this challenge, the need for modeling humans interacting with these complex systems is paramount. As technologies have improved, many of the failure modes originating from the physical systems have been solved. This means the overall proportion of human errors has increased, such that it is not uncommon to be the primary driver of system failure in modern complex systems. Moreover, technologies such as automated systems may introduce emerging contexts that can cause new, unanticipated modes of human error. Therefore, it is now more important than ever to develop models of human behavior to realize overall system error reductions and achieve established safety margins. To support new and novel concepts of operations for the anticipated wave of advanced nuclear reactor deployments, human factors and human reliability analysis researchers need to develop advanced simulation-based approaches. This talk presents a simulation environment suitable to both collect data and then perform Monte Carlo simulations to evaluate human performance and develop better models of human behavior. Specifically, the Rancor Microworld Simulator models a complex energy production system in a simplified manner. Rancor includes computer-based procedures, which serve as a framework to automatically classify human behaviors without manual, subjective experimenter coding during scenarios. This method supports a detailed level of analysis at the task level. It is feasible for collecting large sample sizes required to develop quantitative modelling elements that have historically challenged traditional full-scope simulator study approaches. Additionally, the other portion of this experimental platform, the Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER), is presented to show how the collected data can be used to evaluate novel scenarios based on the contextual factors, or performance shaping factors, derived from Rancor simulations. Rancor-HUNTER is being used to predict operator performance with new procedures, such as results from control room modernization or new-build situations. Rancor-HUNTER is also proving a useful surrogate platform to model human performance for other complex systems.

99 GENERAL AND MISCELLANEOUS↗

Randomized Algorithms for Scientific Computing (RASC)

Randomized algorithms have propelled advances in artificial intelligence (AI) and represent a foundational research area in advancing AI for Science. Future advancements in DOE Office of Science priority areas such as climate science, astrophysics, fusion, advanced materials, combustion, and quantum computing all require randomized algorithms for surmounting challenges of complexity, robustness, and scalability. Advances in data collection and numerical simulation have changed the dynamics of scientific research and motivate the need for randomized algorithms. For instance, advances in imaging technologies such as X-ray ptychography, electron microscopy, electron energy loss spectroscopy, or adaptive optics lattice light-sheet microscopy collect hyperspectral imaging and scattering data in terabytes, at breakneck speed enabled by state-of-the-art detectors. The data collection is exceptionally fast compared with its analysis. Likewise, advances in high-performance architectures have made exascale computing a reality and changed the economies of scientific computing in the process. Floating-point operations that create data are essentially free in comparison with data movement. Thus far, most approaches have focused on creating faster hardware. Ironically, this faster hardware has exacerbated the problem by making data still easier to create. Under such an onslaught, scientists often resort to heuristic deterministic sampling schemes (e.g., low-precision arithmetic, sampling every nth element) and sacrifice potentially valuable accuracy. Dramatically better results can be achieved via randomized algorithms, reducing the data size as much as or more than naive deterministic subsampling can achieve, while retaining the high accuracy of computing on the full data set. By randomized algorithms we mean those algorithms that employ some form of randomness in internal algorithmic decisions to accelerate time to solution, increase scalability, or improve reliability. Examples include matrix sketching for solving large-scale least-squares problems (see Figure 1) and stochastic gradient descent for training machine learning models. We are not recommending heuristic methods but rather randomized algorithms that have certificates of correctness and probabilistic guarantees of optimality and near-optimality. Such approaches can be useful beyond acceleration, for example, in understanding how to avoid measure zero worst-case scenarios that plague methods such as QR matrix factorization.

97 MATHEMATICS AND COMPUTING↗

NewLife Nuclear - An Environmentally and Economically Minded Solution for Fusion Energy Waste Handling

Energy demand is rising as a result of innovative and increasingly more energy intensive processes coming to fruition, particularly through the recent interest in the development of AI data centers as well as manufacturing with the push towards increasing domestic manufacturing interest. Fusion energy can provide virtually limitless energy to support this increase in energy demand. Fusion energy concepts, largely classified as magnetic fusion energy (MFE) and inertial fusion energy (IFE) are being pursued, each having unique challenges to overcome before the successful deployment of electricity to the grid. Achieving fusion ignition on the National Ignition Facility, first in December 2022, and eight times since, has demonstrated the scientific viability of the IFE approach. Meanwhile, MFE test stands continue to improve confinement times, making meaningful strides in progressing towards experimental scientific viability. In each of these approaches, an emphasis is placed on generating more power out of the system than what is required to power the system. An under-researched area applicable to both IFE and MFE is handling activated waste coming out of fusion energy systems, both in the course of normal daily operations, as well as in intermittent periods as structural materials may need to be replaced. In the context of an IFE plant system, commonly discussed plant designs suggest targets are ignited within a chamber at a rate of up to one million targets per day. Between each shot, the chamber housing the ignition event will clear a portion of the chamber – resulting in a mixture of vaporized target gas, target debris, and other materials being expelled from the chamber [source]. Additionally, IFE system concepts typically discuss the modularization of plant designs, which are expected to be replaced periodically as the components degrade over time. This would result in the irradiated chamber structure materials, likely metals and alloys, needing to be removed and safely stored. In MFE plant systems, while targets are not ignited at a repetition rate with the frequent chamber clearing as is expected in IFE plant systems, it is anticipated that portions of the confinement area interfacing with the hot plasma will need to be replaced periodically. In each system, without additional investment and research into alternative processing and recycling methods, the result is storing irradiated materials, and other elements in a safe containment area until they are no longer activated. – resulting in significant waste both economic and environmental.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗