Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “future of computing”

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 865 records · Page 48

General formulas for adiabatic invariants in nearly periodic Hamiltonian systems

While it is well known that every nearly periodic Hamiltonian system possesses an adiabatic invariant, extant methods for computing terms in the adiabatic invariant series are inefficient. The most popular method involves the heavy intermediate calculation of a non-unique near-identity coordinate transformation, even though the adiabatic invariant itself is a uniquely defined scalar. Here, a less well-known method, developed by S. Omohundro, avoids calculating intermediate sequences of coordinate transformations but is also inefficient as it involves its own sequence of complex intermediate calculations. In order to improve the efficiency of future calculations of adiabatic invariants, we derive generally applicable, readily computable formulas for the first several terms in the adiabatic invariant series. To demonstrate the utility of these formulas, we apply them to charged-particle dynamics in a strong magnetic field and magnetic field-line dynamics when the field lines are nearly closed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

2014 HPC Annual Report

Our commitment is to support you through delivery of an IT environment that provides mission value by transforming the way you use, protect, and access information. We approach this through technical innovation, risk management, and relationships with our workforce, Laboratories leadership, and policy makers nationwide. This second edition of our HPC Annual Report continues our commitment to communicate the details and impact of Sandia’s large-scale computing resources that support the programs associated with our diverse mission areas. A key tenet to our approach is to work with our mission partners to understand and anticipate their requirements and formulate an investment strategy that is aligned with those Laboratories priorities. In doing this, our investments include not only expanding the resources available for scientific computing and modeling and simulation, but also acquiring large-scale systems for data analytics, cloud computing, and Emulytics. We are also investigating new computer architectures in our advanced systems test bed to guide future platform designs and prepare for changes in our code development models. Our initial investments in large-scale institutional platforms that are optimized for Informatics and Emulytics work are serving a diverse customer base. We anticipate continued growth and expansion of these resources in the coming years as the use of these analytic techniques expands across our mission space. If your program could benefit from an investment in innovative systems, please work through your Program Management Unit ’s Mission Computing Council representatives to engage our teams.

97 MATHEMATICS AND COMPUTING↗

A Dynamic Waste Isolation Pilot Plant Performance Assessment Tool - 12490

The Waste Isolation Pilot Plant (WIPP) Performance Assessment (PA) methodology comprises a toolbox used to demonstrate regulatory compliance of the repository after facility closure. The PA framework rests upon an extensive suite of computational codes. In some cases, significant alteration of code inputs is a tedious and difficult task. Due to the nature of the application for which they are used, PA codes used in support of WIPP regulatory compliance demonstration must satisfy stringent quality assurance requirements. Consequently, many of the coding practices used during original code development are still implemented today. A more efficient workflow configuration has the potential to alleviate difficulties associated with extensive code input modifications. Here, this potential is assessed via an implementation of a more flexible scientific workflow system for a subset of the codes used in WIPP PA. The scientific workflow approach taken here for a dynamic PA system enables users from disparate backgrounds to dramatically shorten the time between hypothesis and analysis by decreasing the amount of a priori knowledge, from a range of disciplines, needed to execute the code. Having smaller iteration times allows for more ideas to be tested and explored, which leads to safer and more optimized systems. Note that these high-level, dynamic tools are intended only for initial scoping studies on the personal computer of a researcher. Full, regulatory compliance calculations may occur only within a qualified computing environment. However, the WIPP PA tools here may guide future research and indicate regions of the analysis space that are worth further study. This next generation of PA software provides the ability to perform scoping investigations of repository performance quickly and easily, and has an accessible and useful interface to a variety of users, such as fuel cycle systems designers, domain experts such as repository modelers, and policy makers. The purview of this project allows for many opportunities for future work. Foremost among these is the desire to implement the full BRAGFLO suite within the workflow. This will entail porting or wrapping Genmesh, Matset, LHS, and ICSet within Python. Moreover, unifying the two GUIs into a single driver application would be a natural next step. Once the BRAGFLO suite is completed, other portions of WIPP PA could be implemented with corresponding and inter-operable work-flows. Likely first candidates for this are those codes that are similarly computationally intensive, such as the one used to generate complementary cumulative distribution functions used to demonstrate regulatory compliance (code CCDFGF). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reasoning About Computing in the Post-Moore world [Slides]

We are facing new challenges due to technology shifts and evolving mission requirements. Meeting these challenges requires careful reasoning about hardware and our applications. This is important for decisions we’re making for our codes following ATDM and our decisions about future platforms.

97 MATHEMATICS AND COMPUTING↗

Calculation of exact vibrational spectra for P{sub 2}O and CH{sub 2}NH using a phase space wavelet basis

‘‘Exact” quantum dynamics calculations of vibrational spectra are performed for two molecular systems of widely varying dimensionality (P{sub 2}O and CH{sub 2}NH), using a momentum-symmetrized Gaussian basis. This basis has been previously shown to defeat exponential scaling of computational cost with system dimensionality. The calculations were performed using the new “SWITCHBLADE” black-box code, which utilizes both dimensionally independent algorithms and massive parallelization to compute very large numbers of eigenstates for any fourth-order force field potential, in a single calculation. For both molecules considered here, many thousands of vibrationally excited states were computed, to at least an “intermediate” level of accuracy (tens of wavenumbers). Future modifications to increase the accuracy to “spectroscopic” levels, along with other potential future improvements of the new code, are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Information barrier experimental: Toward a trusted and open-source computing platform for nuclear warhead verification

Trusted measurement systems are critical for the verification of future arms-control treaties that involve confirming the authenticity of nuclear warheads on the basis of their radiation signatures. Several research efforts have produced prototype systems, but their designs are typically not published in detail, making it difficult to enable trust in these devices. The Information Barrier Experimental (IBX) is a new prototype built around the Red Pitaya computing platform using passive gamma spectroscopy. It is the first such platform designed to help a broad research community study vulnerabilities and define the required specifications for a common, trusted inspection system. It is low-cost, simple to assemble, and enables comprehensive hardware and software authentication studies. The device follows a digital data acquisition approach, which significantly reduces the number of components between scintillator and spectrum output. Finally, measurements demonstrate that this approach is technically feasible and produces excellent measurement results.

97 MATHEMATICS AND COMPUTING↗

Large-scale Stochastic Optimization for Unit Commitment and Economic Dispatch

Our FY 2012 LDRD research on stochastic optimization with its application in power grid planning and operations has advanced the state-of-the-art of this area in the following aspects: (1) Deploying stochastic unit commitment on large real-world applications, (2) Stochastic optimization coupled with real-time production cost simulation (interleaved), and (3) Real-time stochastic optimization with rolling-horizon look-ahead. In conjunction, a number of computational experiments have been performed to evaluate the capability of the latest stochastic optimization software and explore opportunities for future research. In what follows, we discuss each of these contribution areas, including the results of our computational experiments and their implications on the scalability and solvability of large-scale stochastic optimization models.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Large-Scale Stochastic Optimization for Unit Commitment and Economic Dispatch

Our FY 2012 LDRD research on stochastic optimization with its application in power grid planning and operations has advanced the state-of-the-art of this area in the following aspects: (1) Deploying stochastic unit commitment on large real-world applications, (2) Stochastic optimization coupled with real-time production costsimulation (interleaved), and (3) Real-time stochastic optimization with rolling-horizon look-ahead. In conjunction, a number of computational experiments have been performed to evaluate the capability of the latest stochastic optimization software and explore opportunities for future research. In what follows, we discuss each of these contribution areas, including the results of our computational experiments and their implications on the scalability and solvability of large-scale stochastic optimization models.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Projected Changes in Temperature and Precipitation Over the United States, Central America, and the Caribbean in CMIP6 GCMs

The Coupled Model Intercomparison Project Phase 6 (CMIP6) dataset is used to examine projected changes in temperature and precipitation over the United States (U.S.), Central America and the Caribbean. The changes are computed using an ensemble of 31 models for three future time slices (2021–2040, 2041–2060, and 2080–2099) relative to the reference period (1995–2014) under three Shared Socioeconomic Pathways (SSPs; SSP1-2.6, SSP2-4.5, and SSP5-8.5). The CMIP6 ensemble reproduces the observed annual cycle and distribution of mean annual temperature and precipitation with biases between – 0.93 and 1.27 °C and – 37.90 to 58.45%, respectively, for most of the region. However, modeled precipitation is too large over the western and Midwestern U.S. during winter and spring and over the North American monsoon region in summer, while too small over southern Central America. Temperature is projected to increase over the entire domain under all three SSPs, by as much as 6 °C under SSP5-8.5, and with more pronounced increases in the northern latitudes over the regions that receive snow in the present climate. Annual precipitation projections for the end of the twenty-first century have more uncertainty, as expected, and exhibit a meridional dipole-like pattern, with precipitation increasing by 10–30% over much of the U.S. and decreasing by 10–40% over Central America and the Caribbean, especially over the monsoon region. Seasonally, precipitation over the eastern and central subregions is projected to increase during winter and spring and decrease during summer and autumn. Over the monsoon region and Central America, precipitation is projected to decrease in all seasons except autumn. The analysis was repeated on a subset of 9 models with the best performance in the reference period; however, no significant difference was found, suggesting that model bias is not strongly influencing the projections.

54 ENVIRONMENTAL SCIENCES↗

Virtual prototyping of liquid metal blanket performance in fusion pilot plant

Liquid metal blanket is a dominant design option for the next step fusion devices responsible for harvesting energy from fusion reaction, and simultaneously producing fuel for the same reaction through tritium breeding. Liquid metal blankets introduce additional complexity to the design due to fluid motion, fluid structure interaction, and magnetohydrodynamic (MHD) effects arising from the motion of the conducting fluid through the magnetic field. They are also directly affected by the plasma heat flux and neutronic fluence. PPPL is currently developing a virtual prototyping system for numerical analysis of the liquid metal blankets for future fusion devices. The system has a customized 3D computational fluid dynamics (CFD) code in its core, allowing MHD flow and conjugate heat transfer analysis in blankets fluids and solids. The code was successfully used before for dual coolant blanket analysis [A. Khodak et al., Fusion Eng. and Des. 137 (2018)]. Recently the same code was modified to allow verified simulation of MHD flows at high Hartmann numbers of several thousand typical for blanket applications. CFD code receives volumetric heat source distribution from the neutronic analysis based on MCNP code. In addition, direct tritium breeding simulation will be performed allowing optimization of the blanket performance. 2D axisymmetric version of neutronics code will be used for rapid optimization, with 3D version employed for detailed analysis. The surface heat distribution on the plasma facing wall will be defined by the software HEAT allowing 3D modeling of the heat flux based on the magnetic field distribution including gyro-orbit effects. Results of thermal analysis are imported into structural analysis code also included in the system. Finally, direct import of CAD geometry will be used for analyzing all components and as a result design option can be efficiently optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predictions of delafossite-hosted honeycomb and kagome phases

Delafossites, typically denoted by the formula ABO 2 , are a class of layered materials that exhibit a wide range of electronic and optical properties. Recently, the idea of modifying these delafossites into ordered kagome or honeycomb phases via strategic doping has emerged as a potential way to tailor these properties. In this study, we use high-throughput density functional theory calculations to explore many possible candidate kagome and honeycomb phases by considering dopants selected from the parent compounds of known ternary delafossite oxides from the inorganic crystal structure database. Here our results indicate that while A-site in existing delafossites can host a limited range of elemental specifies, and display a low propensity for mixing or ordering, the oxide sub-units in the BO 2 much more readily admit guest species. Our study identifies four candidate B-site kagome and fifteen candidate B-site honeycombs with a formation energy more than 50 meV f.u. −1 below other competing phases. The ability to predict and control the formation of these unique structures offers exciting opportunities in materials design, where innovative properties can be engineered through the selection of specific dopants. A number of these constitute novel correlated metals, which may be of interest for subsequent efforts in synthesis. These novel correlated metals may have significant implications for quantum computing, spintronics, and high-temperature superconductivity, thus inspiring future experimental synthesis and characterization of these proposed materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The order of three lowest-energy states of the six-electron harmonium at small force constant

The order of low-energy states of six-electron harmonium is uncertain in the case of strong correlation, which is not a desired situation for the model system being considered for future testing of approximate methods of quantum chemistry. The computational study of these states has been carried out at the frequency parameter ω = 0.01, using the variational method with the basis of symmetry-projected, explicitly correlated Gaussian (ECG) lobe functions. It has revealed that the six-electron harmonium at this confinement strength is an octahedral Wigner molecule, whose order of states is different than in the strong confinement regime and does not agree with the earlier predictions. The results obtained for ω = 0.5 and 10 are consistent with the findings based on the Hund’s rules for the s{sup 2}p{sup 4} electron configuration. Substantial part of the computations has been carried out on the graphical processing units and the efficiency of these devices in calculation of the integrals over ECG functions has been compared with traditional processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics of E × B discharges relevant to plasma propulsion and similar technologies

This paper provides perspectives on recent progress in understanding the physics of devices in which the external magnetic field is applied perpendicular to the discharge current. This configuration generates a strong electric field that acts to accelerate ions. The many applications of this set up include generation of thrust for spacecraft propulsion and separation of species in plasma mass separation devices. These “E × B” plasmas are subject to plasma–wall interaction effects and to various micro- and macroinstabilities. In many devices we also observe the emergence of anomalous transport. This perspective presents the current understanding of the physics of these phenomena and state-of-the-art computational results, identifies critical questions, and suggests directions for future research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global estimates of the fitness advantage of SARS-CoV-2 variant Omicron

Abstract New variants of SARS-CoV-2 show remarkable heterogeneity in their relative fitness over both time and space. In this paper we extend the tools available for estimating the selection strength for new SARS-CoV-2 variants to a hierarchical, mixed-effects, renewal equation model. This formulation allows us to estimate selection effects at the global level while incorporating both measured and unmeasured heterogeneity among countries. Applying this model to the spread of Omicron in forty countries, we find evidence for very strong but very heterogeneous selection effects. To test whether this heterogeneity is explained by differences in the immune landscape, we considered several measures of vaccination rates and recent population-level infection as covariates, finding moderately strong, statistically significant effects. We also found a significant positive correlation between the selection advantage of Delta and Omicron at the country level, suggesting that other region-specific explanatory variables of fitness differences do exist. Our method is implemented in the Stan programming language, can be run on standard consumer-grade computing resources, and will be straightforward to apply to future variants.

59 BASIC BIOLOGICAL SCIENCES↗

Hybrid Programming for Near-Term Quantum Computing Systems

Recent computations involving quantum processing units (QPUs) have demonstrated a series of challenges inherent to hybrid classical-quantum programming, compilation, execution, and verification and validation. Despite considerable progress, system-level noise, limited low-level instructions sets, remote access models, and an overall lack of portability and classical integration presents near-term programming challenges that must be overcome in order to enable reliable scientific quantum computing and support robust hardware benchmarking. In this work, we draw on our experience in programming QPUs to identify common concerns and challenges, and detail best practices for mitigating these challenges within the current hybrid classical-quantum computing paradigm. Following this discussion, we introduce the XACC quantum compilation and execution framework as a hardware and language independent solution that addresses many of these hybrid programming challenges. XACC supports extensible methodologies for managing a variety of programming, compilation, and execution concerns across the increasingly diverse set of QPUs. We use recent nuclear physics simulations to illustrate how the framework mitigates programming, compilation, and execution challenges and manages the complex workflow present in QPU-enhanced scientific applications. Finally, we codify the resulting hybrid scientific computing workflow in order to identify key areas requiring future improvement.

Mccaskey, Alex↗

SU-D-17A-04: The Impact of Audiovisual Biofeedback On Image Quality During 4D Functional and Anatomic Imaging: Results of a Prospective Clinical Trial

Purpose: The ability of audiovisual (AV) biofeedback to improve breathing regularity has not previously been investigated for functional imaging studies. The purpose of this study was to investigate the impact of AV biofeedback on 4D-PET and 4D-CT image quality in a prospective clinical trial. We hypothesized that motion blurring in 4D-PET images and the number of artifacts in 4D-CT images are reduced using AV biofeedback. Methods: AV biofeedback is a real-time, interactive and personalized system designed to help a patient self-regulate his/her breathing using a patient-specific representative waveform and musical guides. In an IRB-approved prospective clinical trial, 4D-PET and 4D-CT images of 10 lung cancer patients were acquired with AV biofeedback (AV) and free breathing (FB). The 4D-PET images in 6 respiratory bins were analyzed for motion blurring by: (1) decrease of GTVPET and (2) increase of SUVmax in 4-DPET compared to 3D-PET. The 4D-CT images were analyzed for artifacts by: (1) comparing normalized cross correlation-based scores (NCCS); and (2) quantifying a visual assessment score (VAS). A two-tailed paired t-test was used to test the hypotheses. Results: The impact of AV biofeedback on 4D-PET and 4D-CT images varied widely between patients, suggesting inconsistent patient comprehension and capability. Overall, the 4D-PET decrease of GTVPET was 2.0±3.0cm3 with AV and 2.3±3.9cm{sup 3} for FB (p=0.61). The 4D-PET increase of SUVmax was 1.6±1.0 with AV and 1.1±0.8 with FB (p=0.002). The 4D-CT NCCS were 0.65±0.27 with AV and 0.60±0.32 for FB (p=0.32). The 4D-CT VAS was 0.0±2.7 (p=ns). Conclusion: A 10-patient study demonstrated a statistically significant reduction of motion blurring of AV over FB for 1/2 functional 4D-PET imaging metrics. No difference between AV and FB was found for 2 anatomic 4D-CT imaging metrics. Future studies will focus on optimizing the human-computer interface and including patient training sessions for improved comprehension and capability. Supported by NIH/NCI R01 CA 093626, Stanford BioX Interdisciplinary Initiatives Program, NHMRC Australia Fellowship, and Kwanjeong Educational Foundation. GE Healthcare provided the Respiratory Gating Toolbox for 4D-PET image reconstruction. Stanford University owns US patent #E7955270 which is unlicensed to any commercial entity.

60 APPLIED LIFE SCIENCES↗

Insight from Fukushima Daiichi Unit 3 Investigations using MELCOR

During the emergency response period of the accidents that took place at Fukushima Daiichi in March of 2011, researchers at Oak Ridge National Laboratory (ORNL) conducted a number of studies using the MELCOR code to help understand what was occurring and what had occurred. During the post-accident period, the Department of Energy (DOE) and the US Nuclear Regulatory Commission (NRC) jointly sponsored a study of the Fukushima Daiichi accident with collaboration among Oak Ridge, Sandia, and Idaho national laboratories. The purpose of the study was to compile relevant data, reconstruct the accident progression using computer codes, assess the codes predictive capabilities, and identify future data needs. The current paper summarizes some of the early MELCOR simulations and analyses conducted at ORNL of the Fukushima Daiichi Unit 3 accident. Extended analysis and discussion of the Unit 3 accident is also presented taking into account new knowledge and modeling refinements made since the joint DOE/NRC study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗