Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Clock”

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 1,243 records · Page 69

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗

TCLK Must Stay! CAMAC Must Go! How Does Fermilab Move Forward

The current Timing System at Fermilab has been around for 40 years and currently relies on 7 CAMAC crates and over 100 CAMAC cards to produce the Tevatron Clock (TCLK). Thanks to the ingenuity of those before us, this has allowed Fermilab the flexibility to change the timing and Events for its accelerator as beamlines and projects have changed over the years. With the advent of the Proton Improvement Plan-II (PIP-II), the Timing System at Fermilab is being reimagined into a single chassis with even greater flexibility and functionality for decades to come while tackling the ever-challenging task of maintaining backwards compatibility.

Austin, M. R. [Fermilab]↗

Exposure to ionizing radiation disrupts normal epigenetic aging in Japanese medaka

Alterations to the epigenome are a hallmark of biological aging and age-dependent patterning of the DNA methylome (“epigenetic aging”) can be modeled to produce epigenetic age predictors. Rates of epigenetic aging vary amongst individuals and correlate to the onset of age-related disease and all-cause mortality. Yet, the origins of epigenetic-to-chronological age discordance are not empirically resolved. Here, we investigate the relationship between aging, DNA methylation, and environmental exposures in Japanese medaka (Oryzias latipes). We find age-associated DNA methylation patterning enriched in genomic regions of low CpG density and that, similar to mammals, most age-related changes occur during early life. We construct an epigenetic clock capable of predicting chronological age with a mean error of 61.1 days (~8.4% of average lifespan). To test the role of environmental factors in driving epigenetic age variation, we exposed medaka to chronic, environmentally relevant doses of ionizing radiation. Because most organisms share an evolutionary history with ionizing radiation, we hypothesized that exposure would reveal fundamental insights into environment-by-epigenetic aging interactions. Radiation exposure disrupted epigenetic aging by accelerating and decelerating normal age-associated patterning and was most pronounced in cytosines that were moderately associated with age. These findings empirically demonstrate the role of DNA methylation in integrating environmental factors into aging trajectories.

60 APPLIED LIFE SCIENCES↗

Effective fractonic behavior in a two-dimensional exactly solvable spin liquid

In this work we propose a \mathbb{Z}_N ℤ N clock model which is exactly solvable on the lattice. We find exotic properties for the low-energy physics, such as UV/IR mixing and excitations with restricted mobility, that resemble fractonic physics from higher dimensional models. We then study the continuum descriptions for the lattice system in two distinct regimes and find two qualitative distinct field theories for each one of them. A characteristic time scale that grows exponentially fast with N^2 N 2 (and diverges rapidly as function of system parameters) separates these two regimes. For times below this scale, the system is described by an effective fractonic Chern-Simons-like action, where higher-form symmetries prevent quasiparticles from hopping. In this regime, the system behaves effectively as a fracton as isolated particles, in practice, never leave their original position. Beyond the large characteristic time scale, the excitations are mobile and the effective field theory is given by a pure mutual Chern-Simons action. In this regime, the UV/IR properties of the system is captured by a peculiar realization of the translation group.

Delfino, Guilherme↗

Parallel-in-Time Multigrid Methods for Hyperbolic Problems, with a Focus on the Shallow Water Equations

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations, which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic problems of programmatic interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic problems on future architectures. The difficulty lies in the fact that hyperbolic problems are well-known to be difficult for parallel-in-time methods, with the most common parallel-in-time method, Parareal, diverging in many cases. Here, the chosen methodology for scalably solving these hyperbolic space-time equations is multigrid-reductionin-time (MGRIT), because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. To further research in this area, this project will explore a model hyperbolic problem (the shallow water equations) in the context of recent advances (e.g., by Wingate and Haut) in constructing improved coarse time-grid time-propagators by using the slow asymptotic structure of the equations. In particular, we will investigate if such coarse time-propagators offer significant advantages in an MGRIT setting, and explore any broader insights gained into parallel-in-time for hyperbolic problems.

97 MATHEMATICS AND COMPUTING↗

Antifoam Development for Eliminating Flammability Hazards and Decreasing Cycle Time in the Defense Waste Processing Facility

The Savannah River National Laboratory (SRNL) was requested to develop a new antifoam control method for the Defense Waste Processing Facility’s (DWPF) Chemical Process Cell (CPC). SRNL completed testing of both chemical and nonchemical foam controls. The nonchemical foam controls were either ineffective (or worse, created more foam) or impractical (a water spray can control foam, but excessive water is needed). As a result, the focus of this study was on finding a superwetter or commercial antifoam for controlling foam. Thirty potential antifoams were tested as part of this study. A series of tests were developed to help screen out ineffective alternatives including: 1. Spreading testing of superspreaders, 2. Foam column testing with physical simulants, 3. Boiling testing with physical and chemical simulants, 4. Days-only Sludge Receipt and Adjustment Tank (SRAT) process simulations with sludge (containing noble metals and mercury), Precipitate Reactor Feed Tank (PRFT), and Slurry Mix Evaporator Feed Tank (SEFT) simulants in the RC1 Reaction Calorimeter (purchased for antifoam testing), and 5. Around-the-clock SRAT and Slurry Mix Evaporator (SME) process simulations with sludge(containing noble metals and mercury), PRFT, and SEFT simulants in the RC1 Reaction Calorimeter. Evonik Surfynol® MD20, a commercially available defoamer, was relatively effective in controlling foam, while remaining chemically stable in SRAT and SME processing across the pH range of 4 to 13. No degradation products were detected in the offgas, in the condensate or in the SRAT and SME products. In nitric-glycolic acid flowsheet testing, 250 mg/kg Evonik Surfynol® MD20 was needed for foam control compared to 1,625 mg/kg for Antifoam 747, DWPF’s current antifoam. In nitric-formic acid flowsheet testing, 1,125 mg/kg of Evonik Surfynol® MD20 was needed to control foam throughout the SRAT and SME cycles. The commercially available superspreader Momentive™ Y-17112 was even more effective than Evonik Surfynol® MD20 as both a defoamer and an antifoam. Not only was the foam destroyed upon addition but also was less persistent between additions. It was the most effective antifoam in testing using both the nitric-glycolic acid flowsheet and the nitric-formic acid flowsheet. In nitric-glycolic acid flowsheet testing, only 100 mg/kg Momentive™ Y-17112 was needed to control foam throughout the SRAT and SME cycles. In nitric-formic acid flowsheet testing, 300 mg/kg Momentive™ Y-17112 was needed to control foam throughout the SRAT and SME cycles. Momentive™ Y-17112 is also resistant to hydrolysis as demonstrated by its chemical stability in SRAT and SME processing across the pH range of 4 to 13 and lack of degradation products in offgas or condensate. Both candidates were effective as potential replacements for Antifoam 747, with Y-17112 demonstrating superior foam control. During nitric-glycolic flowsheet testing 50% less antifoam was needed when using Momentive™ Y-17112 compared to MD20. During nitric-formic flowsheet testing 75% less antifoam was needed when using Momentive™ Y-17112 compared to MD20. Foam remediated with Momentive™ Y-17112 was less persistent throughout testing. In addition, no degradation products were detected in the offgas, in the condensate or in the SRAT and SME products. Based on this testing, Momentive™ Y-17112 is clearly superior to Evonik Surfynol® MD20 and Antifoam 747, especially for the nitric-formic acid flowsheet processing; it is recommended that Momentive™ Y-17112 replace Antifoam 747 in DWPF. An antifoam addition strategy is recommended for both the nitric-glycolic acid flowsheet and the nitric-formic acid flowsheet. Implementation of Momentive™ Y-17112 is expected to decrease SRAT and SME boiling times by up to 50%, eliminate the issues resulting from Antifoam 747 degradation products, and minimize foamovers. To validate the effectiveness of these defoaming agents, SRNL recommends irradiation of a SRAT or SME product simulant containing fresh antifoam. The goal of this testing is to determine whether the irradiation causes decomposition of the antifoam that would make it less effective or produce new species in the offgas or slurry. This testing began in April 2020. An evaluation should be completed to determine the thermolytic hydrogen and methane generation rate in downstream equipment, including the High-Level Waste evaporators.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Digital Application Specific Integrated Circuit Design

Single radiation effects can be transient and static. The charge collected from radiation has a prompt (drift) and a slow (diffusion) component. This current spike causes disparate results: SEU, SET, SEB, MBU, and more. A system level SEE mitigation technique, used to overcome, is the triple module redundancy (TMR). TMR for system Verilog (TMR4SV), written in Python, is the first step towards design IP triplication. It allows various triplication methods for flip flops, clocks, and logic.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Battery Monitoring System

The component that is powered by the battery pack being monitored is a valuable asset and must be in working condition at all times. Battery chemistry and characteristics have a major role in how to evaluate the state of the battery. The battery monitoring system has many parts that lead to an accurate battery reading. The components consist of a coulomb counting device, end of life voltage detection, a consideration of use for a real-time clock (RTC), temperature sensor, and non-volatile random-access memory (NVRAM). The combination of these elements allows the monitoring system to be highly reliant. Moving forward a better implementation of the ideas in this paper and further testing should ensure a high-quality battery monitoring system.

25 ENERGY STORAGE↗

Parallel Time Integration for Constrained Optimization

The number of transistors in an average processor continues to increase, but individual clock speeds have plateaued. Those transistors are instead going into additional cores, increasing the number of different things that a processor can do at once and placing an emphasis on parallel computation. Many problems in scientific computing follow a time-evolution model, and it can be difficult to solve such problems in parallel across the temporal domain. The Multi-Grid Reduction In Time (MGRIT) algorithm, developed at Lawrence Livermore National Laboratory (LLNL), solves differential equations with a method designed specifically to take advantage of extreme numbers of processors by parallelizing across time. The Tri-diagonal MGRIT (TriMGRIT) algorithm, also developed at LLNL, is a generalization of MGRIT which enables parallel-in-time solving of a greater number of problems. Constrained optimization problems, in particular, may be solved in parallel using TriMGRIT. These consist of choosing a control function such that an objective functional is minimized, constrained by a differential-equation. We consider two such problems: applying torque to a pendulum to bring it to a gentle stop and moving a crowd of people from one distribution into another. We also perform some miscellaneous theoretical and practical research, including investigating the use of a line-search subroutine to refine intermediate TriMGRIT results and preliminary work on strategies for choosing operators for TriMGRIT to use.

97 MATHEMATICS AND COMPUTING↗

Accomplishments of Sandia and Kitware CMake/CTest/CDash Contract for (FY2017-2020)

We describe the accomplishments jointly achieved by Kitware and Sandia over the fiscal years 2016 through 2020 to benefit the Advanced Scientific Computed (ASC) Advanced Technology Development and Mitigation (ATDM) project. As a result of our collaboration, we have improved the Trilinos and ATDM application developer experience by decreasing the time to build, making it easier to identify and resolve build and test defects, and addressing other issues . We have also reduced the turnaround time for continuous integration (CI) results. For example, the combined improvements likely cut the wall clock time to run automated builds of Trilinos posting to CDash by approximately 6x or more in many cases. We primarily achieved these benefits by contributing changes to the Kitware CMake/CTest/CDash suite of open source software development support tools. As a result, ASC developers can now spend more time improving code and less time chasing bugs. And, without this work, one can argue that the stabilization of Trilinos for the ATDM platforms would not have been feasible which would have had a large negative impact on an important internal FY20 L1 milestone.

97 MATHEMATICS AND COMPUTING↗

Constitutive Model Development for Aging Polymer Encapsulants (ASC P&EM FY2021 L2 Milestone 7836)

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

36 MATERIALS SCIENCE↗

Parallel Multigrid in Time and Space for Extreme-Scale Computational Science

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations, which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of programmatic interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures.

97 MATHEMATICS AND COMPUTING↗

Assessment of the Electrical Substation-Grid Testbed with Inside/Outside Devices and Distributed Ledger Technology

The electrical substation-grid testbed was created to integrate the GOOSE and/or DNP (Distributed Network Protocol) messages with time synchronized sources and Distributed Ledger Technology (DLT). The objective was to study the impact of faults and cyber-events at an electrical substation with inside (protective relays) and outside (power meters) substation devices. The electrical substation-grid testbed was based on the design of a 34.5/ 12.47 kV electrical substation (sectionalized bus configuration) with two power transformers, connected to radial power lines and load feeders. The electrical substation-grid testbed was installed at 252 lab space (Advanced Power System Protection), Grid Research Integration and Deployment Center (GRID-C), Oak Ridge National Laboratory. This testbed was created for Task 5, DarkNet project. The electrical substation-grid testbed was created to simulate fault and/or cyber events that could potentially result in damage to the electrical infrastructure. In addition, tests were run that are usually not allowed to be performed in an operational electrical power grid, because these test scenarios could trip breakers and/or generate fault situations that could potentially damage equipment. The number of tests performed in the electrical substation-grid testbed were executed in a better way than in a real electrical substation and/or power grid, because multiple tests could be run in a short period of time, and complex permits, and safety/ schedule restrictions like in a real electrical substation environment were not needed. The electrical substation-grid testbed was created using real measurement, communication, and protection devices that are used by electrical utilities, to have same conditions that we could observe in a real power grid or electrical substation. The electrical substation-grid testbed was based on using a real time simulator and expansion box with amplifiers that were wired to electrical substation-grid devices. This hardware-in-the-loop (HIL) was provided by protective relays, power meters, ethernet switches, remote terminal units, synchronized timing network clock, DLT devices, workstations, and servers. This report includes the design, installation, and assessment of the electrical substation-grid testbed that was similar to an operational electrical substation, integrating the power system protection, communication, and control systems. The results for the electrical substation-grid testbed were based on:• verifying the analog signals for protective relays and power meters, • observing the synchronized time source frame at devices, • authenticating the GOOSE (IEC 61850) and DNP messages from power meters and protective relays, and • verifying the trip conditions of protective relays at fault tests with the power system fault event detection, using DLT devices. For future work, the electrical substation-grid testbed with protective relays and power meters, using DLT and synchronized time source from DarkNet, will be used to study the impact of cyber-events at inside and outside substation devices. Advanced algorithms for detecting cyber-events produced by non-desired protective relay settings will be studied, to improve the detection and reliability of protection, control, and communication systems at power grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Calibrating the SPECTACULAR constitutive model using legacy Sandia data for two filled epoxy systems: 828/CTBN/DEA/GMB and 828/DEA/GMB

The SPECTACULAR model is a development extension of the Simplified Potential Energy Clock (SPEC) model. Both models are nonlinear viscoelastic constitutive models used to predict a wide range of time-dependent behaviors in epoxies and other glass-forming materials. This report documents the procedures used to generate SPECTACULAR calibrations for two particulate-filled epoxy systems, 828/CTBN/DEA/GMB and 828/DEA/GMB. No previous SPECTACULAR or SPEC calibration exists for 828/CTBN/DEA/GMB, while a legacy SPEC calibration exists for 828/DEA/GMB. To generate the SPECTACULAR calibrations, a step-by-step procedure was executed to determine parameters in groups with minimal coupling between parameter groups. This procedure has often been deployed to calibrate SPEC, therefore the resulting SPECTACULAR calibration is backwards compatible with SPEC (i.e. none of the extensions specific to SPECTACULAR are used). The calibration procedure used legacy Sandia experimental data stored on the Polymer Properties Database website. The experiments used for calibration included shear master curves, isofrequency temperature sweeps under oscillatory shear, the bulk modulus at room temperature, the thermal strain during a temperature sweep, and compression through yield at multiple temperatures below the glass transition temperature. Overall, the calibrated models fit the experimental data remarkably well. However, the glassy shear modulus varies depending on the experiment used to calibrate it. For instance, the shear master curve, isofrequency temperature sweep under oscillatory shear, and the Young's modulus in glassy compression yield values for the glassy shear modulus at the reference temperature that vary by as much as 15 %. Also, for 828/CTBN/DEA/GMB, the temperature dependence of the glassy shear modulus when fit to the Young's modulus at different temperatures is approximately four times larger than when it is determined from the isofrequency temperature sweep under oscillatory shear. For 828/DEA/GMB, the temperature dependence of the shear modulus determined from the isofrequency temperature sweep under oscillatory shear accurately predicts the Young's modulus at different temperatures. When choosing values for the shear modulus, fitting the glassy compression data was prioritized. The new and legacy calibrations for 828/DEA/GMB are similar and appear to have been calibrated from the same data. However, the new calibration improves the fit to the thermal strain data. In addition to the standard calibrations, development calibrations were produced that take advantage of development features of SPECTACULAR , including an updated equilibrium Helmholtz free energy that eliminates undesirable behavior found in previous work. In addition to the previously mentioned experimental data, the development calibrations require data for the heat capacity during a stress-free temperature sweep to calibrate thermal terms.

36 MATERIALS SCIENCE↗

Parallel Multigrid in Time and Space for Extreme-Scale Computational Science: Chaotic and Hyperbolic Problems

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallelin-time methods, with the most common method, parareal, diverging in many cases. The current stateof-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, we take inspiration from the Least Squares Shadowing (LSS by Wang) approach for solving chaotic problems. Here, an optimization approach is able to find “well-conditioned” shadow trajectories/solutions to the original “ill-conditioned” chaotic problem. Thus, the new multigrid methods researched here also arise in an optimization context.

97 MATHEMATICS AND COMPUTING↗

Molecules Functionalized with Cycling Centers for Quantum Information Science

The ability for a single molecule to produce enough laser-induced fluorescence to control and detect its quantum state is known as optical cycling. This property has been used to great effect in atoms, where the relatively small number of states promotes their ability to repeatedly fluoresce, for applications such as laser cooling, precision measurement, atomic clocks, and quantum information processing. The extension of this capability to molecular species, however, has been difficult, primarily due to their propensity to become vibrationally excited during illumination. This project brought together experts from atomic physics, physical chemistry, and computational quantum chemistry to develop means to create molecules capable of optical cycling. In particular, this program sought to produce an optical cycling center (OCC), a quantum functional group capable of chemical attachment to a wide variety of molecular hosts while retaining (and thereby furnishing to its host) the ability to optically cycle. Following our team’s theoretical identification of promising species, including alkaline-earth phenoxides, we then created these species and in many cases confirmed their optical cycle closure. We now not only know of the existence of a handful of large, molecular species that can be deployed for applications in quantum technologies and measurement, but have improved our knowledge of the molecular mechanisms leading to this discovery, which will enable generalization in the future. The extension of these ideas for informed chemical design of cycling molecules with bespoke properties is now being pursued by multiple groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The mapping & avoidance of high order cavity modes

In this paper a technique developed in 2016 is presented for mapping the multidimensional, multi-mode high order mode landscape of a beam driven super conducting RF cavity. The initial use case for this technique was intended to allow longitudinal gradient control of a beam driven super conducting cavity while simultaneously avoiding many dangerous high order modes and excessive component stresses. This paper will focus on the stability and component stress considerations of the 56MHz superconducting cavity, located in the 4 o’clock sector interaction region of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory. The 56MHz cavity is intended to work in conjunction with the Stochastic Cooling system to improve the longitudinal confinement in the central 197MHz RF bucket and thus increase the integrated luminosity during heavy ion beam runs at RHIC.

43 PARTICLE ACCELERATORS↗

Intelligent Proppants for Hydraulic Fracture Diagnostics

This Phase 1 Report details our Phase 1 progress in proving the feasibility of a highly integrated, highly miniaturized sensor ball. The ball contains sensors, analog to digital converters, memory, a clock, and an acoustic communications module all integrated onto a single silicon chip. Phase 2 is to focus on the application specific integrated circuit (ASIC). Phase 1 piecewise proves the feasibility of this innovation by examining a few critical subsystems including the communications module and inertial sensor performance. Acoustic simulations of an acoustically active and traceable proppant are also performed and summarized in this report.

01 COAL, LIGNITE, AND PEAT↗